Multiphase flow modeling with general boundary conditions and automatic phase-configuration changes using a fractional-flow approach

The multiphase flow simulator moving particle semi-implicit (MPS) method is developed based on the fractional-flow approach, originated in the petroleum engineering literature, considering the fully three-phase flow with general boundary conditions. The fractional flow approach employs water saturation, total liquid saturation, and total pressure as primary variables. Most existing models based upon fractional flow are limited to two-phase flow and specific boundary conditions. Although there appear a number of three-phase flow models, they were mostly developed using pressure-based approaches, which require variable-switch techniques to deal with phase appearance and disappearance. The use of fractional flow-based approaches in MPS makes it unnecessary to use variable-switching to handle the change of phase configurations because the water saturation, total liquid saturation, and total pressure exist throughout the solution domain regardless of whether certain phases are present or not. Furthermore, most existing fractional flow-based models consider only specific boundary conditions, usually Dirichlet-type pressure for water phase and flux-type boundary for nonaqueous phase liquid or particular combinations for individual phase. However, the present model considers general boundary conditions of ten most possible and plausible cases. The first eight cases are the combinations of the phase pressure or the phase flux of each of the three individual phases. The other two cases are the variable boundary conditions: one for water-medium interface and the other for the air-medium interface when the directions of fluxes are not known a priori. Thus, the model’s capabilities of handling general boundary conditions extend the simulators’ usefulness in the field system.     

水位动态变化下水体驳岸渗流场及稳定性研究

目前水体驳岸设计只注重景观性与生态性,缺乏水位动态变化对驳岸渗流场及稳定性影响的研究。以澧水河2007-2012年水文地质资料为基础,从驳岸水位动态变化速度、岸坡安全系数两个角度,研究不同水位升降速度的驳岸稳定性变化规律。研究表明:水位动态变化对驳岸渗流场有明显影响,水位变化速度快,驳岸自由水面变化滞后越明显;驳岸水位上升时,其安全系数由大变小再到稳定;随着水位逐渐下降,驳岸安全系数先减小,当水位达到一定高度(此时安全系数最小)时,安全系数逐渐回升。研究结论用以指导水体驳岸设计,对水体驳岸在水位动态变化条件下安全运行具有现实意义。     

利用有机包裹体探讨塔里木盆地依奇克里克构造带下侏罗统油气运移与油气藏的存储条件

本文主要是对塔里木盆地库车坳陷东部的依奇克里克构造、克拉苏及吐格尔明三个构造带砂岩储层中包裹体进行了研究。利用包裹体的均一温度、有机包裹体的显微镜荧光特征、气相色谱、激光拉曼的成分分析,重点讨论了依奇克里克构造带成藏条件。通过对研究区储层中包裹体研究结合油气藏的分布特征认为在库车坳陷东部储层物性不是油气成藏的关键因素,封存条件才是控制该区油气能否成藏的决定因素。断层断在地层内部,对成藏有利,而地层产状较直立,断层直通到近地表,油气沿断层向上运移逃逸,包裹体是最好的证据。由储层包裹体温度、荧光特征和成分上的差异显示出在依奇克里克构造带至少有二次油气运移,第一期有机质低成熟时有过少量的油气运移,大量的油气运移是第二期。爆裂法进行群体包裹体气相色谱分析结果,显示出两个高峰区域,两个高峰区分别反应了轻烃与重烃“C”数上的差异,利用了514nm 和785nm两个激光波长,进行了激光拉曼分析,其测试结果显示出,断层开启的样品包裹体中 H_2O 都很高,CO_2含量很低,断层封闭的样品 H_2O 部很低,CO_2含量高,以次来判别储层的封闭条件。     

川西北马脑壳金矿床流体包裹体研究及对成矿条件的制约

川西北的马脑壳金矿床为一赋存于三叠系碎屑岩、碳酸盐岩系的大型微细浸梁型金矿床,通过详细野外地质调查和室内综合研究,可将矿床的热液成矿作用划分为四个阶段,应用均一法测得：Ｉ、ＩＩ、ＩＩＩ、ＩＶ阶段的成矿温度分别为１７５～２３０℃,１７９～２５５℃,１５７～２０７℃,１２０～２００℃,峰值分别为１９０～２１０℃,２００～２３０℃,１７０～１９０℃,１３０～１７０℃,盐度变化范围,Ｉ、ＩＩ、ＩＩＩ、ＩＶ     

P-T-X conditions of hydrothermal fluids and precipitation mechanism of stibnite-gold mineralization at the Wiluna lode-gold deposits,Western Australia: conventional and infrared microthermometric constraints

The Wiluna lode-gold deposits are located in the Archean Wiluna greenstone belt, in the northern sector of the Norseman-Wiluna belt in the Yilgarn Craton of Western Australia. They are hosted in subgreenschist facies meta-basalts, and controlled by the Wiluna strike-slip fault system and associated shear veins and breccias. The 13 individual lode-gold deposits have produced around 115 t Au from 1901 to 1946 and 1986 to today. Historically, they also produced 38.3 t As and 3.5 t Sb. Gold formed in two stages: stage 1 gold-pyrite-arsenopyrite is finely disseminated in the wallrock and breccia fragments, whereas stage 2 gold-stibnite is located in massive shear veins and breccia matrix, as fracture-fill and in banded-colloform textured veins. Stibnite-gold orebodies only occur in some of the deposits (e.g., Moonlight and northern part of the West Lode) and also display a restricted vertical extent, being preserved only in the uppermost 200 m of stibnite-bearing lodes.Petrographic, conventional, and infrared microthermometric and laser-Raman analysis on stibnite-bearing quartz veins and breccias reveal that the antimony- and gold-rich hydrothermal fluid was of mixed H₂O-NaCl-CO₂±CH₄ type. Microthermometric measurements reveal maximum homogenization temperatures of 340 °C (average 290±25 °C), and a wide range of salinities between 0.2 and 23 eq. wt% NaCl. Aqueous-carbonic fluid inclusions contain variable X_CO2+CH4 (0.03 to 0.82), with the carbonic phase containing a maximum X_CH4 of 0.21.Combined petrographic and microthermometric evidence suggests that the fluid inclusion properties reflect fluid immiscibility of a low-salinity, medium X_CO2+CH4, homogeneous parent fluid at about 290 °C and pressures between 700 and 1,700 bar. Fluid immiscibility was triggered by cyclic pressure release during fault-zone movement. The decompression (adiabatic cooling) of the hydrothermal fluids shifted the ore fluid to lower temperatures, significantly reduced the degree of stibnite undersaturation, and caused stibnite to precipitate. The deposition of stibnite reduced the ore-fluid H₂S concentration, thereby destabilized gold bisulfide complexes in solution, and caused gold precipitation locally. This mechanism explains the intimate spatial association of stibnite and gold in quartz veins and breccias in the stibnite-gold orebodies at Wiluna.Editorial handling: B. Lehmann     

Cross-bedding related anisotropy and its interplay with various boundary conditions in the formation and orientation of joints in an aeolian sandstone

Previous research revealed that the cross-bedding related anisotropy in Jurassic aeolian Aztec Sandstone cropping out in the Valley of Fire State Park, Nevada, affects the orientation of compaction bands, also known as anti-cracks or closing mode structures. We hypothesize that cross-bedding should have a similar influence on the orientation of the opening mode joints within the same rock at the same location. To test this hypothesis, we investigated the relationship between the orientation of cross-beds and the orientation of different categories of joint sets including cross-bed package confined joints and joint zones in the Aztec Sandstone. The field data show that the cross-bed package confined joints occur at high-angle to bedding and trend roughly parallel to the dip direction of the cross-beds. In comparison, the roughly N–S trending joint zones appear not to be influenced by the cross-beds in any significant way but frequently truncate against the dune boundaries.To characterize the anisotropy due to cross-bedding in the Aztec Sandstone, we measured the P-wave velocities parallel and perpendicular to bedding from 11 samples and determined an average P-wave anisotropy to be slightly larger than 13%. From these results, a model based on the generalized Hooke's law for anisotropic materials is used to analyze deformation of cross-bedded sandstone as a transversely isotropic material. In the analysis, the dip angle of cross-beds is assumed to be constant and the strike orientation varying from 0° to 359° in the east (x), north (y), and up (z) coordinate system. We find qualitative agreement between most of the model results and the observed field relations between cross-beds and the corresponding joint sets. The results also suggest that uniaxial extension (ε_zz > ε_xx = ε_yy = 0) and axisymmetric extension (ε_xx = ε_yy < ε_zz and ε_xx = ε_yy > ε_zz) would amplify the influence of cross-bedding associated anisotropy on the joint orientation whereas a triaxial extension (ε_xx > ε_yy > ε_zz) would mitigate this influence. We suggest that the potential implication of different categories of joint sets (i.e., cross-bed package confined joints and joint zones) forming in response to the variation of the boundary conditions (axisymmetric extension and triaxial extension, respectively) and the interplay with the rock anisotropy is significant. These results have important implications for fluid flow through aeolian sandstones in reservoirs and aquifers.     

陕西镇安石炭系--二叠系界线剖面层序地层和相对海平面变化研究

在陕西镇安西口上石炭统逍遥阶至下二叠统隆林阶的剖面上,区分出含生物碎屑灰泥灰岩、生物碎屑粒泥灰岩、生物碎屑泥粒灰岩、生物碎屑颗粒灰岩、团块颗粒灰岩和核形石颗粒灰岩6种主要的碳酸盐岩微相类型。通过详细的微相特征分析,确定了碳酸盐岩沉积过程中的水动力状况和相对水深变化趋势,识别出低能潮下型碳酸盐旋回A1、A2及高能浅滩型碳酸盐旋回B1、B2两类4种高频旋回。根据这些高频旋回在剖面上的有序叠置型式,将镇安西口地区的上石炭统逍遥阶至下二叠统隆林阶划分出12个四级层序和5个II型三级层序地层单元,反映了这一地区晚石炭世逍遥期至早二叠世隆林期曾发生过12次四级相对海平面波动和5次明显的三级相对海平面升降。     

陕西镇安石炭系一二叠系界线剖面层序地层和相对海平面变化研究

在陕西镇安西口上石炭统逍遥阶至下二叠统隆林阶的剖面上，区分出含生物碎屑灰泥灰岩、生物碎屑粒泥灰岩、生物碎屑泥粒灰岩、生物碎屑颗粒灰岩、团块颗粒灰岩和核形石颗粒灰岩6种主要的碳酸盐岩微相类型。通过详细的微相特征分析，确定了碳酸盐岩沉积过程中的水动力状况和相对水深变化趋势，识别出低能潮下型碳酸盐旋回A1、A2及高能浅滩型碳酸盐旋回B1、B2两类4种高频旋回。根据这些高频旋回在剖面上的有序叠置型式，将镇安西口地区的上石炭统逍遥阶至下二叠统隆林阶划分出12个四级层序和5个Ⅱ型三级层序地层单元，反映了这一地区晚石炭世逍遥期至早二叠世降林期曾发牛过12次四级相对海平面波动和5次明显的三级相对海平面升降。     

Temporal and spatial evolution of surface sediments characteristics in the Dagu River estuary and their dynamic response mechanism

Based on the 39 surface sediment samples collected in the flood season and the dry season in 2012 respectively and the measured hydrological data in October 2012, the sediment grain size characteristics has been analyzed and the response mechanism of surface sediments to estuarine hydrodynamics was revealed by calculating the range of waves and tidal currents. The results show that: (1) The grain size of the surface sediment samples decreased gradually from land to sea in the flood season. The fine sediment was redistributed under marine hydrodynamics in the dry season and the sediments showed coarser tendency ingeneral; (2) tidal current stirring sediment was very obvious in Dagu River estuary area, and wave stirring sediments mainly occurred in the tidal flat area and estuary sand bar area; (3) in the flood season, surface sediment sat the estuary were transported towards south and southeast. In the dry season, surface sediments were transported towards southwest at the north area of Jiaozhou Bay Bridge, and sediments were transported towards northeast area at the south of Jiaozhou Bay Bridge.     

On the mechanics of contact between a flexible transducer diaphragm located at a rigid boundary and an elastic material

In this paper we examine the mechanics of contact between a flexible elastic diaphragm and an isotropic elastic halfspace region where contact is established through the application of a normal stress. It is shown that a convenient solution for the maximum displacement and maximum flexural moment in the flexible diaphragm can be obtained by employing a variational technique. Copyright © 2006 John Wiley & Sons, Ltd.     

广西菜子岩剖面中上泥盆统界线附近沉积相与海平面变化*——兼论华南吉维期中期的海退及时间

广西菜子岩剖面处于台地边缘地带,含有底栖生物及浮游生物化石。综合岩相、生物相等资料,按牙形石带时限讨论了该地区吉维期—弗拉斯期早期主要的海平面变化历史。吉维期早期(相当于东岗岭组)经历了持续的海侵,大致下varcus带末期可能出现短暂的海退;“谷闭组”下部岩性复杂,海平面波动强烈,其中较大规模的包括了发生于上varcus带的海退(即吉维期中期海退)及下hermanni-cristatus带初期的海侵。吉维期晚期地层含有丰富的床板珊瑚以及其他礁体生物碎屑,产层孔虫、海百合等,意味着这个时期研究区可能处于礁前的上斜坡部位。“谷闭组”上部（牙形石falsiovalis带至transitans带下部）的岩性为中—薄层与厚层的灰岩互层,以正常的底栖生物为主,反映了由原来的礁前到滩前的变化,但水体深度在中—晚泥盆世界线前后可能总体相似。融县组(底部)以浅灰、灰白色的中—薄层至厚层—块状灰岩为特征,其沉积始于牙形石transitans带中部,证据表明此时海水可能变深。菜子岩剖面的白云岩位于上varcus带的下部,主要为细—中晶结构,见有灰岩残余以及正常海相生物碎屑。横向上与其他剖面对比表明,广西地区上varcus带可能出现了区域性的海退事件。     

安徽五河、河北南皮Propontocypris的发现——兼论2kaB.P.的突然变化事件

继在微山湖西江苏丰县、沛县及安徽萧县上更新亚统见广盐海相介形虫以来，首次在皖苏交界的安徽省五河县２ｋａＢ．Ｐ．以来的河漫滩相粘土层中见广盐海相介形虫ＰｒｏｐｏｎｔｏｃｙｐｒｉｓｅａｒｙｈａｌｉｎａＺｈａｏ，在河北省南皮县２ｋａＢ．Ｐ．西汉遗址文化层中亦见此种。在不能排除鸟类、风暴潮搬运海相介形虫卵的可能时，本文亦考虑其产出层位水化学状况应有较大变化。文内总结了目前对渤海湾两岸２ｋａＢ．Ｐ．“海侵”的研究状况，提出低洼滨海平原当时有地下水位上升或风暴潮加剧的可能，基本确认当时发生过突然变化事件，但不足以称     

A high‐resolution record of mire development and climatic change spanning the Late‐glacial–Holocene boundary at Church Moss,Davenham (Cheshire,England)

Excavations of deposits filling a closed basin within glacial drift at Church Moss, Davenham, near Northwich (Cheshire, England) revealed a sequence of Late‐glacial and Early Holocene sediments. Analyses of pollen and plant and invertebrate macrofossils were undertaken, together with loss‐on‐ignition analyses and a programme of AMS radiocarbon dating, to provide a record of changing biostratigraphy and climatic and ecological regimes. The infilling of features identified as frost‐cracks in the till flooring the basin gave remains that reflected conditions of extreme cold towards the end of the Devensian. The pollen record from a 3.5 m sequence of peat towards the deepest part of the basin, supported by radiocarbon dates, shows that organic deposition was initiated during the Late‐glacial Interstadial and continued into the early part of the Holocene. There was some evidence for a cool episode during the interstadial, with amelioration prior to the rapid onset of the tundra conditions of the Loch Lomond Stadial. Following the stadial, amelioration was rapid. There was evidence from both central and marginal sequences for a mosaic of fen dominated by sedges and often also mosses, with short‐lived small pools through much of the succession. Change to terrestrial conditions proceeded intermittently, probably as a result of subsidence caused by solution of underlying salt‐bearing strata. Copyright © 2000 John Wiley & Sons, Ltd.     

Columnar calcites as testimony of diagenetic overprinting at the boundary between Upper Tournaisian dolomites and limestones (Belgium): multiple origins for apparently similar features

In topographic flat areas, sedimentary settings may vary from one outcrop to another. In these settings, calcite precipitates may yield macroscopically similar columnar features, although they are products of different sedimentary or diagenetic processes. Three columnar calcite crystal fabrics, i.e. rosettes, palisade crusts and macro-columnar crystal fans, have been differentiated near and at the contact between Upper Tournaisian dolomites and limestones along the southern margin of the Brabant-Wales Palaeohigh. Their petrographic characteristics, and geochemical and fluid inclusion data provide information on the (dia)genetic processes involved. Rosettes composed of non-luminescent columnar calcite crystal fans (1–5 cm in diameter) developed on top of one another, forming discrete horizons in repetitive sedimentary cycles. The cycles consist of three horizons: (I) a basal horizon with fragments from the underlying horizon, (II) a micrite/microspar horizon with incipient glaebules, (III) an upper horizon consisting of calcite rosettes, with desiccation features. The petrographical features and δ¹⁸O signatures of −10·0 to −5·5‰ and δ¹³C values of −5·5 to −3·2‰ support either evaporative growth, an evaporative pedogenic origin, or overprinting of marine precipitates. Palisade crusts, composed of a few to 10 mm long non-luminescent calcite crystals, coat palaeokarst cavities. Successive palisade growth-stages occur which are separated by thin laminae of micrite or detrital quartz, displaying a geopetal arrangement. Palisade crusts are interpreted as intra-Mississippian speleothems. This interpretation is supported by their petrographic characteristics and isotopic signature (δ¹⁸O = −8·7 to −6·5‰ and δ¹³C = −4·8 to −2·5‰). Macro-columnar crystals, 1–50 cm long, developed mainly perpendicular to cavity walls and dolomite clasts. Crystal growth stages in the macro-columnar crystals are missing. δ¹⁸O values vary between −16·4 and −6·8‰ and δ¹³C values between −5·2 and −0·9‰. These features possibly support a late diagenetic high temperature precipitation in relation to hydrothermal karstification.     

下垫面蒸发和云中凝结分馏对降水稳定同位素影响的数值试验——时间变化的比较(以长沙降水同位素为例)

利用稳定同位素大气水平衡模式（iAWBM）的模拟数据,分析了在不同的下垫面蒸发和不同的凝结分馏条件下降水中δ¹⁸O的时间变化、降水量效应、负温度效应和大气水线。并通过与长沙站5年实测数据的比较以及模拟试验结果之间的相互比较,揭示下垫面蒸发水汽中稳定同位素的季节性变化和云中稳定同位素分馏对降水中稳定同位素变化的可能影响,增进对季风区水稳定同位素效应的理解和认识。iAWBM给出的4个模拟试验均很好地再现了监测站降水中δ¹⁸O的时间变化,模拟出季风区降水中稳定同位素在暖半年被贫化、在冷半年被富集的基本特点。与平衡分馏相比,动力分馏下降水中稳定同位素被贫化的程度加强、季节差和离散程度减小;由下垫面蒸发水汽中稳定同位素δ_e季节性变化所引起的降水中稳定同位素的变化在不同季节完全相反：在长沙,暖半年降水中δ¹⁸O更低,冷半年降水中δ¹⁸O更高,使得降水中稳定同位素季节差和离散程度增大。4个模拟试验均很好地再现了季风区的降水量效应和负温度效应。与平衡分馏相比,动力分馏下模拟的降水量效应和负温度效应的斜率相对较小;δ_e季节性变化导致模拟的降水量效应和负温度效应的斜率增大。利用iAWBM,模拟出季风区湿热气候条件下的MWL。动力分馏以及δ_e季节变化均使模拟得到的MWL的斜率和截距减小。