声波测井的数值模拟课件.ppt

1、科研活动的简要回顾科研活动的简要回顾1.地震波衰减的层析成像(Stanford)2.声波测井的数值模拟(Stanford)3.油田残余油强化开采的地震监测(Texaco)4.油层的地质统计描述(Texaco)5.地质统计软件(SMT)6.断层的自动识别(SMT)7.用生物芯片发现蛋白质生物标志物(LumiCyte)8.差式共振声谱仪(Stanford)9.地震监测地下存储的CO2(Stanford)10.环形换能器超声成像(Stanford/LANL)11.油田油储的综合监测(Stanford)Stanford University,1991-1996频移法是我们提出和发展的一种简单实用的技术

2、。这一方法自1997发表后得到了广泛的应用。比如说，最近一年的 Geophysics 杂志有6篇文章引用该论文。今年的SEG Expended Abstracts有7篇文章引用了这篇论文。Acoustic SourcesAcoustic ReceiversLower AbsorptionHigher AbsorptionHigher frequencyWaveform Lower frequencyWaveform 波形和衰减波形和衰减用中心频率移动法估计用中心频率移动法估计衰减衰减Medium ResponseH(f)Incident WaveS(f)Transmitted WaveR(f)

3、=S(f)H(f)H(f)=exp-faodl-=raysrsoffdl2axxxxxxxxxxxxooooooooooooooSourcesReceiversVf=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 tomo

4、graphy for 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.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

6、 modelzrzrComplex radially symmetric modelFormationCementFluidSteel Modeling of Waves in a Borehole 测井模型测井模型(横切面横切面)使用的方法使用的方法 Modeling of Waves in a Borehole l 广义反射透射系数法广义反射透射系数法(稳定稳定,速度快速度快)l 归一化的数学处理归一化的数学处理(可以模拟高频和可以模拟高频和厚层的问题厚层的问题)layers Modeling of Waves in a Borehole 计算实例010203040500 0.1 0.2

7、0.3 0.4QpRadius(m)Qs10002000300040005000600000.10.20.30.4VpRadius(m)VsVelocities(m/sec)11.522.500.1 0.2 0.3 0.4DensityRadius(m)A simple fluid-filled open boreholeTime(ms)013422.445.44Source-receiver offset(m)Seismograms in this simple borehole简单充水裸井10002000300040005000600000.10.20.30.4Velocities(m/s

8、ec)Radius(m)VpVs11.522.50 0.1 0.2 0.3 0.4DensityRadius(m)010203040500 0.1 0.2 0.3 0.4QpRadius(m)QsTime(ms)01234Source-receiver offset(m)2.445.44Seismograms in this damaged borehole A damaged fluid-filled open borehole有损伤的裸井oSourceReceivers100 m150 m10 mFormation IFormation IIBoreholexxxxxxxxxx从远方柱面上

9、的反射从远方柱面上的反射(模型模型)ModelVp=3 km/sVs=1.8 km/sVp=1.9 km/sVs=1.4 km/sSeismograms(fo=800 Hz)SP-SS-PP-P-P-PS-SPTube waveP-P101501050100150Time(ms)source-receiver offset(m)从远方柱面上的反射从远方柱面上的反射(地震波地震波)井间地震例子井间地震例子ooooooooooooooo30 mReceiversSourcesCased boreholeVp=5 km/sVs=2.9 km/sVp=5.8 km/sVs=3.3 km/s100 mF

10、ormation(a)Model.There is a fault in the formation(b)A common receiver gatherTube WaveS-waveP-wavexxxxxxxxxxxxxxx30 m102030400Time(ms)0S-Roffset(m)100一种新测井数据的模拟一种新测井数据的模拟(a)Layered borehole model(b)The pressure field recorded at top,middle and bottom(a)(b)Time-lapse CO2 Monitoring D DVp(P Wave Veloc

11、ity Change Vp2-Vp1)Tomogram3400390044004900540001000500150005001000150020002500Depth(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/V

12、p(%)Top of Lower San Andres1.50Distance-100200500800Subsea Depth-1320-1270-1220-1170-1120Porosity Based on Well LogsDistance-100200500800Subsea Depth-1320-1270-1220-1170-1120WELL10WELL8WELL3WELL6 WELL4Z3Z3Z3Z3Z3Z4Z4Z4Z4Z4GRIDSTAT Sep18 97用测井数据计算孔隙度剖面用测井数据计算孔隙度剖面Distance-100200500800Subsea Depth-1320

13、-1270-1220-1170-1120Impedance Based on Seismic AmplitudeDistance-100200500800Subsea Depth-1320-1270-1220-1170-1120WELL10WELL8WELL3WELL6 WELL4GRIDSTAT Sep18 97用地震数据计算的阻抗剖面用地震数据计算的阻抗剖面Distance(m)-100010003000500070009000Subsea Depth(m)-1500-1400-1300-1200-1100-1000Porosity Section Integrating Seismic

14、and Logs.Distance(m)-100010003000500070009000Subsea Depth(m)-1500-1400-1300-1200-1100-10000.1.15.2.25.3.35.4Well_10 Well_8Well_3Well_6Well_4Nm3Nm4Ng1GRIDSTAT Jul 6 98综合地震和测井数据得到的阻抗剖面综合地震和测井数据得到的阻抗剖面-0.1-0.04-0.02-0.01.01.02.04.086GRIDSTAT 3-16-98油层的三维空隙度模型油层的三维空隙度模型 地质统计研究空间上起伏变化的现象地质统计研究空间上起伏变化的现象.

15、地质统计提供一些用于建立能描述空间变化规律地质统计提供一些用于建立能描述空间变化规律的工具的工具.地质统计地质统计 这些工具包括这些工具包括 histogram,covariance,variogram,kriging,等等.质谱仪两组质谱数据以及能把它们分开的几个蛋白质 Differential Resonance Acoustic Spectroscopy(DARS)圆筒共振腔圆筒共振腔Elastic constants,Q-value and permeability may be estimated from these spectra DARS 共振谱共振谱把二氧化碳注入地下把二氧化

16、碳注入地下CO2SealFaultInjection Well x(mm)y(mm)05010015020005010015020014401460148015001520154015601580 x(mm)y(mm)05010015020005010015020014401460148015001520154015601580Given ModelReceiverTransmitter501001502002505010015020025000.20.40.60.811.2Picked Travel Timex(mm)y(mm)50100150200501001502001440146014

17、8015001520154015601580Reconstructed Velocity ModelReconstructed Reflectivity ModelField DataNoYesDoes synthetic data match field data?Inversion Acoustic properties of rocks Final ModelJerry Harris:Borehole Waves,Attenuation Tomography,DARSXiaofei Chen and Dan Moos:Borehole WavesChuntang Xu:DARSLianjie Huang:Ultrasonic tomography