用改进的传输线算法计算水平层状横向同性地层中海洋可控源电磁响应

利用二维Fourier变换与电磁场分解技术将层状横向同性地层中Maxwell方程转化成两个独立的关于横磁(TM)波和横电(TE)波的传输线方程; 借助传输线理论与叠加原理, 仅利用电流源传输线Green函数得到TM波和TE波的解, 改进传输线算法, 建立横向同性地层中频率-波数域电流源电场和磁场并矢Green函数的新算法与新的解析表达式, 提高海洋可控源电磁响应数值模拟效率. 在此基础上, 利用传输线Green函数的基本解以及边界条件, 推导出广义反射系数与振幅递推公式, 得到各个地层中传输线Green函数的解析解; 然后利用Fourier逆变换与Bessel公式将海洋可控源电磁响应表示为Sommerfeld形式的积分, 借助三次样条插值与Lommel积分公式快速计算其数值解. 通过数值模拟结果考察工作频率以及地层各向异性电阻率变化等对海洋电磁响应的影响. 关键词： 传输线法 横向同性地层 海洋可控源电磁 Sommerfeld积分     

耦合势有限体积法高效模拟各向异性地层中海洋可控源的三维电磁响应

本文基于电场矢势与标势分解的耦合势有限体积法研究建立一套各向异性地层中海洋可控源电磁法的三维响应的高效数值模拟技术.首先引入电场的矢势和标势,将电场分解为无散场和无旋场之和,Maxwell方程转换为关于矢势与标势的混合Helmholtz方程,克服低感应数问题.在此基础上,借助Yee氏交错网格和有限体积法以及非均质单元中等效电导率公式,建立混合Helmholtz方程的离散方程.并采用直接法求解器PARDISO求解离散方程,有效保证在大的求解空间中仍然能够获得电磁场稳定可靠的数值解.此外,在数值模拟中利用差异场技术,克服源的奇异性问题,尽可能提高近场的计算精度.与解析解的对比证明了该算法的有效性.数值模拟结果表明,海洋可控源电磁法沿测线方向的电场,对油气藏的纵向电阻率敏感,对横向电阻率不敏感;对油气藏上方的覆盖层的纵向电阻率和横向电阻率都敏感.     

水平层状各向异性介质中电磁场并矢Green函数的一种高效算法

利用高阶窗函数结合连分式展开等技术研究并建立一种水平层状各向异性介质中电磁场并矢Green函数的快速有效算法. 首先借助于高阶窗函数将构成并矢Green函数的Sommerfeld积分转化成广义快速下降路径上积分,并给出高阶窗函数Hankel变换的一种新的更高阶幂级数展开式以及严格的Lommel函数表达式,以满足在全空间上高精度计算并矢Green函数的要求. 在此基础上,用Bessel函数的零点将积分路径划分成一系列小区间并通过改进的自适应Gauss求积公式确定各个小区间上的积分值,然后引入连分式展开法对各个区间上的积分值求和,从而使整个积分的收敛效率得到大大提高. 最后通过数值结果验证本方法的有效性. 关键词： 高阶窗函数 连分式展开 并矢Green函数 层状各向异性介质     

海洋可控源三维电磁响应显式灵敏度矩阵的快速算法

借助电场耦合势三维有限体积法与直接求解技术,研究建立了一套海洋可控源三维电磁响应显式灵敏度矩阵(或称为Fréchet导数)高效算法.首先,利用Yee氏交错网格和有限体积法对电场混合势Helmholtz方程进行离散处理,建立与移动源电磁场正演模拟相对应的大型代数方程组,再应用直接法得到的逆矩阵和三维线性插值技术事先确定插值算子和投影算子,并利用投影算子与各个发射源离散向量的乘积计算多发射源电磁响应,极大地提高了多发射源电磁场正演模拟效率.在此基础上,根据块状模型和像素模型中异常体电导率分片常数分布特征,将电导率摄动产生的一次散射电流场表示成Yee氏剖分网格上散射电流元的叠加,由投影算子与散射电流元的离散向量的乘积快速计算出电场强度与磁场强度的显式灵敏度矩阵.最后,通过数值计算检验算法的有效性,并通过块状模型与像素模型分别研究海洋可控源电磁响应特征.     

非均质各向异性地层中方位随钻电磁测井响应三维有限体积法数值模拟算法

方位随钻电磁测井是一种能够实时探测地层边界、实现地质导向与井眼成像的新型测井技术. 本文根据方位随钻电磁测井仪器的典型线圈系结构, 首先引入柱坐标系下非均质完全各向异性地层中电流源并矢Green函数, 并利用电磁场叠加原理给出倾斜发射线圈激发的电场以及倾斜接收线圈上感应电动势的计算公式; 然后应用电流源电场并矢Green函数的混合势克服非均质地层中电磁数值模拟的低感应数问题, 通过ρ 和z方向上Lebedev网格设法降低网格节点个数, 并且利用标准化算法确定柱坐标系下非均质单元上的等效电导率. 在此基础上, 用三维有限体积法建立柱坐标系电场混合势的离散方法, 得到一个交错网格上电场矢势和标势大型代数方程, 并用不完全LU分解以及稳定双共轭梯度法确定数值解. 最后, 通过数据模拟结果对算法的有效性进行检验, 并考察钻铤、线圈倾斜角度以及地层各向异性等参数对仪器响应的影响. 数值结果表明: 在柱坐标系下用三维有限体积法的数值模拟算法处理非均质各向异性层中方位随钻电磁测井响应可以得到很好的结果. 钻铤、电导率各向异性、层边界均对方位随钻电磁波测井响应产生较大的影响; 在电阻率较大的地层, 幅度比和相位差响应越小; 发射线圈和接收线圈同时倾斜时, 幅度比和相位差响应受地层的影响更灵敏.     

勒让德方程本征值的确定

大部分数学物理方法教材直接给出了勒让德方程本征值的表达式,比较突兀,学生也难以理解.为消除学生在此处的学习障碍,直接从勒让德方程的一般形式入手,通过常点邻域的级数解法,推导出勒让德方程本征值的表示形式.     

水平层状横向同性地层中频率测深资料全参数快速迭代反演

提出水平层状横向同性地层中频率测深资料的全参数快速迭代反演算法,以便从测量资料中同时确定各个地层的横向、 纵向电导率及层界面深度.首先,利用水平层状介质中电磁场并矢Green函数在频率波数域中的解析解和Sommerfeld积分的 快速计算技术确定正演响应.然后,利用摄动理论和Fourier逆变换公式,研究建立一套快速求解全参数Fréchet导数的 有效算法,并利用规范化处理和奇异值分解技术给出迭代反演过程.最后,给出数值试验证明反演理论的有效性以及反演算法 的抗噪声能力.