本文選題：等效理論　切入點：各向異性　出處：《中國石油大學(華東)》2015年博士論文　論文類型：學位論文

【摘要】：裂縫是特殊巖性儲層(碳酸鹽,頁巖等)的重要儲集空間和運移通道,圍繞裂隙介質的正反演研究對于深入了解地球介質結構、勘探開發(fā)復雜油氣藏和監(jiān)測預報地質災害等具有重要理論意義和現(xiàn)實意義。本文基于裂縫等效理論深入研究了裂隙介質正反演方法,推導了裂縫介質地震波傳播特征(相速度、群速度與偏振方向)的精確與近似表達式,建立了裂縫介質三維qP波波動方程有限差分數(shù)值模擬算法,研究了裂隙參數(shù)的反演策略和方法�；诰�性滑動等效理論的模型參數(shù)化是實現(xiàn)裂隙介質正反演問題的理論基礎。本文采用等效理論將裂隙介質等效為各向異性介質的研究思路,并考慮地下壓力的存在使得裂隙介質中的裂縫通常為高陡傾角裂縫的特點,利用線性滑動與Hudson相結合的等效理論將單組裂縫等效為HTI介質,將兩組正交裂縫等效為正交各向異性介質,將兩組斜交裂縫等效為單斜各向異性介質,推導了三類不同裂隙介質的等效各向異性參數(shù)定量表達式,利用Thomsen參數(shù)對等效后的各向異性參數(shù)進行弱各向異性近似,從而建立起裂隙介質正反演的理論基礎。地震波的傳播特征(速度與偏振方向)是開展裂隙介質正反問題研究的核心方面。地震波傳播速度作為研究地震波傳播規(guī)律和描述介質特征的重要參數(shù),主要包括相速度與群速度。本文利用彈性動力學三個基本方程即廣義虎克定律、運動微分平衡方程及幾何方程推導了三類裂隙介質的Christoffel方程,通過求解Christoffel方程的特征值得到了裂隙介質的相速度表達式,以裂隙介質相速度為基礎,利用Berryman方程與Crampin方程推導了裂隙介質的群速度表達式,通過求解Christoffel方程的特征向量得到了裂隙介質的偏振方向表達式,并利用Thomsen參數(shù)實現(xiàn)了相速度與偏振方向的弱各向異性近似。qP波波動方程地震正演模擬是研究裂隙介質中地震波傳播規(guī)律與地震響應特征的重要方法。本文根據彈性動力學三個基本方程即廣義虎克定律、運動微分平衡方程及幾何方程推導了裂隙介質的彈性波方程,結合Alkhalifah聲學近似假設原理(令橫波速度為零)求取了裂隙介質的qP波波動方程,利用有限差分原理將qP波微分方程轉化為qP波差分方程,采用完全匹配層邊界條件推導了裂隙介質帶有PML邊界的有限差分方程。選用層狀模型和凹陷模型進行了正演模擬驗證了該數(shù)值模擬精度和準確度,分析研究了不同模式裂隙介質的地震響應特征,討論了不同裂縫參數(shù)(裂縫密度、裂縫長度、裂縫寬度與裂縫充填物)與地震波場特征的關系,并進一步利用逆時偏移結果驗證了正演模擬方法的正確性。疊前地震反演是獲取裂隙參數(shù)和認知地下裂縫的有效手段。本文利用裂隙介質彈性波相速度與偏振方向,建立了入射波、反射波與透射波的波函數(shù),結合snell定律與介質分界面處的應力連續(xù)和位移連續(xù)邊界條件,構建了qP波入射裂隙介質彈性波分界面的擬Zoeppritz方程,根據弱各向異性近似與彈性界面的相似近似推導了qP波入射的近似反射系數(shù),分析研究了裂縫參數(shù)的AVOZ特征。在此基礎上利用最小二乘的反演策略實現(xiàn)了基于兩個方位地震數(shù)據反演裂縫參數(shù)的方法,并進一步針對模型數(shù)據通過波形反演方法完成了背景縱波速度的反演研究工作。
[Abstract]:Crack is a special lithologic reservoir (carbonate, shale) in space and the important migration channels around the reservoir, inversion of fractured media research is to understand the earth medium structure, has important theoretical significance and practical significance in the exploration and development of complex oil and gas reservoirs and geological disaster monitoring and forecasting. In this paper, in-depth study of the theory of equivalent crack fracture the medium is derived based on inversion method, seismic wave propagation characteristics of fractured medium (phase velocity, group velocity and polarization direction) of the exact and approximate expressions, establish fractured medium 3D qP wave equation finite difference numerical simulation algorithm, inversion strategy and method to study the fracture parameters. The model parameters of the equivalent linear sliding theory based on the basic theory of the realization of fractured medium and inversion problem. This paper uses the theory of equivalent fracture medium equivalent to the research of anisotropic media, and to consider The presence of pressure makes the fracture medium usually characteristics of high steep crack, using the equivalent linear sliding theory combined with Hudson single group crack equivalent to HTI medium, two groups of orthogonal cracks equivalent to orthotropic medium, two sets of oblique cracks using monoclinic anisotropic media. The quantitative formula of equivalent anisotropic parameters are deduced with three different types of fracture medium, the anisotropic parameter equivalent of weak anisotropy approximation using Thomsen parameters, so as to establish the theoretical basis of fractured medium inversion. The characteristics of seismic wave propagation (velocity and polarization direction) is the core of the research problem and medium cracks. As an important parameter of propagation velocity of seismic wave propagation and describing the media characteristics of seismic waves, including the phase velocity and group velocity. This paper using the elastic dynamics of three The basic equations of the generalized Hooke's law, the differential equilibrium equations and geometric equations deduced the Christoffel equation of the three kind of fractured media, the characteristics of the solution of the Christoffel equation to the fracture medium phase velocity expressions in fractured medium phase velocity as the basis, using Berryman equation and Crampin equation deduced the group velocity expression of fractured media, the the feature vector for solving Christoffel equations has been the polarization direction of the fracture media expression, and the use of Thomsen parameters to achieve weak anisotropic phase velocity and polarization direction of the approximate.QP wave equation seismic forward modeling is an important method for the response characteristics of seismic wave propagation and seismic study of fractured medium. According to the three basic equations of elastic dynamics the generalized Hooke's law, the differential equilibrium equations and geometric equations derived fracture medium elastic wave. The process, combined with the Alkhalifah approximation principle (the acoustic shear wave velocity is zero) for qP wave equation for fractured media, using finite difference qP wave differential equations into difference equations of qP wave principle, the perfectly matched layer boundary conditions are deduced by finite difference fracture medium with PML boundary integral equation is selected. Layered model and depression model for forward simulation to verify the precision and accuracy of the numerical simulation and analysis of response characteristics of different models of fractured medium earthquake, different crack parameters were discussed (crack density, crack length, crack width and crack filling) and seismic wave field characteristics, and further use of the inverse when the migration results verified the correctness of the forward simulation method. The prestack seismic inversion is an effective method to obtain the fracture parameters and fracture cognition. This paper uses the elastic fracture medium wave phase velocity and polarization The establishment of the incident wave direction, wave functions, wave and transmission wave reflection, combined with the stress of Snell's law and the continuous medium interface and continuous displacement boundary conditions, construct the quasi Zoeppritz equation between wave fracture medium elastic wave qP, according to the weak anisotropy approximation and approximation derived approximate reflection coefficient of qP wave incident similar elastic interface, AVOZ to study the characteristics of fracture parameter analysis. On the basis of using the least squares inversion strategy realization method of two azimuth seismic data based on the inversion of fracture parameters, and further the data of model by waveform inversion method to complete the work of inversion of P-wave velocity background.

【學位授予單位】：中國石油大學(華東)
【學位級別】：博士
【學位授予年份】：2015
【分類號】：P631.4

【相似文獻】

相關期刊論文 前10條

1 張茂鋒;劉文明;;裂隙介質中多場耦合理論現(xiàn)狀研究[J];西部探礦工程;2009年10期

2 李寧,陳文玲,張平;動荷作用下非貫通裂隙介質的強度性質[J];自然科學進展;2000年11期

3 桂志先,賀振華,黃德濟;兩種類型裂隙介質的彈性常數(shù)分析[J];物探化探計算技術;2000年02期

4 桂志先,賀振華,黃德濟;裂隙介質彈性性質研究[J];礦物巖石;2001年03期

5 馮增朝;趙陽升;呂兆興;;強隨機分布裂隙介質的二維逾滲規(guī)律研究[J];巖石力學與工程學報;2006年S2期

6 李樹文;高軻;王若男;;基于非連續(xù)裂隙介質網絡流的線素法模型研究[J];河北工程大學學報(自然科學版);2013年04期

7 陳蘊生;韓信;李寧;馬立新;蒲毅彬;廖全榮;;非貫通裂隙介質單軸受力條件下的損傷本構關系探討[J];巖石力學與工程學報;2005年S2期

8 陳蘊生;劉晟鋒;李寧;韓信;;裂隙對非貫通裂隙介質強度與變形特性影響的研究[J];西北農林科技大學學報(自然科學版);2007年07期

9 劉金英,楊天行,徐紅敏,黃繼國;放射性核素在裂隙介質中遷移模型研究綜述[J];世界地質;2000年01期

10 畢貴權;李寧;李國玉;;非貫通裂隙介質中波傳播特性試驗研究[J];巖石力學與工程學報;2009年S1期

相關會議論文 前3條

1 陳蘊生;張平;呂強;;非貫通裂隙介質變形與強度效應的試驗研究[A];巖石力學新進展與西部開發(fā)中的巖土工程問題——中國巖石力學與工程學會第七次學術大會論文集[C];2002年

2 林劍鋒;謝金川;黨海軍;李梅;張繼紅;張海濤;;裂隙介質中α-FeOOH對钚遷移行為的影響研究[A];第十一屆全國核化學與放射化學學術討論會論文摘要集[C];2012年

3 田冬冬;溫聯(lián)星;;利用有限差分方法模擬地震波在含裂隙介質中的傳播[A];中國地球物理2013——第十三專題論文集[C];2013年

相關博士學位論文 前4條

1 秦海旭;基于各向異性理論的裂隙介質正反演方法研究[D];中國石油大學(華東);2015年

2 張鵬;裂隙表面幾何形態(tài)對裂隙介質力學、水力學特性的影響規(guī)律研究[D];西安理工大學;2007年

3 黃勇;多尺度裂隙介質中的水流和溶質運移隨機模擬研究[D];河海大學;2005年

4 蘇銳;低滲透裂隙介質滲透特征評價技術及其應用研究[D];核工業(yè)北京地質研究院;2008年

相關碩士學位論文 前8條

1 鄭東柱;基于孔隙—裂隙滲透效應的注漿擴散分析方法[D];山東大學;2016年

2 錢玉萍;基于唐的孔裂隙介質彈性波動統(tǒng)一理論反演裂隙參數(shù)[D];中國石油大學(華東);2014年

3 畢貴權;裂隙介質中波傳播特性試驗研究[D];西安理工大學;2004年

4 高軻;地下水裂隙介質網絡流的計算方法及其應用研究[D];河北工程大學;2013年

5 韓信;裂隙介質靜、動剪切特性試驗研究[D];西安理工大學;2006年

6 陳蘊生;單軸壓縮條件下非貫通裂隙介質損傷演化特征的試驗研究[D];西安理工大學;2002年

7 劉君;放射性核素在花崗巖裂隙介質中遷移數(shù)值模擬[D];哈爾濱工程大學;2013年

8 謝洪毅;弱滲透裂隙介質深埋長隧洞水文地質模型及其涌水量預測研究[D];河海大學;2006年

，

本文編號：1560899