基于空间约束的裂缝密度反演方法

doi:10.13810/j.cnki.issn.1000-7210.2024.02.012

Abstract
Figure/Table
References
Related Citation (15)

Download: PDF (8704 KB) HTML (1 KB)
Export: BibTeX | EndNote (RIS)

Abstract Due to the large error in conventional indirect prediction of fracture density using normal weakness and tangential weakness， this paper proposes a direct inversion method for fracture density in carbonate rocks （HTI medium）. Based on the inversion of amplitudes varying with offset and azimuth （AVAZ）， the paper introduces an inverse distance weighted method to optimize the transverse difference operator of traditional anisotropic total variational（ATV） multichannel inversion. The method fully utilizes the correlation between adjacent seismic tracks to further improve the horizontal continuity and stability of the inversion algorithm. Lateral constraints， low-frequency constraints， and sparse constraints are applied to construct the inversion objective function， and the alternate direction multiplier method （ADMM） is adopted to optimize the inversion objective function. With the two advantages of direct inversion of fracture density and inverse distance weighted interpolation， model tests are conducted on the Marmousi II model to verify the effectiveness and noise resistance of the proposed method. Then， the paper further verifies the practical feasibility of the proposed method by applying it to the actual data of carbonate rocks in YM block， Tarim Basin. The proposed method is more accurate and reliable for predicting fracture density and can be widely applied in similar areas.

Key words： anisotropic fracture density inverse distance weighted transverse constraint ADMM

Received: 10 March 2023

PACS:

P631

	Service

	E-mail this article
	Add to my bookshelf
	Add to citation manager
	E-mail Alert
	RSS
	Articles by authors
	YUAN Jingyi
	CAI Zhenzhong
	ZHANG Yintao
	XIE Zhou
	SUN Chong
	ZAHNG Xiaohong

Cite this article:

YUAN Jingyi,CAI Zhenzhong,ZHANG Yintao等. Inversion method of fracture density based on spatial constraints[J]. OGP, 2024, 59(2): 299-310.

URL:

http://www.ogp-cn.com/EN/10.13810/j.cnki.issn.1000-7210.2024.02.012 OR http://www.ogp-cn.com/EN/Y2024/V59/I2/299

[1]	GOODWAY B, VARSEK J, ABACO C. Anisotropic 3D amplitude variation with azimuth (AVAZ) methods to detect fracture prone zones in tight gas resource plays[C].2007 CSPG CSEG Convention, 2007, 590-596.
[2]	LIU E, MARTINEZ A. Seismic Fracture Characterization:Concepts and Practical Applications[M]. EAGE Publications, Houten, Netherlands, 2012.
[3]	陈怀震, 印兴耀, 高建虎, 等. 基于等效各向异性和流体替换的地下裂缝地震预测方法[J].中国科学(地球科学), 2015, 45(5):589-600.CHEN Huaizhen, YIN Xingyao, GAO Jianhu, et al. Seismic inversion for underground fractures detection based on effective anisotropy and fluid substitution[J].Scientia Sinica(Terrae), 2015, 45(5):589-600.
[4]	LI L, ZHANG G, LIU J, et al. Estimation of fracture density and orientation from azimuthal elastic impedance difference through singular value decomposition[J].Petroleum Science, 2021, 18(6):1675-1688.
[5]	XUE J, GU H, CAI C. Model-based amplitude versus offset and azimuth inversion for estimating fracture parameters and fluid content[J].Geophysics, 2017, 82(2):M1-M17.
[6]	PAN X P, LI L, ZHANG G Z. Multiscale frequency-domain seismic inversion for fracture weakness[J].Journal of Petroleum Science and Engineering, 2020, 195(3):107845.
[7]	WANG Y, WANG C, LI J, et al. Frequency-dependent S-wave splitting parameters analysis:a case study from azimuthal PS data, Sanhu area of the Qaidam Basin, China[J].Geophysics, 2019, 84(6):B375-B386.
[8]	刘开元. 三维方位各向异性分析及裂缝检测[D]. 四川成都:成都理工大学, 2011.LIU Kaiyuan. Three-dimensional Research of Azimuthal Anisotropy and Fracture Detection Methods [D].Chengdu University of Technology, Chengdu, Sichuan, 2011.
[9]	陈怀震, 印兴耀, 高成国, 等. 基于各向异性岩石物理的缝隙流体因子AVAZ反演[J].地球物理学报, 2014, 57(3):968-978.CHEN Huaizhen, YIN Xingyao, GAO Chengguo, et al. AVAZ inversion for fluid factor based on fracture anisotropic rock physics theory[J].Chinese Journal of Geophysics, 2014, 57(3):968-978.
[10]	BAKULIN A, GRECHKA V, TSVANKIN I. Estimation of fracture parameters from reflection seismic data- Part I:HTI model due to a single fracture set[J].Geophysics, 2000, 65(6):1788-1802.
[11]	ZHANG F, LI X. Generalized approximations of reflection coefficients in orthorhombic media[J].Journal of Geophysics and Engineering, 2013, 10(5):054004.
[12]	DOWNTON J E, ROURE B. Interpreting azimuthal Fourier coefficients for anisotropic and fracture parameters[J].Interpretation, 2015, 3(3):ST9-ST27.
[13]	ISABEL V G. Fracture Studies from Amplitude Versus Offset and Azimuth and Vertical Seismic Profile Data[D].The University of Edinburgh, Edinburgh, Britain, 2009.
[14]	刘财, 刘宇巍, 冯晅, 等. 基于方位相交的纵波AVA数据运用SVD反演HTI介质裂缝密度[J].吉林大学学报(地球科学版), 2013, 43(5):1655-1662.LIU Cai, LIU Yuwei, FENG Xuan, et al. Invert crack density of HTI media by using SVD based on PP-wave AVA data from crossing seismic survey lines[J].Journal of Jilin University(Earth Science Edition), 2013, 43(5):1655-1662.
[15]	BAKULIN A, GRECHKA V, TSVANKIN I. Estimation of fracture parameters from reflection seismic data-Part II:fractured models with orthorhombic symmetry[J].Geophysics, 2000, 65(6):1803-1817.
[16]	AL-MARZOUG A M, NEVES F A, KIM J J, et al. P-wave anisotropy from azimuthal AVO and velocity estimates using 3D seismic data from Saudi Arabia[J].Geophysics, 2006, 71(2):E7-E11.
[17]	秦海旭, 吴国忱. 基于各向异性梯度的裂缝密度反演[J].地震学报, 2014, 36(6):1062-1074, 1155.QIN Haixu, WU Guochen. Fracture density inversion based on the anisotropic gradient[J].Acta Seismologica Sinica, 2014, 36(6):1062-1074, 1155.
[18]	魏欣伟, 薛姣, 罗霞. 基于OVT域地震数据的叠前AVOA裂缝密度反演[J].石油物探, 2021, 60(5):816-825.WEI Xinwei, XUE Jiao, LUO Xia. Fracture density estimation using an amplitude-versus-offset-and-azimuth inversion based on prestack seismic data in the offset vector tile domain[J].Geophysical Prospecting for Petroleum, 2021, 60(5):816-825.
[19]	CONNOLLY P. Elastic impedance[J].The Leading Edge, 1999, 18(4):438-452.
[20]	MARTIN S, JORGE L. Elastic impedance in weakly anisotropic media[J].Geophysics, 2006, 71(3):D73-D83.
[21]	陈超, 印兴耀, 刘晓晶, 等. 基于方位各向异性的裂缝密度反演方法及应用[J].地球物理学报, 2022, 65(1):371-383.CHEN Chao,YIN Xingyao, LIU Xiaojing,et al. Fracture density inversion based on azimuthal aniso- tropy and its application[J].Chinese Journal of Geophysics, 2022, 65(1):371-383.
[22]	BOYD S, PARIKH N, CHU E, et al. Distributed optimization and statistical learning via the alternating direction method of multipliers[J].Foundations and Trends in Machine Learning, 2011, 3(1):1-122.
[23]	LIU Z Q, XU B, CHENG B, et al. Interpolation parameters in inverse distance-weighted interpolation algorithm on DEM interpolation error[J].Journal of Sensors, 2021, 2021:3535195.
[24]	SCHOENBERG M, DOUMA J. Elastic wave propagation in media with parallel fractures and aligned cracks[J].Geophysical Prospecting, 1988, 36(6):571-590.
[25]	TENG L. Seismic and Rock-Physics Characterization of Fractured Reservoirs[D].Stanford University, Palo Alto, USA, 1998.
[26]	周城, 杨宇勇, 周怀来,等. 基于HTI介质走时反演的各向异性参数建模与AVAZ反演[J]. 石油地球物理勘探, 2023, 58(4):949-960,969.ZHOU Cheng, YANG Yuyong, ZHOU Huailai, et al. Modeling of anisotropic parameters based on HTI medium travel-time inversion and AVAZ inversion[J]. Oil Geophysical Prospecting, 2023, 58(4): 949-960,969.
[27]	周晓越, 甘利灯, 杨昊,等. 利用叠前振幅和速度各向异性的联合反演方法[J]. 石油地球物理勘探, 2020, 55(5):1084-1091.ZHOU Xiaoyue, GAN Lideng, YANG Hao, et al. A joint inversion method using amplitude and velocity anisotropy[J].Oil Geophysical Prospecting,2020,55(5):1084-1091.
[28]	刘晓晶, 陈祖庆, 陈超,等. 方位弹性阻抗傅里叶级数裂缝预测方法[J]. 石油地球物理勘探, 2022, 57(2):423-433.LIU Xiaojing, CHEN Zuqing, CHEN Chao, et al. Fracture prediction method based on Fourier series of azimuthal elastic impedance. Oil Geophysical Prospecting, 2022, 57(2):423-433.
[29]	于晓东, 桂志先, 汪勇,等. 叠前AVAZ裂缝预测技术在车排子凸起的应用[J]. 石油地球物理勘探, 2019, 54(3):624-633.YU Xiaodong, GUI Zhixian, WANG Yong, et al. Prestack AVAZ fracture prediction applied in Che-paizi Uplift[J].Oil Geophysical Prospecting, 2019, 54(3):624-633.
[30]	EUHANNA G, ANDRE T, IMAN S, et al. Optimal parameter selection for the alternating direction method of multipliers (ADMM):Quadratic problems[J]. IEEE Transactions on Automatic Control, 2015, 60(3):644-658.
[31]	ZHAO L, LIN K, WEN X T, et al. Anisotropic total variation pre-stack multitrace inversion based on Lp norm constraint[J]. Journal of Petroleum Science and Engineering. 2023, 220:111212.
[32]	张雨强, 文晓涛, 吴昊,等. 基于Lp拟范数稀疏约束和交替方向乘子算法的波阻抗反演[J]. 石油物探, 2022, 61(5):856-864.ZHANG Yuqiang, WEN Xiaotao, WU Hao, et al. Seismic acoustic impedance inversion using Lp quasi-norm sparse constraint and alternating direction multiplier algorithm[J]. Geophysical Prospecting for Petroleum, 2022, 61(5):856-864.