声波测井的数值模拟课件

1、全幼黎Stanford UniversityOctober, 2007科研活动的简要回顾内容提要地震波衰减的层析成像 (Stanford)声波测井的数值模拟 (Stanford)油田残余油强化开采的地震监测 (Texaco)油层的地质统计描述 (Texaco) 地质统计软件 (SMT)断层的自动识别 (SMT)用生物芯片发现蛋白质生物标志物 (LumiCyte)差式共振声谱仪 (Stanford)地震监测地下存储的CO2 (Stanford)环形换能器超声成像 (Stanford/LANL)油田油储的综合监测 (Stanford)1. 地震波衰减的层析成像 Stanford Universit

2、y, 1991-1996频移法是我们提出和发展的一种简单实用的技术。这一方法自1997发表后得到了广泛的应用。比如说，最近一年的 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理论数据的例子野外数据实例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. 声波测井的数值模拟 Stanford University, 1991-1996我们提出的一种很有效的计算油井中波场的方法。这一方法发表后许多人来电子邮件索

5、取程序。有几个大石油公司在使用我们的程序。 去年我把这一方法应用到一种新测井数据的解释，并申请到一项美国专利。与井有关的地震观测手段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 广义反射透射系数法 (稳定, 速度快)l 归一化的数学处理 (可以模拟高频和厚层的问题)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一种新测井数据的模拟Layered borehole model(b) The pressure field recorded at top, middle and bottom(a)(b)3. 油田残余油强化开采的地震监测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. 油层的地质统计描

10、述 (QHD 32-6, 1997)用测井数据计算孔隙度剖面用地震数据计算的阻抗剖面综合地震和测井数据得到的阻抗剖面油层的三维空隙度模型5. 地质统计软件 地质统计研究空间上起伏变化的现象. 地质统计提供一些用于建立能描述空间变化规律的工具.地质统计 这些工具包括 histogram, covariance, variogram, kriging, 等等.6. 断层的自动识别 7. 用生物芯片发现蛋白质生物标志物 质谱仪两组质谱数据以及能把它们分开的几个蛋白质8. 差式共振声谱仪 Differential Resonance Acoustic Spectroscopy (DARS)圆筒共振腔E

11、lastic constants, Q-value and permeability may be estimated from these spectra DARS 共振谱9. 地震监测注入地下的CO2把二氧化碳注入地下CO2SealFaultInjection Well为什么要进行监测探测CO2移动和泄漏的路径以保证公众的安全了解被注入的CO2时间和空间分布以便更有效地进行注入操作10. 环形换能器超声成像 环形换能器和仿真模型仿真模型的重建Given ModelPicked Travel TimeReconstructed Velocity ModelReconstructed Refl

12、ectivity Model11. 油田油储的综合监测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