首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 515 毫秒
1.
线性方程组的正则化解法及其应用   总被引:3,自引:0,他引:3  
在重力资料反演解释中,正则化法解线性方程组是一个行之有效的方法。这种方法的特点是引入正则化因子,求泛函极小和构组一族正则化近似并从中寻求最佳解。在求解过程中,要求对所求的解答预先有一个初步的推断解释。为此,要求对所解释的问题先进行一番研究,并对综合利用各种方法和手段所获得的先验资料进行对比,提出比较合适的基望值。本文给出了一个参数为已知的理论模型体的横剖面,与正则化法和其它线性方程解法的结果作对比,说明正则化法的有效性。文章也给出了利用正则化法对横穿吕宋岛、马尼拉海沟的一段地质地球物理剖面以及中国南黄海海域内─地质地球物理调查剖面进行反演解释的应用效果。  相似文献   

2.
羌塘盆地综合地球物理剖面的模拟解释   总被引:4,自引:3,他引:1  
根据羌塘盆地重磁震数据资料,使用GM SYS软件系统对重力、磁法及地震位图三种地球物理剖面进行处理解释,共完成4条地球物理剖面的模拟解释,并由其构组羌塘盆地地质-地球物理综合解释栅栏剖面图.地球物理剖面"栅状图"反映了研究区内深层勘探目标的变化趋势和构造特征,为使用重、磁、震资料进行非震地球物理反演提供了约束骨架.若结合研究区内钻孔及其他地质资料,可为盆地进行三维地球物理精细反演提供约束和初始模型.  相似文献   

3.
为探明郑州市及其周边地区隐伏构造的分布情况,在郑州北郊老鸦陈断层东侧布设了一条EW向音频大地电磁(AMT)测线,并对野外数据进行了去噪、主轴旋转、阻抗分析、静态校正等处理。用正则化反演理论方法——非线性共轭梯度法(NLCG)进行二维反演,利用与AMT测线接近平行的地震剖面分层结果对模型进行约束反演,降低了反演的多解性以及NLCG对初始模型的依赖性,最终取得了较好的电性结果。结果表明测线下方分为3层,并有一隐伏构造。结合地质、地球物理资料推断解释,此隐伏构造有可能是一条发育在晚新近系的正断层,亦或为第四系地层与新近系地层的不整合接触面。本文的研究工作为进一步探明城市隐伏构造提供了地球物理及技术手段依据。  相似文献   

4.
半航空时间域电磁法集成了地面和航空电磁法各自的优势,不仅容易实现大发射磁矩,增大勘探深度,同时还可以提高勘探效率。但是由于半航空时间域电磁法在测量过程中,收、发距会发生变化,因此与地面和航空电磁法相比,其数据处理与反演解释更加困难。以长导线电性源半航空电磁法为例,将自适应正则化反演算法应用于半航空时间域电磁数据反演中。首先给出最平缓模型约束条件下的半航空时间域电磁数据自适应正则化反演算法,并引入并行技术实现了该反演算法的并行计算;然后对两种三层地电模型和一种六层地电模型的理论电磁响应加高斯白噪声后的数据进行反演计算;最后对其反演结果和并行计算的效率进行了详细地分析与讨论。实验结果分析表明:自适应正则化反演算法是一种有效的半航空时间域电磁数据反演解释方法,该方法具有很好的收敛性和稳定性;自适应正则化反演算法并行计算能大大减少反演计算时间,提高反演效率。  相似文献   

5.
给出电阻率、激发极化法数据的二维有限元反演方法和原理 ,并将此方法用于高龙金矿电法勘探资料解释中 ,解释结果与已知地质资料基本吻合。由此证明 ,利用双线性模型进行二维反演的有效性及其独特的优点。  相似文献   

6.
大地电磁测深资料常规反演方法的应用现状及进展   总被引:1,自引:0,他引:1  
文中介绍了大地电磁测深(MT)资料处理解释中几种常规地球物理反演方法的基本原理及该领域内的最新进展,通过应用实例分析说明在进行资料定量解释的过程中应该多种反演手段相互佐证.可以用BOSTICK反演结果建立初始模型,通过一维、二维反演,在结合地质、物探资料并通过正演模拟计算验证的基础上,建立合理的地质-地球物理模型.利用综合信息建模及反演成像技术结合现代地质理论,能有效地降低解的非唯一性,提高解释结果的可靠程度及综合地质解释结果的合理性.  相似文献   

7.
解析延拓是一种广泛应用的位场处理方法,向下延拓可以压制深部地质体的影响,突出浅部异常。但是,向下延拓滤波因子是一个高通滤波器,会造成下延结果震荡,从而限制了该方法在实际资料中的应用。文中详细介绍并实现了迭代Tikhonov正则化向下延拓方法,在理论模型上将该方法与传统频率域延拓方法进行对比,表明迭代Tikhonov正则化向下延拓方法的有效性;并将该方法应用于青海尕林格铁矿区磁测资料的处理解释中,下延结果与钻探情况相符,说明在厚覆盖层的勘查区中,运用迭代Tikhonov正则化向下延拓方法能够有效地提高资料处理解释的效果。  相似文献   

8.
为了减小地球物理反演的多解性,通常采用累加型目标函数,即数据拟合差项加上正则化项。基于累加型目标函数的反演存在如何优选正则化因子的问题,这个过程通常需要做大量的反演计算。本研究系统分析和实现了基于乘积型目标函数(数据拟合差项和正则化项相乘)的电阻率法三维反演,乘积型目标函数的反演不存在优选正则化因子的问题。三维正演采用非结构化有限单元法,三维反演采用有限内存拟牛顿方法。使用理论模型合成数据进行了三维反演,检验了基于乘积型目标函数的电阻率法三维反演的可行性与有效性。反演算例结果表明:基于乘积型目标函数的反演方法能够可靠恢复异常体的电阻率值、形态和位置。  相似文献   

9.
利用直流激发极化法对称四极测深装置对印尼塔里阿布岛铅锌多金属矿化地段进行勘探,结合该区地质特征、成矿规律及岩层地球物理特征,获得该矿区铅锌矿体(化)的产状及规模信息。研究表明,该矿区矿体的低阻高极化异常特征与围岩差异显著。采用RES2DINV电法二维反演软件对测深剖面进行反演处理,在反演时加入剖面地形数据,可有效制约地形起伏对异常体形态及埋深的干扰,提高数据处理的精确度。依据激电测深数据反演结果并结合物性工作成果,对测深剖面进行地质解译推断。利用I线剖面布置的槽探及钻探工程对解译推断成果进行验证,结果表明,见矿位置及产状规模与物探解译推断成果吻合度较高,激电测深法在该区多金属勘查中取得了较好的应用效果。  相似文献   

10.
重磁电联合反演在银—额盆地勘探中的应用   总被引:1,自引:0,他引:1  
依据大地电磁测深连续二维反演结果初步确定构造与地层,通过综合分析区内物性特征给定地层与岩石的物性参数,建立正演地质—地球物理初始模型,在地质、钻井、地震资料约束下进行重磁电联合反演,不断修改地质—地球物理模型和大地电磁测深连续二维反演结果,使重磁正演数据和实测数据拟合度达到最佳。在银—额盆地构造与地层解释中的应用结果表明,重磁电联合反演可以有效克服单方法解释的多解性,准确识别地质层位,大大提高定量解释的精度。  相似文献   

11.
大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。  相似文献   

12.
二维大地电磁正则化共轭梯度法反演算法   总被引:2,自引:0,他引:2  
针对大地电磁二维反演中目标函数收敛速度慢而且解的稳定性较差等问题,提出了大地电磁数据的正则化共轭梯度法反演算法(Regularized Conjugate Gradient Algorithm,RCGA)。此算法在构建目标函数时引入正则化的思想,利用共轭梯度法求解最优化问题。在每次迭代过程中根据目标的收敛情况更新正则化因子,有效地解决了迭代时目标函数发散的问题。最后对两个典型模型进行了试算,对其收敛速度及解的稳定性进行了分析,并与传统反演方法的计算结果进行了对比,结果都表明RCGA具有收敛速度快、精度高、结果稳定等优点。  相似文献   

13.
High-resolution deconvolution can mathematically be viewed as a regularized inverse problem. Besides, the result of the high-resolution deconvolution is generally accepted as reflectivity series of the layered media. On the other hand, lateral continuity is frequently poorer than vertical resolution on post-stack seismic section after application of any high-resolution deconvolution. However, because of the ill-posed inherent of the deconvolution problem, the Cauchy norm regularization term, a non-quadratic prior-information is widely used to provide the stability and uniqueness of the problem. But, it does not provide adequate quality of deconvolution if the noise in the data is strong. In this study, a stable and high-resolution deconvolution of post-stack seismic data was accomplished by an iterative inversion algorithm incorporating the Cauchy norm regularization with FX filter weighting. Cauchy norm regularization was utilized to force the solution to a spikiness structure, while the effective random noise reduction was performed by using the FX filter. Applications to synthetic and real post-stack data showed that the resolution in the vertical direction and continuity in the lateral direction are better improved. Thus, we think that this process makes seismic sections obtained especially from thin layered sedimentary basins more interpretable.  相似文献   

14.
LSQR法具有计算效率高、对计算机内存要求低的优点,适合于大规模问题的求解。为探讨其应用于位场反演的稳定性和可靠性,笔者以加入不同噪声的两个合成模型数据为实验对象,比较分析了Tikhonov正则化与LSQR法求解结果,显示直接利用LSQR法求解位场反问题能够得到满意的正则化解,其解模型相对Tikhonov正则化,最大相对误差仅为0.36%,说明直接利用LSQR法求解位场反问题是可行的。将其应用于四川盆地雅安地区重力三维反演,极大地降低计算成本,获取了区内沉积盆及主要断裂分布情况,为页岩气靶区优选提供了有力支撑。  相似文献   

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

16.
The dissipation test evaluation method presented here employs: (i) a point-symmetric, linear, coupled consolidation model with a new boundary condition and a new time factor, (ii) an automatic and mathematically precise, non-linear, inverse problem solver which includes some reliability testing methods and a kind of regularization technique, and (iii) a method to identify the initial condition.The evaluation method is equally applicable to pore water pressure data displaying monotonic or non-monotonic time variation, with the only difference that the initial condition is identified differently during the inverse problem solution.The necessary testing time is very short provided that the pore water pressure is measured well above the tip. This is attributed to the ability of the one-dimensional, linear consolidation model to account for the geometry and the unloading effects that occur when steady penetration of the static cone penetrometer ceases.  相似文献   

17.
Classical 3D/4D variation fusion is based on the theory that error follows Gaussian distribution. When using minimization iteration, the gradient of objective function is involved, and the solution of which requires the continuity of data. This paper adopted the extended classical 3D/4D variation fusion method, and explicitly applied the prior knowledge, which was based on L1-norm, as regularization constraint to the classical variation fusion method. Original data was firstly projected into the wavelet domain during the implementation process, and new fusion model was adopted for data fusion in wavelet space, then inverse wavelet transform was used to project the result to the observation space. Ideal experiment was carried out by using linear advection-diffusion equation as four-dimensional prediction model, which made a hypothesis of the discontinuity with the data between background and observation, and that meant the derivatives between left and right were not equal on some points. The result of the experiment showed that the method adopted here was practicable. A further research was also done for multi-source precipitation fusion. Firstly, CMORPH inversion precipitation data were corrected through PDF (Probability Density Function, PDF) matching method based on GAMMA fitting function. Then corrected data was fused with the observation one. By comparison with the reference field, the result showed that this method can keep some outliers better, which might represent certain weather phenomenon. The L1-norm regularization variation fusion in this paper provided a possible way to deal with discrete data, especially for jump point.  相似文献   

18.
The present study is an attempt to upgrade the previously used semi-analytical inverse slope method. The upgrade method was applied to the interpretation of VES’s measured with different electrode arrays in two and three layers laboratory tank model. Then the method was applied on VES’s measured at different sites of different geologic conditions. The study showed that the inverse slope method is a quick method for obtaining depths to interfaces. The method can be implemented on a hand-held calculator or more efficiently on a small microcomputer where the data can be processed at a faster rate than that of conventional data acquisition systems. This allows the method to be potent tools for infield data analysis and more cost effective to assess the success of a sounding on-the-spot and modify the experimental set-up as appropriate. The comparison between the results obtained from this method and the traditional digital computerized methods indicated that the geoelectrical models derived from inverse slope method have a good correlation with the actual model. Also, when applying the method on VES’s measured by using Schlumberger and Wenner arrays at areas of different geologic setting, gave results in comfortable with the actual lithology obtained from the lithologic data of the neighboring wells. The results obtained from the lab experiment and field survey indicated that the inverse slope method has the potential to be applied on VES data collected with any electrode array, the advantage that is not available with all the traditional interpretation methods available nowadays.  相似文献   

19.
Doklady Earth Sciences - Existing methods for solving inverse problems, such as the regularization method, look mostly for a quasi-solution that may not be a solution to the original problem, but...  相似文献   

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

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