方位梯度电极系的正反演

针对环井周的非均质地层模型,利用三维有限元方法模拟了方位梯度电极系的测井响应．结果表明,在非均质程度比较小的情况下,各个方位上的响应之间差别不大,与真实的电阻率分布差别也不大;反之,则差别很大．其特点是所测得的响应远远低于较高的真电阻率值,而接近较低的真电阻率值．因此,有必要对测井响应进行校正．考虑到传统的校正方法的弊端,本文利用非线性反演方法求取真实的电阻率分布,数值模拟的结果表明了此方法的可行性．     

利用共轭梯度算法的电阻率三维有限元正演

引入对称超松弛预条件共轭梯度(SSOR PCG)迭代算法求解电阻率三维有限元计算形成的大型线性方程组，并有机结合系数矩阵的稀疏存储模式，使得三维有限元正演计算的速度大大提高而内存需求则大大减少. 该算法可方便地应用于求三维异常电位，在保持快速计算的基础上，正演计算的精度显著提高.     

基于Matlab的一维电测深正反演可视化软件设计

基于Matlab设计一款一维电测深正反演可视化软件,主要包括正演拟合、直接反演和影响系数计算等功能。该软件界面简洁友好,反演速度快,结果误差小。该软件的使用,可为科研工作提供更加可靠的电性结构模型。     

裂缝性地层方位侧向测井响应数值模拟

为研究井周裂缝发育特征,本文提出一种新型方位侧向测井方法,利用三维有限元法,模拟裂缝的方位侧向测井响应.结果显示,深浅侧向电阻率幅度差异受裂缝倾角的控制,低角度缝为负差异,高角度缝为正差异;倾斜裂缝张开度的增大使测井响应值减小,方位电阻率差异增大;井周方位电阻率可反映裂缝方位产状,单一缝或裂缝密度较小时,沿裂缝走向的方位电阻率小,沿裂缝倾向的方位电阻率大;裂缝发育地层的测井响应显示宏观各向异性特征,但方位电阻率的差异显示发生反转现象,即沿裂缝走向/层理方向的方位电阻率大,沿裂缝倾向/垂直层理方向的方位电阻率小;对方位电阻率测井响应进行井周成像,直观显示了裂缝的产状和发育特征.

     

积分方程法复杂介质可控源电磁勘探快速正反演研究现状及发展方向

复杂介质可控源电磁勘探数值模拟及反演算法的研究一直是国内外地球物理学者研究的热点。本文对复杂介质可控源电磁勘探快速正反演算法研究进行综述,重点对复杂介质快速正反演算法及应用进行分析,指出高效并行、特殊边界条件或将是其真正实用化的关键,当前仍然是极具挑战的研究方向。着重对未受关注的可控源电磁法复杂介质积分方程法正反演算法及其应用研究,如二维、2.5维快速正反演算法;地面、井筒电磁勘探实例、起伏地形异常场模拟等进行讨论。指出国内积分方程法的研究相对滞后,但应用前景较可观;特别是大尺度隐伏资源勘探领域,高精度、高效电磁勘探正反演需求较迫切。通过体积分方程法快速正反演算例分析,表明该方法可适用于大尺度勘探生产,具有较好实用性。复杂地形模拟,高效正反演算法等是积分方程法实用化的关键。      

可控源电磁场三维自适应矢量有限元正演模拟

本文实现了可控源电磁（CSEM）场三维自适应矢量有限元正演算法,该算法采用非结构四面体单元进行三维网格剖分,能够真实模拟地形起伏和复杂电性异常体.采用一次场和二次场分离的方式计算电磁场响应,能够有效解决有限元模拟中的源点奇异性,提高场源附近电磁场数值精度,其中一次场利用CSEM一维正演算法解析求得,二次场采用矢量有限元方法求得.并利用基于后验误差估计的自适应网格细化算法指导网格细化,以减少人为设计网格导致的误差.通过一维和三维模型的数值模拟,验证了本文算法的有效性：一维模型有限元解与解析解吻合得很好,电磁场振幅相对误差在1%左右,相位差整体小于1°;三维模型有限元解与有限体积解吻合得也很好.模拟了一个含三维倾斜板状异常体的可控源电磁响应,表明了本文算法模拟复杂地电结构电磁场的能力和有效性.

     

基于二次场方法的并行三维大地电磁正反演研究

快速且高精度的三维大地电磁法正反演是目前研究的热点.由于大地电磁法场源的平面波特性,以往的正演方法大多采用直接求解总场的方法,在边界强加二维边界条件.本文提出了一种基于二次场方法的三维大地电磁法正演算法,将平面波在层状背景模型中的响应作为场源项,得到二次场满足的偏微分方程,并利用交错网格有限差分法求取二次场.与其他学者的基于总场方法的结果的对比证明了本文采用方法的正确性.在基于二次场的正演算法基础上,实现了基于L-BFGS的三维反演方法,并对公开的数据集进行了反演.另外,针对大地电磁法的多频率观测特性,采用了基于MPI的分频并行策略对程序进行并行化,可达到接近线性的加速比.

     

关于三电极系侧向测井纵向几何因子

作为《三电极系侧向测井径向几何因子》一文的继续,本文讨沦了三电极系侧向测井仪的纵向特性。除了现在考虑的介质是层状结构,并且略去了井眼和侵入带的影响之外,所有的前提都和参考资料[1]相同。给出了近似电位场和纵向几何因子。     

新型贴井壁式阵列方位侧向测井数值模拟与响应特性

在现有侧向测井仪器基础上,提出了一种新型贴井壁式阵列方位侧向测井电极系,该仪器能提供径向和周向方位电阻率测量.该电极系采用贴井壁测量方式,提供5种径向探测深度.利用三维有限元方法模拟了阵列方位侧向测井电极系的井眼影响特性、径向探测特性、纵向分层能力、方位分辨能力,并模拟其对井周地层、水平井非对称泥浆侵入和倾斜地层的响应.在导电泥浆中最大探测深度为1.23m,纵向分辨率为0.3m,可以识别出0.1m薄层,方位分辨率为20°.贴井壁测量时,纵向分辨率不受泥浆和围岩电阻率的影响,能够准确测量井周方位电阻率,较不贴井壁测量具有很大优势,同时利用12条方位电阻率曲线能够反映出水平井泥浆非对称侵入特性,倾斜地层倾角和倾斜方向.     

用梯度和感应测井联合确定地层电阻率各向异性

实际地层可能出现2种电阻率的宏观各向异性:微观的统计平均效应和裂缝走向等因素的构造效应.三分量感应测井方法可能是确定地层电阻率各向异性的最好方法,但是,该仪器目前在中国还没有应用.已有的数据一般是梯度和双感应曲线,梯度和双感应曲线单独使用都难以揭示各向异性,但是,将二者联合使用则有可能.在反演模型中考虑了电阻率各向异性的存在,给出了梯度和双感应曲线联合求地层电阻率各向异性的方法.通过直井水平层的二维人工模型表明,用感应方法可确定地层的水平电阻率和厚度.而梯度方法的视电阻率可以近似看作水平电阻率和纵向电阻率的几何平均值,梯度方法的视厚度是各向异性参数和真实厚度的乘积.因而不能用梯度方法单独确定各向异性参数和真实厚度2个参数中的任何1个.然而,联合梯度和双感应的方法则可能确定3个参数:各向异性参数,水平电阻率和地层厚度.人工数据的例子表明了这种联合的可行性.     

用梯度和感应测井联合确定地层电阻率各向异性(英文)

实际地层可能出现 2种电阻率的宏观各向异性 :微观的统计平均效应和裂缝走向等因素的构造效应。三分量感应测井方法可能是确定地层电阻率各向异性的最好方法 ,但是 ,该仪器目前在中国还没有应用。已有的数据一般是梯度和双感应曲线 ,梯度和双感应曲线单独使用都难以揭示各向异性 ,但是 ,将二者联合使用则有可能。在反演模型中考虑了电阻率各向异性的存在 ,给出了梯度和双感应曲线联合求地层电阻率各向异性的方法。通过直井水平层的二维人工模型表明 ,用感应方法可确定地层的水平电阻率和厚度。而梯度方法的视电阻率可以近似看作水平电阻率和纵向电阻率的几何平均值 ,梯度方法的视厚度是各向异性参数和真实厚度的乘积。因而不能用梯度方法单独确定各向异性参数和真实厚度 2个参数中的任何 1个。然而 ,联合梯度和双感应的方法则可能确定 3个参数 :各向异性参数 ,水平电阻率和地层厚度。人工数据的例子表明了这种联合的可行性     

地震勘探中广义弹性阻抗的正反演

常规的地震道反演方法建立在反射P波垂直入射假设 的基础上,而实际地震资料采集时多数是非零炮检距的,反射振幅是共中心点道集叠加的结 果 . 因此,利用常规地震道反演方法就不能得到可靠的波阻抗或其他岩性信息. 本文利用Patr ick Connolly弹性阻抗的思想,通过对Zoeppritz方程的进一步简化,推导出适合常规叠后 资料的、非零炮检距条件下纵波反射系数递推公式,提出了广义弹性阻抗的概念,解决了非 零炮检距条件下,常规叠后地震道正反演的关键问题. 广义弹性阻抗不仅包含波阻抗,还包 含了纵横波速度等岩性信息,具有很好的实用价值. 进行广义弹性阻抗的反演,能较常规地 震道反演获得更多、更可靠的流体、孔隙度、砂泥含量等信息,有助于解释常规地震道反演 和道积分剖面中的假象,降低反演的多解性,提高储层预测的精度.     

三维井间电磁场的正反演计算

The method of modified local nonlinear iteration （MLNI） for volume integral equations is introduced and used to calculate forward and inverse problems of 3 D cross hole electromagnetic fields. The large scale scatters are divided into two kinds, the near region scatters and the far region scatters, whose locations and sizes are determined by the position of the field point. The influence of near region scatters is computed by the method of local nonlinear approximation, while that of the far region scatters is treated as external exciting source and is calculated reiteratively. The method has the advantages, such as faster calculating speed, less computer memory requirements and more efficient convergency. The nonlinear inversion algorithm based on MLNI is applied to inversion problems which can concentrate the imaging region to a particular range instead of whole 3 D domain. Since the algorithm considers the effect of the nonlinear term in each element of the Jacobi matrix, it is more accurate and stable in comparison with the convenitional iterative Born inversion method. The double cell decomposition method is applied during the first rough imaging process considering the limitation of the information and computer memory, then the imaging range is reduced and the second more accurate imaging process can be carried out. Numerical results indicate that MLNI is an efficient method to compute the scattering fields of large scale scatters and that the application of the method to inversion can give high resolution 3 D cross hole conductivity images.     

地壳三维构造反演和速度层析成像

本文给出利用空间深地震测深反射资料重建地壳三维构造和速度分布的方法，适用于任意分层并且每层可由若干断块组成的地壳模型。通过采用适当的反演方法获得地壳界面的三维构造，并在此结果基础之上，进一步将剩余走时残差归因于地壳速度的不均匀性，采用东分块的反演方法，重建地壳三维速度分布，数值模拟的结果表明算法是有效的。     

三维井间电磁场的正反演计算

提出了一种计算体积分方程的方法－改进型局域非线性迭代（MLNI），并用其对三维井间电磁场进行反演计算，该方法将井间大尺度散射体分为近场和远场区域两部分，它们的位置和尺寸均随场点位置的变化而改变，采用局域非线性近似计算近场区域的影响，将远场区域的影响作为外部激励源，采用迭代方法计算，该方法具有计算速度快，所需内存量少，收敛性好的优点，在反演中采用基于MLNI的非线性反演方法，利用这种方法可以将成像区域集中于一定范围内而不是整个三维空间，该反演方法由于考虑到了雅可比矩阵元素的非线性项，因而与传统的Born迭代反演方法相比更为精确和稳定，考虑到信息量和计算机内存的限制，第1次成像采用双重体元分割法进行较粗略的成像，然后缩小成像范围进行第2次较精确的成像，数值计算结果表明，MLNI是一种有效的计算井间大尺度异常体散射场的方法，将该方法用于反演过程能够得到较高分辨率的三维井间电导率图像。     

APPLYING 3D INVERSION OF SINGLE-PROFILE MAGNETOTELLURIC DATA TO IDENTIFY THE SHADE AND YUNONGXI FAULTS

Many synthetic model studies suggested that the best way to obtain good 3D interpretation results is to distribute the MT sites at a 2D grid array with regular site spacing over the target area. However, MT 3D inversion was very difficult about 10 years ago. A lot of MT data were collected along one profile and then interpreted with 2D inversion. How to apply the state-of-the-art 3D inversion technique to interpret the accumulated mass MT profiles data is an important topic. Some studies on 3D inversion of measured MT profile data suggested that 2D inversions usually had higher resolution for the subsurface than 3D inversions. Meanwhile, they often made their interpretation based on 2D inversion results, and 3D inversion results were only used to evaluate whether the overall resistivity structures were correct. Some researchers thought that 3D inversions could not resolute the local structure well, while 2D inversion results could agree with the surface geologic features much well and interpret the geologic structures easily. But in the present paper, we find that the result of 3D inversion is better than that of 2D inversion in identifying the location of the two local faults, the Shade Fault(SDF)and the Yunongxi Fault(YNXF), and the deep structures. In this paper, we first studied the electrical structure of SDF and YNXF based on a measured magnetotelluric(MT) profile data. Besides, from the point of identifying active faults, we compared the capacity of identifying deep existing faults between 2D inversion models and 3D models with different inversion parameters. The results show that both 2D and 3D inversion of the single-profile data could obtain reasonable and reliable electrical structures on a regional scale. Combining 2D and 3D models, and according to our present data, we find that both SDF and YNXF probably have cut completely the high resistivity layer in the upper crust and extended to the high conductivity layer in the middle crust. In terms of the deep geometry of the faults, at the profile's location, the SDF dips nearly vertically or dips southeast with high dip angle, and the YNXF dips southeast at depth. In addition, according to the results from our measured MT profile, we find that the 3D inversion of single-profile MT data has the capacity of identifying the location and deep geometry of local faults under present computing ability. Finally, this research suggests that appropriate cell size and reasonable smoothing parameters are important factors for the 3D inversion of single-profile MT data, more specifically, too coarse meshes or too large smoothing parameters on horizontal direction of 3D inversion may result in low resolution of 3D inversions that cannot identify the structure of faults. While, for vertical mesh size and data error thresholds, they have limited effect on identifying shallow tectonics as long as their changes are within a reasonable range. 3D inversion results also indicate that, to some extent, adding tippers to the 3D inversion of a MT profile can improve the model's constraint on the deep geometry of the outcropped faults.     

重磁异常反演的拟BP神经网络方法及其应用

把神经网络与重磁异常反演理论相结合,提出了用于重磁反演的一种拟BP神经网络方法．基于3层神经网络结构,把隐含层神经元设定为三维空间物性（磁化强度或密度）单元．对实测与理论重磁异常经S型函数变换,采用自动修改物性单元物性值的拟BP算法,反演三维空间的物性分布．利用该网络对理论模型数据和内蒙古某花岗岩体上的航磁资料进行了反演计算,取得了满意的反演效果．     

新丰江水库三维速度结构和震源参数的联合反演

本文对ACH方法做简化处理,先以一个地震为单元形成子矩阵,并将慢度扰动用一平均量来代替;在修正参数时,震源部分按常规进行,对慢度部分则利用代数重建技术将其分配到各个块体上.数值模拟结果表明,该方法能很好地反演震源位置和速度结构.运用新丰江遥测台网的数据进行研究,表明地震强活动区位于大坝附近,特别是在人字石断裂和高寨断裂的交汇处,震源集中沿罗坑—双下、葫芦凹—燕岩和黄竹蒿3条地震带,呈密集的高倾角分布,地震活动主要与北北西向构造和东西向构造有关,而与地表甚为发育的北东东向构造关系不大.速度图像中的低速区对应地表的破碎区,即地震带区域,也是该区重磁异常急剧变化的区域;而高速区内地震较少.上述特征同水的渗透作用有密切关系.