基于序列二次规划方法的高密度电阻率反演

针对高密度电阻率成像法反演计算中存在的困难，探讨了基于序列二次规划（SQP）方法建立可行的电阻率反演算法的问题。在对三维点电源二维地电体电位场模型有限元法正演计算的基础上．建立了基于SQP方法的高密度电阻率成像法反演算法。根据勘测中测量电极间距的不同，提出按不同的电极间距分别建立优化模型进行参数优化。通过对模拟电阻率模型和实际观测资料的反演计算．表明该算法进行高密度电阻率反演是可行的，其具有对初始模型无特殊要求、收敛速度快的特点。这一反演算法可望在高密度电阻率成像的反演计算中得到应用。     

利用改进模拟退火算法反演甘肃兰州地区地幔转换带导电模型

地磁测深一维反演可以获得地磁台站下方地球深部电导率结构.将改进模拟退火算法应用于地磁测深数据的反演,并将基于反演过程中的所有优质解进行统计分析得到的期望值作为模型最优解.地磁测深合成数据的一维反演测试结果表明:对于无噪声的理论数据,传统模拟退火算法可以恢复出真实的电导率剖面;存在噪声时,相较传统模拟退火算法,改进模拟退...     

重力数据约束的二维地震层析成像在台湾海峡南部的应用

地震-重力联合反演可以降低多解性,针对大勘探范围内岩石波速-密度关系存在较大散布的问题,引入交叉梯度结构约束构建统一的目标函数是一种有效的解决方案。利用台湾海峡南部HX-13测线进行验证:由于火成岩侵入体的存在,部分区域难以拾取可靠的地震初至,常规的走时反演无法准确恢复基底面的形态。搜集研究区的船测重力资料,根据地震地质条件和数据特点设置合适的反演参数,实现基于物性和结构双重约束的重力-地震联合反演。由对反演结果的分析和解释可以看出,该方法可在很大程度上弥补地震数据不完备的缺陷,使反演过程稳定,并提高模型的可靠性。     

概率积分法沉陷预计与参数反演优化算法及实现

任意多边形工作面沉陷预计及精确的求取概率积分法反演参数是开采沉陷研究的重点。提出了使用Delaunay三角剖分将任意多边形开采区域划分为若干小三角形,然后进行沉陷预计的方法,解决了任意多边形沉陷预计的难题;综合模矢法和遗传算法优点,提出了组合算法,提高了概率积分法参数反演精度和效率;开发了集实测数据处理、预计模型参数反演和移动变形预计模块于一体的地表移动变形数据处理软件。优化算法的提出及系统开发,可以更好的服务于"三下采煤"、矿区建筑物稳定性评估、矿区土地复垦、采煤塌陷地治理规划设计等工作。     

Application of constrained sparse spike inversion in reservoir predication:a case study of Oriente Basin in South America

In the conditions of low Signal-to-Noise Ratio(SNR) of seismic data and a small quality of log information,the consequences of seismic interpretation through the impedance inversion of seismic data could be more precise. Constrained sparse spike inversion(CSSI) has advantage in oil and gas reservoir predication because it does not rely on the original model. By analyzing the specific algorithm of CSSI,the accuracy of inversion is controlled. Oriente Basin in South America has the low amplitude in geological...     

地震勘探全波形反演的应用与发展分析

本文首先对20世纪80年代发展起来的全波形反演应用及其在勘探地球物理领域的发展进行了分析;其次,面对定量化、精细化的地震勘探要求,提出了将地震勘探全波形反演与其他数据处理环节或处理技术相结合的研究设想,并展望了全波形反演的发展趋势;最后,论述了全波形反演研究中地震波场数值模拟、反演初始速度模型获取、目标函数形式选择、寻优算法启用及各向异性介质中的应用等关键问题,并总结了通过Laplace域的全波形反演获取反演初始速度模型、结合射线追踪并充分发挥并行计算之于波动方程方法来模拟地震波场的巨大优势,及灵活选用反演目标函数形式和寻优算法更新速度模型参数来加快全波形反演方法的实用化进程。     

Study of frequency domain full waveform inversion based on Huber norm and L-BFGS algorithm

Full waveform inversion( FWI)is a high resolution inversion method,which can reveal detailed in-formation of the structure and lithology under comPlex geological background. It is limited by many kinds...     

Application of constrained sparse spike inversion in reservoir predication: a case study of Oriente Basin in South America

In the conditions of low Signal-to-Noise Ratio(SNR) of seismic data and a small quality of log information,the consequences of seismic interpretation through the impedance inversion of seismic data could be more precise. Constrained sparse spike inversion(CSSI) has advantage in oil and gas reservoir predication because it does not rely on the original model. By analyzing the specific algorithm of CSSI,the accuracy of inversion is controlled. Oriente Basin in South America has the low amplitude in geological structure and complex lithologic trap. The well predication is obtained by the application of CSSI.     

Joint inversion of magnetotelluric and fullwaveform seismic data based on alternating cross-gradient structural constraints

Magnetotelluric(MT)inversion and seismic inversion are important methods for the interpretation of subsurface exploration data,but separate inversion of MT and seismic produces ambiguous and non-unique results due to various factors.In order to achieve accurate results,the authors propose a joint inversion method of two-dimensional MT and seismic data in the frequency domain.The finite element method is used for numerical simulation of electromagnetic data in the forward modelling,and the Gauss-Newton method is used for the inversion.The 9-point-finite-difference method is used to solve the seismic wave field in the acoustic wave equation,and the inverse problem of seismic data is solved by full waveform inversion with a conjugate gradient,a simple and fast method.Cross gradient functions are used to provide constraint structure between resistivity and velocity parameters to carry out the joint inversion.The joint inversion algorithm is tested by double-rectangular model synthesis data,and the accuracy of the algorithm is verified.The results show that the joint inversion results are better than those from separate inversion.The algorithm is applied to a geophysical model of a metalliferous deposit in Jinchuan and is compared with the separate inversion results.It shows that the results obtained with joint inversion are much closer to the real model.     

3D fast inversion of gravity data based on GPU

Large-scale gravity 3D interpretation depends on efficient and high-resolution 3D inversion processing of massive data.The authors applied the conjugate gradient method with minimum support function and prior model constraint to reduce multi-solutions in gravity inversion.Based on the parallel programming and compu-ting platform NVIDIA CUDA with C++language, we achieve fast 3D gravity inversion by adopting GPU parallel technique into the most time consuming part relating to sensitivity matrix.The results of theoretical model show that the abnormal body can be clearly located and the inversion speed is improved greatly.Comparing with in-version speed of the Matlab program, speed of inversion with GPU parallel technique has improved more than 100 times under the hardware condition of Geforce GTX 1060 graphics card.     

位错模式反演的算法研究

综述了大地测量反演算法的发展状况，详细介绍了模拟退火，随机耗费和区间算法等3种优化方法，并基于位错模式，采用模拟的重力测量观测数据，比较了3种算法的反演效果，结果表明：模拟退火法优于随机耗费法，而区间算法的可靠性又优于模拟退火法，区间算法是最可靠的反演方法。     

小波变换和改进Burg算法的GNSS-IR海面高度反演模型

针对传统的GNSS-IR海面高度监测方法信号分离不佳且精度有待提高的问题,提出结合小波变换和改进Burg算法的新型GNSS-IR海面高度反演模型。相比于传统的多项式拟合法,小波变换得到的SNR振荡项更加完整、精确。改进的Burg算法能有效抑制峰值偏移或谱分裂现象,提高谱分析精度。基于瑞典Onsala空间天文台提供的GNSS数据和验潮仪数据的实验结果表明,优化后的海平面测高模型的反演结果与验潮仪数据具有较高的一致性,相比于传统的GNSS-IR海面测高模型精度提高约20%。     

基于压缩感知原理的L1范数AVO反演算法

基于压缩感知稀疏信号采样与重构理论,利用AVO反演方法将传统的L2范数改变为L1范数,反演地下地层在L1范数下的稀疏脉冲反射系数。反演得到的稀疏尖峰将局部地下结构通过有限数量的层状结构的叠加来表示,能够提高纵向精度,与传统的AVO反演算法相比提高了薄层的反演效果且具有一定的抗噪性。数值模型及实际数据的结果表明,基于压缩感知原理的L1范数AVO反演方法更加准确、分辨率更高。     

Application of artificial bee colony algorithm in Rayleigh wave inversion

In order to solve the problems of multi-parameter, multi-extreme and multi-solution in the nonlinear iterative optimization process of Rayleigh wave inversion, the artificial bee colony(ABC)algorithm is selected for global nonlinear inversion. The global nonlinear inversion method does not rely on a strict initial model and does not need to calculate the derivative of the objective function. The ABC algorithm uses the local optimization behavior of each individual artificial bee to finally highl...     

����IPSO-BPģ�͵Ĵ�Ӷ�Դ�������Ԥ���뷴��

???????????????????IPSO??????????????????????м??????????????????????????????????????????????????????????????????????????????BP?????????γ?IPSO_BP??????????????????μ????????????????????????IPSO_BP???????????????Ч??????????????????????????????     

An algorithm for multi-exponential inversion of T2spectrum in nuclear magnetic resonance

NMＲ logging can provide the permeability parameter and abundant stratigraphical information such as total porosity,oil,gas and water saturation,oil viscosity,etc． And these physical parameters can be ob- tained by T2spectrum inversion． NMＲ inversion is an important part in logging interpretation． The authors de- scribe a multi-exponential inversion algorithm,solid iteration redress technique (SIＲT),and apply the algo- rithm in real data and compare the results with those based on singular value decomposition (SVD)． It shows that SIＲT algorithm is easier to be understood and implemented,and the time spent in SIＲT is much shorter than that of SVD algorithm． And the non-negative property of T2spectrum is much easier to be implemented． It can match the results based on SVD very well． SIＲT algorithm can be used in T2spectrum inversion for NMＲ analysis．     

3D Joint inversion of gravity and gravity tensor data

The processing and interpretation of gravity and gradient data plays an important role in geophysics. The cross gradient joint inversion is usually used for achieving structure coupling of multiple geophysical mod-els.In order to realize the coupling of gravity and gravity tensor data, the authors analyzed each component. The results show that different types of data contain different direction information, and derived the joint inver-sion based on cross gradient function and applied it to model data.The theoretical model results show that the cross gradient method can reduce the multi solution and significantly improve the resolution of the inversion. The method was also applied to inverse the gravity tensor data in Vinton salt dome, showing that this method can get higher resolution results than the separate linear inversion, and be closer to the real density from drilling data.     

基于GNSS-IR双系统组合的土壤湿度多星线性回归反演模型

提出一种GPS/BDS双系统组合的土壤湿度多星线性回归反演模型,并以GNSS接收机实测数据为例,对比分析不同GPS和BDS卫星组合反演土壤湿度的效果。实验表明：1)GPS和BDS双系统组合相对于单系统在短观测时间内可以提高有效卫星数,通过多元线性回归原理可实现双系统多卫星的有效融合,提高土壤湿度反演的精度;2)当GPS和BDS组合卫星数达到6颗以上时,反演效果趋于稳定,反演结果与土壤湿度的相关系数均优于0.90,RMSE相对于单星至少提高25.8%。     

青藏东缘三维Moho界面的位场遗传算法反演

针对波数域中Parker-Oldenburg位场迭代反演方法中存在向下延拓的问题,提出了利用改进的遗传算法进行三维Moho界面反演的位场遗传反演方法。通过模型试算证明了该方法的有效性和稳定性,并利用该方法反演了青藏高原东缘三维Moho界面的分布。     

附加区域电离层约束的PPP定位收敛性分析

为了改善PPP的定位精度和收敛速度,提出一种有效的解决思路:在非差非组合模型PPP定位的基础上,附加高精度区域电离层先验信息约束,利用美国西海岸CORS网内183个参考站的观测数据进行区域电离层建模,获取高精度电离层延迟信息;通过实验对比分析无电离层组合模型、非差非组合模型和附加电离层约束的非差非组合模型3种算法的PP...