用户名: 密码: 验证码:
中、下地壳深地震反射地震资料的AVO岩性识别技术
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
随着深反射地震技术从解决构造问题转向解决深部的岩性问题,深反射地震的AVO(Amplitude versus Offset)技术已经成为深部地球物理的重要发展方向。因此,在我的研究中,基于Hudson和Thomsen裂隙介质理论,将中、下地壳岩体等效为裂隙介质,考虑裂隙介质各向异性的特点,采用R(u|¨)ger(1998)的HTI(包括VTI)介质中的AVO公式,考虑裂缝密度、裂缝方位、孔隙度、流体充填物、饱和状态,甚至地层压力、温度等因素分析中、下地壳介质的AVO响应,并采用实际的WLR(Wollaston Lake Reflector)反射剖面数据进行验证。通过我们的研究,希望能够利用反射振幅的AVO属性,来建立中下地壳岩石物性的识别方法与准则,从而进一步深化对中下地壳岩石物性、形成机制及其动力学过程的认识。
     结果表明,WLR层的真实振幅值的总体趋势是随着炮检距逐渐增大的,AVO斜率很小。将真实振幅值与四个各向异性模型数据拟合之后,发现模型一与真实数据拟合的最好,而模型三偏离真实数据较远。因此,在WLR层反射亮点的四个解释成因中,模型一(辉绿岩,Christensen,1996)是最佳成因。且随着泊松比的减小,拟合误差也逐渐降低,这意味着WLR的强振幅可能是由低泊松比的地层引起的。运用各向异性理论AVO交绘图,我们发现泊松比和纵波速度是检测岩性的最佳因子。在反射系数交绘图中,AVO曲线的斜率和截距交绘图也能很好的区分岩性。从实际地震剖面中分别提取了含流体层段和不含流体层段的AVO斜率和截距,发现斜率和截距交绘图能很好的区分含流体与不含流体层段。而且从含流体到不含流体层段的变化趋势和理论模型的变化趋势是一样的。这说明我们的各向异性AVO理论是正确的。通过对裂隙和孔隙并存介质的有效压力敏感性分析,发现随着有效压力的增大,介质中的裂隙随之闭合,各向异性现象逐渐消失;速度和反射系数也增大,但增量逐渐减小;比较考虑与不考虑等径孔隙的模型表明,在相同的弹性参数下,有效压力从0MPa增加到30MPa的过程中,其速度和反射系数的差值是逐渐增大的,在20MPa—30MPa左右,其差值达到最大。而有效压力继续增大时,其两者的速度差值逐渐减小,最终重合。因为在0—30MPa的范围内裂隙孔隙度和等径孔隙度的变化也最大。随着有效压力的继续增大,介质的弹性参数也趋于常数。
With the deep seismic reflection techniques changing from solving structural problems to lithology problems, Amplitude versus offset (AVO) technique has become an important research direction in deep seismic sounding. In my research, based on Hudson and Thomsen's fracture media theory, I made fracture media rock physics model in middle and lower crust. Considering the anisotropic characterization of fracture media and the factor of cracks density, cracks directions, porosity, filled fluid, pressure and temperature, I used Ruger's AVO formula for HTI (include VTI) media, and applied it to the real WLR (Wollaston Lake Reflector) seismic reflection profile. I hope to establish the identification method and rule for discriminating lithology in middle or lower crust by using the AVO attributes . Furthermore, my study will improve our understanding of the petrophysics, form-mechanism and the dynamics process in middle or lower crust.
     My study shows that, the whole trend of AVO for WLR1 is increasing. The AVO gradient is small. We calculated four anisotropic models to fit WLR1 data, and found model one is the best fit model. Model three is the worst which was considered the best fittest model by using isotropic theory. So I think that in the present the four candidates of WLR reflection bright, the Model one is the best candidate. With the decreasing of Possion's ratio, the misfit error is gradually decreasing in Model one. Therefore, I can infer that the high amplitude of WLR may be caused by the lower Possion's ratio formation. By using the anisotropic AVO theory, I found that Possion's ratio and the P-wave velocity is the best factor to discriminate lithlogy for middle and lower crust metamorphic rock. Reflection coefficient crossplot, AVO gradient and intercept crossplot can also be used to discriminate lithology. AVO gradient and intercept crossplot of WLR reflection had prove this. It means that our anisotropic AVO theory are correctly to be used here as well.
     By analyzing the sensitivity of fracture-pore media, I found that with the increasing of effective pressure, the fracture in media is closed, the anisotropic phenomena is gradually disappeared, the velocity and reflection coefficient also increased. However, the increment of velocity and reflection coefficient is gradually decreasing. Comparsion of no-equant variation modle with the same elastic parameters shows that reflectivity or velocity difference is increasing and reach to 30MPa, and continue decreases and finally tend to constant with the increasing of effecive pressure. The reason is that variation of crack- porosity andpore-porosity is the largest from 0-30Mpa, and gradually tend to constant with the increasingof pressure.
引文
[1] Ake, J. P., and A. R. Sanford. New evidence for the existence and internal structure of a thin layer of magma at mid-crustal depths near Socorro[J], New Mexico, Bull. Seismol.Soc. Am., 2003, 78 (3), 1335-1359
    
    [2] Aki, K., and P. G. Richards. Quantitative Seismology, Second Edition[J], University Science Books, Sausalito, CA,2002
    [3] Bakulin, A., Grechka,V., & ,I. Estimation of fracture parameters from reflection seismic data[J], Part II: Fractured models with orthorhombic symmetry. Geophysics,2000b,65,1803-1817
    [4] Barnes, A. E., and T. J. Reston.A study of two mid-crustal bright spots from Southeast Georgia (USA)[J], Geophys. J. Int., 1998,108(2), 683-691
    [5] Cameli, G. M., I. Dini, and D. Liotta.Upper crustal structure of the Larderello geothermal field as a feature of post-collisional extensional tectonics (southern Tuscany, Italy)[J],Tectonophysics, 1993, 224(4), 413-423
    [6] Castagna, J., Swan, H., and Foster, D., Framework for AVO gradient and intercept interpretation[J], Geophysics, 1998,63, 948-956
    [7] DEKOPR Research Group. Result of the DEKOPR 4/KTB Oberpfalz deep seismic reflection investigations[J]. Geophysics, 62:69
    [8] Domenico S N, Effect of Brine-gas mixture on velocity in an unconsolidated sand reservoir[J]. Geophysics, 1967,41(5):882-894
    [9] Domenico S N.Elastic properties of unconsolidated porous sand reservoirs [J].Geophysics, 1977,42(7): 1339-1368
    
    [10] Hudson, J.A., Effective-medium theories for fluid-saturated materials with aligned cracks[J]: Geophysical Prospecting,2001, 49,509-522
    
    [11] Hudson,J.A.,Wave speeds and attenuation of elastic waves in material containing cracks[J], Geophys.J.R.astr.Soc,1981,64,133-150
    
    [12] Hudson, J.A.,The effect of fluid pressure on wave speeds in cracked solid[J],Geophysical Journal International,2000,143,302-310
    [13] Hilterman, F.J., Seismic amplitude interpretation[M],SEG/EAGE Dist. Instr.Series No.4,Soc. Expl. Geoph.2001
    
    [14] Hobbs, R. W., and B. J. Drummond. Can we accurately estimate meaningful reflection coefficients from 2-D seismic data?[J] 11th International Symposium on Deep Structure of the Continents and their Margins: Results from Reflection, Refraction and Teleseismic Seismology, 26 Sept - 1 Oct 2004, Mont-Tremblant, Quebec, Canada, A29
    [15] Hussenoeder, S. A., J. A. Collins, G. M. Kent, R. S. Detrick, and TERA-Group, Seismic analysis of the axial magma chamber reflector along the southern East PacificRise from conventional reflection profiling[J],J. Geophys. Res., 101, 22,087-22,105
    [16] Hyndman, R. D., and P. M. Shearer, Water in the lower continental crust: Modeling magnetotelluric and seismic reflection results[J], Geophys.J. R. Astron. Soc, 1989,98(2),343-365
    [17] Juhlin, C, Interpretation of the reflections in the Siljan Ring area based on results from the Gravberg-1 borehole[J].Tectonophysics,1990,173:345-360
    [18] Jin F M, Igor B. M. Structure of the Wollaston Lake Reflector (Trans-Hudson Orogen,Canada) from reflection AVO analysis: Fractured diabase intrusion, fluids, or silicified shear zone?[J] Tectonophysics, 2007, 441:97-114
    [19] Jarchow, M. J., G. A. Thompson, R. D. Catchings, and W. D. Mooney, Seismic evidence for active magmatic underplating beneath the Basin and Range province, western United States[J],J. Geophys. Res., 1993,98, 22,095-22,108
    [20] Knott,C.G., Reflection and refraction of elastic waves with seismological applications[J].Phil.Mag.,1899,48,64-97
    [21] Koefoed, O.,On the effect of Poisson's ratios of rock strata on the reflection coefficients of plane waves[J]. Geophys.Prosp.,1955,3,381-387
    [22] Makovsky, Y., Klemperer, S.L., Measureing the seismic properties of Tibetan bright spots:Evidence for free aqueous fluids in the Tibetan middle crust [J]. Journal of Geophysical Research, 1955 ,104: 10795 -10825
    [23] Mandler, H. A. F., and R. M. Clowes,The Wollaston Lake reflector: a 160 km long "bright spot" in the Precambrian crust of the Trans-Hudson orogen, Canada, in Deep Seismic Reflection Profiling of the Continents[J], II: A Global Survey, edited by S. L.Klemperer and W. D. Mooney, Tectonophysics, 1998,288, 71-81
    [24] Mandler, H. A. F. and Clowes, R. M. Evidence for extensive tabular intrusions in the Precambrian shield of western Canada: A 160-km-long sequence of bright reflections[J].Geology, 1997, 25: 271-274
    [25] Mandler, H. A., F., and Clowes, R. M., The HSI bright reflector: further evidence for extensive magmatism in the Precambrian in western Canada[J], Tectonophysics,1998,288: 71-81
    [26] Matsumoto, S., and A. Hasegawa, Distinct S-wave reflector in the midcrust beneath Nikko-Shirane Volcano in the northeastern Japan Arc[J], J. Geophys. Res., 1996,101, 3067-3083,
    [27]Simon.M,AVO analysis by offset-limited prestack migrations of crustal seismic data[J],Tectonophy -sics,1998,286:143-153
    [28]Ostrander W J,Plane waves reflection coefficients for gas sands at normal angles of incidence[J].Geophysics,1984,49(10)1 1637-1648.
    [29]O'Connell,R.G.,and Budiansky,B.,Viscoelastic properties of fluid-saturated cracked solids[J]:Journal of Geophysical Research,1997,82,5719-5735.
    [30]Philip T.C.Hammer,Ron M.Clowes,and Kumar Ramachandran.High-resolution seismic reflection imaging of a thin,diamondiferous kimberlite dyke[J].Geophysics,Vol.69(5),1143-1154,2004
    [31]Pratt,T.L.,Mondary,J.E,Brown,L.D.,Christensen,N.I.,and Danbom,S.H.,Crustal structure and deep reflector pproperties:Wide angle shear and compressional wave studies of the midcrust Surrency bright spot beneath southeastern Georgia[J].Journal of Geophysical Research,98:17723-17735,1993
    [32]Ross,G.M.,and Eaton,D.W.,The Winagami reflector sequence:Seismic evidence for post-collisional magmatism in the Proterozoic of western Canada[J]:Geology,v.25,p.199-202.1997
    [33]R(u|¨)ger A.,Variation of P-wave reflectivity with offset and azimuth in anisotropic media[J]:Geophysics,63,935-947.1997
    [34]Ryberg,T.,and G.S.Fuis,The San Gabriel bright reflective zone:Possible evidence of young mid-crustal thrust faulting in southern California,in Deep Seismic Reflection Profiling of the Continents,I:General Results and New Methods[J],edited by S.L.Klemperer and W.D.Mooney,Tectonophysics,286,31-46,1998
    [35]Rutherford S R,Williams R H,Amplitude versus offset variations in gas sands [J].Geophysics,1989,54(6)1680-688v[36]Shuey R T,1985,Amplification of the zoeppritz's equations[J].Geophysics,50(4):609-614
    [37]Smith,G.C.,and Gidlow,P.M.,Weighted stacking for rock property estimation and detection of gas[J]:Geophys.Prosp.,1987,35,993-1014
    [38]Thomsen L.,Elastic anisotropy due to aligned cracks in porous rock:Geophysical Prospecting,1995,43,805-829
    [39]Thomsen L.,Weak elastic anisotropy[J]:Geophysics,1986a,51,1954-1966
    [40]Tod,S.R.,The effects of stress and fluid pressure on the anisotropy of interconnected cracks[J]:Geophysical Journal International,2002,149,149-156
    [41]Unsworth,M.J.,Jones,A.G.,Wei,W.,Marquis,G.,Gokam,S.G.,Spratt,J.E.& the INDEPTH-MT team.Crustal rheology of the Himalaya and Southern Tibet inferred from magnetotelluric data[J].Nature,2005,438,78-81(2005)
    [42]Vladimir G,Multiple cracks in VTI rocks:Effective properties and fracture characterization[J]:Geophysics,2002,72(5),81-91
    [43]Vernik,L.,S.Hickman,D.Lockner,and M.Rusanov.Ultrasonic velocities in cores from the Kola Superdeep well and the nature of subhorizontal seismic reflections[J].J.Geophys.Res.,2004,99,24209 -24219
    [44]Warner,M.,Basalts,water or shear zones in the lower continental crust?[J]Tectonophysics,2002,173:163-174
    [45]Xin-Gong Li Han De-Hua,/nversion of Sw and porosity from seismic AVO[J]SEG/Houston 75~(th) SEG Annual Meeting,2005
    [46]Yanghua Wang.Approximations to the Zeoppritz equations and their use in AVO analysis[J].Geophysics,1999,Vol.64(6):1920-1927
    [47]吕庆田,侯增谦,杨竹森,等.长江中下游地区的底侵作用及动力学演化模式:来自地球物理资料的约束[J],中国科学,D辑,2004,34(9),783-794
    [48]杨宝俊,李勤学,唐建人,等.松辽盆地反射地震莫霍面的形态、三瞬处理结果及其地质解释[J],地球物理学报,2003,46(3),398-402
    [49]杨宝俊,刘财,杨平华,等.用近垂直地震反射方法研究莫霍面的特征与成因[J],地球物理学报,1999,42(5),617-619
    [50]赖晓玲,李松林,张先康,等.玛沁-兰州-靖边剖面壳幔复杂性的研究[J],地球物理学进展,2001,16(2),66-72
    [51]赵文津,K.D.纳尔逊,车敬凯,等.喜马拉雅地区深反射地震-揭示印度大陆北缘岩石圈的复杂结构[J].地球学报,1996,17(2),138-152
    [52]王海燕,高锐,卢占武,等.地球深部探测的先锋-深地震反射方法的发展与应用[J],勘探地球物理进展,2006,29(1),7-14
    [53]王海燕,高锐,马永生,等.若尔盖盆地-西秦岭造山带结合部位深反射资料的静校正和去噪技术[J].地球物理学进展.2007,22(3),743-749
    [54]刘葵,赵文津.地震相分析在深反射地震勘探资料解释中的应用[J],地质力学学报。2005,11(2),129-134
    [55]刘永祥,伍旭光.高温高压下岩石波速与地壳深部的物质组成[J].地球物理学报,1997, 40(2):211-219
    [56]刘斌,王宝善,席道瑛.水饱和裂纹对地壳岩样中地震波速及各向异性的影响[J].地球物理学报.1999,42(5):702-710
    [56]郑晓东,Zoeppritz方程的近似及其应用[J].石油地球物理勘探,1991,26(2):126-144
    [57]杨绍国,周熙襄.Zoeppritz方程的级数表达式及近似[J].石油地球物理勘探,1994,29(4):399-412
    [58]董树文,高锐,李秋生,等,.大别山造山带前陆深地震反射剖面[J].地质学报.2005,79(5),595-603.
    [59]王海燕,高锐,薛爱民,等,.适用于造山带深地震反射资料的动校正方法[J].吉林大学学报,.2006,36(4),622-627.
    [60]王海燕,高锐,马永生,等,.青藏高原东北缘盆山结合部位深反射资料处理方法与初步地质认识[J].中国地质,2007,34(1),142-149
    [61]王紫娟.AVO交绘图分析[J],天然气工业,2003,23(增刊),49-51
    [62]孙鹏远,孙建国,卢秀丽,.PP波AVO多属性交绘图分析[J],天然气工业,2004,24(2),44-47
    [62]孙鹏远,孙建国,卢秀丽,P-SV波反射系数近似及其AVO属性特征[J],地球学报,2006,27(2),85-89

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700