简介：坐标变换法通过将物理空间的曲网格映射为计算空间的矩形网格,将起伏地表转化为水平地表,同时将物理空间的波动方程转化为计算空间的波动方程,在计算空间完成数值模拟,坐标变换的方法对处理起伏自由边界具有较好的适应性和应用效果。本文在传统坐标变换方法的基础上,根据计算区域速度差异采用不同的网格大小和采样时间步长,提出了一种基于时空双变网格的起伏地表坐标变换正演模拟方法。在编程实现算法的基础上,通过典型模型波场模拟试算结果分析可知：（1）变网格方法与常规方法波场模拟误差在0.5%左右;（2）变网格方法计算效率视不同的变网格区域面积及变网格大小可提高几倍量级,在本文模型和计算参数下提高约5倍。（3）在满足模拟精度及频散条件要求下,变网格方法较全局细网格算法能显著节约计算内存。为此,针对起伏地表数值模拟,本文方法具有较高的模拟计算精度和一定的适应性。

简介：阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以Zxx,Zxy,Zyx,和Zyy为元素的2×2矩阵。在本次研究中,6个元素的阻抗张量的计算使用了一个含有Zxx,Zxy,Zyx,Zyy,Zxz和Zyz分量的2x3矩阵。对上述两类阻抗张量元素的属性进行了分析。利用由印度古吉拉特邦卡奇沉积盆地采集的5个分量大地电磁数据测试了文中的方法。从视电阻率和相位的计算中我们是观察到在大部分的频带范围内4个元素阻抗和6个元素阻抗Zxy和Zyx两类元素区别不大。然而,较长周期时间的数据,如超过100秒,观察到视电阻率的增加和相位的减少。我们还注意到,倾子幅度在大部分时间几乎是零,但较长周期（超过100秒）,逐渐呈增加的趋势。卡奇沉积盆地的地电断面表明在较长的周期内浅层近水平层和深层异常高电导性的不均质层都可能是引起大的Hz分量的原因。这表明,磁场垂直分量Hz对在大的2D/3D结构区域内的电场参数估计发挥的重要作用。

简介：逆时偏移被认为是处理陡倾角断层、回转波和逆掩推覆构造最有效的成像方法。考虑到当今勘探对地下复杂构造体、薄互层和隐蔽储层的高分辨率成像需求，井间、井中地震资料成为地面地震资料的一个很重要的补充。但由于地面地震、井中地震和井间地震存在很大的频率差异，本文将相应的变网格算法应用到偏移中。另外，也将Lanczos滤波引入到波场外推计算中，从而可以很好地消除变网格引起的人为反射。试算表明，对比地面、井中和井间的偏移，可以看到微构造、陡反射层都得到了更好地刻画和描述。基于变网格的逆时偏移在储层预测和岩性解释方面有广阔的应用前景。

简介：地震波场数值模拟方法对理解和分析地震波的传播规律具有着重要的意义。弹性波动方程能够模拟地下介质的实际情况，为偏移和成像提供有效的依据。在弹性波波场数值模拟中，旋转交错网格数值模拟（RSM）修改了标准交错网格数值模拟（SSM）方法，将同类的参数定义在同样的节点上，拓宽了稳定性条件的约束，但在低速区会出现较严重的频散。变阶数差分方法是自适应空间算子长度方法的一种变化和推广。它以理论频散误差研究为基础，结合实际波场传播的情况进行误差计算，对不同速度匹配不同的差分阶数。本文研究了变阶数旋转交错网格数值模拟（VRSM），即是籽变阶数方法应用到RSM中，它可以很好地解决RSM在低速区域的数值频散问题，以及减少不必要的时间损耗；同时讨论了旋转交错网格的理论频散特性，并基于波场分离的方法分析了实际波场传播的频散误差，将原方法的应用范围由声波推广到剪切波，由理论值推广到时变值。在数值模拟试验中，VRSM将被应用于水平层模型和Overthrust模型。通过阶数分配以及相应波场传播效果和计算时间的分析，验证了该方法应用于复杂介质波场模拟中的实用性和有效性。实验的结果表明VRSM能够合理分配不同速度所对应的差分阶数，能保证计算的精确性，并合理控制计算的时间。

简介：通过地震数据获取裂缝储藏中流体的性质并对流体类型进行识别,是地震勘探岩性反演的重要问题之一。由于地震波的速度、储层的密度等弹性参数对某些流体不具有很强的敏感性,使只依赖振幅信息进行流体识别的传统AVO方法面临困境。作为传统叠前振幅反演的一个拓展,频变AVO（FDAVO）技术进一步考虑了振幅对频率的依赖关系,将这种依赖关系与地下裂缝结构、流体填充对应起来,能带来更丰富的流体信息。利用该技术,本文提出了一种基于地震数据参数化Chapman模型的贝叶斯反演新方法（BIDCMP）,它包含两步算法,即,FDAVO反演储层的非弹性属性和贝叶斯框架下的流体识别。首先,通过匹配观测数据和模型数据,构造差函数反演裂缝储层非弹性参数。随后,在贝叶斯框架下,使用马尔科夫随机场（MRF）作为先验模型,联合多参数场识别流体。本方法在计算过程中,除综合考虑了弹性参数场、测井资料等常规信息外,还特别地加人了第一步中反演得的非弹性参数的约束,从而充分利用了流体粘性差异,最后在最大后验概率（MAP）准则下输出最佳岩性一流体识别结果。分别对合成地震记录和模拟岩性—流体剖面验证本文方法的有效性,结果证明本文方法获得的流体识别结果准确可信。

简介：标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论：采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。

简介：在直流电法有限元数值模拟中，针对常规结构化网格源点附近网格节点数值精度低的问题，并考虑计算效率和反演成像的需求，本文在二维点源问题中提出一种新的网格加密一收缩方法。其核心思想是在结构化网格计算区域中先后引入网格节点加密和网格节点收缩两个环节。通过在计算区域水平方向上均匀加密网格节点密度，降低源点奇异性的影响范围，并提升对地形的模拟能力；通过在计算区域垂向方向上将多个网格收缩到一个网格中，降低网格节点的规模，进而提高数值计算的效率。理论模型检验表明，网格加密-收缩方法能有效地解决源的奇异性问题，与均匀加密网格相比，网格节点总数降低了约80％。

简介：标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。

简介：本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。

简介：由于CRS叠加考虑了反射层的局部特征和第一菲涅耳带内的全部反射，从而更充分地利用了多次覆盖反射数据的信息。就目前的地震资料处理技术而言，它是最佳的零偏移距成像方式。本论文利用改进型的参数优化技术，得到高质量的CRS运动学参数剖面，并利用参数剖面计算出叠加孔径，实现了基于最优孔径的CRS叠加，使CRS参数的用途得到了充分利用。模型数据和实际资料的试算表明，基于最优孔径的CRS叠加的成像剖面与传统CRS叠加剖面相比，有着较高的信噪比和同相轴的连续性。

简介：含洞穴的碳酸盐岩地层具有强烈的非均质性及储集空间预测难度大的特点，利用随钻电阻率测井方法对井眼环境含洞穴的储层进行准确识别和划分，是当前研究的一个焦点问题。本文使用一种新型的高效和高精度自适应有限元方法（hp—FEM）模拟和分析了含洞穴地层随钻电阻率测井仪器响应。本文所提的hp—FEM与传统h-FEM相比，其结果具有网格自适应的特点，并且计算能够以指数速率收敛于较高的精度。数值实例使用自适应有限元方法研究地层中洞穴的大小、洞穴距离井眼的远近和仪器发射频率改变对测井响应的影响，并提供了识别含洞穴地层的方法。研究结果可以为实际测井中遇到的各种地层洞穴的准确识别和定量评价提供理论依据。

简介：Anewmulti-moderesistivityimagingsonde,withtoroidalcoilsassource,canconductthreeresistivitymeasurements:azimuthalresistivity,lateralresistivity,andbitresistivitymeasurements.Thus,theloggingtimeandcostaregreatlysaved.Thetoroidalcoilsaresimplifiedasanextendedvoltagedipoleandtheresponseequationsarederivedforahomogenousformation.Basedon3DFEM,thedepthofinvestigation(DOI),verticalresolution,circumferentialazimuthalcapacity,boreholediameter,mudresistivity,thicknessoftargetformation,andtheresistivityofthesurroundingformationandmudinvasionaresimulated.TheresultssuggestthatthethreemeasurementmodesofthenewsondearedifferentinverticalresolutionsandDOIs.Thecircumferentialdetectionabilityoftheazimuthbuttondependsonthecontrastbetweentheanomalyandformationresistivityandtheopenangleoftheanomaly.Whethertheboreholeistruncatedatthebitornothasagreatinfluenceonthesimulationresults.Theboreholeandmudinvasionaffecttheapparentresistivityinallmodes,buttheeffectsofresistivityofsurroundingformationandthicknessofthetargetformationareonlycorrectedforlateralresistivitymeasurement.