不规则井眼的斯通利波反射数值模拟

讨论了不规则井眼内低频斯通利波反射.利用简化的一维模型研究低频斯通利波的物理性质.用质量平衡边界和传播矩阵表示井眼不规则度与斯通利波的相互关系.模拟结果表明,井径变化比岩性变化产生更强的斯通利波反射.这对实际测井资料解释有重要意义,即斯通利波解释裂缝时,应注意井径变化的影响.将合成模拟的波形与实测波形比较,可以识别裂缝的存在.     

基于非结构化网格的2.5维时间域激发极化反演

针对我国西部山区起伏地形,提出通过Delaunay三角剖分生成非结构化三角网格,可以根据实际地形情况,较好地控制节点分布,从而降低初始模型的几何离散误差。在反演中对模型参数施加最小光滑和背景约束信息,从而提高反演的分辨率。通过在正演求解稀疏矩阵时采用稳定型双共轭梯度法和计算系数矩阵时利用电位互换等手段,可有效地提高反演计算速度。最后,为验证反演方法的有效性,编制了跨平台反演软件;并对2例地电模型和实测数据进行反演试算。计算结果表明:反演耗费时间少,仅迭代3次拟合差便趋于稳定,反演结果能较好地刻画异常体形态。     

含地形起伏的MT模型非规则四边形FEM正演与反演

从大地电磁(MT)边值问题满足的变分方程出发,采用非规则四边形网格、双线性插值有限单元法(FEM)开展复杂起伏地形MT模型正演,探讨二维Jacobian变换行列式的计算方法,推导任意非规则四边形单元的插值方式及单元系数矩阵表达式,实现起伏地形MT模型的高精度正演。然后,介绍光滑约束的Tikhonov正则化反演算法,针对反演中正则化参数选取困难的问题,将广泛应用的L-curve法引入反演的正则化参数选取中。研究结果表明:L-curve法的曲线中曲率最大的拐点准确地指示了最优正则化参数;L-curve法选取的最优正则化参数对应的反演结果与原模型所示结果吻合度最高,反演效果最好。     

基于有限单元法的三维地电断面电阻率反演

在三维电阻率的有限单元正演模拟中,简化边界条件以减少计算时间,推导了用异常电位法计算偏导数矩阵的方法,克服了总电位法中的电源点附近节点电位模拟误差的影响,使偏导数计算精度大幅度提高;在最小二乘反演中,加入体积因子以解决多解性问题,并提高了反演精度;通过对典型模型的试算,反演结果与真实模型较吻合.     

利用合成地震波振幅预测煤层厚度及其方法研究

为了更精确预测煤层厚度变化,分析了波阻抗反演与煤厚的相关性,对测井资料和反演后的合成地震波振幅进行线性回归,建立了原生煤的正演模型.将该模型应用到非井点的地震属性上,利用同位协同克里金方法实现了对淮北桃园矿区煤层厚度的预测,结果与已知钻孔探测资料十分吻合,证明了利用地震振幅预测煤厚以及分类正演模型是可行的.     

基于三次B样条模型的油水相对渗透率曲线数值反演方法

引入三次B样条相渗表征模型,通过半迭代集合卡尔曼滤波(En KF)算法对见水前、后的生产动态数据进行自动历史拟合,提出了一种新的油水相对渗透率曲线数值反演方法。以某径向流数值岩心模型为例,对比分析了三次B样条模型与幂律模型的优劣。结果表明,相比于传统的幂律模型,三次B样条模型可对油水相对渗透率曲线进行局部逼近,灵活性更强,反演精度更高。通过对径向流岩心驱替实验数据进行自动历史拟合,反演估算了油水相对渗透率曲线。对驱替压差、累产油及累产水等观测数据拟合效果好,相对误差均小于5%,说明本文方法可靠性强,计算结果能够满足工程实际的需要。     

大地电磁主轴各向异性二维理论模型反演分析

地下介质普遍存在各向异性,为反应地下真实地质情况,需建立各向异性模型进行研究。本文主要进行主轴各向异性大地电磁二维反演研究,采用有限单元法作为正演模拟方法,将正演响应结果与前人计算结果进行对比分析,验证了算法的正确性,并对不同方向上电阻率响应灵敏度进行分析。采用非线性共轭梯度法进行反演研究,该方法不需要直接计算雅克比矩阵,相对于其他反演方法节省了计算时间,计算效率高;目标函数的建立对于不同方向上的电阻率采用不同的正则化参数来进行约束。通过复杂的各向异性体模型进行反演,结果显示对于不同方向上电阻率的恢复以及异常体位置的圈定都取得了较好的效果,说明了反演算法的有效性。     

基于射线追踪的微地震信号正演模型

为了更准确地收集微地震信号,采用射线追踪的方法对多震源微地震数据进行正演模拟。首先根据射线 追踪试射法原理,不断修正初始出射方向,通过反复迭代获得最优出射角,最终得到精确的射线路径以及多道 合成微地震响应记录,实现正演过程。基于射线追踪的正演结果能准确地反映地层结构的运动学特征,是一种 快速、准确的微地震正演方法。通过对地震波传播规律的分析,反演结果表明,建立的正演模型能为微地震信 号处理、反演以及震源定位研究提供帮助。     

基于二维高密度电阻率勘探数据的三维反演及应用

针对目前高密度电阻率法所采用的数据处理方法主要是把地质结构体视为二度体进行二维处理,从而二维数据资料处理结果只是一种近似解释,其计算精度与反演效果达不到精确反演的要求,设计一种典型的地质体模型,利用有限单元法进行正演计算,对含有高斯随机误差的各测线二维正演数据分别进行高密度电阻率法二维、三维反演,对比二维和三维的反演结果.最后在广东某场地采用高密度电阻率法对岩层划分,利用二维勘探的数据进行二、三维反演.研究结果表明:三维反演受高斯随机误差的影响更小,反演结果与实际地质模型更接近;三维反演在岩层划分中更能凸显岩层分界面区域,岩层的空间分布及岩层、覆盖层的整体连续性能更真实地体现,与钻探结果吻合程度更好.     

桩基动力学模型参数反演识别方法

针对桩基模型参数存在不确定性的问题,提出一种桩基模型参数反演识别方法.以测点和敲击点作为噪声因素,综合考虑模型中不同测点和敲击工况下加速度响应对反演结果的影响,利用多信息融合方法将多个评分值融合成唯一的参数反演识别目标响应值.研究结果表明:利用本研究方法得到的反演曲线与仿真曲线拟合良好.本研究为桩基模型参数反演识别提供了一定的借鉴.     

各向异性孔隙地层随钻声波测井理论模拟及应用

储层中普遍存在的各向异性和渗透性均会对弹性波传播产生影响。为满足油气资源勘探需求，对复杂储层中各向异性与渗透率的测量评价已逐渐从常规的电缆测井发展为随钻测井。本文建立了随钻声波测井模型，以Biot理论为基础研究了弹性波在横观各向同性孔隙地层中的传播机制，并对含有钻铤的充液井孔内各部分声场进行分析，导出了横观各向同性孔隙地层随钻声波测井的声场表达式。研究了随钻多极子声场模式波的频散特征。结合时间—慢度相关法分析了阵列波形中地层波和钻铤波的到时及慢度特征。计算结果表明随钻斯通利波和螺旋波速度频散曲线对地层各向异性和渗透率有不同程度的灵敏度，可用于实际资料的反演，并在此基础上处理了实钻井声波数据，对地层渗透率和各向异性进行了综合评价。本文的计算结果可作为复杂介质地层中随钻声波测井解释的理论基础。     

用有限元法计算井中水平裂缝的反射斯通利波

用有限元法研究了当井中存在水平裂缝时的声波全波列 ,并对裂缝宽度对低频斯通利波的影响进行了数值分析。当井中存在水平裂缝时 ,在接收的声波全波列波形中 ,既有直达的斯通利波 ,也存在来自裂缝的反射斯通利波。数值计算结果与实验测量结果具有很好的一致性。计算结果表明 ,裂缝宽度越大 ,反射斯通利波的幅度也越大 ,而直达波幅度变小     

斯通利波在测井应用中的研究进展和现状

斯通利波是一种沿井壁传播的声波,当声波脉冲与井壁和井内流体的界面相遇时就会产生斯通利波。近十几年来,由于长源距和偶极子阵列声波仪在测井中的广泛应用和全波形阵列声波测井资料中各组分波的处理分析方法的进步,可将斯通利波分离出来单独研究并应用于测井资料分析。斯通利波技术是评价裂缝及其渗透性的有效方法。     

采用混合法和递推矩阵算法模拟层状介质中随钻电磁波电阻率测量仪器的响应

将混合法与递推矩阵算法相结合计算轴对称条件下二维层状介质中随钻电磁波电阻率测量仪器的响应。根据层界面处电磁场的连续性条件得到确定所有待定系数的矩阵方程组并通过递推方法快速求解,以摒弃繁琐的上行波和下行波模式并避免指数增加项的出现。利用上述算法分析井眼、侵入和钻铤凹槽对随钻电磁波电阻率测井响应的影响。结果表明,由于金属钻铤的存在导致井眼内钻井液所占空间相对变小,从而使钻井液对接收线圈响应的影响变小;钻铤凹槽对随钻电磁波电阻率测量仪器相位差和幅度比的影响可忽略不计,但对单个接收线圈处电磁信号的相位和幅度产生明显影响,并且对幅度的影响大于对相位的影响。     

弹性半空间-理想流体层中Stoneley波的动力特性

为了解弹性半空一理想流体层中Stoneley波的动力特性，并为进一步开发利用海洋资源提供调查方法，根据弹性固体与理想流体动力学方程，推导了弹性半空间一理想流体层中Stoneley波的特征方程、基于该特征方程，分析了固液两介质的物理特性对流体动压力的影响，以及固体和流体介质的位移分布情况．结果表明，该体系中Stoneley波在固体中的位移分布类似于半空间Rayleigh波，而在流体中的位移分布有明显的规律性；弹性半空一理想流体层中Stoneley波具有频散性，且其波速小于弹性半空间中Rayleigh波的波速．     

各向异性人工材料的电磁波吸收特性

利用传输矩阵法研究了含各向异性人工材料的一维结构的电磁波吸收特性,分析了入射波频率、极化方向、入射角以及材料厚度对电磁波吸收率的影响.结果表明,电磁波垂直入射到单层各向异性人工材料中,尽管s波和p波的吸收率极大值出现在不同的频率,但两种波的吸收率随层厚增大的变化规律一致.斜入射时,不管入射角如何变化,s波和p波的吸收率都随着层厚的增大而趋于定值,该极限值由入射角决定.改变入射波频率,当某一频率对应磁导率张量的矩阵元的实部为负值,对应极化波的吸收率随层厚的增大而趋于定值;当磁导率张量的矩阵元的实部为正值,对应极化波的吸收率随层厚的增大而振荡,且峰值出现在层厚为λ/2的整数倍时.含各向异性人工材料的三层结构的吸收率优于单层结构.以上结论可为使用各向异性人工材料制作电磁波吸收器提供理论支持.     

地震勘探中平均速度的求取方法及应用

针对地震数据处理和资料解释过程中经常用到地层的平均速度,而平均速度的求取不易,尤其在测井资料很少的情况下,平均速度的求取更难等问题.依据地震勘探原理,在中、浅层地震勘探中,可将沉积地层近似为连续介质的地质模型,在连续介质中速度随深度的变化存在一定规律,不同地区有不同的经验公式,应用均方根速度在连续介质情况下的近似公式及数据处理中所得到的速度谱资料能够计算出平均速度的经验公式,得出的平均速度能够满足中、浅层地震勘探中时深转换的需要,另外,应用平均速度及钻孔资料可以对地震时间剖面进行相位标定.     

粗糙裂缝对井眼声波传播的影响

由于天然裂缝并非光滑裂缝,在通过声、电成像对裂缝进行识别过程中,如果将裂缝视为光滑裂缝,将会产生误差。斯通利波反射系数对裂缝宽度具有较高的灵敏度,常用斯通利波识别裂缝。本文分别采用物理实验和数值模拟的方法来研究粗糙裂缝对声波响应的影响,在实验方面,进行了模型井裂缝粗糙度实验,得到了粗糙裂缝的裂缝宽度和衰减系数之间的关系。对水平裂缝井眼模型进行数值模拟,研究正弦形状的粗糙裂缝对斯通利波的影响。采用迂曲度对裂缝表面的粗糙程度进行描述。实验结果显示随着裂缝表面刻槽的增加,斯通利波的衰减系数也会增加。模拟结果显示粗糙裂缝的宽度越小,和同样宽度的光滑的裂缝衰减系数之间的差距越大。通过拟合建立了斯通利波衰减系数与裂缝宽度、迂曲度的关系式。     

寺家庄矿煤层顶板岩性有色反演预测

煤层的顶板岩石的封堵性对煤层气的储集非常关键。有色反演可以获得相对波阻抗的变化,对测井的依赖度低、平面分辨率高,可用于少井区煤层顶板岩性的分布进行预测。在原始地震数据谱蓝化处理、提高分辨率的基础上,对寺家庄矿三维地震数据进行了有色反演,预测了9号、15号煤层顶板岩性。结果表明:相对波阻抗体异常与钻孔岩性较为吻合,并通过相对波阻抗平面分布图预测了寺家庄地区9号、15号煤层顶板岩性,且不同时间地层切片的反演结果与区域钻孔资料相吻合,说明有色反演结果可靠,可以对煤层顶板岩性分布特征进行有效预测。