聲波測井的數(shù)值模擬課件

1、全幼黎Stanford UniversityOctober, 2007科研活動的簡要回顧內(nèi)容提要地震波衰減的層析成像 (Stanford)聲波測井的數(shù)值模擬 (Stanford)油田殘余油強化開采的地震監(jiān)測 (Texaco)油層的地質(zhì)統(tǒng)計描述 (Texaco) 地質(zhì)統(tǒng)計軟件 (SMT)斷層的自動識別 (SMT)用生物芯片發(fā)現(xiàn)蛋白質(zhì)生物標志物 (LumiCyte)差式共振聲譜儀 (Stanford)地震監(jiān)測地下存儲的CO2 (Stanford)環(huán)形換能器超聲成像 (Stanford/LANL)油田油儲的綜合監(jiān)測 (Stanford)1. 地震波衰減的層析成像 Stanford Universit

2、y, 1991-1996頻移法是我們提出和發(fā)展的一種簡單實用的技術(shù)。這一方法自1997發(fā)表后得到了廣泛的應(yīng)用。比如說，最近一年的 Geophysics 雜志有6篇文章引用該論文。今年的SEG Expended Abstracts有7篇文章引用了這篇論文。Acoustic SourcesAcoustic ReceiversLower AbsorptionHigher AbsorptionHigher frequencyWaveform Lower frequencyWaveform 波形和衰減用中心頻率移動法估計衰減Medium ResponseH(f)Incident WaveS(f)Tran

3、smitted WaveR(f)=S(f)H(f)H(f)=exp-faodlxxxxxxxxxxxxooooooooooooooSourcesReceiversVf=5 kft/sQf=2070 ft100 ftVp=11.8 kft/sVs=6.9 kft/sQp=30Vp=12 kft/sVs=7 kft/sQp=602565(a) Original model(b) ReconstructionCrosswell geometry, RT method for modeling理論數(shù)據(jù)的例子野外數(shù)據(jù)實例The attenuation and velocity tomography fo

4、r a west Texas seismic survey. The geological structure in this area is complex. The main features are two low-velocity/low-Q zones indicated by A and B, respectively. Area B is interpreted as a carbonate mound or reef.2. 聲波測井的數(shù)值模擬 Stanford University, 1991-1996我們提出的一種很有效的計算油井中波場的方法。這一方法發(fā)表后許多人來電子郵件索

5、取程序。有幾個大石油公司在使用我們的程序。 去年我把這一方法應(yīng)用到一種新測井數(shù)據(jù)的解釋，并申請到一項美國專利。與井有關(guān)的地震觀測手段ooVertical seismic profilingoCross-boreholeprofilingSingleboreholeprofilingo SonicloggingFluid-filledboreholexSourceFluid-filledborehole Modeling of Waves in a Borehole Modeling of Waves in a Borehole 測井模型 (縱切面)Radially layered modelz

6、rzrComplex radially symmetric model Modeling of Waves in a Borehole 測井模型 (橫切面)使用的方法 Modeling of Waves in a Borehole l 廣義反射透射系數(shù)法 (穩(wěn)定, 速度快)l 歸一化的數(shù)學處理 (可以模擬高頻和厚層的問題)layers Modeling of Waves in a Borehole 計算實例簡單充水裸井有損傷的裸井oSourceReceivers100 m150 m10 mFormation IFormation IIBoreholexxxxxxxxxx從遠方柱面上的反射 (模

7、型)ModelVp=3 km/sVs=1.8 km/sVp=1.9 km/sVs=1.4 km/sSeismograms (fo = 800 Hz)101501050100150Time (ms)source - receiver offset (m)從遠方柱面上的反射 (地震波)井間地震例子ooooooooooooooo30 mReceiversSourcesCased boreholeVp=5 km/sVs=2.9 km/sVp=5.8 km/sVs=3.3 km/s100 mFormation (a) Model. There is a fault in the formation(b)

8、 A common receiver gatherxxxxxxxxxxxxxxx30 m102030400Time (ms)0S-Roffset(m)100一種新測井數(shù)據(jù)的模擬Layered borehole model(b) The pressure field recorded at top, middle and bottom(a)(b)3. 油田殘余油強化開采的地震監(jiān)測Time-lapse CO2 Monitoring DVp (P Wave Velocity Change Vp2-Vp1) Tomogram340039004400490054000100050015000500100

9、0150020002500Depth (ft)North (ft)East (ft)94 - Source well. 106, 81 - Receiver wells. Well separations - 106-94:1216 ft, 81-94: 1479 ft.71, 72, 82 - CO2 injection wells35, 90, 101, 78, 290 - Production wells.GR logs are shown at source/receiver wells.-DVp/Vp (%)Top of Lower San Andres1.504. 油層的地質(zhì)統(tǒng)計描

10、述 (QHD 32-6, 1997)用測井數(shù)據(jù)計算孔隙度剖面用地震數(shù)據(jù)計算的阻抗剖面綜合地震和測井數(shù)據(jù)得到的阻抗剖面油層的三維空隙度模型5. 地質(zhì)統(tǒng)計軟件 地質(zhì)統(tǒng)計研究空間上起伏變化的現(xiàn)象. 地質(zhì)統(tǒng)計提供一些用于建立能描述空間變化規(guī)律的工具.地質(zhì)統(tǒng)計 這些工具包括 histogram, covariance, variogram, kriging, 等等.6. 斷層的自動識別 7. 用生物芯片發(fā)現(xiàn)蛋白質(zhì)生物標志物 質(zhì)譜儀兩組質(zhì)譜數(shù)據(jù)以及能把它們分開的幾個蛋白質(zhì)8. 差式共振聲譜儀 Differential Resonance Acoustic Spectroscopy (DARS)圓筒共振腔E

11、lastic constants, Q-value and permeability may be estimated from these spectra DARS 共振譜9. 地震監(jiān)測注入地下的CO2把二氧化碳注入地下CO2SealFaultInjection Well為什么要進行監(jiān)測探測CO2移動和泄漏的路徑以保證公眾的安全了解被注入的CO2時間和空間分布以便更有效地進行注入操作10. 環(huán)形換能器超聲成像 環(huán)形換能器和仿真模型仿真模型的重建Given ModelPicked Travel TimeReconstructed Velocity ModelReconstructed Refl

12、ectivity Model11. 油田油儲的綜合監(jiān)測An integrated monitoring study: (a) Creation of a geology model for the reservoir (b) Flow simulation (c) Conversion from flow properties to seismic models (d) Seismic simulation for a given monitoring survey (e) Seismic imaging (Migration, Tomography) (f) Comparison with field data (g) Updating of reservoir model with inversion results Work FlowField DataNoYesDoes synthetic data match field data? Inversion Acoustic pr