用等效介质理论模拟裂缝溶洞型孔隙对速度和弹性模量的影响

针对缝洞型储层流体性质识别的困难,采用Kuster-Toks(o)z(K-T)模型、Berryman 自恰模型(SCM)和差分等效介质(DEM)三类近似模型,用不同纵横比的孔隙模拟了缝洞型储层的地震波速度和弹性参数响应,定量分析了溶蚀孔洞和裂缝二类次生孔隙对地下储层弹性性质的影响.明确了孔隙几何结构对岩石弹性参数的影响趋势和幅度大小,为利用地震信息评价孔隙结构提供了理论基础.根据物理模型的裂缝介质参数,分别采用K-T模型、SCM自恰模型和DEM差分等效介质模型,计算了裂缝介质的纵横波速度,考察了三个模型对缝洞型储层弹性响应计算的适用性,得到了有意义的认识.     

可控源海洋电磁勘探中空气波影响的理论分析和数值模拟

可控源海洋电磁勘探(MCSEM)中空气波对海底电磁响应的影响已为业界所熟知,且已提出了多种压制方法.然而,由于空气波与海水层中其它信号相互作用的复杂性,至今其作用机理及对有效信号的影响方式尚未完全清楚,这也阻碍了浅水域MCSEM油气勘探的应用.本文在前人研究的基础上, 基于一维层状介质的电磁位和电磁场分析,采用电磁场的模式分解理论导出了空气、海水层和海底地层之间的相互作用关系.分析了MCSEM电磁响应的影响因素和空气波的作用机理.以无限水深假想模型、三层介质模型和四层介质模型为基础,导出了水平电偶极子(HED)的横电(TE)和横磁(TM)模式的电磁场关系,分析显示了前人提出的上下行波分离进行空气波压制所存在的缺陷.通过理论和数值模拟知道, 利用横电(TE)和横磁(TM)模式受空气相互作用影响的差异和电磁场水平分量与垂直分量受空气层相互作用影响程度的不同可减弱空气波对有效信号的影响,这将有助于实现浅水域海洋电磁勘探数据的有效利用.     

地震相控混沌反演在岩性油气藏储层预测中的应用——以RY坳陷M井区为例

RY坳陷M井区位于浅水三角洲前缘斜坡带,多发育水下分流河道沉积,砂岩储层厚度薄、横向变化快,且同一时期的砂体在平面上表现为多套油水系统,容易形成岩性油气藏,利用常规手段进行储层预测的难度较大。本文通过引入地震相控制,利用纵向及横向上的构造起伏、岩性变化、沉积环境变迁等先验信息降低了反演的多解性;同时将混沌理论与反演算法结合,根据反演迭代的实时状态自适应地修改正则化算子,调节反演稳定性和分辨率达到最佳平衡,既避免了混沌现象的出现,又最大限度地揭示薄层信息。M井区的实际应用显示,相控混沌反演可以突出地震波形的细节变化,准确刻画储层尖灭和分叉点位置,达到精准预测岩性油气藏的目的,对后续岩性油气藏的勘探开发具有较大指导意义。     

2.5维电磁响应的有限元模拟与波数取值研究

利用有限元法实现了任意方向偶极子源在二维介质中频率域电磁响应的数值模拟,研究了波数取值对模拟结果的影响.通过对构造走向的Fourier变换,将全三维电磁问题,转化为一系列二维问题,并在波数域求解,极大地减小了计算工作量,导出了波数域耦合适用于二维电性介质中任何方向电或磁偶极子响应计算的电磁场方程.针对每个给定的波数,上述耦合电磁场方程用等参有限元方法在x-z平面内求解.采用Fourier逆变换,将波数域解积分,得到空间域电磁场.针对电磁模拟计算中,源点的奇异性,采用具有一定面积的伪δ函数表达源电流分布,使数值解精度得以提高.另外,采用等参有限元,使地下复杂地质体得到准确表达.利用不同波数值对均匀介质与层状介质的模拟结果与解析解的对比,验证了算法的正确性与精度.利用层状介质模型的解析解与数值计算结果的对比,分析了波数的优化取值范围及取值点数对数值模拟结果的影响,考察了算法对非均匀介质的适应性.     

可控源海洋电磁勘探中空气波影响的理论分析和数值模拟

可控源海洋电磁勘探(MCSEM)中空气波对海底电磁响应的影响已为业界所熟知,且已提出了多种压制方法.然而,由于空气波与海水层中其它信号相互作用的复杂性,至今其作用机理及对有效信号的影响方式尚未完全清楚,这也阻碍了浅水域MCSEM油气勘探的应用.本文在前人研究的基础上, 基于一维层状介质的电磁位和电磁场分析,采用电磁场的模式分解理论导出了空气、海水层和海底地层之间的相互作用关系.分析了MCSEM电磁响应的影响因素和空气波的作用机理.以无限水深假想模型、三层介质模型和四层介质模型为基础,导出了水平电偶极子(HED)的横电(TE)和横磁(TM)模式的电磁场关系,分析显示了前人提出的上下行波分离进行空气波压制所存在的缺陷.通过理论和数值模拟知道, 利用横电(TE)和横磁(TM)模式受空气相互作用影响的差异和电磁场水平分量与垂直分量受空气层相互作用影响程度的不同可减弱空气波对有效信号的影响,这将有助于实现浅水域海洋电磁勘探数据的有效利用.     

浅水域可控源海洋电磁测量中空气波的渐近表达与波场分离方法

可控源海洋电磁(MCSEM)勘探中空气波对海底电磁响应的影响已为业界所重视, 它是由水平电偶极子源发射的电磁信号沿空气-海水界面传播与来自海底地层的有效信号相互作用产生的。在浅水域勘探时, 空气波淹没来自地层的有效信号, 使浅水域MCSEM实现油气层识别产生困难。基于层状介质模型, 采用电磁场的模式分解理论导出了半空间电阻率模型的空气波表达式, 利用该式将空气波在海水层的传播近似用海水-空气界面与海底地层之间的多次鸣震表达。结合源和接收器两边电磁场的传播特征, 导出了有限水深时空气波近似关系, 用于近似模拟空气波响应。为了对比不同方法压制空气波的效果, 基于上行波场和下行波场的分解方法, 获得了含油气高阻储层上上行波的异常幅度增大数倍的结果, 显示了波场分解方法压制空气波的良好效果。最后, 利用不含油气层的背景模型和含油气储层模型电磁响应的数值模拟结果, 比较空气波渐近表达和波场分离2种空气波压制方法可知:对于水平层状模型后者效果更好;前者可适用于崎岖海底地层的数据处理, 后者只适用于水平海底地层。     

2.5维井间电磁反演在中国东部孤岛油田的应用

In this study, we present a practical technique of transforming cross-hole EM data into the inter-well resistivity distribution. The a priori information constraint is incorporated into an iterative regularized inversion procedure and a variable roughness is added into the inversion process. Finite element approximation based on a two and a half-dimensional （2.5D） model has been developed for the forward problem and the ＂pseudo-forward＂ problem needed for constructing the sensitivity matrix and synthetic data set. The regularized least-squares inversion scheme, constrained with the a priori information obtained from well logs, was adopted to reconstruct the inter-well resistivity profile from two synthetic electromagnetic data sets and field data acquired in the Gudao Oil Field, East China. The partial derivatives of the sensitivity matrix were computed by the adjoint equation based on the reciprocity principle. Inversion results of the synthetic and field data examples suggest that our method is robust and stable in the presence of random noise in the field data and can be used for cross-hole EM field data interpretation.     

电性各向异性地层频率域电磁响应模拟

在现有计算机内存和运算速度条件下,电性各向异性地层中频率域电磁响应的数值模拟是一个难点. 为探索各向异性介质模型电磁响应的有效算法,同时为进一步研究各向异性地层测量响应的校正方法奠定基础,从Maxwell电磁响应方程出发,使用由三个彼此垂直的发射线圈和三个彼此垂直的接收线圈组成的线圈结构模型,应用交错网格有限差分法推导了三维频率域电磁响应的差分计算格式,采用预条件双共轭梯度法对所形成的复系数矩阵进行了迭代求解,从而得到电性各向异性地层模型下的线圈测量响应. 应用这些方法进行了多种条件下电磁响应的计算,计算结果表明,所形成的计算方法是准确可靠的,预条件双共轭梯度法的运算效率也较高;电性各向异性地层的电磁响应与地层倾角（或井斜角）、仪器测量方位、地层的水平电阻率和垂直电阻率有关;根据水平方向磁场强度分量可以确定地层界面;在电磁响应计算过程中,背景电导率的取值要尽量接近整个计算区域的平均电导率.     

用快速汉克尔变换计算各向异性层状介质中点电流源的电位响应

针对任意各向异性地层,利用快速汉克尔（Hankel）变换导出了点电流源在地下介质中的电位响应关系,计算了各向异性地层的直流电视电阻率,重点分析了电阻率测深方法对地下各向异性介质的探测能力。这里采用二维傅立叶变换,首先导出了均匀各向异性层的电位一般解函数,利用各层界面电位和电流密度垂向分量的连续性及地下电流源的耦合条件,推导出不同层之间电位和电流密度垂向分量的传递函数,建立了空间电位和电流密度的递归计算关系。针对源点附近电位函数的奇异性,采用柱坐标变换,将二维傅立叶变换转化为汉克尔（Hankel）变换,并用快速汉克尔（Hankel）变换计算所有的位函数值。采用均匀各向异性半空间的解析解,验证了文中算法的正确性,给出了多层各向异性地层模型的视电阻率响应曲线,分析了直流电法探测裂缝性地层,估计裂缝分布性状的可能性。