苏里格气田盒8气藏气水层测井识别研究

鄂尔多斯盆地苏里格气田盒8气藏为弹性驱动的河流相低效气藏,储层非均质性强,在开发中普遍见水。根据试采情况将气藏储层分为气层、水层、气水层三类,开展气、水层的测井识别研究。以区内26口井50个射孔层段各层的试油结果作为依据,选出了700个样本点作为判识建模的标本,以声波时差和中子孔隙度等8个原始测井值作为判识变量组合,采用模糊判别分析,建立了各层的判识模型。回判结果表明,模型对气层精度达到97％,水层为89％,气水层为94％。利用未参加建模的出水井射孔段对模型进行验证,判识结果基本符合实际情况,模型具有较高精度。     

The dynamics and thermodynamics of large ash flows

 Ash flow deposits, containing up to 1000 km³ of material, have been produced by some of the largest volcanic eruptions known. Ash flows propagate several tens of kilometres from their source vents, produce extensive blankets of ash and are able to surmount topographic barriers hundreds of metres high. We present and test a new model of the motion of such flows as they propagate over a near horizontal surface from a collapsing fountain above a volcanic vent. The model predicts that for a given eruption rate, either a slow (10–100 m/s) and deep (1000–3000 m) subcritical flow or a fast (100–200 m/s) and shallow (500–1000 m) supercritical flow may develop. Subcritical ash flows propagate with a nearly constant volume flux, whereas supercritical flows entrain air and become progressively more voluminous. The run-out distance of such ash flows is controlled largely by the mass of air mixed into the collapsing fountain, the degree of fragmentation and the associated rate of loss of material into an underlying concentrated depositional system, and the mass eruption rate. However, in supercritical flows, the continued entrainment of air exerts a further important control on the flow evolution. Model predictions show that the run-out distance decreases with the mass of air entrained into the flow. Also, the mass of ash which may ascend from the flow into a buoyant coignimbrite cloud increases as more air is entrained into the flow. As a result, supercritical ash flows typically have shorter runout distances and more ash is elutriated into the associated coignimbrite eruption columns. We also show that one-dimensional, channellized ash flows typically propagate further than their radially spreading counterparts. As a Plinian eruption proceeds, the erupted mass flux often increases, leading to column collapse and the formation of pumiceous ash flows. Near the critical conditions for eruption column collapse, the flows are shed from high fountains which entrain large quantities of air per unit mass. Our model suggests that this will lead to relatively short ash flows with much of the erupted material being elutriated into the coignimbrite column. However, if the mass flux subseqently increases, then less air per unit mass is entrained into the collapsing fountain, and progressively larger flows, which propagate further from the vent, will develop. Our model is consistent with observations of a number of pyroclastic flow deposits, including the 1912 eruption of Katmai and the 1991 eruption of Pinatubo. The model suggests that many extensive flow sheets were emplaced from eruptions with mass fluxes of 10⁹–10¹⁰ kg/s over periods of 10³–10⁵ s, and that some indicators of flow "mobility" may need to be reinterpreted. Furthermore, in accordance with observations, the model predicts that the coignimbrite eruption columns produced from such ash flows rose between 20 and 40 km. Received: 25 August 1995 / Accepted: 3 April 1996     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

华北典型箕状断陷盆地层序地层学模式及其与油气赋存关系

本文根据层序地层学的基本原理,立足东部箕状断陷盆地的特点,以冀中拗陷和二连盆地为重点研究对象,提炼出陆相断陷盆地层序地层学模式,在缓坡带,以风暴浪基面作为“陆架边角”的参照物,层序地层特征相似于海相被动大陆边缘背景;在陡坡带则截然不同,并以湖面与边界断裂的交点取代“陆架边角”,各体系域所含沉积体系与海相模式不尽相同,在盆地形成初期的层序,低水位体系域发育,湖泊消亡以后,主要发育高水位体系域。此外,本文还探索了层序地层与油气赋存的关系,指出富砂的“盆底扇”是今后勘探开发的重要方向。     

珠江口盆地西江主洼油气差异分布机制

西江主洼是珠江口盆地一个低勘探程度洼陷，油气差异聚集特征明显，但其机理不清。本文在断陷盆地油气成藏理论指导下，利用研究区地质、地球物理和地球化学等资料，开展了烃源岩、断盖组合、储集体系和运聚模式等方面的研究。结果表明，裂陷期沉降、沉积中心有序迁移，造成主力烃源岩自东向西由文四段迁移至文三段、文一+二段，控制东、西部油气差异分布。恩平组区域性泥岩发育且晚期断裂缺乏，导致油气纵向上更易聚集于下构造层；区域泥岩减薄尖灭或晚期断裂切开盖层的区域，上构造层有一定油气分布。"源-汇"类型从宏观上控制储层优劣，影响下构造层油气富集程度；"仓储"运移是控制上构造层油气规模聚集的主要模式。下构造层古近系应作为重点部署方向，东部围绕文四段、西部围绕文三段、文一+二段烃源岩构成的含油气系统展开；上构造层勘探需关注珠海组，在隆起周边寻找具备"仓储"运移模式的有利区带。     

地铁施工对济南白泉泉群流量的影响

为研究地铁建设对济南白泉泉群的影响，在综合分析白泉泉域地质、水文地质条件的基础上，假定研究区岩溶强径流带位置及水力性质，利用FEFLOW软件建立地下水流数值模型。以规划地铁M1号线为研究对象，分析了济南东站、梁王站、梁王东站分别施工及3个站同时施工4种情景下，采用施工降水或施工降水+人工回灌两种施工方式对白泉泉群流量的影响。结果表明：单独采用施工降水的施工方式使得白泉泉群流量衰减，其中3个站同时施工对泉流量的影响最大，泉流量最大衰减达5.48%；各站分别施工时，济南东站对泉流量影响最大，泉流量较未施工时减少了0.043×10⁴ m³/d。采用施工降水+人工回灌的施工方式，能够有效缓解泉流量的衰减，各车站施工时的泉流量衰减由仅施工降水时的2.26%~5.48%降低至0.08%~1.21%。岩溶强径流带有利于地下水形成优势径流，促进白泉泉群补给，一定程度上缓解因地铁施工引起的泉流量衰减。     

A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.     

福建省第三产业发展的区域比较分析

运用可拓工程方法和主成分分析方法建立了福建省第三产业发展水平的综合评判物元模型。模型的输出结果表明：福建省各地区第三产业发展水平差距悬殊；综合关联度将第三产业发展水平划分成4个等级，单指标关联度则揭示了各地区第三产业发展的优、劣势。     

Large eddy simulation of hydrocyclone—prediction of air-core diameter and shape

Hydrocyclones are widely used in the mining and chemical industries. An attempt has been made in this study, to develop a CFD (computational fluid dynamics) model, which is capable of predicting the flow patterns inside the hydrocyclone, including accurate prediction of flow split as well as the size of the air-core. The flow velocities and air-core diameters are predicted by DRSM (differential Reynolds stress model) and LES (large eddy simulations) models were compared to experimental results. The predicted water splits and air-core diameter with LES and RSM turbulence models along with VOF (volume of fluid) model for the air phase, through the outlets for various inlet pressures were also analyzed. The LES turbulence model led to an improved turbulence field prediction and thereby to more accurate prediction of pressure and velocity fields. This improvement was distinctive for the axial profile of pressure, indicating that air-core development is principally a transport effect rather than a pressure effect.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.