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

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

2、y, 1991-1996頻移法是我們提出和發(fā)展的一種簡(jiǎn)單實(shí)用的技術(shù)。這一方法自1997發(fā)表后得到了廣泛的應(yīng)用。比如說，最近一年的 Geophysics 雜志有6篇文章引用該論文。今年的SEG Expended Abstracts有7篇文章引用了這篇論文。Acoustic SourcesAcoustic ReceiversLower AbsorptionHigher AbsorptionHigher frequencyWaveform Lower frequencyWaveform 波形和衰減用中心頻率移動(dòng)法估計(jì)衰減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ù)實(shí)例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. 聲波測(cè)井的數(shù)值模擬 Stanford University, 1991-1996我們提出的一種很有效的計(jì)算油井中波場(chǎng)的方法。這一方法發(fā)表后許多人來電子郵件索

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

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

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)從遠(yuǎn)方柱面上的反射 (地震波)井間地震例子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一種新測(cè)井?dāng)?shù)據(jù)的模擬Layered borehole model(b) The pressure field recorded at top, middle and bottom(a)(b)3. 油田殘余油強(qiáng)化開采的地震監(jiān)測(cè)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)計(jì)描

10、述 (QHD 32-6, 1997)用測(cè)井?dāng)?shù)據(jù)計(jì)算孔隙度剖面用地震數(shù)據(jù)計(jì)算的阻抗剖面綜合地震和測(cè)井?dāng)?shù)據(jù)得到的阻抗剖面油層的三維空隙度模型5. 地質(zhì)統(tǒng)計(jì)軟件 地質(zhì)統(tǒng)計(jì)研究空間上起伏變化的現(xiàn)象. 地質(zhì)統(tǒng)計(jì)提供一些用于建立能描述空間變化規(guī)律的工具.地質(zhì)統(tǒng)計(jì) 這些工具包括 histogram, covariance, variogram, kriging, 等等.6. 斷層的自動(dòng)識(shí)別 7. 用生物芯片發(fā)現(xiàn)蛋白質(zhì)生物標(biāo)志物 質(zhì)譜儀兩組質(zhì)譜數(shù)據(jù)以及能把它們分開的幾個(gè)蛋白質(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)測(cè)注入地下的CO2把二氧化碳注入地下CO2SealFaultInjection Well為什么要進(jìn)行監(jiān)測(cè)探測(cè)CO2移動(dòng)和泄漏的路徑以保證公眾的安全了解被注入的CO2時(shí)間和空間分布以便更有效地進(jìn)行注入操作10. 環(huán)形換能器超聲成像 環(huán)形換能器和仿真模型仿真模型的重建Given ModelPicked Travel TimeReconstructed Velocity ModelReconstructed Refl

12、ectivity Model11. 油田油儲(chǔ)的綜合監(jiān)測(cè)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