声电效应测井的有限差分模拟

本文研究声电效应测井波场的有限差分模拟算法.忽略井外地层中诱导电磁场对孔隙弹性波的影响,将求解动电耦合波方程组的问题解耦,先计算孔隙弹性波,再计算其诱导电磁场.基于轴对称柱坐标系下的速度-应力交错网格,采用时域有限差分计算井孔流体声波和井外地层孔隙弹性波.将电磁场近似看作似稳场,基于轴对称柱坐标系下的5点式有限差分网格,求解不同时刻的电位Poisson方程,计算诱导电场.结果表明:本文算法可准确模拟频率6.0 kHz的声电效应测井全波;在声波测井频率范围内,电导率、动电耦合系数和动态渗透率的低频近似对伴随电磁场的计算影响不大;地层水平界面导致伴随反射斯通利波的电场和显著的界面电磁波,后者对于探测地层界面具有潜在的应用价值.     

孔隙岩样动电特性的实验研究

本文针对流体饱和孔隙介质的动电效应,在实验室内进行了小岩样的动电实验(流动电势实验)测量,获得了10块岩样在6种不同饱和浓度下的流动电势系数,进而分析了饱和溶液浓度、岩样渗透率和泥质含量对流动电势系数的影响,探讨了流动电势系数对上述参数的敏感性,并进一步给出了流动电势系数随这些参数的变化规律.实验结果表明:前人流动电势系数频响曲线中的凹点现象是由实验装置的共振引起,并非岩样动电效应的自身特性.实验还发现:孔隙介质的动电耦合能力与溶液浓度有关,流动电势系数随溶液浓度增大而减小,但在很高浓度溶液中,动电现象依然存在,流动电势系数趋于常数.此外,流动电势系数对地层参数(渗透率和泥质含量)的敏感性受溶液浓度影响较大,随着溶液浓度增大,敏感性降低.因此,可在低浓度饱和情况下(低于0.4mol/L),利用流动电势系数的幅度对地层渗透率进行直接评价.本文结果对动电测井的可行性及应用条件给出了实验说明.     

动电测井实验研究Ⅰ:渗透率的评价

本文针对流体饱和孔隙介质的动电效应,在砂岩模型井中开展了动电测井实验研究.记录到渗透率不同模型中单极源动电测井的全波波形,清晰地观测到了伴随纵波、横波和斯通利波的动电转换信号,给出了本文实验条件下,动电测井全波中各分波的幅度,并通过电声比的大小说明了各分波的动电转换能力.实验结果表明:孔隙介质动电效应与渗透率密切相关,动电信号的幅度随着地层渗透率的增大而增加,在高渗透率地层中记录动电信号的幅度较大,这一特性可用于地层渗透率的井下动电评估.本文实验结果还验证了动电测井的可行性,同时指出:与声波测井相比,动电测井信号对地层渗透率更加敏感,这为渗透率等地层参数的井下动电测量奠定了实验基础.     

基于动电效应的岩芯渗透率实验测量

根据孔隙介质的动电耦合理论设计了一套岩芯渗透率测量系统.实验采用交流锁相放大技术,在低频12~42 Hz范围内完成了砂岩岩样流动电势和电渗实验,得到了流动电势系数K_S和电渗压力系数K_E,进而计算出岩样动电渗透率,对于中、高渗透率岩样,测量得出的动电渗透率与常规气测渗透率差异较小,两者具有很好的相关性.实验表明,动电测量可作为岩样渗透率测量的一种方法,同时揭示了利用地层动电测井信号反演地层参数的可能性,实验结果对于分析天然地震动电效应也有参考意义.     

动电测井实验研究Ⅱ:伴随动电场和界面动电场

动电效应指孔隙介质中与固/液界面的双电层和孔隙流体渗流有关的弹性-电磁耦合现象,探索基于动电效应的勘探和测井新方法则是石油工业重点关注的研究方向之一.本文对流体饱和孔隙介质中的动电效应进行了实验测量研究,记录到了不同模型井中伴随声波的动电转换信号和界面动电转换信号,对比分析了这两种动电信号的产生条件及传播特性,验证了理论分析结果,并进一步研究了动电信号的分波成分及其与声波信号的关系,探讨了声源激发模式、电极接收方式及数据处理方法对动电信号的影响,为动电测井仪器设计奠定实验基础.     

声电效应测井电声比及其与地层渗透率的关系

基于流体饱和孔隙介质中声波-电磁场耦合效应的测井方法具有潜在的应用价值.本文从Pride动电耦合波方程组入手,推导了伴随斯通利波的井孔电场与声压比值(电声比)的频率域表达式.结果表明,在较低频率条件下,电声比幅角的正切值与渗透率呈反比.在此基础上,提出了利用低频的声电效应测井全波反演地层渗透率的方法.针对砂岩地层,从计...     

储层岩石流动电位频散特性的数学模拟

利用储层岩石流动电位的频散特性评价复杂储层已经成为勘探地球物理领域关注的热点,但是目前还没有形成基于储层岩石储渗特性及电化学性质的具有普遍指导意义的理论方法和数学模型.本文利用微观毛管理论,通过随时间谐变条件下渗流场和电流场的耦合模型,建立了描述储层岩石流动电位频散特性的数学方法,定量分析了频率域储层岩石动态渗透率、动电耦合系数和流动电位耦合系数随储层岩石孔隙度、溶液浓度和阳离子交换量的变化规律.研究结果表明:储层岩石流动电位频散特性是储层流体惯性力与流体黏滞力相互作用的结果.储层岩石孔隙度越大,储层维持流体原有运动状态的能力越大,临界频率越小;储层岩石的溶液浓度和阳离子交换量对临界频率没有影响.储层岩石的孔隙度越大,流体流动能力越强,流动电位各耦合系数的数值越大;溶液浓度越小或阳离子交换量越大,孔隙固液界面的双电层作用越强,各耦合系数的数值越大.     

多极源随钻动电测井波场特性实验研究

探索基于孔隙介质动电效应的随钻测井新方法是当前测井勘探领域的热点研究之一.本文针对随钻测井情况下的动(声)-电转换机理,在砂岩模型井中进行了多极源(单/偶/四极源)随钻动电测井实验研究,记录了井中声波激发动电转换波场的时域波形,给出了各分波的时间-速度图,并与多极源随钻声波测井数据进行对比,结果表明:钻铤波诱导的动电转换信号在随钻动电测井全波中的相对幅度显著降低,从而可准确提取到地层声波的传播速度,实现了利用随钻动电测井方法测量地层波速的目的,验证了随钻动电测井技术可有效抑制钻铤波的思想.本文还进一步分析了钻铤波的动电耦合特性,并结合有机玻璃模型井中随钻动电测井和无井孔模型下的随钻声波测量,研究了伴随钻铤波电场信号的传播特性和衰减规律,深入探讨了随钻动电测井全波中钻铤波被削弱的原因,并给出了定量的实验数据分析,为随钻动电测井理论模拟的解释提供有效实验数据支撑.     

孔隙介质声学理论中的动态渗透率

解释了Johnson动态渗透率k(ω)的概念，指出k(ω)包含了Biot理论中质量耦合系数和粘滞修正系数的意义，对k(ω)与Biot理论中相应表达式k_B(ω)进行了比较，当结构因子β≈64α_∞/m时k(ω)与k_B(ω)的差异很小.计算表明，当声源频率接近临界频率时声测井波形中导波的衰减随质量耦合系数增大而减小.     

油水饱和泥质砂岩流动电位的频散特性

油水饱和泥质砂岩中流动电位的研究对于揭示含油储层震电勘探和动电测井的机理有着重要的意义.本文首先从岩石孔隙的微观结构出发,构造了描述水润湿条件下油水饱和泥质砂岩储层的毛管模型.在模型中依据油水流动遵守的Navier-Stokes方程和电化学传质动力学理论,建立了描述油水饱和泥质砂岩流动电位的数学方程,并数学模拟了岩石储渗参数对流动电位频散特性的影响规律.研究结果表明:储层孔隙内流体受到的粘滞力与惯性力控制着水相和油相的流动,从而决定了流动电位的频散特性.随着孔隙度的增大,油水两相各自的有效渗透率均增大;而含水饱和度的升高使得水相有效渗透率增大,油相有效渗透率减小.在水润湿条件下,流动电位耦合系数随含水饱和度升高而增大,随束缚水饱和度的升高而减小.另外,流动电位相对耦合系数也随含水饱和度的升高而增大,但无频散现象.     

An equivalent medium model for wave simulation in fractured porous rocks

Seismic wave propagation in reservoir rocks is often strongly affected by fractures and micropores. Elastic properties of fractured reservoirs are studied using a fractured porous rock model, in which fractures are considered to be embedded in a homogeneous porous background. The paper presents an equivalent media model for fractured porous rocks. Fractures are described in a stress‐strain relationship in terms of fracture‐induced anisotropy. The equations of poroelasticity are used to describe the background porous matrix and the contents of the fractures are inserted into a matrix. Based on the fractured equivalent‐medium theory and Biot's equations of poroelasticity, two sets of porosity are considered in a constitutive equation. The porous matrix permeability and fracture permeability are analysed by using the continuum media seepage theory in equations of motion. We then design a fractured porous equivalent medium and derive the modified effective constants for low‐frequency elastic constants due to the presence of fractures. The expressions of elastic constants are concise and are directly related to the properties of the main porous matrix, the inserted fractures and the pore fluid. The phase velocity and attenuation of the fractured porous equivalent media are investigated based on this model. Numerical simulations are performed. We show that the fractures and pores strongly influence wave propagation, induce anisotropy and cause poroelastic behaviour in the wavefields. We observe that the presence of fractures gives rise to changes in phase velocity and attenuation, especially for the slow P‐wave in the direction parallel to the fracture plane.     

Developing a new form of permeability and Kozeny–Carman constant for homogeneous porous media by means of fractal geometry

The semi-empirical Kozeny–Carman (KC) equation is the most famous permeability–porosity relation, which is widely used in the field of flow in porous media and is the starting point for many other permeability models. However, this relation has many limitations from its inception, and the KC constant is an empirical parameter which was proved to be not a constant. In this paper, we briefly reviewed the KC equation, its modifications and various models for the KC constant. We then derived an analytical expression for the permeability in homogeneous porous media based on the fractal characters of porous media and capillary model. The proposed model is expressed as a function of fractal dimensions, porosity and maximum pore size. The analytical KC constant with no empirical constant is obtained from the assumption of square geometrical model. Furthermore, a distinct linear scaling law between the dimensionless permeability and porosity is found. It is also shown that our analytical permeability is more closely related to the microstructures (fractal dimensions, porosity and maximum pore size), compared to those obtained from conventional methods and models.     

The pore structure of tight limestone—Jurassic Ziliujing Formation,Central Sichuan Basin,China

The pore structure of the tight limestone in the Daanzhai Member of the Ziliujing Formation, Jurassic System, in central Sichuan Basin, China, is complex but essential to the exploration and development of tight oil. The pore structure of the tight limestone is studied by using scanning electron microscopy (SEM), nitrogen adsorption, high-pressure mercury intrusion, and nuclear magnetic resonance (NMR). The experimental results suggest that the pores are mainly slit pores and mesopores and macropores contribute to the pore volume and specific surface. The displacement pressure, average pore size, and homogeneity coefficient correlate with porosity and permeability and can be used to evaluate the pore structure. The full pore-size distribution was obtained by combining nitrogen adsorption and high-pressure mercury intrusion. We find that the limestone mainly contains mesopores with diameter of 2?50 nm. The T₂ distribution was converted into pore-size distribution, well matching the full pore-size distribution. The relation between T₂ and pore size obeys a power law and the geometric mean of T₂ correlates with the pore structure and can be used in the pore structure evaluation.     

A transport phase diagram for pore-level correlated porous media

Transport in porous media is often characterized by the advection–dispersion equation, with the dispersion coefficient as the most important parameter that links the hydrodynamics to the transport processes. Morphological properties of any porous medium, such as pore size distribution, network topology, and correlation length control transport. In this study we explore the impact of correlation length on transport regime using pore-network modelling. Earlier direct simulation studies of dispersion in carbonate and sandstone rocks showed larger dispersion compared to granular homogenous sandpacks. However, in these studies, isolation of the impact of correlation length on transport regime was not possible due to the fundamentally different pore morphologies and pore-size distributions. Against this limitation, we simulate advection–dispersion transport for a wide range of Péclet numbers in unstructured irregular networks with “different” correlation lengths but “identical” pore size distributions and pore morphologies. Our simulation results show an increase in the magnitudes of the estimated dispersion coefficients in correlated networks compared to uncorrelated ones in the advection-controlled regime. The range of the Péclet numbers which dictate mixed advection–diffusion regime considerably reduces in the correlated networks. The findings emphasize the critical role of correlation length which is depicted in a conceptual transport phase diagram and the importance of accounting for the micro-scale correlation lengths into predictive stochastic pore-scale modelling.     

Effect of mineral reactions on the hydraulic properties of unsaturated soils: Model development and application

The selective radius shift model was used to relate changes in mineral volume due to precipitation/dissolution reactions to changes in hydraulic properties affecting flow in porous media. The model accounts for (i) precipitation/dissolution taking place only in the water-filled part of the pore space and further that (ii) the amount of mineral precipitation/dissolution within a pore depends on the local pore volume. The pore bundle concept was used to connect pore-scale changes to macroscopic soil hydraulic properties. Precipitation/dissolution induces changes in the pore radii of water-filled pores and, consequently, in the effective porosity. In a time step of the numerical model, mineral reactions lead to a discontinuous pore-size distribution because only the water-filled pores are affected. The pore-size distribution is converted back to a soil moisture characteristic function to which a new water retention curve is fitted under physically plausible constraints. The model equations were derived for the commonly used van Genuchten/Mualem hydraulic properties. Together with a mixed-form solution of Richards’ equation for aqueous phase flow, the model was implemented into the geochemical modelling framework PHREEQC, thereby making available PHREEQC’s comprehensive geochemical reactions. Example applications include kinetic halite dissolution and calcite precipitation as a consequence of cation exchange. These applications showed marked changes in the soil’s hydraulic properties due to mineral precipitation/dissolution and the dependency of these changes on water contents. The simulations also revealed the strong influence of the degree of saturation on the development of the saturated hydraulic conductivity through its quadratic dependency on the van Genuchten parameter α. Furthermore, it was shown that the unsaturated hydraulic conductivity at fixed reduced water content can even increase during precipitation due to changes in the pore-size distribution.     

An analytical interpretation of liquid injection induced microseismicity in porous reservoirs

The paper presents a rigorously derived analytical method to describe and interpret the low-magnitude earthquakes caused by injection of the borehole fluids into surrounding porous reservoirs. Microseismicity is induced due to changes in the pore pressure, which, in turn, is influenced mainly by the low-frequency slow Biot (P2) wave. The classical Biot model is used to obtain the distribution of pore pressure in a reservoir. The constructed solution to the Biot system of equations and the spatio-temporal cloud of microseismic events allow one to assess the critical value of the pore pressure, sufficient for the generation of a microearthquake, and the values of hydromechanical parameters (e.g. permeability) of a saturated porous rock.     

基于Kelvin-Voigt黏弹性骨架的含非饱和流体孔隙介质BISQ模型(英文)

本文综合考虑了在波传播过程中孔隙介质的三种重要力学机制——"Biot流动机制一squirt流动机制-固体骨架黏弹性机制",借鉴等效介质思想,将含水饱和度引入波动力学控制方程,并考虑了不同波频率下孔隙流体分布模式对其等效体积模量的影响,给出了能处理含粘滞性非饱和流体孔隙介质中波传播问题的黏弹性Biot/squirt(BISQ)模型。推导了时间-空间域的波动力学方程组,由一组平面谐波解假设,给出频率-波数域黏弹性BISQ模型的相速度和衰减系数表达式。基于数值算例分析了含水饱和度、渗透率与频率对纵波速度和衰减的影响,并结合致密砂岩和碳酸盐岩的实测数据,对非饱和情况下的储层纵波速度进行了外推,碳酸盐岩储层中纵波速度对含气饱和度的敏感性明显低于砂岩储层。     

裂缝-孔隙介质储层渗透率表征及其影响因素分析

含裂缝多孔介质渗透率预测是非常规油气资源勘探开发的一个紧迫问题.现有多孔介质岩石物理模型通常利用圆形孔管模拟宏观岩石孔隙空间,难以定量描述软孔隙/裂缝在压力作用下的闭合情况,缺乏裂缝/孔隙间流量交换的连通机制.本文提出含三维裂缝/软孔隙网络多孔介质模型,将储层岩石裂缝/软孔隙表示为椭圆截面微管,建立了周期性压力作用下微观裂缝流量表达式,通过网络模型和流量守恒条件,得到含有三维裂缝/软孔隙网络的多孔介质渗透率计算方法.数值算例表明,预测结果与实验数据分布范围吻合很好,能够给出不同类型岩心对应孔隙纵横比的分布图.三维裂缝/软孔隙网络模型建立了宏观可观测量与裂缝参数之间关系,能够定量分析岩石渗透率随裂缝体密度、纵横比、孔隙流体类型和围压等因素的变化规律,为复杂条件下储层渗透率预测提供了一种有效方法.     

波浪荷载作用下海底隧道-孔隙海床的动力分析

孔隙海床在波浪荷载作用下其有效应力减小、孔隙水压力增加,这将影响孔隙海床中隧道的稳定性,因此,研究在波浪荷载作用下孔隙海床与隧道的动力响应具有重要的工程意义。本文基于Biot的动力固结理论和弹性动力学理论建立波浪荷载作用下孔隙海床与海底隧道的动力分析模型,并同时考虑海底隧道与海床之间的接触效应、边界效应等对海底隧道内力的影响。最后通过变换海床的变形模量、渗透系数和海底隧道的半径、埋深等,观察其对海床土孔隙水压力和海底隧道内力的影响,为海底隧道的设计提供依据。