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.     

基于双精度与四精度的重力场解算精度分析

基于动力学方法比较分析了双精度与四精度模式下重力场模型的解算精度,主要包括缔合勒让德函数计算、数值积分器及重力场反演结果。结果显示,在勒让德函数计算方面,部分角度在双精度模式下计算至1 900阶以后会出现溢出问题,而在四精度模式下任何角度都满足精度要求,并且计算结果比双精度模式高8个量级。数值积分器Adams预测校正法积分1 d的位置和速度误差,在四精度模式下比在双精度模式下高4个量级。在精密轨道反演重力场计算方面,动力学方法在双精度及四精度模式下反演结果一致,统计其计算至60阶的累计大地水准面误差为1.29×10~(-5 )m,这是因为动力学方法的线性误差相对计算误差而言是主要误差;非线性动力学方法在四精度模式下比在双精度模式下高7个量级,其大地水准面误差分别为8.92×10~(-15) m和8.16×10~(-8) m。     

Bostick反演法误差修正

对于水平2层均匀介质模型,为提高Bostick反演法精度,推导希尔伯特变换近似所造成的误差公式,给出误差的数值计算方法,得到修正后的Bostick反演公式。对比2种方法反演结果发现,修正后的Bostick反演法能够更好地反映分界面及第2层介质的电阻率。     

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

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

近地层大气折光系数的新反演法

以近地层大气垂直折光系数的传统反演方法和理论计算模型为基础,提出了一种可以直接利用对向EDM三角高程测量结果计算往测与返测折光系数的新反演方法。运用该方法时,不需要增加外业工作量,只需通过联合求解传统反演式和理论计算式,便可求得相对可靠的折光系数值,从而有效提高EDM三角高程测量的精度。     

Inversion analysis of the shear strength parameters for a high loess slope in the limit state

It is difficult to obtain reliable shear strength parameters for the stability analysis and evaluation of high loess slopes.Hence,this paper determines slope elements and physical parameter of 79 slopes with heights of[40,120]m based on the measured loess slopes in the Ganquan and Tonghuang subregions of northern Shaanxi Province,China.In the limit state of the loess slope(stability factor Fs is 1.0),a fitting equation for the slope height and width is established.Then,the model is developed by stability analysis software-SLOPE/W,and the comprehensive shear strength parameters corresponding to different slope heights of the high slope in the study area are obtained by inversion using the Morgenstern-Price method.The results show that when the height of the slope increases,the cohesion c increases in the soil,and the internal friction angle j decreases.This change is consistent with the characteristics reflected in the composition and physical properties of the slopes,and the comprehensive strength parametric curves are very similar between the Ganquan subregion and the Tonghuang subregion.A landslide that occurred in Miaodian-zaitou of Jingyang County,Shaanxi Province,is selected to verify this inversion method,and the results show that the calculated shear strength parameters of the inversion are consistent with the measured value of the actual slope.     

一种改进的平方根滤波在卫星钟差实时解算中的应用

目前普遍采用卡尔曼滤波方法来实现非差法卫星钟差的实时解算。平方根滤波可以增强数值计算的稳定性,避免滤波发散,但矩阵的求逆会耗费较长的计算时间;序贯算法可以避免对矩阵直接求逆,能提升计算效率。综合两种算法的优点,提出一种基于序贯算法的平方根滤波器,并应用于卫星钟差的实时解算,得到的卫星钟差精度优于0.2 ns,且计算时长缩短40%。     

Multi-component joint inversion of gravity gradient based on fast forward calculation

With the development of gravity gradient full tensor measurement technique , three-dimensional (3D) inversion based on gravity gradient tensor can provide more accurate information .But the forward calcula-tion of 3D full tensor sensitivity matrix is very time-consuming, which restricts its development and application . According to the symmetry of the kernel function , the authors reconstruct the underground source of geological body to avoid repeat computation of the same value , and work out the corresponding relationship between the re-sponse of geological body to the observation point and the response of reconstructed geological body to the obser -vation point .According to the relationship , rapid calculation of full tensor gravity sensitivity matrix can be achieved .The model calculation shows that this method can increase the speed of 30-45 times compared with the traditional calculation method .The sensitivity matrix is applied to the multi-component inversion of gravity gradient .The application of this method on the measured data provides the basis for the promotion of the meth -od.     

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

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

重磁力联合反演在研究复杂盆地基底构造中的应用

本文采用重磁力联合反演方法,计算了民和盆地基底深度和视磁化强度。针对该盆地特点,设计了等效和非均质重力反演模型,在约束条件下分块反演得到准确可靠的基底深度;用基底深度和航磁异常反演基底视磁化强度,根据岩石磁化率参数,判定基底岩性横向分布。该方法不仅能提供盆地基底构造形态,而且能提供基底岩性信息。     

CRInSAR与PSInSAR融合探测地表线性变形速率

提出利用CRInSAR与PSInSAR融合解算地表变形参数的算法。将CR点上获取的线性变形速率与高程改正结果作为PS基线网络的起算数据,依据最小二乘原理求解PS点上待求线性变形速率与高程改正值的最优解。这一算法可融合PSInSAR与CRInSAR两种算法的优势,起到很好的互补作用。实验显示,这一算法获取的线性变形速率精度可达±0.37 mm/a,高程改正值精度可达±0.5 m,证实该算法在实际地表变形监测中具备可行性。     

��ʽԼ�����������Ϸ��ݵĶԱ��о�

???????????????????????????????????????????????????????????????????????????????????壬??????????????ó????Щ????????     

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.     

NOTEA NEW METHOD TO DETERMINE THE START DEPTH IN GRAVITY ANOMALY INVERSION

INTRODUCTIONThevalueofthestartdepthDshouldbegivenfirstinthegraVityanoInalinversion.DisthecmstthicknesswherethecorresPondinggraVityanomaliszerooragivenknownvalue.TheburieddepthofthedensityinterfacecalculatedfrOmthegraVityanomalinver-siondependsonthestaItdePthD-Whnthegivenvalue0fthestartdepthissInall(large)lthewholeburieddepthofthecalculateddensityinteifaceissInall(larg),sointhegravityan0Inalyinversion,thedeterminati0nofthestartdepthDisveryimPortant.Usually,Disgivenbasedonseismic,ear…     

地面倾斜的形变特征及计算方法

对地倾斜的基本概念、形变特征以及计算方法作了系统的阐述,对地倾斜面的几何物理性质作了深入分析,并指出地面倾斜量是个伪矢量。最后,对地倾斜面的正演与反演问题进行了讨论。     

顾及先验信息的模型参数反演——澜沧-耿马地震位错模型参数的求定

采用澜沧－耿马地震前后两期大地测量资料，结合发震区的地质资料和震后地表形变破坏带的调查资料，运用顾及参数先验信息的贝叶斯反演方法，联合求定了澜沧－耿马地震的两个发震断层的位错模型参数。由反演后的位借模型推算的地震应力降同用地震波资料估算的应力降属同一量级     

Common offset ground penetrating radar data inversion based on ray theory

The authors use the common offset ground penetrating radar(GPR) data inversion based on ray theory to estimate interval velocity and to obtain the relative permittivity. In the ray-tracing based inversion, the input data are the offset distance between antennas, the velocity of the first layer, the pick-up amplitude and re-ference amplitude of each reflection layer. The thickness and velocity of each layer are calculated by this recursive method. Firstly, the horizontal homogeneous layered medium model is established, and the ideal inversion results are obtained. Subsequently, Monte Carlo method is used to establish a randomly undulating homogeneous layered medium model. The common offset GPR data for the built geological model is then simulated by finite-difference time-domain(FDTD). It proved that this ray-tracing based inversion method is feasible for the horizontal layered geological model, even the layered geological model with random undulation. Undulation, represented by RMS height and CL(correlation length), influences the inversion results. Finally, a more complex geological model--pinch-out model was established. In the pinch-out model, the pinch-out interface can be clearly identified, though there is a false anomaly, which will not significantly affect the identification of the underground medium structure.     

一种计算补偿最小二乘正则化参数的最优化方法

在半参数模型补偿最小二乘估计中,如何有效地计算正则化参数是求解的关键。在建立一种新的补偿最小二乘模型的基础上,基于广义交叉核实与L曲线准则,采用最优化方法计算正则化参数。仿真计算表明,本文方法是一种计算正则化参数更为有效的方法。     

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.     

����ɢ�ҵ�ľֲ�n��������ϼ������ʼ���

???????????????????????????????????????ξ??????????????????????????????????????????????????n???????????????????????????????????????????????Ρ????Ρ????????????????????????????????????????????????????????????