地面微地震资料震源定位的贝叶斯反演方法
|
摘要
针对地面微地震资料信噪比低、初至拾取不准、速度模型难以准确建立等问题,以及地面微地震资料多条测线测量和浅地表地层速度变化复杂特点,研究了地面微地震资料震源定位的贝叶斯反演方法,把所有测线反演结果设定为一个全概率事件,每条测线反演问题设定为一个划分,讨论利用贝叶斯最大后验方法反演震源位置。在反演时浅部采用横向变速模型,中深部采用水平横向均匀速度模型模型。对目标函数的后验概率密度函数、加权函数后验密度函数、速度参数方差的后验概率密度函数进行理论模型拟合,并取拟合后结果作为估计概率密度。采用极快速模拟退火方法加网格法的混合算法作为搜索方法,以网格算法为先导使搜索落入最优解所在的凸区间,再利用极快速模拟退火算法搜索最优解,这样既可以防止算法收敛于局部极值点,又极大地提高了算法的收敛速度。通过理论模型和实际资料验证了该方法的应用效果,即对随机跳动误差较大初至反演能够保证反演结果的精度。
The surface monitoring microseismic data is characterized by low S/N,inaccurate first-break picking-up,difficult to establish accurate velocity model.Meanwhile,surface monitoring microseismic data is measured by multi survey lines and its shallow surface layers have complex strata velocity.In order to solve the above problems,the Bayesian inversion method for the source positioning of surface monitoring microseismic data was probed.We assume the inversion results of all lines as a total probability and the inversion of every line as a classification,and then discuss the Bayesian maximum posterior method for source position inversion.During inversion,lateral velocity-variable model is adopted in shallow layers and horizon lateral even velocity model in middle-deep layers.The theoretical models are fitted on posterior probability density function of objective function,weighted function and velocity parameter variance.The fitting results are regarded as estimation probability density.Taking the hybrid algorithm of extremely fast simulated annealing grid method as searching method,and grid algorithm as guide to make searching fall into the convex interval of the optimal solution.Then,extremely fast simulated annealing algorithm is used to search the optimal solution.This is to prevent the algorithm converges to a local extreme point,and greatly improve the speed of convergence of the algorithm.The application results on theoretical model and actual data indicate the validation of the method that is to guarantee the accuracy of first-break inversion with large random error.
The surface monitoring microseismic data is characterized by low S/N,inaccurate first-break picking-up,difficult to establish accurate velocity model.Meanwhile,surface monitoring microseismic data is measured by multi survey lines and its shallow surface layers have complex strata velocity.In order to solve the above problems,the Bayesian inversion method for the source positioning of surface monitoring microseismic data was probed.We assume the inversion results of all lines as a total probability and the inversion of every line as a classification,and then discuss the Bayesian maximum posterior method for source position inversion.During inversion,lateral velocity-variable model is adopted in shallow layers and horizon lateral even velocity model in middle-deep layers.The theoretical models are fitted on posterior probability density function of objective function,weighted function and velocity parameter variance.The fitting results are regarded as estimation probability density.Taking the hybrid algorithm of extremely fast simulated annealing grid method as searching method,and grid algorithm as guide to make searching fall into the convex interval of the optimal solution.Then,extremely fast simulated annealing algorithm is used to search the optimal solution.This is to prevent the algorithm converges to a local extreme point,and greatly improve the speed of convergence of the algorithm.The application results on theoretical model and actual data indicate the validation of the method that is to guarantee the accuracy of first-break inversion with large random error.
引文
[1]Stoffa P L,Sen M K.Nonlinear multiparameter opti-mization using genetic algorithms:inversion of plane-wave seismograms[J].Geophysics,1991,56(11):1794-1810
[2]Pei D H,Quirein J A,Cornish B E,et al.Velocity cali-bration for microseismic monitoring:a very fast simula-ted annealing(VFSA)approach for joint-objective opti-mization[J].Geophysics,2009,74(6):47-55
[3]Lisbeth E S.Inversion of Arreal times of microearth-quake sources in the noth sea using a 3-D velocitystructure and prio information:part I[J].MethodBulletin of the Seismological Society of America,1991,81(4):1183-1194
[4]Willis R B,Fontaine J,Paugh L,et al.Geology andgeometry:a review of factors affecting the effective-ness of hydraulic fracturing[R].SPE Eastern Re-gional Meeting,97993-MS,14-16
[5]Chunduru R K,Sen M K,Stoffa P L.Hybrid optimi-zation methods for geophysical inversion[J].Geo-physics,1997,62(4):1196-1207
[6]Phillip W S,House L S.Micro-seismic mapping of aCotton Valley Hydraulic Fracture using decimateddownhole arrays[J].International Exposition andSixty Eighth Annual Meeting,1998,9:13-18
[7]Dumay J,Fournier F.Multivariate statistical analy-ses applied to seismic facies recognition[J].Geo-physics,1988,53(9):1151-1159
[8]De Meersman K,Kendall J M,van der Baan M.The1998 Valhall microseismic data set:an integratedstudy of relocated sources,seismic multiplets,and S-wave splitting[J].Geophysics,2009,74(6):183-195
[9]詹海刚,施平,陈楚群.基于贝叶斯反演理论的海水固有光学特性准分析算法[J].科学通报,2006,51(2):204-210Zhan H G,Shi P,Chen C Q.Based on Bayesian in-version theory of inherent optical properties of sea-water anaysis algorithm[J].Chinese Science Bulle-tion,2006,512(2):204-210
[10]王华忠,方正茂,徐兆涛,等.地震波旅行时计算[J].石油地球物理勘探,1999,34(2):155-163Wang H Z,Fang Z M,Xu Z T,et al.Computation ofseismic travel time[J].Oil Geophysical Prospecting,1999,34(2):155-163
[11]万永革,李鸿吉.遗传算法在确定震源位置中的应用[J].地震地磁观测与研究,1995,16(6):1-7Wan Y G,Li H J.The preliminary study on the seis-mic hypocenter location using genetic algorithms[J].Seismological and Geomagnetic Observation andResearch,1995,16(6):1-7
[12]屈永华,王锦地,刘素红,等.贝叶斯网络支持的地表参数混合反演模式研究[J].遥感学报,2006,10(1):6-13Qu Y H,Wang J D,Liu S H,et al.Study on hybridinversion scheme under Bayesian network[J].Jour-nal of Remote Sensing,2006,10(1):6-13
[13]宋维琪,刘军,陈伟.改进射线追踪算法的微震源反演[J].物探与化探,2008,32(3):275-277Song W Q,Liu J,Chen W.Microearthquake sourceinversion of an improved ray tracing algorithm[J].Geophysical and Geochemical Exploration,2008,32(3):275-277
[2]Pei D H,Quirein J A,Cornish B E,et al.Velocity cali-bration for microseismic monitoring:a very fast simula-ted annealing(VFSA)approach for joint-objective opti-mization[J].Geophysics,2009,74(6):47-55
[3]Lisbeth E S.Inversion of Arreal times of microearth-quake sources in the noth sea using a 3-D velocitystructure and prio information:part I[J].MethodBulletin of the Seismological Society of America,1991,81(4):1183-1194
[4]Willis R B,Fontaine J,Paugh L,et al.Geology andgeometry:a review of factors affecting the effective-ness of hydraulic fracturing[R].SPE Eastern Re-gional Meeting,97993-MS,14-16
[5]Chunduru R K,Sen M K,Stoffa P L.Hybrid optimi-zation methods for geophysical inversion[J].Geo-physics,1997,62(4):1196-1207
[6]Phillip W S,House L S.Micro-seismic mapping of aCotton Valley Hydraulic Fracture using decimateddownhole arrays[J].International Exposition andSixty Eighth Annual Meeting,1998,9:13-18
[7]Dumay J,Fournier F.Multivariate statistical analy-ses applied to seismic facies recognition[J].Geo-physics,1988,53(9):1151-1159
[8]De Meersman K,Kendall J M,van der Baan M.The1998 Valhall microseismic data set:an integratedstudy of relocated sources,seismic multiplets,and S-wave splitting[J].Geophysics,2009,74(6):183-195
[9]詹海刚,施平,陈楚群.基于贝叶斯反演理论的海水固有光学特性准分析算法[J].科学通报,2006,51(2):204-210Zhan H G,Shi P,Chen C Q.Based on Bayesian in-version theory of inherent optical properties of sea-water anaysis algorithm[J].Chinese Science Bulle-tion,2006,512(2):204-210
[10]王华忠,方正茂,徐兆涛,等.地震波旅行时计算[J].石油地球物理勘探,1999,34(2):155-163Wang H Z,Fang Z M,Xu Z T,et al.Computation ofseismic travel time[J].Oil Geophysical Prospecting,1999,34(2):155-163
[11]万永革,李鸿吉.遗传算法在确定震源位置中的应用[J].地震地磁观测与研究,1995,16(6):1-7Wan Y G,Li H J.The preliminary study on the seis-mic hypocenter location using genetic algorithms[J].Seismological and Geomagnetic Observation andResearch,1995,16(6):1-7
[12]屈永华,王锦地,刘素红,等.贝叶斯网络支持的地表参数混合反演模式研究[J].遥感学报,2006,10(1):6-13Qu Y H,Wang J D,Liu S H,et al.Study on hybridinversion scheme under Bayesian network[J].Jour-nal of Remote Sensing,2006,10(1):6-13
[13]宋维琪,刘军,陈伟.改进射线追踪算法的微震源反演[J].物探与化探,2008,32(3):275-277Song W Q,Liu J,Chen W.Microearthquake sourceinversion of an improved ray tracing algorithm[J].Geophysical and Geochemical Exploration,2008,32(3):275-277
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心