首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Surface wave methods are becoming increasingly popular in many geotechnical applications and in earthquake seismology due to their noninvasive characteristics.Inverse surface wave dispersion curves are a crucial step in most surface wave methods.Many inversion methods have been applied to surface wave dispersion curve inversion,including linearized inversion and nonlinearized inversion methods.In this study,a hybrid inversion method of Damped Least Squares(DLS) with Very Fast Simulated Annealing(VFSA) is developed for multi-mode Rayleigh wave dispersion curve inversion.Both synthetic and in situ fi eld data were used to verify the validity of the proposed method.The results show that the proposed method is superior to the conventional VFSA method in aiming at global minimum,especially when parameter searching space is adjacent to real values of the parameters.The advantage of the new method is that it retains both the merits of VFSA for global search and DLS for local search.At high temperatures,the global search dominates the runs,while at a low temperatures,the local search dominates the runs.Thus,at low temperatures,the proposed method can almost directly approach the actual model.  相似文献   

2.
正则参数控制下的波阻抗约束反演   总被引:16,自引:4,他引:16       下载免费PDF全文
通过势函数方式将波阻抗反演的病态问题转为良态问题,并且给出了边界保护势函数所具备的条件. 在反演过程中,通过改变正则参数数值以及合理地选择正则参数的初值,改善反演结果,提高反演收敛速度. 同时,在具体反演中使用快速模拟退火算法,可以克服目标函数局部极值的限制,从而获得全局最优解. 通过理论模型试算和实际资料处理,说明本文方法具有精度高、实用性强的特点.  相似文献   

3.
本文基于Metropolis接收原则改进快速模拟退火算法同时进行地下介质层速度和层厚度的联合反演.通过数值模拟验证,该方法在调整好起始温度、温度下降方案等参数的情况下,通过非线性反演可以得到接近理论模型的结果.通过对沿86.E经线南北跨越印度板块、青藏高原和塔里木盆地的垂直剖面的纯频散进行反演,得到S波速度结构,由此识别其Moho 面深度和岩石圈底界面等层面位置.  相似文献   

4.
为了精确得到大陆地震震源位置,利用模拟退火法对大陆地震参数进行反演,并对改变算法中相关参数后的反演结果进行分析,表明反演精度受降温速度和收敛阈值的影响较大,初始温度只影响迭代次数,模拟退火法在对模型参数进行全局搜索情况下能获得质量较好的最优解。  相似文献   

5.
量子退火最优化与地球物理反演方法   总被引:2,自引:1,他引:1       下载免费PDF全文
量子退火最优化方法的思想是基于模拟退火方法发展而来的,但利用的是量子跃迁隧道效应的机制,这和模拟退火基于热力学的退火原理是不一样的.本文首先概述了量子退火和模拟退火的区别所在;其次主要论述了量子退火在最优化方面的应用发展情况,比如伊辛自旋玻璃问题,旅行商问题,布尔可满足性问题等,讨论了具体测试问题中量子退火方法的应用情况以及和模拟退火方法的比较;随后对基于量子退火方法的地球物理反演原理和方法做了相应的阐述,讨论了不同参数比如温度、横向场的设置,比较了处于不同的横向场和势能状态的概率大小问题;最后就量子退火方法在地球物理反问体中的应用前景以及存在的问题做了简要概述.  相似文献   

6.
瞬变电磁反演存在高度的非线性特征,常用的最小二乘等线性反演方法往往对初始模型高度依赖,并且极易陷入局部最优解.本文基于观测数据与模拟数据的L1范数建立目标函数,采用模拟退火非线性全局最优化方法实现瞬变电磁一维反演.初始模型完全随机产生,通过指数函数退温机制模拟系统能量最小实现迭代,通过接收概率函数评价当前模型,实现局部最优解的跳出,最终实现全局最优化求解.通过数值算例发现,无论给定的反演层数等于还是大于设计模型,都可以获得较好的反演效果,因而可以在反演初始就设计较多的层数,实现反演模型的自动拟合;同时,利用含噪声数据反演进一步验证算法的稳定性.最后,对实测数据进行了反演测试,结果与钻孔编录基本一致,表明提出的基于L1范数的模拟退火反演可用于实测数据处理.  相似文献   

7.
量子退火反演的原理和实现   总被引:22,自引:5,他引:17       下载免费PDF全文
现有的非线性反演方法,大都存在着求解的局部收敛性和计算速度慢的问题. 在与模拟退火方法进行比较的基础上,本文研究了基于量子跃迁过程中量子隧道效应的量子退火反演算法. 数值试验结果表明,量子退火方法应用于地球物理反问题的求解是成功的,它适合于非线性、多极值的地球物理反演问题,与模拟退火反演相比,量子退火反演在退火收敛速度和避免陷入局部极小等方面有着一定的优势,该方法也适用于其他领域非线性最优化问题的求解,具有较强的普适性.  相似文献   

8.
套管井井壁附近地层横波速度径向分布反演   总被引:2,自引:0,他引:2       下载免费PDF全文
王兵  马明明  刘鹤  刘志军 《地球物理学报》2016,59(12):4782-4790
套管外地层受异常地应力、油气开采的影响,在径向上表现出非均质性;采用声波测井可以对该非均质性进行探测,利用偶极子横波测井数据可以对横波速度的径向分布进行反演.本文建立了套管井外地层横波速度径向分层参考模型,采用修正的微扰法计算了该模型的偶极弯曲波频散曲线,建立了横波速度径向分布反演目标函数,采用高斯牛顿法和快速模拟退火法对目标函数进行了求解,得到了套管井外地层横波速度的径向分布.分析了偶极弯曲波频段、套管横波速度对反演结果的影响,对比了高斯牛顿法和快速模拟退火法对反演过程的影响.分析对比结果表明,采用偶极弯曲波激发强度较高的频段与采用全频段的反演结果相近;套管的横波速度准确度越高,反演结果越准确;高斯牛顿法和快速模拟退火法计算精度相同,都可以得到高精度的横波速度径向分布;快速模拟退火法的计算效率略低于高斯牛顿法,但其收敛性对初始值依赖更小,实际处理中应选择快速模拟退火法.  相似文献   

9.
In this paper, a global inversion method is developed for seismic moment tensor inversion by using the body wave forms. The algorithm depends on neither the selection of starting model nor the forms of objective function and constraints. When the error function, measure of the difference between the observed and synthetic waveforms, is chosen as the objective function, the best fitting source model is found; when a certain combination in seismic moment tensor elements is selected as the objective function and the values of error function are constrained in a suitable bound, the extreme source models can be produced by minimizing or maximizing this combination. By changing the form of the combination of moment tensor elements, a variety of different source characteristics can be considered. Therefore the extreme solution provides an estimation of the uncertainty in the best fitting source model. The seismic waveform data was used to evaluate the effectiveness of this algorithm. This research was supported by the National Natural Science Foundation of China.  相似文献   

10.
基于精确震源函数的解调包络多尺度全波形反演   总被引:3,自引:3,他引:0       下载免费PDF全文
本文提出解调包络方法来重构地震记录中缺失的低频信号,同时该方法能够降低全波形反演的非线性程度;提出伴随状态震源函数反演方法来得到精确的震源函数,并推导了梯度计算公式;解调包络方法结合低通滤波技术,实现了从低频到高频的多尺度反演策略,有效缓解了全波形反演的周波跳跃问题.数值算例证明了解调包络、伴随状态震源函数反演方法和低通滤波多尺度反演策略的可行性及优越性.震源函数反演精度测试结果表明:即使观测记录在缺失低频信息的情况下,也能反演得到精确的震源函数.缺失低频测试和抗噪能力测试结果表明:即使地震数据中缺失9Hz以下的低频信号或者信噪比极低的情况下,利用反演得到的精确震源函数进行解调包络多尺度全波形反演,同样可以得到高精度的全波形反演结果.与Hilbert包络全波形反演对比结果表明:解调包络在重构低频和降低伴随震源主频方面具有一定优势.  相似文献   

11.
快速模拟退火算法中,要求扰动后新的模型参数必须在取值范围内,对此,人们通常采用两种方法:(1)反复扰动直到新的模型参数满足要求;(2)当新的模型参数超过最大值时,参数等于最大值;新的模型参数小于最小值时,参数等于最小值.研究发现,方法1在计算过程中增加了寻找新模型参数的时间开销,计算效率明显降低;方法2当目标函数在边界或边界的小邻域存在多个极小值时,在有限的时间内,其陷入局部极小的概率大大增加,优化效果明显下降.本文在前人研究的基础上提出了改进的模拟退火算法,同时利用改进的算法反演了地层泊松比、纵波速度和密度.结果表明,改进的算法加快了模拟退火算法的收敛速度,提高了模拟退火算法的效率和精度.  相似文献   

12.
Inversion for seismic impedance is an ill-posed and band-limited problem so that inversion results are non-unique and unstable and low and high frequency components of inversion results are missed. Combining regularization with constraints of sonic log data and geological structure information can help to alleviate these problems. To achieve this, we developed an inversion method by constructing a new objective function which includes edge-preserving regularization and soft constraint based on Markov random field (MRF). The method directly introduces absolute constraints with prior impedance and sonic log data in the objective function and indirectly achieves relative constraints with geologic structures of layer interfaces and faults by adjusting the regularization parameter which is the scaling parameter δ. Moreover, we improved the inversion result using anisotropic diffusion smoothing method. Optimization approach utilized in inversion is a fast simulated annealing (FSA). We test the method on both synthetic and field data examples. Tests on 2-D synthetic data indicate that aspects of the discontinuity in the inversion results are significantly improved by adding δ values in faults and layer interfaces. We obtained better results by combining the first-order neighborhood and the third-order neighborhood of MRF. The inversion results of the field data provide more detailed information of the layers. The results of nearby faults were improved by introducing the geological structure constraints.  相似文献   

13.
基于遗传算法的CSAMT最小构造反演   总被引:15,自引:4,他引:11       下载免费PDF全文
利用遗传算法进行不考虑近场校正的全场资料CSAMT反演研究.遗传算法属于全局最优化方法,具有对初始模型依赖小,不易陷入局部极值的优点,然而,当未知数较多时,多解性仍是该方法的瓶颈.为了减小多层反演的多解性,在反演中引入最小构造约束,针对CSAMT的遗传算法反演问题定义了最小构造目标函数,经过模型试验找到了其具体表达式,并找到了适合CSAMT资料反演的拉格朗日乘子的最佳取值μ=0.5,实现了基于遗传算法的CSAMT最小构造反演.利用H、A、K、Q和HKH、KHA模型对方法进行了数值试验,在无噪和加入10%噪声情况下,反演结果与模型一致;加入20%噪声后,反演仍取得良好结果,与理论模型基本吻合.将该方法用于水平层状地层和横向变化地层的实测资料反演,结果与地质资料吻合.不同的计算实例表明了该方法的有效性.  相似文献   

14.
近年来,海洋可控源电磁法(MCSEM)被引入油气勘探领域以降低勘探风险.在海洋环境中,受沉积因素所造成的电阻率各向异性的影响,地电模型往往会非常复杂.为更好地反映地下电性结构,本文实现了基于VTI各向异性介质的频率域海洋可控源电磁三维反演.其中,正演采用基于Yee氏交错网格的三维有限差分算法,所形成的离散线性系统通过大规模并行矩阵直接求解器(MUMPS)进行求解.反演采用基于不等式约束的有限内存BFGS(L-BFGS)算法.最后,利用VTI各向异性介质合成数据,分别进行了电阻率各向异性覆盖层和电阻率各向异性高阻层的三维反演,结果表明:(1)基于并行直接法的MCSEM非常适用于海洋电磁所特有的多场源问题;(2)针对各向异性覆盖层模型进行三维各向异性约束反演,提高了解的可靠性;(3)针对电阻率各向异性高阻层,Inline和broadside数据覆盖的反演结果对异常体位置有很好的反映.  相似文献   

15.
为充分利用微动信号中基阶和高阶模式瑞雷波,本文研究了基于多阶瑞雷波SPAC系数直接反演的方法.该方法首先基于地层介质响应计算多阶瑞雷波的能量占比,考虑实际观测台阵有限台站个数对SPAC系数影响,正演计算多阶瑞雷波SPAC系数,再采用快速模拟退火算法对其反演以获得地下介质横波速度结构.在此基础上,本文通过数值模拟验证该方法的可靠性,分别选取三种典型地质模型,基于模式叠加算法合成理论微动信号,采用本文方法计算其理论多阶瑞雷波SPAC系数并反演,给出反演结果与真实模型对比.我们将该方法应用于上海中心城区的地质调查中,通过与钻探结果对比,进一步验证该方法的有效性.本文理论与实际应用研究表明,基于多阶瑞雷波SPAC系数直接反演的微动探测方法有助于提高反演结果的可靠性,尤其对含软硬夹层的复杂地层介质,可提高探测精度.  相似文献   

16.
单相介质AVO反演的精度分析   总被引:8,自引:8,他引:0       下载免费PDF全文
振幅随偏移距变化(AVO)反演是一个非线性的组合最优化过程,理论上可先将该非线性问题线性化,然后求解线性问题;或者直接利用非线性的模拟退火、遗传算法等方法求解.但无论哪种反演思路,实际中影响其精度的因素很多,因此分析AVO反演中的误差来源对提高反演精度和评价反演方法的可靠性非常重要.本文对能造成反演误差的主要因素,噪音、薄层调谐、地震数据处理中的误差、入射角范围等进行了分析,讨论了这些因素可能对反演结果造成的影响,发现在AVO反演过程中可以从优化参数选择和针对性处理方面来减小这些误差,提高反演精度.  相似文献   

17.
青藏高原是研究地球动力学问题天然的试验场,也是研究全球变化的关键地区之一,中国科学院青藏高原研究所于2004年在西藏那木错建立了圈层相互作用综合观测研究站.为了解那木错站下方的地质构造,于2005年8月在本站布设了宽频带地震仪(记录器为Reftek130,摆为STS2),并于2006年5月取得首批数据.本文利用宽频带地震仪提供的三分量地震波形记录,应用转换函数及快速模拟退火算法对那木错站下的地壳横波速度结构进行了反演.反演结果表明,那木错站Moho面深度在85公里左右,地壳结构复杂,尤其在中上地壳,明显呈高低速互层结构,反映了该地区构造活动、物质交换活跃,表明这些地区还未达到均衡;同时,为了便于研究该地区地壳结构,还对IRIS提供的临近地区的地震宽频带进行了处理,为研究那木错地区圈层结构提供地球物理依据.  相似文献   

18.
Summary A method of determining simultaneously the moment tensor and source time function in the point source approximation is presented. For trial values of the moment tensor components and of the source time function, parametrized by the sum of overlapping triangles delayed in time, theoretical seismograms can be synthetized and compared with the recorded ones. The iterative procedure determines the adjustment of source parameters until a good correlation of both synthetic and observed records is reached. The Green functions in a horizontally stratified medium are constructed with the use of a modal summation method.The limits of applicability of the algorithm are illustrated by the inversion of four synthetic seismograms constructed for two horizontally stratified models of the structure in Friuli area, Italy. The records constructed for the same structural model as for which the Green functions were computed can be inverted even in the high-frequency range. In the opposite case, when the records and Green functions used corresponded to different structural models, a good correlation of the input records with the final synthetics was obtained for low - pass filtered data only.Additional tests performed with input seismograms contaminated with random noise yielded good resolution of the moment tensor and the duration of the source time function even for a high noise to signal ratio.  相似文献   

19.
In the context of wide-angle seismic profiling, the determination of the physical properties of the Earth crust, such as the elastic layer depth and seismic velocity, is often performed by inversion of P- and/or S-phases propagation data supplying the geometry of the medium (reflector depths) or any other structural parameter (P- or S-wave velocity, density...). Moreover, the inversion for velocity structure and interfaces is commonly performed using only seismic reflection travel times and/or crustal phase amplitudes in isotropic media. But it is very important to utilize more available information to constrain the non-uniqueness of the solution. In this paper, we present a simultaneous inversion method of seismic reflection travel times and polarizations data of transient elastic waves in stratified media to reconstruct not only layer depth and vertical P-wave velocity but also the anisotropy feature of the crust based on the estimation of the Thomsen’s parameters. We carry out a checking with synthetic data, comparing the inversion results obtained by anisotropic travel-time inversion to the results derived by joint inversion of seismic reflection travel times and polarizations data. The comparison proves that the first procedure leads to biased anisotropic models, while the second one fits nearly the real model. This makes the joint inversion method feasible. Finally, we investigate the geometry, P-wave velocity structure and anisotropy of the crust beneath Southeastern China by applying the proposed inversion method to previously acquired wide-angle seismic data. In this case, the anisotropy signature provides clear evidence that the Jiangshan-Shaoxing fault is the natural boundary between the Yangtze and Cathaysia blocks.  相似文献   

20.
重力与地震资料的模拟退火约束联合反演   总被引:25,自引:8,他引:17       下载免费PDF全文
联合反演是综合地球物理研究的重要定量解释手段.本文在总结和分析重力与地震资料联合反演的研究现状基础上,利用改进的全局寻优的快速模拟退火算法,实现了重力和地震资料的约束同步联合反演.针对性地设计了密度和速度界面不完全一致的模型,理论模型的试验说明了方法的效果和适用性.结合最近完成的广东徐闻地区实际资料的处理和解释,表明该方法可准确确定复杂构造物性界面的密度和速度结构,在该地区的油气勘探中发挥了作用.在先验信息约束下,该联合反演方法要明显优于单独的重力反演.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号