重力及其梯度异常正演的Moving-footprint大尺度模型分解方法

重力及其梯度异常正演计算效率决定了反演的可行性,也是高效构建足量、多样深度学习样本数据的基础。为了进一步提高重力及其梯度异常正演的计算速度,受航空电磁领域"Moving-footprint"快速正演技术的启发,本文在基于网格点几何格架函数空间域快速正演的基础上,提出了一种"Moving-footprint"大尺度模型分解的重力及其梯度异常正演计算方法。此方法在地下半空间内选择观测点正下方一定有效范围的子空间,该观测点异常近似为其对应子空间内长方体单元的异常和,而忽略子空间外长方体单元产生的异常;当观测点移动,其对应的子空间随之移动,从而可以将地下半空间长方体模型进行大尺度模型分解,为每一个计算点所对应的子空间进行正演计算。模型实验表明,在256x256x15个长方体模型的地下半空间内选取32x32x15的子空间进行计算,重力异常及部分梯度异常的相对平均误差小于10%,计算速度提高19倍;文中方法针对1024x1024x15个长方体模型计算时间约为32min,相比已有算法中存在的超常规计算量的瓶颈问题具有显著的计算优势。     

重力归一化总梯度法勘探浅层油气藏浅析

随着仪器采集精度的提高,分析、反演解释技术的长足进步以及油气藏重力勘探的需要,出现了重力归一化总梯度法。结合测井资料,通过重力模型正反演计算,得出高精度重力勘探浅层油气藏是可行的结论。重力归一化总梯度法自提出以来获得了不断的完善发展,并在某气田的实践中被证明是勘探浅层油气藏的有效方法。     

井中重力矢量数据联合反演研究

井中重力矢量测量可以近距离感知地下目标体,获得地质体不同方向的重力异常特征,提高地下介质的纵向分辨率.这里首先基于点元法实现了井中重力矢量正演,通过模型正演计算,分析了井中三分量重力异常响应特征,然后基于相关性搜索的黄金分割算法实现了井中重力矢量联合反演,通过模型反演分析了井中重力不同分量数据联合反演的效果,验证了反演...     

起伏地形频率域可控源电磁二维快速正反演

起伏地形频率域可控源电磁正反演研究受到众多学者的关注和重视,目前普遍实用的为未考虑地形的二维正反演算法程序。本文基于积分方程法正演和多场源、多频率、对比源反演算法,开展起伏地形条件下频率域可控源电磁二维快速正反演研究。通过引入层状参考模型、将地形与目标体整体作为异常场剖分区域,实现起伏地形可控源电磁正反演计算。针对包含地形在内的大尺度剖分区域引起的大型计算代价问题,采用快速傅里叶算法提高正反演计算效率。通过与现有正反演算法进行模型算例对比,验证了本文采用方法的可行性与有效性。     

球坐标系密度界面反演方法及在华南大陆的应用

密度界面反演方法在油气勘探、推断区域构造、获取地壳结晶基底面、莫霍面起伏形态等方面具有重要意义。现有的密度界面反演方法大多基于笛卡尔坐标系统,当涉及到大区域乃至全球尺度的密度界面反演时,地球曲率的影响将不可忽视,需考虑基于球坐标系Tesseroid模型的密度界面反演方法。然而,受计算精度和效率制约,已有的基于Tesseroid模型的密度界面反演方法并不能很好地适用于地表重力观测数据的反演计算。笔者基于前人研究,给出了一种适用于地表观测数据的球坐标系密度界面反演方法。该方法首先将常规的球坐标系高斯—勒让德重力积分公式进行简化,提高了重力正演计算效率。随后,引入并改进了前人的自适应剖分方案,提高了重力正演计算精度。在此基础上,采用Cordell迭代优化算法,得到了适用于地表观测数据的球坐标系密度界面反演方法。通过模型数据试验,对本文球坐标系密度界面反演方法进行了验证。结果表明,笔者对高斯—勒让德积分公式加以改进和引入改进的自适应剖分方案后,很好地克服了计算精度和效率对地表观测重力计算的掣肘,并且,基于球坐标系的密度界面反演结果优于基于笛卡尔坐标系的密度界面反演结果。在华南大陆实际数据试验中...     

基于起伏地层的曲面重力场快速高精度正演

Parker公式在起伏地层重力场正演方面以其简洁、快速等优点被广泛应用于地形改正以及界面反演等,但传统的Parker公式正演方法存在一定的缺陷,在起伏地层模型重力场的正演计算中,由于正演式中e的指数项数值稳定性较差,导致正演结果精度不高。在前人研究的基础上,对Parker公式进行了改进,通过增加上、下界面平均值的方法,提高正演的数值稳定性,改进后的正演精度得到大幅提升。在此基础上,提出了一种曲面观测算法,将观测面由传统的水平面推广到任意起伏面,实现了上下界面起伏地层在起伏观测面上的快速高精度重力正演计算和带地形的密度界面起伏模型的快速高精度重力正演模拟。为了保障数值精度,采用了Gauss-FFT算法,模型实验均取得很好效果。     

重磁位场综合优化解释方法研究及其应用

1　基于改进的模拟退火法的重磁位场反演方法在深入分析重磁三维位场场源位置函数特点的基础上,利用改进的模拟退火法实现了重磁三维位场反演,并研制开发了相应的软件。采取纵横向不等间距对地下场源空间进行剖分。这种剖分方式可起到沿深度加权的作用,极大地提高了解的稳定性。利用快速正演计算和变步长搜索方法,极大地压缩了模拟退火法与局部最优化结合中正、反演的计算量,计算表明两者的结合明显优于各自独立反演,极大地提高了计算速度。在对场源施加一定约束条件的基础上,采用基于位场特点的加权的目标函数,取得了较为精确的反演结果。…     

有限单元法三维电阻率最小二乘反演中存在问题的研究

针对有限单元法在三维电阻率最小二乘反演中关于深部单元的反演能力、反演稳定性和速度问题进行了讨论与改进.首先构造了全新的"体积因子"作为先验信息加入到电阻率反演中,使地质体赋存区域的电阻率变化更符合实际情况,提高了深度较大的地质体的分辨能力.反演计算中对实测电阻率数据和模型正演数据分别取对数,并改进迭代修正量的计算方法,将以前的加减修正量改为与修正量指数之积,这样,不但可以改进电阻率模型出现负值的缺陷,提高反演计算的稳定性,而且还可减少迭代次数;针对反演时间过长的缺陷,将三维全域反演改为三维局域反演,使每次迭代的时间由4 h减少到目前的15 min;对模型的反演计算表明,以上改进是非常有效的.     

基于分形先验信息的非线性反演

基于分形先验信息的非线性反演方法能综合利用测井数据和地震信息,在贝叶斯框架下,通过分形高斯噪音算法得到基于分形理论的先验信息,然后根据地震资料构建似然函数,最终利用基于快速模拟退火算法(Very Fast Simulated Annealing,VFSA)改进的粒子群优化(Particle Swarm Optimization,PSO)算法(VFSA-PSO)实现后验概率密度的抽样。与确定性反演结果相比,该方法能够有效地融合测井资料中的高频信息,提高反演结果的分辨率,并且目标函数的建立融合了确定性反演中的低频约束,从而得到宽频带的反演结果。数值模拟试验表明:基于分形先验的非线性反演结果与理论模型吻合较好,实际资料的应用效果也证明了该反演方法的有效性。     

核磁共振地下水勘查的模拟退火反演

简要地介绍了NMR找水的理论方法,分析了核磁共振找水对于不同含水层的正演计算方法,并用模拟退火反演方法对一维含水模型进行反演,给出了三个不同含水层模型的模拟退火反演效果图。经对比验证,该方法在技术上是可行的,在理论上给出了对NMR反演进行改进的可行性方法。     

一维层状介质大地电磁模拟退火反演法

大地电磁模拟退火反演法是一种最优化的非线性反演方法,与传统的线性反演方法相比该方法具有：（１）不依赖于初始模型的选择;（２）能寻找全局最小点而不陷入局部极小;（３）在反演过程中不用计算雅可比偏数矩阵等优点;通过对各种类型的大地电磁测深理论曲线试算,结果表明模拟退火法能准确地自动反演地电参数（地层电阻率,厚度）最后对实际资料进行了处理,取得了较好的效果。     

基于模拟退火粒子群优化算法的裂缝型储层各向异性参数地震反演

发育垂直定向排列裂缝的地下岩石可等效为具有水平对称轴的横向各向同性（horizontal transverse isotropic,HTI）介质。针对HTI介质模型,本文研究了裂缝型储层的各向异性参数地震振幅随方位角变化（amplitude variations with azimuth,AVAZ）的反演方法。首先,在地震AVAZ反演流程中,提出采用模拟退火粒子群优化算法实现裂缝型储层各向异性参数反演。之后,通过理论模型测试,验证了基于模拟退火粒子群优化算法的地震AVAZ反演的有效性。最后,将反演方法应用于四川盆地龙马溪组页岩气储层实际方位地震数据;在反演之前先利用傅里叶级数方法估计裂缝方位并对实际数据进行方位校正,以提供更准确的输入数据;通过计算得到的P波、S波各向异性参数可用于评价裂缝发育程度。反演结果表明,研究区域构造顶部裂缝较发育,与地质基本理论一致,验证了反演方法的合理性。     

重力异常AlexNet深度神经网络反演

针对传统反演方法存在的初始模型依赖、计算时间较长等问题,提出了一种新的基于AlexNet深度神经网络的重力异常反演方法。该方法首先借鉴经典的深度神经网络AlexNet设计了一种用于重力异常反演的Alex反演网络(AlexInvNet),接着设计大量密度异常体模型并通过正演计算得到带标签的数据集,然后用该数据集训练AlexInvNet网络,最后将重力异常数据输入训练好的AlexInvNet网络直接得到反演结果。理论模型反演结果表明,该方法相较于全连接网络深度学习反演方法,能够更好地反演出异常体的位置和密度,具有较好的泛化能力和抗噪声能力。实测数据反演结果表明,该方法能有效解决重力异常反演问题。     

勘查油气田化探新方法——土壤相态烃测量法

烃气测量是油气化探最主要的方法。目前应用的烃气测量方法有两类:一类是测量土壤中游离烃气如K—V指纹法、吸附丝法、物上气法及直接取土壤气测量法等。此类方法由于需特殊设备及壤气测量方法本身固有的缺陷,目前还未得到普遍应用,另一类为土壤吸附烃气法即国内外目前普遍使用的酸解烃气法、此法在多数油气田上取得好效果,但也有相当多的地区如黄土、荒漠及水网发育等特殊景观区效果不好或不显著。本文介绍一种新的烃气测量法──土壤相态烃测量法。这是用一种特殊的技术,只选择性提取由油气藏渗漏过程中保存在运移通道围岩及近地表土壤各种次生矿物结构中以"穴居"态存在的轻烃组分。"穴居"态烃为累积性的指标,不受外界温度、修水等的影响,保存较牢固,测定结果稳定。相态烃法避免了酸解烃法存在的生物成因烃、矿物同生烃及碳酸盐的干扰影响,提高了测量指标的信噪比,突出了应有的异常。本方法已在黄土厚覆盖区两个已知油气田和一个风成沙覆盖的断裂发育油田上试验获得成功,结果表明相态烃法优于酸解烃法。     

长江源生态环境问题与可持续发展对策

长江源自然环境严酷,以草原退化、森林萎缩、水土流失、气候暖旱化、冰川退缩、雪灾频繁为代表的生态环境退化正在进行。宜采用保护、恢复和改善林草植被系统来对其进行综合整治。     

一种大坝渗透系数分区反演新方法研究

基于人工神经网络的非线性映射特征,在渗流有限元计算的基础上,结合水头和渗流量等实测资料提出了大坝渗透系数的反演方法。为了克服经典神经网络存在的缺陷,提出了模拟退火的交替迭代算法神经网络新方法。在相同的初始条件下,用该新方法和经典网络进行了比较,得出前者的优越性和有效性。同时将该方法用于大坝的渗流反分析,利用反演出的渗透系数进行渗流场计算。数值计算结果表明,这种方法对大坝渗透系数反演问题具有较高的识别精度,反演结果可靠,可以用于实际工程。     

波阻抗约束反演的几种方法

介绍三种理论与实际效果比较好的波阻抗反演方法,即宽带约束反演(BCI)、快速模拟退火算法(FSA)和正则化方法.BCI是一种线性迭代方法,FSA是一种具有全局寻优特点的随机搜索方法,正则化方法是通过对目标函数施加先验约束,来解决反演的病态问题--多解性和不稳定性.最后,使用正则化方法并结合快速模拟退火算法,进行了理论模型试算和实际资料处理,实用效果比较好.     

改进快速模拟退火方法进行AVO岩性参数反演

以广义波乐兹曼的吉布斯（Ｂｏｌｔｚｍａｎｎ－Ｇｉｂｂｓ）统计理论为基础,将数认（ｎｕｍｃｒａｌ）和依赖于温度的似柯西（Ｇａｕｃｈｙ）分布相结合建立一种新的扰动模型,改进接受概率,提出了一种非常快速模拟退火法（ｖｅｒｙ ｆａｓｔ ｓｉｍｕｌａｔｅｄ ａｎｎｅａｌｉｎｇ）。通过多参数的振幅随炮检距变化（ＡＶＯ）的合成数据和实际投资试算结果表明,该方法能在高温下进行大范围内的搜寻,解的稳定性好和抗随机干     

波阻抗反演在南海北部神狐海域天然气水合物勘探中的应用

针对神狐海域的地质构造和天然气水合物的赋存特征,以重点测线三维地震数据为基础,分析讨论了基于宽带约束的模拟退火波阻抗反演方法、流程和关键技术问题,定量获得了含天然气水合物沉积物的波阻抗特征。结果表明:基于宽带约束的模拟退火波阻抗反演数据具有较高的有效垂向分辨率和较好的横向连续性;神狐海域高波阻抗异常反映了含天然气水合物沉积层,而不连续异常低波阻抗层是水合物层之下游离气的表现,这与钻探结果吻合。由此可见,基于宽带约束的模拟退火波阻抗反演可为天然气水合物层识别和预测、勘探目标圈定、钻探井位选择提供重要依据。