Abstract:In view of the shortcomings of available simultaneous joint inversion methods,this paper pre-sents a simultaneous inversion method of gravity and seismic data based on variable density-velocity relation,derives the formulas for getting joint solutions of density model and velocity model in details,and presents the updating algorithm of density-velocity relation using higher-order statistics.During inversion,the initial relation is updated according to previous and present density and velocity models.Model applications show that the me-thod reduces the adverse effects resulted by an inaccurate initial relation.Practical applications show that the velocity model established by this method can improve seismic imaging quality,and has a great application potential in the area less explored.
相鹏, 王金铎, 谭绍泉, 陈学国. 一种变密度—速度关系的重力与地震同步联合反演方法[J]. 石油地球物理勘探, 2020, 55(3): 686-693.
XIANG Peng, WANG Jinduo, TAN Shaoquan, CHEN Xueguo. Simultaneous and joint inversion of gravity and seismic data based on variable density-velocity relation. Oil Geophysical Prospecting, 2020, 55(3): 686-693.
徐桂芬,何展翔,石艳玲,等.电磁-地震联合研究深层火山岩储层——以辽东凹陷火山岩勘探为例[J].石油地球物理勘探,2019,54(4):937-946.XU Guifen,HE Zhanxiang,SHI Yanling,et al.Deep volcanic reservoir exploration with jointed electromagnetic and seismic data:an example of volcanic investigation in Liaodong Sag[J].Oil Geophysical Prospecting,2019,54(4):937-946.
[2]
潘豪杰,张研,李红兵.基于贝叶斯理论的天然气水合物储层弹性-电性数据联合反演[J].石油地球物理勘探,2018,53(3):568-577.PAN Haojie,ZHANG Yan,LI Hongbing.Joint inversion of elastic-electrical data for gas hydrate reservoirs based on Bayesian theory[J].Oil Geophysical Prospecting,2018,53(3):568-577.
[3]
李华东,于鹏,刘振友,等.基于随机分布共网格模型的重磁电震联合反演技术及应用[J].石油地球物理勘探,2015,50(4):742-748.LI Huadong,YU Peng,LIU Zhenyou,et al.Joint inversion of gravity,magnetotelluric and seismic data based on common gridded model with random distri-butions[J].Oil Geophysical Prospecting,2015,50(4):742-748.
[4]
李德春,杨书江,胡祖志,等.三维重磁电震资料的联合解释——以库车大北地区山前砾石层为例[J].石油地球物理勘探,2012,47(2):353-359.LI Dechun,YANG Shujiang,HU Zuzhi,et al.Integrated interpretation of 3D gravity,magnetic,electromagnetic and seismic data:a case study of conglome-rate mass investigation in piedmont area of Kuche Depression[J].Oil Geophysical Prospecting,2012,47(2):353-359.
[5]
Vozoff K,Jupp D L B.Joint inversion of geophysical data[J].Geophysical Journal of the Royal Astronomical Society,1975,42(3):977-991.
[6]
Gomez-Trevino E,Edwards R N.Electromagnetic sounding in the sedimentary basin of southern Ontario:a case history[J].Geophysics,1983,48(3):311-330.
[7]
Raiche A P,Jupp D L B,Rutter H,et al.The joint use of coincident loop transient electromagnetic and Schlumberger sounding to resolve layered structures[J].Geophysics,1985,50(10):1618-1627.
[8]
Lines L,Schultz A K,Treitel S.Cooperative inversion of geophysical data[J].Geophysics,1988,53(1):8-2.
[9]
Dobroka,Gyulai A,Ormos T,et al.Joint inversion of seismic and geo-electric data recorded in an underground coal mine[J].Geophysical Prospecting,1991,39(5):643-665.
[10]
Hering A,Misiek R.A joint inversion algorithm to process geoelectric and surface wave seismic data,Part 1:basis idea[J].Geophysical Prospecting,1995,43(2):135-156.
[11]
Abubakar A,van den Berg P M.Joint inversion of electrode and induction logging data[C].SEG Technical Program Expanded Abstracts,2000:2166-2168.
[12]
Bosch M,McGaughey J.Joint inversion of gravity and magnetic data under lithologic constraints[J].The Leading Edge,2001,20(8):877-881.
[13]
阎汉杰,阎红,李云平.重磁震统计推断联合建模反演技术研究与应用[J].石油大学学报(自然科学版),2003,27(3):36-39.YAN Hanjie,YAN Hong,LI Yunping.Application of gravity-magnetic-seismic combination inversion mode-ling technology based on statistical inference[J].Journal of University of Petroleum (Edition of Natural Science),2003,27(3):36-39.
[14]
Colombo D,de Stefano M.Geophysical modeling via simultaneous joint inversion of seismic,gravity,and electromagnetic data:Application to prestack depth imaging[J].The Leading Edge,2007,26(3):326-331.
[15]
Meju M A.Joint inversion of TEM and distorted MT soundings:some effective practical considerations[J].Geophysics,1996,61(1):56-65.
[16]
Meju M A,Fontes S L,Oliveira M F B.Regional aquifer mapping using combined VES-TEM-AMT/EMAP methods in the semiarid eastern margin of Parnaiba Basin,Brazil[J].Geophysics,1999,64(2):337-356.
[17]
Zhang G L,Xiang P,Wang J D,et al.Geologic mode-ling study of foothill area of the Junggar Basin rim[J].Interpretation,2018,6(4):SM19-SM26.
[18]
Birch F.The velocity of compressional waves in rocks to 10 Kilobars[J].Journal of Geophysical Research Atmospheres,1961,66(7):2199-2224.
[19]
Archie G E.The electrical resistivity log as an aid in determining some reservoir characteristics[J].Tran-sactions of AIME,1942,146(1):54-62.
[20]
Gallardo L A,Meju M A.Characterization of heterogeneous near-surface materials by joint 2D inversion of DC resistivity and seismic data[J].Geophysical Research Letters,2003,30(13):1658-1661.
[21]
Gallardo L A,Meju M A.Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints[J].Journal of Geophysical Research Atmospheres,2004,109(B3):B03311.
[22]
相鹏.一种改进的二维MT预条件非线性共轭梯度反演方法[J].中国石油大学学报(自然科学版),2014,38(4):42-49.XIANG Peng.Two-dimensional MT inversion method based on an improved preconditioned nonlinear conjugate gradient algorithm[J].Journal of China University of Petroleum(Edition of Natural Science),2014,38(4):42-49.
[23]
Delyon B,Lavielle M,Moulines E.Convergence of a stochastic approximation version of the EM algorithm[J].Annals of Statistics,1999,27(1):94-128.
[24]
Zhdanov M S.Geophysical Inverse Theory and Regularization Problem[M].Elsevier,Amsterdam,Netherlands,2002.