共查询到20条相似文献，搜索用时 171 毫秒
1.
研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性. 相似文献
2.
在频率域弹性波有限元正演方程的基础上,依据匹配函数（也就是观测数据和正演数据残差的二次范数）最小的准则,用矩阵压缩存储与LU分解技术来存储和求解频率域正演方程中的大型稀疏复系数矩阵、用可调阻尼因子的Levenberg Marquard方法求解反演方程组,直接求取地下介质的弹性波速度,导出了频率域弹性波有限元最小二乘反演算法. 为了利用地下地质体的分布规律,减少反演所求的未知数个数,本文又提出了规则地质块体建模方法引入到反演中来. 经数值模型验证,在噪声干扰很大（噪声达到50髎）或初始模型与真实模型相差很大的情况下,反演也能取得很满意的效果,证明本方法具有很好的抗噪性与“强壮性”. 相似文献
3.
从各向异性弹性波的有限元正演方程出发 ,导出了反问题中时间域雅可比矩阵求解的计算公式 .它具有与时间域有限元正演方程相同的表达形式 ,故可通过有限元正演计算来获得雅可比矩阵 .研究了有限元正演算法的效率和精度、吸收边界条件等方面的问题 ,以提高反演系统的效率和精度 .在此基础上 ,实现了叠前全波场各向异性弹性参数反演 .计算表明 ,在初始模型偏离真实模型较大的情况下 ,层状模型和横向不均匀模型的反演结果均能准确地收敛到真实模型上 . 相似文献
4.
5.
推导出频率域有限元声波正演方程,为了消除边界反射,将Clayton-Engquist旁轴波动方程吸收边界条件引入频率域,并对有限元刚度矩阵和质量矩阵进行压缩存储,利用广义共轭梯度法求解有限元方程获得正演解.在此基础上,推导出在某一频率下波场数据残差δU与单元物性参数修改量δλ之间关系的Jacobi矩阵,反演方法允许利用地面二维炮集全波场资料与给出初始模型参数的正演值的差值δU,迭代求得δλ.由于计算机内存的限制,方法计算不允许有过多数目的未知数个数,因此还提出了对同一介质物性单元的Jacobi矩阵元素进行压缩组装的措施,从而使反演的未知量个数减少,结合采用共轭梯度迭代法,使得只需利用有效波频段的少数一些频率即可进行迭代反演.正演和反演理论模型的数值模拟结果表明方法是有效的. 相似文献
6.
7.
8.
《应用地球物理》2019,(4)
阵列侧向测井具有分辨率高、探测深度多等优势,但受围岩、侵入等因素影响会导致油层漏失和错判,在阵列侧向测井资料的反演时需要考虑这些因素,以使反演结果反映原状地层真实情况。阵列侧向测井资料的反演遇到的问题是:反演参数多,导致雅克比矩阵计算困难;正演速度慢,反演初值选取困难。为此,本文研发了一种阵列侧向测井资料快速处理方法。首先,分析层厚、围岩及侵入等对阵列侧向测井响应的主控因素;然后,基于深度开窗技术,将全井段的阵列侧向测井资料的处理问题转化为一系列连续深度窗口内地层的2D反演;针对深度窗口内每个地层,结合1D等效快速正演算法,实现地层径向电阻率剖面的快速提取;最后,以1D反演结果为初值,进一步消除围岩和层厚对视电阻率响应的影响。数值模拟结果表明:阵列侧向测井响应受侵入影响严重,而层厚和围岩影响次之;开窗技术极大降低了反演参数数量,提高反演的速度。实际井资料处理结果表明,2D反演实现了对地层纵横向电阻率分布的快速重构,提升了储层饱和度计算的精度。 相似文献
9.
在大斜度井/水平井环境下,阵列侧向测井受钻井液侵入、地层倾角和各向异性等多种因素影响,导致测井响应复杂,需借助反演手段提取地层真实电阻率.然而,阵列侧向测井三维正演效率低,难以满足测井资料快速反演和油气藏快速评价的需求.为此,本文基于深度学习并联合混合多种群遗传(MPGA)与列文伯格马奎特(LM)算法建立了一种快速反演方法.首先从地层参数敏感性出发,基于严格的三维有限元正演算法,依次确定侵入、各向异性和地层倾角的敏感性大小;其次,引入深度学习和模型可视化技术,实现斜井各向异性地层阵列侧向测井响应的快速正演;最后,基于数据集分解技术和混合MPGA-LM算法,实现斜井各向异性地层电阻率剖面快速精确重构.数值模拟结果表明:斜井各向异性地层中,阵列侧向测井响应对侵入深度敏感性最高,各向异性和地层倾角次之;相较于反向传播神经网络(BPNN),二维卷积神经网络(2D-CNN)能够实现阵列侧向测井响应的快速精确计算,计算一个测井点仅需0.36 ms,精度可达99％左右;基于三层反演模型的MPGA-LM混合算法稳定性强,电阻率参数反演精度高的优点,可用于阵列侧向测井资料的快速处理. 相似文献
10.
三维陆地可控源电磁法有限元快速正演的主要瓶颈在于多频率大型稀疏方程组求解问题.本文引入一种基于模型降阶的Krylov子空间投影算法,推导了有限元刚度矩阵的模型降阶形式,构建了频率域传递函数;采用标准正交向量序列,构建一个远远小于有限元刚度矩阵维度的矩阵,该矩阵与频率无关,通过一次模型降阶即可实现多频点有限元方程快速求解.采用基于电场的变分方程,加入散度校正条件,以消除伪解;引入伪δ函数,消除了源点的奇异性,可适用于复杂背景模型三维有限元数值模拟,并为多源的求解奠定了基础;以层状介质模型解析解为标准,通过和基于Pardiso直接求解器的有限元算法（3DFEM）进行比较,模型降阶法计算时间小于前者的1/10,平均相对误差在1.72%,在满足精度要求下,实现了高效率三维有限元数值求解;分别设计了横向高低阻模型和纵向高低阻模型,分析了从近区到远区电场和卡尼亚视电阻率的变化规律,假极值的表现特征,阴影效应的影响等,从而也验证了该算法的正确性.最后,建立了一个地层陷落柱模型,通过模型降阶有限元正演模拟,发现视电阻率断面图在陷落柱上方出现"凹陷",与模型设计吻合,表明该算法对复杂地层模拟具有同样的适用性. 相似文献
11.
电阻率测井反演问题解的局部收敛性和非唯一性增加了测井解释的难度，影响了实际的应用，把测量原理迥异的感应型和电流型测井结果联合起来反演，是减小这些影响的有效方法，利用经测井理论和反演技术，建立了适用感应－横向组合测井响应的联合演算法，研究了单一的横向测井反演的感应－横向联合反演的局部收敛性、非唯一性及解的偏差，结果表明，联合反演更有利于同时确定地层电阻率、侵入半径和侵入带电阻率，当初值地较大范围内变化时，反演结果均能在真值或真值附近收敛，改善了反演结果的局部收敛性，增加了反演结果的可信度，避免了差异较大乃至相互矛盾的解，现场实例证实了这一结论。 相似文献
12.
13.
本文将应用2.5维数值模式匹配算法研究建立柱状横向同性地层中偏心条件下磁流源并矢Green函数的高效算法,并利用该算法高效计算多分量阵列感应井眼校正库.首先,通过Fourier展开技术将偏心条件下磁流源并矢Green函数的数值模拟转化为一系列的轴对称问题;然后利用模式匹配算法求解轴对称问题,得到柱状介质中偏心条件下磁流源并矢Green函数的半解析解;最后,结合中国石油集团测井有限公司开发的三维感应测井仪器(TDIT)的结构参数,推导出复合线圈系视电导率张量的计算公式.在此基础上,根据井场实际地层资料归纳出理论地层模型相应参数的变化范围,高效高精度建立一套直接应用于工业生产的多分量阵列感应井眼校正库,并利用井眼校正库详细地分析考察各模型参数对TDIT响应的影响. 相似文献
14.
原状地层水电阻率是重要的储层参数,也是进行精细储层评价的基础.基于泥浆侵入数值模拟与侵入过程中井周岩石物理特征分析,确定了薄膜电位的产生位置,针对储层高、低侵等不同侵入特征,提出了可适用于包括存在"低阻环"等不同侵入特征时储层电阻率分布的"五参数"地层模型,基于几何因子理论与有限元方法,建立了阵列感应与自然电位测井联合反演算法,实现了地层电阻率参数反演,重构了地层径向电阻率剖面,进而精确求取了地层水电阻率.通过对实际资料处理表明:反演算法稳定可靠,"五参数"模型能很好地表征储层电阻率分布形态,重构储层电阻率剖面,确定薄膜电位产生位置;基于阵列感应与自然电位的联合反演,能精确计算原状地层水电阻率,为储层评价与流体性质识别提供依据. 相似文献
15.
为研究地层电性变化时不同区域对电磁波测井响应的贡献分布,从响应信号对地层参数求偏导的角度,给出一种新的电阻率敏感性函数定义,引入模式匹配法对纵向成层、径向非均匀介质敏感性分布进行快速模拟;通过对敏感性函数纵、横向积分给出了单发双收线圈系纵、横向探测范围,研究了井眼、频率、背景地层电阻率等对探测特性的影响.结果表明:敏感性函数能够定量表征响应对地层纵向与径向微观与宏观敏感性,幅度比与相位差敏感性分布形态类似,幅度比较广较深,而相位差分辨率高,敏感范围小;敏感性函数进行径向积分后可表征仪器的探测深度,与伪几何因子对比达到了同样的效果;背景地层电阻率在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,远小于地层电阻率的影响. 相似文献
16.
过套管电阻率测井是通过测量套管壁电势实现测量地层的视电阻率,基于传输线方程理论,针对层状地层,给出了套管壁电势、电流对地层横向电阻导数的微分方程(称Jacobi矩阵微分方程)及边界条件;利用Jacobi矩阵微分方程边值问题导出了过套管电阻率测井反演地层参数的Jacobi矩阵系数的解析表示,利用Marquardt方法实现了过套管测井的地层电阻率反演;通过计算对Jacobi矩阵的特性进行了探讨,并获得了较快的计算速度(因为Jacobi矩阵是用解析解表示的),反演结果与地层模型取得了较好的逼近.本文实现了过套管电阻率测井地层参数的Jacobi系数矩阵的快速计算及地层电阻率反演,为进一步开展电阻率测井数据处理提供了理论依据和快速反演算法. 相似文献
17.
One of the best indicators of the potential erosion risks is the rainfall–runoff erosivity factor (R) of the revised universal soil loss equation (RUSLE). Frequently, however, there is not enough data available to compute the R value, and other parameters, such as the modified Fournier index (Fmod), are used instead. But RUSLE is less effective if only the alternative procedures exist. One of the major discrepancies between R and the alternative parameters is time resolution: individual storms are used to calculate R while monthly averages over the year are used to calculate Fmod.
In this study, a multiple linear regression (r2=0.89) involving monthly EI30, monthly rainfall for days with ≥10.0 mm and monthly number of days with rainfall ≥10.0 mm, for the Algarve region, is presented. Twenty-seven years of monthly rainfall erosivity values were computed for the 32 standard daily-read raingauge stations of the Algarve region. 相似文献
18.
霍山地震区位于大别造山带北缘华北板块与扬子板块接触带上,是大别造山带及周边地震活动最频繁、最集中的地区.83个大地电磁测点组成的大地电磁三维阵列覆盖了整个霍山地震区.用多重网格法、印模迭代重构法和非线性共轭梯度法对阵列数据进行三维带地形反演,获得了地震区深部三维电性结构.电性结构显示,北大别、北淮阳区的中上地壳为电阻率1000Ωm以上的高阻区,中下地壳为电阻率数十欧姆米的相对低阻区;六安盆地电阻率整体较低,中地壳存在显著的电阻率为几欧姆米的壳内高导层.北西向的晓天—磨子潭断裂分隔了北大别高阻层和北淮阳高阻层,在浅部向NE倾,深部向SW倾;北东向的落儿岭—土地岭断裂切穿北大别上地壳高阻层.小震双差定位结果表明,地震主要发生在NE向延伸的落儿岭—土地岭断裂附近的北大别、北淮阳中上地壳的高阻区,并集中于NW向的晓天—磨子潭断裂运动所造成的构造薄弱带中;2014年M S4.3霍山地震震源深度较深,位于北大别高阻区内部的电性梯度较大的区域.综合上述结果我们认为,霍山地震区的主要发震断裂为落儿岭—土地岭断裂,断裂的运动变形充分利用了晓天—磨子潭断裂早先活动所形成的构造薄弱带,断裂下方壳源高导体中的流体沿断层传播使断层强度弱化,使得这些薄弱带区易于发生小地震.由于北大别、北淮阳构造区显著高阻层的存在,我们认为霍山地震区存在发生6级以上中强震的深部孕震环境. 相似文献
19.
随钻电磁波传播方位电阻率仪器在钻井过程中可以提供地层边界的方位及距离信息,因此在地质导向应用中发挥着重要作用,是提高油气资源开采率的重要手段.该技术是国外近十年来发展起来的前沿技术,在国内尚属空白.本文首次构思并实现了一种应用"交联天线"的随钻电磁波传播方位电阻率仪器,并从天线结构,测量原理,数据处理以及地层电阻率成像及解释等几个方面详细阐述.交联天线由极化方向正交的线圈串绕组成,同时具有方位探测和地层背景电阻率测量的功能.通过分析交联天线的电压在仪器旋转过程中的变化规律,可以计算出仪器所在地层的电阻率,判断地层边界的方位以及估算地层边界相距仪器的距离.该随钻电磁波传播方位电阻率仪器还包括用于常规电磁波传播电阻率测试的多频多测距补偿天线结构.结合地层电阻率测量及其在方位上的相对变化,可以实现对地层电阻率的全方位成像.该仪器在国内油田进行了多次实井测试,测试结果证明仪器能够在水平井地质导向中准确提供仪器所在地层的电阻率,以及地层边界方位和距离信息.已有实井测试结果表明仪器在油层(100Ωm)和泥岩(5Ωm)中的地层边界探测深度分别为2.2m和1.6m,平均误差在0.2m以内. 相似文献
20.
二维核磁共振能从"弛豫-扩散"两个维度上展现流体性质,在识别稠油储层方面具有理论优势,是当前核磁共振测井技术的研究热点.本文深入研究二维核磁共振测井原理,系统分析CPMG-DE脉冲序列测量扩散系数与弛豫时间的方法,结合核磁共振二维谱数理模型,提出一种基于SVD和BRD的正则化反演算法.该算法通过SVD压缩数据,采用带非负约束的Tikhonov正则化方法求解流体"弛豫-扩散"分布,并基于BRD算法迭代确定最佳正则化因子.模拟实验与数值分析表明,该算法无需先验信息、运算效率高、相对误差小,在原始数据信噪比低至50时,仍可有效获取流体(T2,D)二维分布.在二维核磁共振测井数据实时解释应用中,该方法较传统反演算法(如TSVD)具有较大优势.同时,在自主研发的核磁共振测井仪测量CuSO4溶液(T2,D)分布的实验显示,本文设计算法对弛豫时间和扩散系数的反演误差分别仅为2%和4%,较TSVD算法有较大改善. 相似文献