Undersea Buried Pipeline Reconstruction Based on the Level Set and Inverse Multiquadric Regularization Method

The electric inversion technique reconstructs the subsurface medium distribution from acquired data. On the basis of electric inversion, objects buried under the earth or seabed, such as pipelines and unexploded ordnance, are detected and located in a contactless manner. However, the process of accurately reconstructing the shape of the target object is challenging because electric inversion is a nonlinear and ill-posed problem. In this work, we present an inverse multiquadric (IMQ) regularization method based on the level set function for reconstructing buried pipelines. In the case of locating underwater objects, the unknown inversion area is split into two parts, the background and the pipeline with known conductivity. The geometry of the pipeline is represented based on the level set function for achieving a noiseless inversion image. To obtain a binary image, the IMQ is used as the regularization term, which 'pushes' the level set function away from 0. We also provide an appropriate method to select the bandwidth and regularization pa-rameters for the IMQ regularization term, resulting in reconstructed images with sharp edges. The simulation results and analysis show that the proposed method performs better than classical inversion methods.     

用U曲线法确定地震同震滑动分布反演正则化参数

针对地震滑动分布反演中正则化参数选取问题,提出利用U曲线法确定地震滑动分布反演正则化参数。利用U曲线法、L曲线法设计模拟实验,并将两种方法应用到芦山实际震例反演中。模拟实验以及芦山实际震例反演结果表明,利用U曲线法确定正则化参数反演地震滑动分布结果与L曲线法相比具有反演精度高、无需依赖数据拟合精度等优势。     

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...     

基于归一化磁源强度的聚焦反演方法

针对剩磁条件下铁磁物质反演中存在的问题,提出基于归一化磁源强度的聚焦反演方法。首先,利用归一化磁源强度作为实测数据对磁性目标进行反演,减弱剩余磁化对反演结果的影响;然后,利用深度加权矩阵和最小支撑矩阵对经典Tikhonov正则化理论框架下的反演模型进行约束得到目标函数,并有效解决了核函数随深度增大而快速衰减的问题;最后,通过对目标函数进行迭代奇异值分解获得最佳物性参数,并根据Morozov偏差原则自适应地确定目标函数在迭代过程中的正则化参数,提高了迭代速度和求解精度。     

��λƫ����Ƽ�����PPPģ���ȹ̶��е�Ӧ�÷���

???????λ??????????λ??????????????????????????????????????????????????????λ?????????????????????λ????????????????????????λ????????????????????????????????λ???????????????λ???????????????????????????3/4??????????????Ratio?????δ?????????3??????λ????????????????     

MODIS_LST与AMSR-E_BT的相关性及地表温度反演

本文以2007年和2008年MODIS每日地表温度（LST）数据及AMSR-E地表亮温（BT）数据为研究对象,结合土地覆盖类型数据,统计分析MODIS_LST与AMSR-E_BT在不同土地覆盖类型、频率和极化方式条件下的相关性。结果表明,频率在18.7、23.8和36.5 GHz的AMSR-E-BT与MODIS_LST的相关性较大,且在垂直极化通道上的相关性较在水平极化上大;不同土地覆盖类型,与MODIS_LST相关性较大所对应的AMSR-E微波通道不同。同时,考虑混合像元问题对相关性的影响,对25种不同地物类型组合下MODIS_LST与AMSR-E-BT的相关性进行统计分析,发现混合像元中地物类型越多,则二者相关性越小。最后,采用多元线性回归分析法,根据不同土地覆盖类型建立反演回归模型,对部分研究区域MODIS-LST进行反演,误差平均在±3.15 K以内,与不考虑下垫面覆盖的模型比较,反演MODIS_LST精度平均提高了1.5 K。     

Application of homotopy parameter inversion method in Miyun Reservoir

The large-scale convergence of homotopy parametric inversion method on the water quality model parameters calculated was used, with application in parametric inversion calculation of total phosphorus of Beijing Miyun Reservoir. Through calculated and compared the error of sedimentation rate by homotopy parametric inversion method and genetic inversion calculation method, the results indicate that homotopy parametric inversion method has good stability, calculating speed, and even if the initial selection away from the objective function, the solution still has a good convergence.     

Inversion Study on Pollutant Discharges in the Bohai Sea with the Adjoint Method

The adjoint method is presented which could be used to estimate the time-varying contamination concentration(CC) from pollution source(PS). Then the pollutant flux is calculated. In order to guarantee the continuity of pollutant distribution and make the calculated results more accurate, the independent point(IP) scheme is proposed. The contamination concentrations(CCs) at some time steps are selected as the independent points(IPs), and only CCs at these IPs are optimized while CCs at other points are calculated through linear interpolation of the independent CCs. In twin numerical experiments, all the given distributions are successfully inverted with the adjoint method. The cost functions and the mean absolute errors(MAEs) in concentrations and pollutant fluxes decrease greatly after assimilation, and the cost functions are reduced by about 5 orders of magnitude compared with their initial values. The results indicate that the adjoint method is computationally efficient to recover CCs from PS. It is easier to invert the given distribution which is less complex. The inversion efficiency with IP scheme is raised compared to that without this scheme. The IP scheme is significant for the inversion result, in which appropriate IP number could yield better inversion results. More work will be done to apply this method to real experiment.     

���ʶ���������С�������������ͬ�׷�����Ӧ��

???????????????????????Jacobian???????????????????????????????С?????????????????????????С?????????????????????????в?????????С??ж?????????????????????2??????????????????仮??????????????????С??????????н???????????÷???????????????????????????????????????????????????????????????????????????????????????????????????????????С??????????????С??????     

利用大地测量数据反演断层滑移分布的MCMC方法

断层滑移的反演可以表达为一个不等式约束下的正则化问题。本文基于改进的MCMC方法,利用GPS位移数据反演2015年尼泊尔地震的同震滑移分布。结果表明,主要滑移集中在震中东侧靠近加德满都的区域,最大滑移约7 m,最大滑移误差约为滑移量的10%,地震矩为7.67×1020 Nm,相当于MW7.86,与其他手段得到的结果一致。利用反演模型预测的GPS位移与观测相吻合。震例结果表明,改进的MCMC方法适用于断层滑移分布的反演,并能够从统计角度给出滑移误差。     

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 highlight the global optimal value in the colony,and the convergence speed is faster.While searching for the global optimal solution,an effective local search can also be performed to ensure the reliability of the inversion results.This paper uses the ABC algorithm to perform Rayleigh wave dispersion inversion on the actual seismic data to obtain a clear undergrounding of shear wave velocity profile and accurately identify the location of the high-velocity interlayer.It is verified that the ABC algorithm used in the inversion of the Rayleigh wave dispersion curve is stable and converges quickly.     

基于重力垂直梯度的欧拉反褶积改进措施

欧拉反褶积反演解的发散与不稳定性一直制约着其实用效果。为此,提出利用垂直梯度反演和筛选准则并在分辨信号极值区域基础上进行重力反演等改进方法。通过模型试算,全面分析它们相对常规重力数据反演方法的效果。结果表明,改进方法能明显提高反演解水平位置和深度的精度和收敛性。     

重力-地震联合反演的改进及在遂宁-阿坝剖面的应用

针对传统顺序联合反演的缺陷,通过改变地下空间的网格剖分方式和约束条件的约束方式来进行相应的改进,并将拟BP神经网络算法和小波多尺度分析应用到反演中,以提高反演结果的分辨率和收敛速度。龙门山地区遂宁-阿坝剖面的反演结果表明,本文方法具有更好的反演效果。     

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

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

Modal wave number tomography and bottom parameter dependence

A shallow water tomography scheme based on the modal wave number inversion technique is considered in this paper. The scheme is based on the assumption that modal wave number for transpped mode can be measured in a suitable way. The tomographic inversion is accomplished into steps; firstly, the bottom parameter are inverted by using the bottom reflection phase shift with the known sound speed profile; secondly, the variation of sound speed profile at different time is inverted provided the bottom parameters are known. A numerical simulation shows that the proposed scheme works well, and the sensitivity analysis of sound speed profile inversion is performed, for shallow water environmental parameters: sound speed, density and attenuation coefficient of the bottom.     

Application of homotopy parameter inversion method in Miyun Reservoir

The large-scale convergence of homotopy parametric inversion method on the water quality model parameters calculated was used, with application in parametric inversion calculation of total phosphorus of Beijing Miyun Reservoir. Through calculated and compared the error of sedimentation rate by homotopy parametric in- version method and genetic inversion calculation method, the results indicate that homotopy parametric inversion method has good stability, calculating speed, and even if the initial selection away from the objective function, the solution still has a good convergence.     

Full waveform inversion based on improved MLQN method

Full waveform inversion( FWI) is a challenging data-fitting procedure between model wave field value and theoretical wave field value. The essence of FWI is an optimization problem,and therefore,it is important to study optimization method. The study is based on conventional Memoryless quasi-Newton( MLQN)method. Because the Conjugate Gradient method has ultra linear convergence,the authors propose a method by using Fletcher-Reeves( FR) conjugate gradient information to improve the search direction of the conventional MLQN method. The improved MLQN method not only includes the gradient information and model information,but also contains conjugate gradient information. And it does not increase the amount of calculation during every iterative process. Numerical experiment shows that compared with conventional MLQN method,the improved MLQN method can guarantee the computational efficiency and improve the inversion precision.     

重磁线性反演的正则化方法

重磁线性反演最后都要归结于解数据方程。由于重磁数据方程是奇异方程,通常大多数解方程组的方法都失效。因此,本文采用了正则化方法,通过模型试验,证明了该方法具有解奇异方程和压制高频干扰的功能。为了保证反演结果的可靠,本文提出了一套比较可行的处理资料过程。最后,反演计算了山西太原盆地的莫霍面和居里等温而,其结果与已有的地质、地球物理资料相吻合。     

联合迭代重构算法在对流层水汽三维重构中的应用研究

通过数值模拟仿真将SIRT算法应用于水汽层析过程,结果显示该算法可以较好地重构水汽的空间分布,并且不受观测方程性质的影响,收敛较快,较之求逆的层析算法更易于实现。     

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

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