简介：无网格法形函数构造不依赖预定义的单元，具有计算精度高、处理复杂模型便利等优点。本文介绍了无单元Galerkin法（EFGM）、点插值法（PIM）与径向基点插值法（RPIM）三种全域弱式无网格法的近似原理及特点；以二维泊松方程为例研究了支持域无量纲尺寸、场节点与背景网格设置对无网格法计算精度的影响。将RPIM与EFGM应用于频率域线源二维正演，给出了RPIM形状参数的推荐值；分析了均匀介质模型大地电磁（MT）二维正演无网格法边界条件直接加载与罚函数法加载的精度差异，结合PIM与RPIM边界条件加载便利及EFGM计算复杂模型精度高的优势，提出了EFG—PIM及EFG．RPIM耦合算法，数值计算结果验证了耦合算法的有效性。研究发现：无网格法及其耦合方法适用于电磁法数值模拟；支持域无量纲尺寸取1．0时无网格法精度与效率高，场节点与背景网格重合时计算效果佳；泊松方程求解PIM及RPIM精度较EFGM低，计算均匀介质MT响应精度较EFGM高；RPIM改善了PIM计算涉及的奇异性问题，对应支持域无量纲尺寸选择空间大。

简介：通过引入岩石骨架参数表达式到Gassmann近似流体方程中，实现了基于该方程的地震孔隙度反演。但是岩石骨架参数表达式有很多种，使得地震孔隙度反演公式也各不相同，实际应用不方便，反演结果也难以评价。作者在研究现有常用岩石骨架参数表达式如Esheby－Walsh，Pride，Geertsma，Nur，Keys－Xu以及Krief等的基础上，提出了含有调节参数的岩石骨架模型统一表示式，发展了基于Gassmann方程的地震孔隙度反演方法。该方法应用范围宽，参数调节方便、灵活。为验证所提公式和方法的实用性，作者们结合ZJ地区钻井岩心的岩石物理样品测试数据和测井数据，进行了地震孔隙度反演。由于该地区油气储层的存在与中等程度孔隙度的分布有关，因此地震孔隙度反演对预测和识别油气层、干层和含水层有重要作用。反演结果与工区内孔隙度测井瞌线吻合度高，说明本文所提公式与方法是有效的，可靠的。

简介：地层品质因子Q反映了地层对地震波吸收衰减的强弱,是油气检测的重要指示因子,准确提取Q值十分重要。本文根据地震子波能量衰减特征,提出基于子波衰减前后能量比的Q值提取方法。该方法是综合利用多点信息估算子波能量,规避了常规Q值方法（如质心频移法）对地震子波谱的假设,适用于任意波谱形状,但对资料的保真度要求较高。VSP正演模型测试发现,在非固有衰减（如几何扩散、反射透射损失等）得到准确补偿的前提下,能比法能较好地反演地层Q值。实际零偏VSP资料应用表明,能比法提取的Q值趋势与谱比法基本一致,具有一定的可信度。

简介：震源定位是微地震监测关键技术之一。本文提出用于微地震定位的弹性波和加权弹性波（WEW）干涉成像方法。该方法在保证定位精度的同时，还可避免震源假象。通过各向同性水平层介质状模型的数值试验，初步表明该方法可适应低信噪比微震信号、速度随机扰动、较稀疏的检波器分布等情况，并在速度模型存在一定的系统误差时也仍保持较高的定位精度。由于干涉成像方法不需要进行初至拾取，定位效率相对传统走时方法也得到了提高。采用二维断层模型试算Nnumericalresultsofusingatwo-dimensionalfaultmodel，表明方法还能实现多震源定位，且比逆时成像有更高的定位精度。

简介：本文介绍了在B隧道勘查中应用高分辨电磁成像技术探测隐伏构造的一个成功实例.数据采集采用高频电磁成像系统(STRATAGEMEH4,频率范围从1Hz到90kHz).野外采集正交的电磁场分量,从地磁成像剖面中提取目标体有关电磁学信息.为获取高质量的采集数据,野外采用包含天然场源和全张量可控源的混合场源.B隧道坐落在位于中国中部的湖北省的西部,隧道埋深不到200m,但区域地质作用导致它的地质机构十分复杂.第一次勘查过程中,勘探人员误把一个脱落体的露头当成基岩面,施工过程中出现了冒顶.第二次勘探时采用高分辨电磁法和折射地震法,这次勘探找到了隐伏的基岩面和一个隐伏断裂.勘探结果与后来隧道挖掘揭示的构造吻合.

简介：常规的重力梯度法对地质体位置划分很模糊,当深度越大,地质体越小时,对Vzz和Vzx求导的误差也越大,导致结果与实际偏差很大。本文提出状态判别因子对角点位置进行优化确定,比常规高阶重力导数法更精确,使重力梯度法在对小地质体和断裂上有更好的分辨率,理论模型试算和实际资料处理效果很好,为找矿或划分地下地质体产状提供了更有利的依据。

简介：具有较高自然频率的高灵敏度检波器在采集数据时可以通过压制低频信号来相对提高高频能量，但这也造成了其低频响应差的问题，如果将高灵敏度数据与常规数据进行优势组合，就可以达到拓宽频带的目的。为此，本文提出分频段匹配滤波的方法，即在保持常规数据低频优势的前提下，对其高频端进行匹配滤波，实现不同频带范围内的优势互补，从而改善地震记录。通过引入不同主频的雷克子波模拟得到具有常规数据和高灵敏度数据特点的理论模型，论证了分频匹配滤波方法的可行性。在对野外单炮地震记录处理中发现，分频匹配滤波方法拓宽了地震记录的有效频带宽度，提高了地震记录的分辨率。

简介：本文利用取自南海神狐海域钻井岩芯进行了核磁共振测量,研究了岩石中所含水合物对地层渗透率的影响,该岩芯样品为一＂孔隙充填型＂水合物样品。核磁共振测量结果表明岩石中所含水合物对地层渗透率影响很大,随着含水合物饱和度的增加,地层渗透率迅速降低。通过对测量结果的进一步分析可知＂孔隙充填型＂水合物对地层渗透率的影响-存在多方面的因素,包括水合物对地层孔隙度、孔隙迂曲度的影响以及对水合物对孔隙喉道的堵塞等。最后本文利用Masuda下降模型建立了适用于南海神狐海域的相对渗透率与含水合物饱和度的经验公式,下降指数N=7.9718。

简介：随着煤矿的深部开采，多层积水采空区探测成为煤田水文物探工作的新内容。中心回线瞬变电磁（简称TEM）法因场结构与地层的耦合关系，对低阻层探测有利。但低阻层的屏蔽作用不仅使得探测同样的深度需要更长的观测时间，而且还会减弱下伏地层的异常响应。本文通过直接时域数值模拟和水平分层大地的模型正演，估算了探测目的层所需要的时间长度，根据噪声对观测数据造成的影响，给出了多层积水采空区可分辨的标准。山西大同达子沟煤矿水文勘探中获得的TEM实测曲线，表现了对多层积水采空区的探测能力。在实测曲线的定量反演解释中，利用电测井资料作为初始参数进行定量反演，约束了煤系薄层的等值性。所推断的三层和两层积水采空区，为钻孔所证实。研究结果表明，当观测时间有足够的长度、下伏地层的异常显示大于观测误差，中心回线TEM法探测多层采空区积水是可行性的。

简介：Theattenuationfactororqualityfactor(Q-factororQ)hasbeenusedtomeasuretheenergyattenuationofseismicwavespropagatinginundergroundmedia.ManymethodsareusedtoestimatetheQ-factor.WeproposeamethodtocalculatetheQ-factorbasedontheprestackQ-factorinversionandthegeneralizedS-transform.TheproposedmethodspecifiesastandardprimarywaveletandcalculatesthecumulativeQ-factors;then,itfindstheinterlaminarQ-factorsusingtherelationbetweenQandoffset(QVO)andtheDixformula.TheproposedmethodisalternativetomethodsthatcalculateinterlaminarQ-factorsafterhorizonpicking.Becausethefrequencyspectrumofeachhorizoncanbeextractedcontinuouslyona2Dtime–frequencyspectrum,themethodiscalledthecontinuousspectralratioslope(CSRS)method.ComparedwiththeotherQ-inversionmethods,themethodoffersnearlyeffortlesscomputationsandstability,andhasmathematicalandphysicalsignificance.WeusenumericalmodelingtoverifythefeasibilityofthemethodandapplyittorealdatafromanoilfieldinAhdeb,Iraq.TheresultssuggestthattheresolutionandspatialstabilityoftheQ-profileareoptimalandcontainabundantinterlaminarinformationthatisextremelyhelpfulinmakinglithologyandfluidpredictions.

简介：为了有效处理人工源的影响,本文开发了带源的CSAMT二维正反演算法,可用于全区（近区、过渡区和远区）资料的反演。引入正则化因子完成磁法二维反演,并且将模型参数调整为磁化率的对数,保证反演过程中磁化率始终为正值。本文基于交叉梯度原理,将CSAMT和磁法进行联合反演,通过搜索交叉梯度项权重的方法,避免了不同异常源引起的两种异常相互干扰的问题。理论模型算例表明基于交叉梯度的联合反演方法优于单独反演。本文开发的带源CSAMT二维正反演算法,有效处理了人工源的影响,保证了最终联合反演算法的可靠性。

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

简介：由于孔隙型碳酸盐岩储层的复杂孔隙结构和强非均质性,岩石导电效率与含水孔隙度之间理想的线性关系相对复杂或并不存在,故前人基于岩石导电效率理论建立含水饱和度的计算公式并不完全适用于碳酸盐岩储层。基于岩石导电效率理论,推导了岩石导电效率的计算公式,阐明了岩石导电效率与含水孔隙度之间线性关系相对复杂或并不存在的根本原因,发现了岩石导电效率与电阻率之间的幂函数关系,分析了所建立含水饱和度计算公式误差的主控因素,得出了岩石导电效率的计算精度是该方法是否能推广应用的关键。与Archie公式相比,基于岩石导电效率理论建立的含水饱和度计算公式能更准确地计算研究靶区孔隙型碳酸盐岩储层的含水饱和度。在伊拉克某油田和印尼某气田3口井碳酸盐岩储层中的应用,表明当计算岩石导电效率的相对误差不大于0.1时,所计算储层含水饱和度的绝对误差不大于0.1,基本满足孔隙型碳酸盐岩储层精细评价的需求。

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

简介：本文提出了一种基于Kirchhoff积分偏移和逆时偏移的联合速度分析及成像方法,采用剩余曲率分析及层剥离策略进行偏移速度建模。本文方法改善了Kirchhoff积分偏移在复杂构造时计算精度不高和逆时偏移计算效率慢的缺点,兼有计算效率高、成像精度高的优点,并将其应用在反射波法隧道超前预报中。通过模型试算,发现隧道中使用逆时偏移的成像结果在多方面优于Kirchhoff积分偏移的成像结果;通过对实测数据的处理,验证本方法计算效率较高,建立的速度模型合理,成像剖面清晰,结合地质调绘资料可对隧道开挖前方地质构造做出较准确预报。

简介：受沉积环境影响,海底地层层理发育,各向异性介质更接近于真实海洋地质条件。对海洋各向异性研究可更好认识海底构造特征,有效地进行海底资源勘查。本文从各向异性麦克斯韦电磁方程出发,根据电磁场的无源特征引入标量位函数;利用电场和磁场的连续性分别向海底深部和海水中延拓,并将延拓后的电磁位函数在海底耦合到发射源上,从而实现海底电场和磁场的递推求解。我们首先探索如何利用海洋电阻率法识别和求解海底电各向异性特征。研究发现海底各向异性可从视电阻率测深曲线和海底视电阻率极性图进行求解。进而,我们通过海底各向异性地层中高阻体(油气藏)模型进行正演模拟,发现海洋电法视电阻率在浅水区对地下高阻薄层有明显的异常反应。相比海洋可控源电磁法,海洋直流电法不受空气波影响,在浅海区油气资源勘查有着较好的优势,而前人大多基于各向同性模型进行研究,本文实现一维海洋直流电阻率法各向异性正演模拟,算法计算精度高,能很好的为二维、三维正演模拟提供理论参考。更多还原

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

简介：MostcurrentprestackAVAjointinversionmethodsarebasedontheexactZoeppritzequationanditsvariousapproximations.However,theseequationsonlyreflecttherelationbetweenreflectioncoefficients,incidenceangles,andelasticparametersoneithersideoftheinterface,whichmeansthatwave-propagationeffects,suchassphericalspreading,attenuation,transmissionloss,multiples,andeventmismatchingofP-andS-waves,arenotconsideredandcannotaccuratelydescribethetruepropagationcharacteristicsofseismicwaves.ConventionalAVAinversionmethodsrequirethatthesewave-propagationeffectshavebeenfullycorrectedorattenuatedbeforeinversionbuttheserequirementscanhardlybesatisfiedinpractice.Usingaone-dimensional(1D)earthmodel,thereflectivitymethodcansimulatethefullwavefieldresponseofseismicwaves.Therefore,weproposeanonlinearmulticomponentprestackAVAjointinversionmethodbasedonthevectorizedreflectivitymethod,whichusesafastnondominatedsortinggeneticalgorithm(NSGAII)tooptimizethenonlinearmultiobjectivefunctiontoestimatemultipleparameters,suchasP-wavevelocity,S-wavevelocity,anddensity.Thisapproachisrobustbecauseitcansimultaneouslycopewithmorethanoneobjectivefunctionwithoutintroducingweightcoefficients.Modeltestsprovetheeffectivenessoftheproposedinversionmethod.Basedontheinversionresults,wefindthatthenonlinearprestackAVAjointinversionusingthereflectivitymethodyieldsmoreaccurateinversionresultsthantheinversionbyusingtheexactZoeppritzequationwhenthewave-propagationeffectsoftransmissionlossandinternalmultiplesarenotcompletelycorrected.

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

简介：Tilt—depth法可以用于快速反演磁源的上顶埋深，但其无法反演磁源底部深度，同时该方法的反演解过于单一。针对于此，本文在有限厚度台阶磁异常导数公式基础上，推导出了可以同时反演地质体上、下界面埋深的改进型Tilt—depth法计算公式，并采用在Tilt梯度图上选取多特征点进行联合反演的计算模式来提高反演解的可靠性。二维及三维组合模型试验均证实了改进型Tilt—depth法可以有效地反演出磁性体的上顶与下底深度，且上顶埋深的反演精度明显高于常规方法的。最后将改进型Tilt—depth法用于松辽盆地长岭断陷航磁数据反演，其上顶埋深的反演结果与钻孔钻遇火山岩的深度基本一致，证实改进型Tilt—depth法相对于常规方法的反演结果更加准确。