基于分形先验信息的非线性反演方法能综合利用测井数据和地震信息,在贝叶斯框架下,通过分形高斯噪音算法得到基于分形理论的先验信息,然后根据地震资料构建似然函数,最终利用基于快速模拟退火算法(Very Fast Simulated Annealing,VFSA)...基于分形先验信息的非线性反演方法能综合利用测井数据和地震信息,在贝叶斯框架下,通过分形高斯噪音算法得到基于分形理论的先验信息,然后根据地震资料构建似然函数,最终利用基于快速模拟退火算法(Very Fast Simulated Annealing,VFSA)改进的粒子群优化(Particle Swarm Optimization,PSO)算法(VFSA-PSO)实现后验概率密度的抽样。与确定性反演结果相比,该方法能够有效地融合测井资料中的高频信息,提高反演结果的分辨率,并且目标函数的建立融合了确定性反演中的低频约束,从而得到宽频带的反演结果。数值模拟试验表明:基于分形先验的非线性反演结果与理论模型吻合较好,实际资料的应用效果也证明了该反演方法的有效性。展开更多
A very fast simulated annealing(VFSA) global optimization is used to interpret residual gravity anomaly.Since,VFSA optimization yields a large number of best-fitted models in a vast model space;the nature of uncerta...A very fast simulated annealing(VFSA) global optimization is used to interpret residual gravity anomaly.Since,VFSA optimization yields a large number of best-fitted models in a vast model space;the nature of uncertainty in the interpretation is also examined simultaneously in the present study.The results of VFSA optimization reveal that various parameters show a number of equivalent solutions when shape of the target body is not known and shape factor 'q' is also optimized together with other model parameters.The study reveals that amplitude coefficient k is strongly dependent on shape factor.This shows that there is a multi-model type uncertainty between these two model parameters derived from the analysis of cross-plots.However,the appraised values of shape factor from various VFSA runs clearly indicate whether the subsurface structure is sphere,horizontal or vertical cylinder type structure.Accordingly,the exact shape factor(1.5 for sphere,1.0 for horizontal cylinder and 0.5 for vertical cylinder)is fixed and optimization process is repeated.After fixing the shape factor,analysis of uncertainty and cross-plots shows a well-defined uni-model characteristic.The mean model computed after fixing the shape factor gives the utmost consistent results.Inversion of noise-free and noisy synthetic data as well as field data demonstrates the efficacy of the approach.展开更多
In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position inf...In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.展开更多
文摘A very fast simulated annealing(VFSA) global optimization is used to interpret residual gravity anomaly.Since,VFSA optimization yields a large number of best-fitted models in a vast model space;the nature of uncertainty in the interpretation is also examined simultaneously in the present study.The results of VFSA optimization reveal that various parameters show a number of equivalent solutions when shape of the target body is not known and shape factor 'q' is also optimized together with other model parameters.The study reveals that amplitude coefficient k is strongly dependent on shape factor.This shows that there is a multi-model type uncertainty between these two model parameters derived from the analysis of cross-plots.However,the appraised values of shape factor from various VFSA runs clearly indicate whether the subsurface structure is sphere,horizontal or vertical cylinder type structure.Accordingly,the exact shape factor(1.5 for sphere,1.0 for horizontal cylinder and 0.5 for vertical cylinder)is fixed and optimization process is repeated.After fixing the shape factor,analysis of uncertainty and cross-plots shows a well-defined uni-model characteristic.The mean model computed after fixing the shape factor gives the utmost consistent results.Inversion of noise-free and noisy synthetic data as well as field data demonstrates the efficacy of the approach.
基金Projects(40804027,41074085) supported by the National Natural Science Foundation of ChinaProject(09JJ3048) supported by the Natural Science Foundation of Hunan Province,ChinaProject(200805331082) supported by the Research Fund for the Doctoral Program of Higher Education,China
文摘In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.