一种高效的基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演策略

doi:10.6038/cjg20170927

摘要
图/表
参考文献
相关文章 (15)

全文: PDF (5555 KB) HTML (1 KB)
输出: BibTeX | EndNote (RIS)

摘要基于交叉梯度结构约束的不同数据的联合反演可以提高地球物理成像的可靠度，但是由于不同观测数据对地下介质的灵敏度不同以及不同反演系统灵敏度矩阵元素的数值可能存在较大的差异，导致构建和求解联合反演系统存在很大的挑战.针对以上问题，本文提出一种新的基于单独反演模型更新量与交叉梯度结构约束相结合的联合反演策略.该策略利用单独反演系统分别确定出各个模型的更新量，然后利用它们约束交叉梯度系统的反演，得出新的模型更新量.通过这样的方式，有效实现了数据拟合与结构约束的平衡，实现了基于交叉梯度联合反演的目标.新的联合反演策略不需要对原来单独反演程序修改即可实现联合反演，减小了联合反演实现的难度，极大地提高了联合反演的易实现性，而且避免了联合反演矩阵存储及结构过于复杂难以求解的问题.基于新的策略，本文首次实现了基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演.合成模型测试表明，与单一成像相比，联合成像减少了地震走时反演中出现的干扰异常并提高了电阻率反演的分辨率.

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	高级
	张海江
	方洪健
	李楠

关键词 ：交叉梯度, 模型更新量, 结构相似, 地震走时成像, 电阻率成像, 联合反演

Abstract：Joint inversion of different geophysical data based on cross-gradient structure constraint can improve the reliability of the imaging result. However, because different geophysical observations have varied sensitivities to the Earth medium properties, which are associated with different inversion systems, it is a big challenge to construct and resolve the joint inversion system. To solve the above problems, we propose a new joint inversion strategy by combining the model updates from separate inversions and cross-gradient structure constraint. This strategy first determines the model updates from separate inversions and then utilizes them to constrain the cross-gradient inversion system to get the new model updates. In this way, we can effectively implement the cross-gradient based joint inversion and realize the balance between fitting separate data types and increasing structure similarity between the two models. The new inversion strategy does not need to modify the individual inversion code to realize the joint inversion, making it much easier to conduct the cross-gradient based joint inversion. In addition, it also overcomes the difficulties in storing and inverting large joint inversion matrix. Based on the new strategy, for the first time we implement 3-D joint inversion of seismic travel-time and DC electrical resistivity data based on cross-gradient structure constraint. The test on a synthetic model shows that compared with separate inversions, the joint inversion method can effectively reduce the artifacts shown in seismic travel time inversion results and significantly improve the resolution of the resistivity inversion.

Key words： Cross-gradient Model update Structural similarity Seismic travel time tomography Resistivity tomography Joint inversion

收稿日期: 2016-11-15

PACS:

P631

基金资助:国家自然科学基金（41474039，41704039，51504250），中国矿业大学煤炭资源与安全开采国家重点实验室开放基金（SKLCRSM16KF09），中国科学技术大学青年创新基金（WK2080000089）联合资助.

通讯作者: 张海江,男,1973年生,教授,博士生导师,主要从事发展新的地震成像算法和联合地球物理成像算法以及对断层、火山、俯冲带结构的成像,油气矿产开发过程诱发的微地震等研究.E-mail:zhang11@ustc.edu.cn E-mail: zhang11@ustc.edu.cn

作者简介: 高级,男,1983年生,博士,主要从事地球物理联合成像研究.E-mail:gaoji617@mail.ustc.edu.cn

链接本文:

http://www.geophy.cn/CN/10.6038/cjg20170927 或 http://www.geophy.cn/CN/Y2017/V60/I9/3628

Aster R C, Borchers B, Thurber C H. 2011. Parameter Estimation and Inverse Problems. Boston:Academic Press.
Bellmunt F, Marcuello A, Ledo J, et al. 2012. Time-lapse cross-hole electrical resistivity tomography monitoring effects of an urban tunnel. Journal of Applied Geophysics, 87:60-70.
Bennington N L, Zhang H J, Thurber C H, et al. 2015. Joint inversion of seismic and magnetotelluric data in the Parkfield region of California using the normalized cross-gradient constraint. Pure and Applied Geophysics, 172(5):1033-1052.
Bosch M. 1999. Lithologic tomography:From plural geophysical data to lithology estimation. Journal of Geophysical Research:Solid Earth, 104(B1):749-766.
Calvetti D, Morigi S, Reichel L, et al. 2000. Tikhonov regularization and the L-curve for large discrete ill-posed problems. Journal of Computational and Applied Mathematics, 123(1-2):423-446.
Cao Q, Xiang W, Jia H L, et al. 2013. Application of cross-hole ultra-density resistivity method to detection of spherically weathered granite. Journal of Engineering Geology (in Chinese), 21(5):730-735.
Carcione J M, Gei D, Picotti S, et al. 2012. Cross-hole electromagnetic and seismic modeling for CO₂ detection and monitoring in a saline aquifer. Journal of Petroleum Science and Engineering, 100:162-172.
Duan C L, Yan C H, Xu B T, et al. 2013. The application of cross-hole seismic CT method in the karst cave exploration of metro engineering construction. Geological Review (in Chinese), 59(6):1242-1248.
Ehosioke S I, Fechner T. 2014. Application of cross-hole seismic tomography in characterization of heterogeneous aquifers.//Near Surface Geoscience 2014-20th European Meeting of Environmental and Engineering Geophysics. EAGE. Extended Abstract, doi:10.3997/2214-4609.20142018.
Fregoso E, Gallardo L A. 2009. Cross-gradients joint 3D inversion with applications to gravity and magnetic data. Geophysics, 74(4):L31-L42.
Gallardo L A, Meju M A. 2004. Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints. Journal of Geophysical Research:Solid Earth, 109(B3):B03311, doi:10.1029/2003JB002716.
Gallardo L A. 2007. Multiple cross-gradient joint inversion for geospectral imaging. Geophysical Research Letters, 34(19):L19301, doi:10.1029/2007GL030409.
Gao J, Zhang H J. 2016. Two-dimensional joint inversion of seismic velocity and electrical resistivity using seismic travel times and full channel electrical measurements based on alternating cross-gradient structural constraint. Chinese Journal of Geophysics (in Chinese), 59(11):4310-4322, doi:10. 6038/cjg20161131.
Greve A K, Acworth R I, Kelly B F J. 2012. 3D Cross-hole resistivity tomography to monitor water percolation during irrigation on cracking soil. Soil Research, 49(8):661-669.
Guo Y F, Hu X J, You J M. 2015. The exploitation of straight and curve tracing method in boulders prospected by seismic cross-hole tomography. Chinese Journal of Engineering Geophysics (in Chinese), 12(3):361-366.
Haber E, Oldenburg D. 1997. Joint inversion:a structural approach. Inverse Problems, 13(1):63-77.
Hu X Y, Yang D K, Liu S H, et al. 2006. The developing trends of environmental and engineering geophysics. Progress in Geophysics (in Chinese), 21(2):598-604, doi:10.3969/j.issn.1004-2903.2006.02.041.
Korenaga J, Holbrook W S, Singh S C, et al. 1997. Natural gas hydrates on the southeast US margin:Constraints from full waveform and travel time inversions of wide-angle seismic data. Journal of Geophysical Research-All Series, 102(B7):15345-15365.
Leontarakis K, Apostolopoulos G V. 2013. Model Stacking (MOST) technique applied in cross-hole ERT field data for the detection of Thessaloniki ancient walls' depth. Journal of Applied Geophysics, 93:101-113.
Li S C, Su M X, Xue Y G, et al. 2014. Study on computed tomography of cross-hole resistivity in urban subway geological prediction. Chinese Journal of Rock Mechanics and Engineering (in Chinese), 33(5):913-920.
Li X C, Liu Z Y, Liu B, et al. 2015. Boulder detection method for metro shield zones based on ross-hole resistivity tomography and its physical model tests. Chinese Journal of Geotechnical Engineering (in Chinese), 37(3):446-457.
Li Z X, Tan H D, Fu S S, et al. 2015. Two-dimensional synchronous inversion of TDIP with cross-gradient constraint. Chinese Journal of Geophysics (in Chinese), 58(12):4718-4726, doi:10.6038/cjg20151232.
Lines L R, Schultz A K, Treitel S. 1988. Cooperative inversion of geophysical data. Geophysics, 53(1):8-20.
McGillivray P R, Oldenburg D W. 1990. Methods for calculating Fréchet derivatives and sensitivities for the non-linear inverse problem:a comparative study. Geophysical Prospecting, 38(5):499-524.
Menu F, Kepic A, Greenwood A. 2015. Crosshole seismic reflection to delineate Rosebery type volcanogenic massive sulphide ore lens.//85th Annual Meeting. Society of Exploration Geophysicists. Expanded Abstracts.
Moorkamp M, Heincke B, Jegen M, et al. 2011. A framework for 3-D joint inversion of MT, gravity and seismic refraction data. Geophysical Journal International, 184(1):477-493.
Morency C, Luo Y, Tromp J. 2011. Acoustic, elastic and poroelastic simulations of CO₂ sequestration crosswell monitoring based on spectral-element and adjoint methods. Geophysical Journal International, 185(2):955-966.
Ogunbo J N, Zhang J. 2014. Joint seismic traveltime and TEM inversion for near surface imaging.//84th Annual Meeting. Society of Exploration Geophysicists. Expanded Abstracts.
Osher S, Fedkiw R. 2006. Level Set Methods and Dynamic Implicit Surfaces. New York:Springer Science & Business Media.
Paige C C, Saunders M A. 1982. LSQR:An algorithm for sparse linear equations and sparse least squares. ACM Transactions on Mathematical Software (TOMS), 8(1):43-71.
Peng M. 2012. Joint inversion of magnetotelluric and teleseismic data[Ph. D. thesis] (in Chinese). Beijing:China University of Geosciences (Beijing).
Pepe P, Martimucci V, Parise M. 2015. Geological and geophysical techniques for the identification of subterranean cavities.//Lollino G, Manconi A, Guzzetti F eds. Engineering Geology for Society and Territory. Cham:Springer International Publishing, 483-487.
Pidlisecky A, Haber E, Knight R. 2007. RESINVM3D:A 3D resistivity inversion package. Geophysics, 72(2):H1-H10.
Qu C, Yang W C, Yu C Q. 2010. An application of high resolution P and S wave crosshole tomography in understanding the reservoir features. Chinese Journal of Geophysics (in Chinese), 53(12):2944-2954, doi:10.3969/j.issn.0001-5733.2010.12.018.
Rovetta D, Colombo D, Curiel E S, et al. 2013. 3D seismic-gravity simultaneous joint inversion for near surface velocity estimation.//75th EAGE Conference & Exhibition Incorporating. SPE.
Sánchez A M L, Delgado L A G. 2015. 2D cross-gradient joint inversion of magnetic and gravity data across the Capricorn Orogen in Western Australia.//85th Annual Meeting. Society of Exploration Geophysicists. Expanded Abstracts.
Sambridge M S. 1990. Non-linear arrival time inversion:constraining velocity anomalies by seeking smooth models in 3-D. Geophysical Journal International, 102(3):653-677.
Schmidt-Hattenberger C, Bergmann P, Labitzke T, et al. 2014. CO₂ Migration monitoring by means of electrical resistivity tomography (ERT)-review on five years of operation of a permanent ERT system at the Ketzin pilot site. Energy Procedia, 63:4366-4373.
Tikhonov A N, Arsenin V I A. 1977. Solutions of Ill-Posed Problems. New York:Winston.
Um E S, Commer M, Newman G A. 2014. A strategy for coupled 3D imaging of large-scale seismic and electromagnetic data sets:Application to subsalt imaging. Geophysics, 79(3):ID1-ID13.
Vidale J. 1988. Finite-difference calculation of travel times. Bulletin of the Seismological Society of America, 78(6):2062-2076.
Vidale J E. 1990. Finite-difference calculation of traveltimes in three dimensions. Geophysics, 55(5):521-526.
Wang H, Li Y G, Chen C. 2015. 3D joint inversion of gravity gradiometry and magnetic data in spherical coordinates with the cross-gradient constraint.//85th Annual Meeting. Society of Exploration Geophysicists. Expanded Abstracts.
Wang J C, Shi X M, Wang F F, et al. 2012. The physical experiment research in laboratory of cross-hole electric resistivity tomography for detecting boulders. Computerized Tomography Theory and Application (in Chinese), 21(4):647-657.
Wu X P, Xu G M. 2000. Study on 3D resistivity inversion using conjugate gradient method. Chinese Journal of Geophysics (in Chinese), 43(3):420-427, doi:10.3321/j.issn:0001-5733.2000.03.016.
Zhang F J, Juhlin C, Cosma C, et al. 2012. Cross-well seismic waveform tomography for monitoring CO2 injection:a case study from the Ketzin Site, Germany. Geophysical Journal International, 189(1):629-646.
Zhang W J, Li X C, Su M X, et al. 2014. Detection method of karst caves in city subway based on the cross-hole resistivity tomography. Journal of Shandong University (Engineering Science) (in Chinese), 44(3):75-82.
Zhou J J, Meng X H, Guo L H, et al. 2015. Three-dimensional cross-gradient joint inversion of gravity and normalized magnetic source strength data in the presence of remanent magnetization. Journal of Applied Geophysics, 119:51-60.
Zhou L F. 2012. Two dimensional joint inversion of MT and seismic data[Master's thesis] (in Chinese). Beijing:China University of Geosciences (Beijing).
曹权, 项伟, 贾海梁等. 2013. 跨孔超高密度电阻率法在球状风化花岗岩体探测中的应用. 工程地质学报, 21(5):730-735.
段成龙, 阎长虹, 许宝田等. 2013. 跨孔地震CT技术在地铁工程施工溶洞探测方面的应用. 地质论评, 59(6):1242-1248.
高级, 张海江. 2016. 基于交叉梯度交替结构约束的二维地震走时与全通道直流电阻率联合反演. 地球物理学报, 59(11):4310-4322, doi:10. 6038/cjg20161131.
郭云峰, 胡晓娟, 游敬密. 2015. 地震波跨孔CT探测孤石中直线与曲线追踪法的应用. 工程地球物理学报, 12(3):361-366.
胡祥云, 杨迪坤, 刘少华等. 2006. 环境与工程地球物理的发展趋势. 地球物理学进展, 21(2):598-604, doi:10.3969/j.issn. 1004-2903.2006.02.041.
李术才, 苏茂鑫, 薛翊国等. 2014. 城市地铁跨孔电阻率CT超前地质预报方法研究. 岩石力学与工程学报, 33(5):913-920.
李术才, 刘征宇, 刘斌等. 2015. 基于跨孔电阻率 CT 的地铁盾构区间孤石探测方法及物理模型试验研究. 岩土工程学报, 37(3):446-457.
李兆祥, 谭捍东, 付少帅等. 2015. 基于交叉梯度约束的时间域激发极化法二维同步反演. 地球物理学报, 58(12):4718-4726, doi:10.6038/cjg20151232.
彭淼. 2012. 大地电磁与天然地震数据联合反演研究[博士论文]. 北京:中国地质大学(北京).
瞿辰, 杨文采, 于常青. 2010. 井间高分辨率纵横波层析成像研究井间油藏. 地球物理学报, 53(12):2944-2954, doi:10.3969/j.issn.0001-5733.2010.12.018.
王俊超, 师学明, 万方方等. 2012. 探测孤石高阻体的跨孔电阻率CT水槽物理模拟实验研究. CT理论与应用研究, 21(4):647-657.
吴小平, 徐果明. 2000. 利用共轭梯度法的电阻率三维反演研究. 地球物理学报, 43(3):420-426, doi:10.3321/j.issn:0001-5733.2000.03.016.
张文俊, 李术才, 苏茂鑫等. 2014. 基于井间电阻率成像的城市地铁溶洞探测方法. 山东大学学报(工学版), 44(3):75-82.
周丽芬. 2012. 大地电磁与地震数据二维联合反演研究[硕士论文]. 北京:中国地质大学(北京).