共查询到17条相似文献，搜索用时 171 毫秒
1.
2.
3.
针对采用单发双收探头同时测量相位差和幅度比,测井响应与地层介电常数和电阻率两者都有关的电磁波测井资料反演问题，构造了其反演目标函数Jacobi矩阵〖WTHX〗J的快速算法.J直接给出了测井响应随地层参数变化的关系，利用J和它转置矩阵J〖WTBZ〗T的有意义乘积〖WTHX〗J·J〖WTBZ〗T和〖WTHX〗J〖WTBZ〗T.〖WTHX〗J〖WTBZ〗可以方便地对在非均质地层中测井响应特征和仪器的探测特性进行考察分析，为测井资料处理方法研究与仪器改进提供依据. 相似文献
4.
5.
过套管电阻率测井是通过测量套管壁电势实现测量地层的视电阻率,基于传输线方程理论,针对层状地层,给出了套管壁电势、电流对地层横向电阻导数的微分方程(称Jacobi矩阵微分方程)及边界条件;利用Jacobi矩阵微分方程边值问题导出了过套管电阻率测井反演地层参数的Jacobi矩阵系数的解析表示,利用Marquardt方法实现了过套管测井的地层电阻率反演;通过计算对Jacobi矩阵的特性进行了探讨,并获得了较快的计算速度(因为Jacobi矩阵是用解析解表示的),反演结果与地层模型取得了较好的逼近.本文实现了过套管电阻率测井地层参数的Jacobi系数矩阵的快速计算及地层电阻率反演,为进一步开展电阻率测井数据处理提供了理论依据和快速反演算法. 相似文献
6.
利用三维时域有限差分技术（FDTD）研究建立倾斜各向异性地层中多分量电磁测井响应的数值模拟算法.首先,利用Yee氏交错差分格式和线性内插公式对各向异性介质中时域Maxwell方程进行离散化,得到各向异性介质中各个交错节点上电磁场的时域递推公式,同时将面积加权平均与旋转矩阵技术相结合确定各个网格节点的等效电导率张量,此外,单轴各向异性吸收边界条件（UPML）被用于降低求解区域外边界反射.在此基础上,利用单频正弦磁偶源计算时域电磁场空间分布以及特殊的双方程算法（2E2U）从数值结果中快速提取电磁场幅度和相位,确定多分量电磁测井仪器的三维响应.最后,通过数值结果证明该算法的有效性,并系统考察不同层厚、地层倾角以及侵入带半径等变化对多分量电磁测井响应的影响. 相似文献
7.
Baker Hughes INTEQ的MPR是应用较广的传播电阻率测井仪器,它有4条相位差和4条幅度比测井曲线,本文以它在垂直井中的测井响应为例开展电磁传播电阻率资料的二维反演方法研究.由于其测井响应对地层纵向边界的探测灵敏度比对其他地层参数的大得多,在其资料反演中必须重视地层纵向边界位置的提取和反演精度.本文首先应用特征识别方法由分辨率最好的测井曲线提取地层纵向边界,接着对浅、深探测的两条相位差曲线分别做纵向一维反演,利用此结果来建立二维反演的初始地层模型,采用改进的阻尼型正则化高斯 牛顿优化算法,实现了MPR资料的二维全参数反演.对含23个纵向分层的两种典型地层模型的数值模拟测井资料的二维反演得到了令人满意的结果,它可显著提高MPR的资料解释水平.发现只用它的5条测井曲线就可得到与8条测井曲线几乎相同的反演效果. 相似文献
8.
为研究地层电性变化时不同区域对电磁波测井响应的贡献分布,从响应信号对地层参数求偏导的角度,给出一种新的电阻率敏感性函数定义,引入模式匹配法对纵向成层、径向非均匀介质敏感性分布进行快速模拟;通过对敏感性函数纵、横向积分给出了单发双收线圈系纵、横向探测范围,研究了井眼、频率、背景地层电阻率等对探测特性的影响.结果表明:敏感性函数能够定量表征响应对地层纵向与径向微观与宏观敏感性,幅度比与相位差敏感性分布形态类似,幅度比较广较深,而相位差分辨率高,敏感范围小;敏感性函数进行径向积分后可表征仪器的探测深度,与伪几何因子对比达到了同样的效果;背景地层电阻率在1~100Ωm变化时,工作频率2 MHz下幅度比探测深度约为0.6~2.3m,相位差为0.3~0.8m,幅度比50%纵向积分敏感性层厚约为0.3~1.6m,相位差约为0.2~0.6m;异常体与背景地层电阻率对比度在1~50变化时,引起的探测深度与敏感性层厚差异约为0.1~0.2m,远小于地层电阻率的影响. 相似文献
9.
为了提高随钻电磁波测井仪器的抗震能力,仪器传感器中所有发射与接收线圈均安装在柱状金属钻铤上的环状刻槽中,同时槽中充填高磁导率的铁氧体,以提高发射和接收效率.本文针对多环状刻槽钻铤中发射与接收线圈分布,研究建立一套随钻电磁响应的混合模拟算法.首先,根据环状刻槽几何尺寸以及轴向边界,将整个钻铤与外部的井眼和地层划分成多个不同电导率和磁导率分布的轴对称非均质各向同性层状模型,并应用混合法建立层状非均质地层中电磁场的半解析解与各个线圈上感应电动势的计算方法.其次,利用摄动原理给出随钻电磁测井中振幅比与相位差的空间灵敏度计算公式.然后,通过数值结果考察铁氧体磁导率和地层电导率变化对发射线圈的发射功率和接收线圈上感应电动势的影响,并确定各个接收线圈上感应电动势的动态变化范围,为仪器设计过程中选择合适的铁氧体提供理论依据.最后,通过空间灵敏度分布特征考察振幅比与相位差在探测特征上的差异. 相似文献
10.
为了实现交错沉积等复杂环境中的电磁场数值模拟,本文在常规横向同性模型的基础上引入了电导率主轴坐标系相对地层坐标系的层理方位角和倾角,建立了交错地层模型.并利用传播矩阵法建立了一维层状交错地层模型中的多分量感应测井仪器响应的正演模拟算法.首先将频率-波数域中的电磁场分解为上行和下行模式波,给出了任意朝向的磁偶极子在无限大地层中模式波的解析解.进一步通过引入地层界面上的透射、局部反射以及广义反射系数矩阵,推导了一维层状地层中的模式波表达式.在此基础上,利用二维Gauss-Legendre积分实现了Fourier逆变换,得到了可用于多分量感应测井模拟的频率-空间域磁场并矢格林函数.最后,通过多个数值模拟结果考察了井眼倾角、层理方位角和倾角变化对多分量感应测井响应的影响. 相似文献
11.
12.
In this paper, a method is proposed to invert permeability from seismoelectric logs in fluid‐saturated porous formations. From the analysis of both the amplitude and the phase of simulated seismoelectric logs, we find that the Stoneley wave amplitude of the ratio of the converted electric field to the pressure (REP) is sensitive to porosity rather than permeability while the tangent of the REP's phase is sensitive to permeability. The REP's phase reflects the phase discrepancy between the electric field and the pressure at the same location in the borehole. We theoretically derive the frequency‐dependent expression of the REP of the low‐frequency Stoneley wave and find that the tangent of the REP's argument is approximately in inverse proportion to permeability. We obtain an inversion formula and present the permeability inversion method by using the tangent of the REP's phase. To test this method, the permeabilities of different sandstones are inverted from the synthetic full‐waveform data of seismoelectric logs. A modified inversion process is proposed based on the analysis of inversion errors, by which the relative errors are controlled below 25% and they are smaller than those of the permeability inversion from the Stoneley wave of acoustic logs. 相似文献
13.
采用改进型逐次逼近解法(MSAM)和Aitken加速技术相结合的迭代算法计算二维积分方程.该算法将纵向成层原状地层作为背景地层,将计算区域限制在井眼和侵入带内,具有未知量数目少、收敛速度快、计算精度高的优点.利用该算法对随钻电磁波电阻率测量仪器在轴对称二维地层中的响应进行数值模拟.模拟结果显示,幅度衰减曲线和相位移曲线受井眼、侵入和围岩的影响程度不同,二者径向探测深度和垂向分辨率也有差异,利用补偿后的幅度衰减电阻率曲线和相位移电阻率曲线交叉点的坐标可精确确定地层层界面位置. 相似文献
14.
Maksim Bano 《Geophysical Prospecting》2004,52(1):11-26
The attenuation of ground‐penetrating radar (GPR) energy in the subsurface decreases and shifts the amplitude spectrum of the radar pulse to lower frequencies (absorption) with increasing traveltime and causes also a distortion of wavelet phase (dispersion). The attenuation is often expressed by the quality factor Q. For GPR studies, Q can be estimated from the ratio of the real part to the imaginary part of the dielectric permittivity. We consider a complex power function of frequency for the dielectric permittivity, and show that this dielectric response corresponds to a frequency‐independent‐Q or simply a constant‐Q model. The phase velocity (dispersion relationship) and the absorption coefficient of electromagnetic waves also obey a frequency power law. This approach is easy to use in the frequency domain and the wave propagation can be described by two parameters only, for example Q and the phase velocity at an arbitrary reference frequency. This simplicity makes it practical for any inversion technique. Furthermore, by using the Hilbert transform relating the velocity and the absorption coefficient (which obeys a frequency power law), we find the same dispersion relationship for the phase velocity. Both approaches are valid for a constant value of Q over a restricted frequency‐bandwidth, and are applicable in a material that is assumed to have no instantaneous dielectric response. Many GPR profiles acquired in a dry aeolian environment have shown a strong reflectivity inside dunes. Changes in water content are believed to be the origin of this reflectivity. We model the radar reflections from the bottom of a dry aeolian dune using the 1D wavelet modelling method. We discuss the choice of the reference wavelet in this modelling approach. A trial‐and‐error match of modelled and observed data was performed to estimate the optimum set of parameters characterizing the materials composing the site. Additionally, by combining the complex refractive index method (CRIM) and/or Topp equations for the bulk permittivity (dielectric constant) of moist sandy soils with a frequency power law for the dielectric response, we introduce them into the expression for the reflection coefficient. Using this method, we can estimate the water content and explain its effect on the reflection coefficient and on wavelet modelling. 相似文献
15.
探地雷达（GPR）时间域全波形反演计算量巨大,内存要求高,在微机上计算难度大.本文中作者基于GPU并行加速的维度提升反演策略,采用优化的共轭梯度法,避免了Hessian矩阵的计算,在普通微机上实现了时间域全波形二维GPR双参数（介电常数和电导率）快速反演.论文首先推导了二维TM波的时域有限差分法（FDTD）的交错网格离散差分格式及波场更新策略.然后,基于Lagrange乘数法,将约束问题转化为无约束最小问题,构建了共轭梯度法反演目标函数,采用Fletcher-Reeves公式与非精确线搜索Wolfe准则,确保了梯度方向修正因子及迭代步长选取的合理性.而GPU并行计算及维度提升反演策略的应用,数倍地提升了反演速度.最后,开展了3个模型的合成数据的反演实验,分别从观测方式、梯度优化及天线频率等方面,分析了这些因素对雷达全波形反演的影响,说明双参数的反演较单一的介电常数反演,能提供更丰富的信息约束,有效提高模型重建的精度. 相似文献
16.