Characterization of reservoir properties and pore structure based on micro-resistivity imaging logging: Porosity spectrum, permeability spectrum, and equivalent capillary pressure curve

According to the capillary theory,an equivalent capillary model of micro-resistivity imaging logging was built.On this basis,the theoretical models of porosity spectrum(Ф_(i)),permeability spectrum(K_(i))and equivalent capillary pressure curve(pe)were established to reflect the reservoir heterogeneity.To promote the application of the theoretical models,the Archie's equation was introduced to establish a general model for quantitatively characterizing bi,K,and pei.Compared with the existing models,it is shown that:(1)the existing porosity spectrum model is the same as the general equation of gi;(2)the Ki model can display the permeability spectrum as compared with Purcell's permeability model;(3)the per model is constructed on a theoretical basis and avoids the limitations of existing models that are built only based on the component of porosity spectrum,as compared with the empirical model of capillary pressure curve.The application in the Permian Maokou Formation of Well TsX in the Central Sichuan paleo-uplift shows that the Ф_(i),K_(i),and p_(ci) models can be effectively applied to the identification of reservoir types,calculation of reservoir properties and pore structure parameters,and evaluation of reservoir heterogeneity.

Rock physics model for velocity–pressure relations and its application to shale pore pressure estimation

Acoustic wave velocity has been commonly utilized to predict subsurface geopressure using empirical relations.Acoustic wave velocity is, however, affected by many factors. To estimate pore pressure accurately, we here propose to use elastic rock physics models to understand and analyze quantitatively the various contributions from these different factors affecting wave velocity. We report a closed-form relationship between the frame flexibility factor(γ) in a rock physics model and differential pressure, which presents the major control of pressure on elastic properties such as bulk modulus and compressional wave velocity. For a gas-bearing shale with abundant micro-cracks and fractures, its bulk modulus is much lower at abnormally high pore pressure(high γ values) where thin cracks and flat pores are open than that at normal hydrostatic pressure(low γ values) where pores are more rounded on average. The developed relations between bulk modulus and differential pressure have been successfully applied to the Upper Ordovician Wufeng and Lower Silurian Longmaxi formations in the Dingshan area of the Sichuan Basin to map the three-dimensional spatial distribution of pore pressure in the shale, integrating core, log and seismic data. The estimated results agree well with field measurements. Pressure coefficient is positively correlated to gas content. The relations and methods reported here could be useful for hydrocarbon exploration, production, and drilling safety in both unconventional and conventional fields.

Influential factors of loess liquefaction and pore pressure development

The liquefaction of loess under dynamic loading is studied experimentally with a dynamic triaxial test system. The effects of over-consolidation ratio (OCR), saturation degree and the frequency of dynamic loading upon loess liquefaction are investigated. The development of pore pressure within loess samples is also discussed. Based on the experimental results, the empirical relationship between pore pressure ratio and loading cycle number ratio is established for normal consolidated saturated loess.

Influence of pore water pressure on upper bound analysis of collapse shape for square tunnel in Hoek-Brown media

To investigate the effective shape of collapsing block in square tunnel subjected to pore water pressure,the analytical solution of detaching curve was derived using upper bound theorem of limit analysis with Hoek-Brown failure criterion. The work rate of pore water pressure,which was regarded as an external rate of work,was taken into account in the framework of limit analysis. Taking advantages of variational calculation,the objective function with respect to detaching curve was optimized to obtain the effective shape of collapsing block for square tunnel. According to the numerical results,it is found that the varying pore water pressure coefficient only affects the height and width of the collapsing block,whereas the shape of collapsing block remains unchanged.

Pre-Drilling Prediction Techniques on the High-Temperature High-Pressure Hydrocarbon Reservoirs Offshore Hainan Island,China

Decreasing the risks and geohazards associated with drilling engineering in high-temperature high-pressure(HTHP) geologic settings begins with the implementation of pre-drilling prediction techniques(PPTs). To improve the accuracy of geopressure prediction in HTHP hydrocarbon reservoirs offshore Hainan Island, we made a comprehensive summary of current PPTs to identify existing problems and challenges by analyzing the global distribution of HTHP hydrocarbon reservoirs, the research status of PPTs, and the geologic setting and its HTHP formation mechanism. Our research results indicate that the HTHP formation mechanism in the study area is caused by multiple factors, including rapid loading, diapir intrusions, hydrocarbon generation, and the thermal expansion of pore fluids. Due to this multi-factor interaction, a cloud of HTHP hydrocarbon reservoirs has developed in the Ying-Qiong Basin, but only traditional PPTs have been implemented, based on the assumption of conditions that do not conform to the actual geologic environment, e.g., Bellotti's law and Eaton's law. In this paper, we focus on these issues, identify some challenges and solutions, and call for further PPT research to address the drawbacks of previous works and meet the challenges associated with the deepwater technology gap. In this way, we hope to contribute to the improved accuracy of geopressure prediction prior to drilling and provide support for future HTHP drilling offshore Hainan Island.

Nanoscale Pore Characteristics of the Upper Permian Mudrocks from a Transitional Environment in and around Eastern Sichuan Basin, China

Nanoscale pore characteristics of the Upper Permian Longtan transitional mudrocks and their equivalent strata Wujiaping Formation marine mudrocks in and around the eastern Sichuan Basin was investigated using field emission scanning electron microscopy(FE-SEM)and low-pressure N2 adsorption experiments.The results indicate that the Upper Permian mudrock is at a mature stage with total organic carbon(TOC)values ranging between 0.47%and 12.3%.The Longtan mudrocks mainly contain vitrinite,and their mineral composition is primarily clay.In contrast,the Wujiaping mudrocks are dominated by sapropelinite and solid bitumen,and their mineral compositions are mainly quartz and a notably high amount of pyrite.The FE-SEM reveals that clay mineral pores and microcracks are the common pore types in the Longtan mudrocks.The specific surface area and pore volume depend on the clay content but are negatively correlated with the TOC.The generation of nanometer pores in the Longtan mudrocks is caused by high clay mineral contents.Meanwhile,the Wujiaping mudrock mainly contains OM pores,and the pore parameters are positively correlated with the TOC.The OM pore development exhibits remarkable differences in the Longtan and Wujiaping mudrocks,which might be related to their sedimentary facies and maceral fractions.Vitrinite and inertinite appear as discrete particles in these mudrocks and cannot generate pores during thermal maturation.Sapropelinite often contains many secondary pores,and solid bitumen with large particles,usually with several pores,is not the major contributor to the pore system of the investigated mudrock.

Classification of microscopic pore-throats and the grading evaluation on shale oil reservoirs

On the basis of the characterization of microscopic pore-throats in shale oil reservoirs by high-pressure mercury intrusion technique, a grading evaluation standard of shale oil reservoirs and a lower limit for reservoir formation were established. Simultaneously, a new method for the classification of shale oil flow units based on logging data was established. A new classification scheme for shale oil reservoirs was proposed according to the inflection points and fractal features of mercury injection curves: microscopic pore-throats(less than 25 nm), small pore-throats(25-100 nm), medium pore-throats(100-1 000 nm) and big pore-throats(greater than 1 000 nm). Correspondingly, the shale reservoirs are divided into four classes, I, II, III and IV according to the number of microscopic pores they contain, and the average pore-throat radii corresponding to the dividing points are 150 nm, 70 nm and 10 nm respectively. By using the correlation between permeability and pore-throat radius, the permeability thresholds for the reservoir classification are determined at 1.00× 10^(-3) μm^2, 0.40×10^(-3) μm^2 and 0.05×10^(-3) μm^2 respectively. By using the exponential relationship between porosity and permeability of the same hydrodynamic flow unit, a new method was set up to evaluate the reservoir flow belt index and to identify shale oil flow units with logging data. The application in the Dongying sag shows that the standard proposed is suitable for grading evaluation of shale oil reservoirs.

New Pore Pressure Evaluation Techniques for LAGIA-8 Well, Sinai, Egypt as a Case Study

Drilling into a geopressured zone will generally cause a change in a number of basic formation/ drilling relationships. This change is usually seen as a reversal of a gradual depth related trend in a lithologically uniform formation. Of all the geophysical methods, the reflection seismic method is essentially the only technique used to predict pore pressures. The seismic method detects changes of interval velocity with depth from velocity analysis of the seismic data. These changes are in turn related to lithology, pore fluid type, rock fracturing and pressure changes within a stratigraphic column. When the factors affecting the velocity are understood for a given area, a successful pressure prediction can be made. For clastic environments such as the Tertiary section of the Gulf of Mexico or the Niger delta, the interval velocity of the rocks increases with depth because of compaction. In these areas, deviations from normal compaction trends are related to abnormally high pore pressures. The adapted methods provide a much easier way to handle normal compaction trend lines. In addition to well log methods, pressure detection can be obtained via drilling parameters by applying Eaton’s DXC methods. Seismic velocities have long been used to estimate pore pressure, indeed both these quantities are influenced by variations in rock properties such as porosity, density, effective stress and so on, and high pore pressure zones are often associated with low seismic velocities. Pressure prediction from seismic data is based on fundamentals of rock physics and seismic attribute analysis. This paper hence tries to assess the use of seismic waves as a viable means to calculate pore pressure, especially in areas where no prior drilling history can be found. Then we applied these methods on LAGIA-8 well, Sinai, Egypt as a case study. Pore pressure prediction from Seismic is a very essential tool to predict pore pressure before drilling operation. This could prevent the well problem as well blowout and to prevent formation damage, especially in areas where no prior drilling history can be found.

New Pore Pressure Evaluation Techniques for LAGIA-8 Well, Sinai, Egypt as a Case Study

Drilling into a geopressured zone will generally cause a change in a number of basic formation/ drilling relationships. This change is usually seen as a reversal of a gradual depth related trend in a lithologically uniform formation. Of all the geophysical methods, the reflection seismic method is essentially the only technique used to predict pore pressures. The seismic method detects changes of interval velocity with depth from velocity analysis of the seismic data. These changes are in turn related to lithology, pore fluid type, rock fracturing and pressure changes within a stratigraphic column. When the factors affecting the velocity are understood for a given area, a successful pressure prediction can be made. For clastic environments such as the Tertiary section of the Gulf of Mexico or the Niger delta, the interval velocity of the rocks increases with depth because of compaction. In these areas, deviations from normal compaction trends are related to abnormally high pore pressures. The adapted methods provide a much easier way to handle normal compaction trend lines. In addition to well log methods, pressure detection can be obtained via drilling parameters by applying Eaton’s DXC methods. Seismic velocities have long been used to estimate pore pressure, indeed both these quantities are influenced by variations in rock properties such as porosity, density, effective stress and so on, and high pore pressure zones are often associated with low seismic velocities. Pressure prediction from seismic data is based on fundamentals of rock physics and seismic attribute analysis. This paper hence tries to assess the use of seismic waves as a viable means to calculate pore pressure, especially in areas where no prior drilling history can be found. Then we applied these methods on LAGIA-8 well, Sinai, Egypt as a case study. Pore pressure prediction from Seismic is a very essential tool to predict pore pressure before drilling operation. This could prevent the well problem as well blowout and to prevent formation damage, especially in areas where no prior drilling history can be found.

Development technology of rigidity-drain pile and numerical analysis of its anti-liquefaction characteristics

Pile foundation is widely used in the offshore engineering. The pile can be seriously destroyed by the soil liquefaction during strong earthquakes. The potentials of liquefaction and damages of pile foundation due to the liquefaction can be reduced by the implementation of the drainage in the liquefiable foundation. A patented pile technology, named rigidity-drain pile, was introduced. The partial section of the pile body was filled by materials with higher penetrability which forms some effective drainage channels in the pile. The principles and construction methods were presented. 3D models for both rigidity-drain pile and ordinary pile were built in FLAC3D code. The dynamic loadings were applied on the bottom of the model. According to the numerical results, in the case of the rigidity-drain pile, the water in the relevant distance range around the pile flows toward the pile drainage, the contour of the pore pressure shows a funnel form. Contrast to the ordinary pile, the rigidity-drain pile can dissipate the accumulated excess pore water, maintain effective stress and obviously reduce the possibility of surrounding soil liquefaction.

Upper bound solutions of stability factor of shallow tunnels in saturated soil based on strength reduction technique

Based on the upper bound theorem of limit analysis,the factor of safety for shallow tunnel in saturated soil is calculated in conjunction with the strength reduction technique.To analyze the influence of the pore pressure on the factor of safety for shallow tunnel,the power of pore pressure is regarded as a power of external force in the energy calculation.Using the rigid multiple-block failure mechanism,the objective function for the factor of safety is constructed and the optimal solutions are derived by employing the sequential quadratic programming.According to the results of optimization calculation,the factor of safety of shallow tunnel for different pore pressure coefficients and variational groundwater tables are obtained.The parameter analysis shows that the pore pressure coefficient and the location of the groundwater table have significant influence on the factor of safety for shallow tunnel.

Numerical analysis of landslide-generated debris flow on July 3, 2021 in Izu Mountain area, Shizuoka County, Japan

With the increasing of extreme rainfall frequency, landslides accompanied by mudslides often lead to serious casualties and property damage. On 3rd July 2021, a debris flow occurred in Izu Mountain area, Shizuoka County, Japan, and then resulted in 26 deaths and 131 houses destroyed, where houses were mainly built along the banks of the creek. In order to analyse the landslide state and distribution, a two-dimensional debris flow dynamic model(Massflow) was selected to simulate the process of the landslide-generated debris flow. When the model results are considered together with remote sensing images, the volume distribution of the unstable landslide is also able to be determined. The results show that(1) the affected areas are mainly concentrated at the outfall of the gully and on both sides of the streets.(2) The pore pressure ratio is an important factor affecting the damage range of this debris flow.(3) The increase of the pore pressure ratio in the landslide make the movement distance of debris flow increase significantly.

Dynamic response analysis of liquefiable ground due to sinusoidal waves of different frequencies of shield construction

Vibration induced by shield construction can lead to liquefaction of saturated sand.Based on FLAC3D software,a numerical model of tunnel excavation is established and sinusoidal velocity loads with different frequencies are applied to the excavation face.The pattern of the excess pore pressure ratio with frequency,as well as the dynamic response of soil mass under different frequency loads before excavation,is analyzed.When the velocity sinusoidal wave acts on the excavation surface of the shield tunnel with a single sand layer,soil liquefaction occurs.However,the ranges and locations of soil liquefaction are different at different frequencies,which proves that the vibration frequency influences the liquefaction location of the stratum.For sand-clay composite strata with liquefiable layers,the influence of frequency on the liquefaction range is different from that of a single stratum.In the frequency range of 5-30 Hz,the liquefaction area and surface subsidence decrease with an increase in vibration frequency.The research results in this study can be used as a reference in engineering practice for tunneling liquefiable strata with a shield tunneling machine.

国产一次性电子输尿管软镜联合负压吸引鞘治疗单侧重复肾重复输尿管畸形合并上尿路结石的疗效分析

目的:探讨国产一次性电子输尿管软镜联合负压吸引鞘治疗单侧重复肾重复输尿管畸形合并上尿路结石的有效性及安全性。方法:回顾性分析2018年5月至2023年9月武汉大学人民医院泌尿外科收治的17例单侧重复肾重复输尿管畸形合并上尿路结石的患者资料,该17例患者均采用输尿管软镜联合负压吸引鞘治疗上尿路结石。分析其结石分布、结石大小、手术时间、术中术后并发症、术后住院时间、结石清除率等,并随访患者的术后情况。结果:所有17例患者手术顺利,手术时间27~65 min,平均(43.4±6.7)min,术中术后均未出现输尿管撕裂、穿孔、感染性休克等严重并发症。术后1~3个月复查腹部平片(kidney ureter bladder,KUB)或泌尿系计算机断层扫描(computed tomography,CT),15例患者结石清除,结石清除率15/17(88.24%)。二期手术患者2例,二期手术结石清除率100%。结论:国产一次性电子输尿管软镜联合负压吸引鞘治疗单侧重复肾重复输尿管畸形合并上尿路结石,结石清除率高,手术并发症少,具有良好的有效性及安全性。

渤海某区块浅层气对工程影响的评估

以渤海某区块项目为例,采用保压取样方法获取待评估区域海底土层样品并对所含气体开展组分测定,定性判断浅层气成因。采用静力触探孔压消散试验方法获取地层中超孔隙水压力,定量评估对含浅层气地层承载力的折减效应。结果表明,内浅层气主要组分为甲烷和二氧化碳,气体含量较低。内地层中的浅层气可能属于生物成因气。含浅层气地层的超孔隙水压力相对上覆有效压力比值约3%,承载力折减因数为1.0%~1.1%。所采取的方法能够有效评估浅层气对海上基础设施和生产安全的影响程度,具有推广价值。

基于合理地层压力的X油田吸水能力评价

地层压力是影响注水井的注入压力和吸水能力的重要因素,保持合理地层压力是实现油田稳产的基础。针对X油田A、B两个区块地层压力差异导致的油井产液降低的问题,采用生产指示曲线法结合油水井生产动态对X油田的地层压力进行了评价,确定了X油田的合理地层压力,构建了不同注水量下吸水指数与地层压力的关系图版,给出了地层压力恢复方案。研究结果表明,目前A区块处于合理地层压力范围内,B区块地层压力偏低,应补充注入井注水量;同时两区块实际注水量均偏低,应进一步提高注水强度。基于合理地层压力的油田吸水能力评价方法对油田水驱开发具有重要指导作用。

考虑土体非均质性和孔隙水压力耦合作用的地下连续墙槽壁稳定性研究

针对非均质和富水地层中地下连续墙成槽施工诱发槽壁坍塌失稳的发生机理和破坏模式研究尚不透彻这一现状,利用空间离散技术构建了地下连续墙槽壁极限状态下的三维离散型破坏机制。基于该破坏机制和极限分析上限定理,通过将土体非均质特性和孔隙水压力引入内外能耗功率计算获得了极限状态下槽壁的安全系数目标函数。利用非线性规划算法对该目标函数进行优化计算,得到了地下连续墙槽壁安全系数的最优上限解。为了验证该方法的有效性,将计算得到的槽壁安全系数分别与已有计算结果以及数值模拟结果进行了对比分析,对比结果表明:该方法计算的安全系数与已有计算结果和数值模拟结果基本一致。在此基础上,考虑土体非均质特性和孔隙水压力的耦合作用对地下连续墙槽壁的安全系数进行了参数分析,分析结果表明:孔隙水压力和土体非均质性对地下连续墙槽壁安全系数有较大影响,地下连续墙槽壁安全系数随地下水位高度Hw的升高和黏聚力非均质系数kc的增大而减小,随重度非均质系数kg的增大而增大。

基于多对数螺旋曲线破坏机构的富水软岩隧道开挖面支护应力分析

为了避免对胡克-布朗准则进行任何线性化处理,将广义胡克-布朗(generalized Hoek-Brown, GHB)强度准则直接引入到开挖面稳定性分析中,基于极限分析理论构建了富水地层浅埋隧道开挖面多对数螺旋线破坏机构。在考虑孔隙水压力的作用下,推导出能量平衡方程,得到开挖面极限支护应力表达式,并且通过程序优化得到支护应力最优上限解。在此基础上,分析了GHB强度参数、孔隙水压力系数、岩土体重度以及隧道埋深比对开挖面稳定性的影响。最后通过归一化处理得到不同情况下的开挖面极限支护应力稳定性图表,通过该图表能够快速便捷地获得富水地层隧道开挖面维持稳定性所需的极限支护应力,为盾构隧道初步设计提供理论依据。

风机荷载作用下频率比对砂土动力响应的影响

海上风机在运行过程中风荷载与波浪荷载频率耦合作用对土体动力特性的影响不容忽略,为了探究这种影响,采用多向循环动单剪系统(VDDCSS)对砂土进行了一系列试验研究,分析了双向剪切频率比f_(r)对砂土动应变、动孔压比以及动强度的影响。试验结果表明:剪应变发展受双向频率耦合及应力组合效应的影响。循环应力比(C_(SR))为0.15时,双向频率耦合效应大于应力组合作用;随着应力水平的增加(C_(SR)=0.20或0.25),频率耦合效应和应力组合作用均得到了加强。在低应力水平时,破坏试样的孔压曲线呈“上凹型”且频率比对试样孔压发展影响不大。而随着应力水平增加(C_(SR)=0.20或0.25),孔压最终值u_(max)随之增大。除f_(r)=1.00与10.00对应试样外,土体动强度随f_(r)与C_(SR)值增大而减少。f_(r)=5.00对应试样应变大、孔压值高、强度低的原因与此比值下土体受双向频率耦合与应力组合的叠加效果增强有关。

竖向荷载对高承台群桩基础侧向位移影响规律的振动台试验

基于大型振动台模型试验,以高承台群桩为研究对象,模拟了地震作用下高桩承台基础的振动响应,通过超静孔压比(简称孔压比)、加速度和地基土位移研究饱和砂土液化场地竖向荷载对桩基侧向位移的影响规律。结果表明:近桩区域受桩基竖向荷载影响最为显著,孔压比与竖向荷载呈负相关关系,桩基两侧孔压比差值与地基土浅层加速度是造成桩基侧向位移的主要原因;地基土对台面输出加速度主要起放大作用,且随着竖向荷载的增加,加速度放大系数显著减小;在振动输入与竖向荷载相同的条件下,地基土在地震动作用下引起的侧向流动及显著变形极大影响了桩基的侧向位移;竖向荷载为9%桩基极限承载力条件下,桩基倾斜角为4.80°,侧向位移较小,结合地基土抗液化能力,桩基抗倾覆能力最优。