构造断层的相干识别技术
|
摘要
在地震勘探资料解释中,探测和识别断层和断裂带是最重要的问题之一。识别断层的技术中,地震相干技术是最具代表性的。它是旨在计算由于构造、岩性、地层、空隙的变化及油气层等因素引起的地震响应横向变化的一种方法。
在C1相干算法中,主要的影响因素有:搜索范围、时窗大小、不同道数的三角形组合模式。通过分析确定了搜索范围、时窗的大小和不同道数的三角形组合模式对低幅度构造的分辨率和追踪效果。
C1相干算法通过波形的相似来对模型进行追踪,其不连续的位置为异常构造,但不能断定是否为断层构造。而小断层在地震剖面上表现出连续的地震反射发生错断或扭动,其波形的相似性发生变化,连续性变差。相干技术可以根据断层的这一地震反射特点,定量地算出在一定时窗内相邻地震道间的相干系数,就可以在相关值分布图上把细小断层反映出来。本文应用相干算法,计算其搜索范围为零相似系数(C(0))值,通过观察C(0)变化来识别断层构造。
建立了二维水平断层、倾斜断层、三维水平断层、三维阶梯构造、三维盆岭构造、含有噪音阶梯构造等模型,应用相干算法对这些模型进行同相轴追踪;针对出现的异常构造求搜索范围τ=0时相似系数(C(0))值,观察其变化,来识别断层构造。并对结果进行分析,结果表明,应用上述方法可以识别效果较好。
Fault identification and interpretation of research has a vital role in oil and gas fields in the process of exploration and development.Fault is a very important for oil and gas reservoirs,one of the geological features.Detection and identification of faults and fault zone is seismic data interpretation in one of the most important issues.in recent years,The interpretation of seismic data have become increasingly demanding when the development of oil and gas exploration work.Conventional seismic interpretation methods and traditional methods of seismic interpretation by a certain amount of constraint.Must seek a more accurate method of seismic data interpretation.With the seismic attributes the rapid development of technology and computer technology made the use of the wide,For example,Coherent.seismic technology is one of them.Coherent seismic technology is calculated changes in response to a method of horizontal for as the structure,lithology,stratigraphy,and changes in oil and gas gap layer and so on.Earthquake-induced changes in response to a method of horizontal is axis with the phase detection.In particular when the in-phase axis and the amplitude value of the smaller hidden in noise
C1 coherent selection algorithm in my article.The preparation of the procedures first of Algorithm for the coherent analysis of the parameters for algorithm.Analysis of search area:When the search area is too small, Coherence algorithm can reduce the computing time.But not suitable for large structures with seismic data.When the search area is too small,Structure is conducive to the identification of But the increase in computing time,the impact of efficiency.Therefore,the scope of the search to obtain an appropriate value. The value of the general search area T~*/2～3 T~*/2(T~* Wavelet as the cycle).Because this article is mainly for low-amplitude structures.Therefore,the value of search to T~*/4～T~*/2.Analysis of time window size:Access time window is too small can increase the resolution,is conducive to the identification of small faults,a clearer portrait of the details,but the impact on computational speed.But when there is noise,the small faults of other noise to cover up information,it is difficult to identify,only suitable for local area high-quality seismic data,if too much time window,the faster the more obvious big fault characterization,but easy to lose sight of the details of information,is not conducive to the identification of small faults.Range of time windows: T~*/2～3 T~*/2(T~* Wavelet as the cycle).A few different combinations of the triangle pattern analysis,The impact of noise in the absence of conditions.The use of C1 Coherence Algorithm.Combination of model space triangle 3 load triangle 5 load triangle 7 load of ladder model to track.The results showed that one of the three highest resolution.Clear not only to track a fault.To determine the location of faults than the 5 and 7 accurately
Coherence parameters of analysis and discussion,Determine the parameters.The geological model was used of cl algorithm on different for tracking and identification.C1 coherent waveform similarity algorithm to track the model.Its not unusual for the location for the structure.But we can not determine whether the fault structure.Small faults show continuous seismic reflection or twist off the wrong took place in Seismic profiles.The similarity of the waveform changes.Coherence tomography technique can be characterized by the seismic reflection.Quantitatively calculated within a certain time window between adjacent coherent seismic coefficient.Value can be related to the distribution of small faults is reflected in.In my paper,coherence algorithm.Calculation of the scope of its search values zero coherence.By observing the C(0) changes to identify the fault structure.
In this paper,Established a two-level faults,tilted fault,three-dimensional level of fault,three-dimensional ladder structure,three-dimensional basin and range structure,ladder structure to contain noise.Coherent application of these model algorithm phase axis tracking.Abnormal structure appeared for the scope of the search forτ=0 when coherence coefficient(C(0)) values.Observe the changes to identify the fault structure.The level of fault for the two-dimensional, three-dimensional level of fault,three-dimensional ladder structure,three-dimensional model of basin and range structure.Phase axis tracking.Calculation of the scope of its search whenτ=0 coherence coefficient(C(0)) values.Two-dimensional model of the level of fault.Coherence coefficient C(0) There are two numerical value and a point mutation Coherence coefficient C(0). Region in which the level of value formation,point mutations for the breakpoint. Two-dimensional tilted fault model,delay coherence coefficient(C(0)) is a uniform change,the two breakpoints in the Department,C(0) mutations,we indicate the location for the fault location.
Discussed the different signal to noise ratio of the number of Triangle Road C1 coherent mix of the impact of algorithm.SNR when less than 1:0.2.Application of the three best.SNR when less than 1:0.4.3 and 5 poor results,7 the effect of general.Contain the noise level in the structural model. SNR when less than 1:0.2.Application of three triangles,to strike a coherent factor.Stratigraphic level of the coherence coefficient(C(0)) for the change in value,so coherent approximation factor.Coherence coefficient C(0) There are two numerical approximation,and a mutation at the level of approximation of regional stratigraphy,point mutation for the fault structure,so can it be said that the level of the structure for the fault structure,when less than 1:0.4.5 applications to strike a coherent triangle coefficient.Stratigraphic level of the coherence coefficient(C(0)) for the changes in the value of.Therefore,the approximation of coherent factor.Coherence coefficient C(0) There are three numerical approximation and two-point mutation.Approximation to the level of the regional stratigraphy,point mutation for the fault structure.So can it be said that the ladder structure for the fault structure.
在C1相干算法中,主要的影响因素有:搜索范围、时窗大小、不同道数的三角形组合模式。通过分析确定了搜索范围、时窗的大小和不同道数的三角形组合模式对低幅度构造的分辨率和追踪效果。
C1相干算法通过波形的相似来对模型进行追踪,其不连续的位置为异常构造,但不能断定是否为断层构造。而小断层在地震剖面上表现出连续的地震反射发生错断或扭动,其波形的相似性发生变化,连续性变差。相干技术可以根据断层的这一地震反射特点,定量地算出在一定时窗内相邻地震道间的相干系数,就可以在相关值分布图上把细小断层反映出来。本文应用相干算法,计算其搜索范围为零相似系数(C(0))值,通过观察C(0)变化来识别断层构造。
建立了二维水平断层、倾斜断层、三维水平断层、三维阶梯构造、三维盆岭构造、含有噪音阶梯构造等模型,应用相干算法对这些模型进行同相轴追踪;针对出现的异常构造求搜索范围τ=0时相似系数(C(0))值,观察其变化,来识别断层构造。并对结果进行分析,结果表明,应用上述方法可以识别效果较好。
Fault identification and interpretation of research has a vital role in oil and gas fields in the process of exploration and development.Fault is a very important for oil and gas reservoirs,one of the geological features.Detection and identification of faults and fault zone is seismic data interpretation in one of the most important issues.in recent years,The interpretation of seismic data have become increasingly demanding when the development of oil and gas exploration work.Conventional seismic interpretation methods and traditional methods of seismic interpretation by a certain amount of constraint.Must seek a more accurate method of seismic data interpretation.With the seismic attributes the rapid development of technology and computer technology made the use of the wide,For example,Coherent.seismic technology is one of them.Coherent seismic technology is calculated changes in response to a method of horizontal for as the structure,lithology,stratigraphy,and changes in oil and gas gap layer and so on.Earthquake-induced changes in response to a method of horizontal is axis with the phase detection.In particular when the in-phase axis and the amplitude value of the smaller hidden in noise
C1 coherent selection algorithm in my article.The preparation of the procedures first of Algorithm for the coherent analysis of the parameters for algorithm.Analysis of search area:When the search area is too small, Coherence algorithm can reduce the computing time.But not suitable for large structures with seismic data.When the search area is too small,Structure is conducive to the identification of But the increase in computing time,the impact of efficiency.Therefore,the scope of the search to obtain an appropriate value. The value of the general search area T~*/2～3 T~*/2(T~* Wavelet as the cycle).Because this article is mainly for low-amplitude structures.Therefore,the value of search to T~*/4～T~*/2.Analysis of time window size:Access time window is too small can increase the resolution,is conducive to the identification of small faults,a clearer portrait of the details,but the impact on computational speed.But when there is noise,the small faults of other noise to cover up information,it is difficult to identify,only suitable for local area high-quality seismic data,if too much time window,the faster the more obvious big fault characterization,but easy to lose sight of the details of information,is not conducive to the identification of small faults.Range of time windows: T~*/2～3 T~*/2(T~* Wavelet as the cycle).A few different combinations of the triangle pattern analysis,The impact of noise in the absence of conditions.The use of C1 Coherence Algorithm.Combination of model space triangle 3 load triangle 5 load triangle 7 load of ladder model to track.The results showed that one of the three highest resolution.Clear not only to track a fault.To determine the location of faults than the 5 and 7 accurately
Coherence parameters of analysis and discussion,Determine the parameters.The geological model was used of cl algorithm on different for tracking and identification.C1 coherent waveform similarity algorithm to track the model.Its not unusual for the location for the structure.But we can not determine whether the fault structure.Small faults show continuous seismic reflection or twist off the wrong took place in Seismic profiles.The similarity of the waveform changes.Coherence tomography technique can be characterized by the seismic reflection.Quantitatively calculated within a certain time window between adjacent coherent seismic coefficient.Value can be related to the distribution of small faults is reflected in.In my paper,coherence algorithm.Calculation of the scope of its search values zero coherence.By observing the C(0) changes to identify the fault structure.
In this paper,Established a two-level faults,tilted fault,three-dimensional level of fault,three-dimensional ladder structure,three-dimensional basin and range structure,ladder structure to contain noise.Coherent application of these model algorithm phase axis tracking.Abnormal structure appeared for the scope of the search forτ=0 when coherence coefficient(C(0)) values.Observe the changes to identify the fault structure.The level of fault for the two-dimensional, three-dimensional level of fault,three-dimensional ladder structure,three-dimensional model of basin and range structure.Phase axis tracking.Calculation of the scope of its search whenτ=0 coherence coefficient(C(0)) values.Two-dimensional model of the level of fault.Coherence coefficient C(0) There are two numerical value and a point mutation Coherence coefficient C(0). Region in which the level of value formation,point mutations for the breakpoint. Two-dimensional tilted fault model,delay coherence coefficient(C(0)) is a uniform change,the two breakpoints in the Department,C(0) mutations,we indicate the location for the fault location.
Discussed the different signal to noise ratio of the number of Triangle Road C1 coherent mix of the impact of algorithm.SNR when less than 1:0.2.Application of the three best.SNR when less than 1:0.4.3 and 5 poor results,7 the effect of general.Contain the noise level in the structural model. SNR when less than 1:0.2.Application of three triangles,to strike a coherent factor.Stratigraphic level of the coherence coefficient(C(0)) for the change in value,so coherent approximation factor.Coherence coefficient C(0) There are two numerical approximation,and a mutation at the level of approximation of regional stratigraphy,point mutation for the fault structure,so can it be said that the level of the structure for the fault structure,when less than 1:0.4.5 applications to strike a coherent triangle coefficient.Stratigraphic level of the coherence coefficient(C(0)) for the changes in the value of.Therefore,the approximation of coherent factor.Coherence coefficient C(0) There are three numerical approximation and two-point mutation.Approximation to the level of the regional stratigraphy,point mutation for the fault structure.So can it be said that the ladder structure for the fault structure.
引文
[1]王志军,黄军斌.利用相干技术和三维可视化技术识别微小断层和砂体[J].石油地球物理勘探,2001,36(3):378-381.
[2]吕小惠.三维地震资料中断层识别方法研究[D].南京:南京理工大学,2004.
[3]李雄志.基于三维可视化的断层解释技术[D].西安:西安科技大学,2004.
[4]郭栋,韩文功.高分辨率地震资料综合解释技术及其应用[J].勘探地球物理进展,2004,27(4):290-296.
[5]张进铎.地震解释技术现状及发展趋势[J].地球物理学进展,2006,21(2):578-587.
[6]Brown A R,Edwards G S,Howard R E.Fault slicing:a new approach to the interpretation of fault detail[J].Geophysics,1987,10:1319-1327.
[7]张延章,韩品龙,池永红,等.地震相干技术的应用及效果分析[J].中国海上油气(地质),2003,17(3):215-217.
[8]张延章,李淑恩,黄国平,等.地震切片的分类及应用价值[J].油气地质与采收率,2002,9(3):67-69.
[9]Bahorich M S,Bridges S R.The seismic sequence at tribute map[J].Expanded Abstracts of the 62nd Annual Internat SEG Meeting,1992:227-230.
[10]Bahorich M S,Farmer S L.3-D seismic coherency for faults andstratigrAphi features[J].The Leading Edge,1995,14(10):1053-1058.
[11]Bahorich M S,Lopez J A.Stratigraphic and structural interpre-tation with 3-D coherence[J].Expanded Abstracts of the 65th Annual Internat SEG Meeting,1995,97-100.
[12]崔鹏飞,赵爱华.利用模式识别方法解释微小断层[J].石油地球物理勘探,1995,20(4):556-566.
[13]崔若飞.煤田地震资料精细构造解释技术[J].物探化探计算技术,1998,13(4):312-315.
[14]王大伟,刘震,陈小宏,等.地震相干技术的进展及其在油气勘探中的应用[J].地质科技情报,2005,24(2):71-76.
[15]张金淼,王春红.相干技术在全三维地震解释中的应用[J].中国海上油气(地质),2000,14(4):277-282.
[16]E.W斯宾塞.地球构造导论[M].北京:地质出版社,1981.
[17]庄培仁,常志忠.断裂构造研究[M].北京:地质出版社,1993.
[18]王桂梁.地质构造图册[M].北京:地质出版社,1992.
[19]JI普拉特,J查利诺.简单的地质构造[M].北京:地质出版社,1959.
[20]G威尔逊,J科斯格罗夫.小型地质构造[M].北京:地质出版社,1989.
[21]旺新文,林建平.地球科学概论[M].北京:地质出版社,2005.
[22]陆基孟.地震勘探原理(下)[M].东营:中国石油大学出版社,2006.
[23]何樵登,熊维纲.应用地球物理教程-地震勘探[M].长春:吉林大学出版社,2003.
[24]R.E.谢里夫.勘探地震学(上)[M].北京:石油工业出版社,1999.
[25]陆基孟.地震勘探原理(上)[M].东营:中国石油大学出版社,2006.
[26]董敏煜.地震勘探[M].北京:石油大学出版社,2000.
[27]李大心.地球物理方法综合应用与解释[M].武汉:中国地质大学出版社,2003.
[28]宗孔德,胡广书.数字信号处理[M].北京:清华大学出版社,1999.
[29]许庆山.信号与系统[M].北京:航空工业出版社,1998.
[30]R.E.谢里夫.勘探地震学(下)[M].北京:石油工业出版社,1999.
[31]丁玉美,高西全.数字信号处理[M].西安:西安电子科技大学出版社,2006.
[32]吕小惠.三维地震资料中断层识别方法研究[D].南京:南京理工大学,2004
[33]付雷,田晓冬,闫庆学,等.相干技术在三维地震资料中的应用[J].世界地质,2000,19(3):279-281.
[34]佘德平,曹辉,王咸斌.相干数据体及其在三维地震解释中的应用[J].石油物探,1998,37(4):75-79.
[35]Marfurt K J,Kirlin R L.3-D seismic attrihotes using a running window setnblance -basedalgorithm[J].Geophysics,1998,63(4):1150-1165.
[36]Lu W,Li Y,Xiao H.Higher-order statisticsbased coherency estimation algorithm [J].Expanded Abstracts of the 73th Annual Internat SEG Meeting,2003:1732-1735.
[37]Gersztenkorn A,Marfurt K J.Eigenstructure based coherence computationns [J].Expanded Abstracts of the 66th Annual Internat SEG Meeting,1996:328-331.
[38]Gersztenkorn A,Marfurt K J.Eigenstructure based coherence computations as an aid to 3-D structural and stratigraphic mapping[J]Geophysics,1999,64(5):1468-1479.
[39]Cohen I,Coifman R R.Local discontinuity measures for 3-D seismicdata[J].Geophysics,2002,67(6):1933-1945.
[40]Marfurt K J,Sndbakar V.Gerszrenkorn A,et al.Coherency calculations in the presence of structural dip[J].Geophysics.1999,64(1):104-111.
[41]汤磊.地质约束相干方法研究及应用[D].北京:中国石油勘探开发研院,2002.
[42]汤祥武,宋中应.三维地震勘探小断层的解释方法及应用[J].煤矿开采,2002,7(3):12-13.
[43]李玲,冯许魁.用地震相干数据体进行断层自动解释[J].石油地球物理勘探.1998,33(1):105-111.
[44]苑书金.地震相干体技术的研究综述[J].勘探地球物理进展,2007,30(1):7-15.
[45]闫德庆,杨飞鹏,郭淑敏,等.应用改进的相干算法提高三维地震资料解释精度[J].物探化探计算技术,2001,23(4):314-317.
[46]张军华,王永刚.相干技术算法改进及其在TJH地区的应用[J].物探与化探,2002,26(1):50-52.
[47]闫德庆,庞留彦,仲其涛,等.C1相干算法及其三维地震勘探解释中的[J].物探化探计算技术,2001,23(3):206-209.
[48]陈凤云,杭远,康建林.相干和方差数据体的算法研究及应用[J].物探与化探,2006,30(3):250-253.
[49]叶增炉,何建军.地震相干体算法的改进及应用[J].物探化探计算方法,2006,28(2):121-123
[50]邓小英.Ricker子波最小二乘支持向量回归机及其在地震勘探资料去噪中的应用研究[D].长春:吉林大学博士论文,2007.
[51]Sheriff R.E.Encyclopedic dictionary of exploration geophy-sics[M].Tulsa:Society of Exploration Geophysicists,1994.
[52]李庆忠.走向精确勘探的道路—多分辨地震勘探系统工程剖析[M].北京:石油工业出版社,1994.
[53]求是科技.Matlab7.0从入门到精通[M].北京:人民邮电出版社,2006.
[54]邓薇.Matlab函数速查手册[M].北京:人民邮电出版社,2006.
[2]吕小惠.三维地震资料中断层识别方法研究[D].南京:南京理工大学,2004.
[3]李雄志.基于三维可视化的断层解释技术[D].西安:西安科技大学,2004.
[4]郭栋,韩文功.高分辨率地震资料综合解释技术及其应用[J].勘探地球物理进展,2004,27(4):290-296.
[5]张进铎.地震解释技术现状及发展趋势[J].地球物理学进展,2006,21(2):578-587.
[6]Brown A R,Edwards G S,Howard R E.Fault slicing:a new approach to the interpretation of fault detail[J].Geophysics,1987,10:1319-1327.
[7]张延章,韩品龙,池永红,等.地震相干技术的应用及效果分析[J].中国海上油气(地质),2003,17(3):215-217.
[8]张延章,李淑恩,黄国平,等.地震切片的分类及应用价值[J].油气地质与采收率,2002,9(3):67-69.
[9]Bahorich M S,Bridges S R.The seismic sequence at tribute map[J].Expanded Abstracts of the 62nd Annual Internat SEG Meeting,1992:227-230.
[10]Bahorich M S,Farmer S L.3-D seismic coherency for faults andstratigrAphi features[J].The Leading Edge,1995,14(10):1053-1058.
[11]Bahorich M S,Lopez J A.Stratigraphic and structural interpre-tation with 3-D coherence[J].Expanded Abstracts of the 65th Annual Internat SEG Meeting,1995,97-100.
[12]崔鹏飞,赵爱华.利用模式识别方法解释微小断层[J].石油地球物理勘探,1995,20(4):556-566.
[13]崔若飞.煤田地震资料精细构造解释技术[J].物探化探计算技术,1998,13(4):312-315.
[14]王大伟,刘震,陈小宏,等.地震相干技术的进展及其在油气勘探中的应用[J].地质科技情报,2005,24(2):71-76.
[15]张金淼,王春红.相干技术在全三维地震解释中的应用[J].中国海上油气(地质),2000,14(4):277-282.
[16]E.W斯宾塞.地球构造导论[M].北京:地质出版社,1981.
[17]庄培仁,常志忠.断裂构造研究[M].北京:地质出版社,1993.
[18]王桂梁.地质构造图册[M].北京:地质出版社,1992.
[19]JI普拉特,J查利诺.简单的地质构造[M].北京:地质出版社,1959.
[20]G威尔逊,J科斯格罗夫.小型地质构造[M].北京:地质出版社,1989.
[21]旺新文,林建平.地球科学概论[M].北京:地质出版社,2005.
[22]陆基孟.地震勘探原理(下)[M].东营:中国石油大学出版社,2006.
[23]何樵登,熊维纲.应用地球物理教程-地震勘探[M].长春:吉林大学出版社,2003.
[24]R.E.谢里夫.勘探地震学(上)[M].北京:石油工业出版社,1999.
[25]陆基孟.地震勘探原理(上)[M].东营:中国石油大学出版社,2006.
[26]董敏煜.地震勘探[M].北京:石油大学出版社,2000.
[27]李大心.地球物理方法综合应用与解释[M].武汉:中国地质大学出版社,2003.
[28]宗孔德,胡广书.数字信号处理[M].北京:清华大学出版社,1999.
[29]许庆山.信号与系统[M].北京:航空工业出版社,1998.
[30]R.E.谢里夫.勘探地震学(下)[M].北京:石油工业出版社,1999.
[31]丁玉美,高西全.数字信号处理[M].西安:西安电子科技大学出版社,2006.
[32]吕小惠.三维地震资料中断层识别方法研究[D].南京:南京理工大学,2004
[33]付雷,田晓冬,闫庆学,等.相干技术在三维地震资料中的应用[J].世界地质,2000,19(3):279-281.
[34]佘德平,曹辉,王咸斌.相干数据体及其在三维地震解释中的应用[J].石油物探,1998,37(4):75-79.
[35]Marfurt K J,Kirlin R L.3-D seismic attrihotes using a running window setnblance -basedalgorithm[J].Geophysics,1998,63(4):1150-1165.
[36]Lu W,Li Y,Xiao H.Higher-order statisticsbased coherency estimation algorithm [J].Expanded Abstracts of the 73th Annual Internat SEG Meeting,2003:1732-1735.
[37]Gersztenkorn A,Marfurt K J.Eigenstructure based coherence computationns [J].Expanded Abstracts of the 66th Annual Internat SEG Meeting,1996:328-331.
[38]Gersztenkorn A,Marfurt K J.Eigenstructure based coherence computations as an aid to 3-D structural and stratigraphic mapping[J]Geophysics,1999,64(5):1468-1479.
[39]Cohen I,Coifman R R.Local discontinuity measures for 3-D seismicdata[J].Geophysics,2002,67(6):1933-1945.
[40]Marfurt K J,Sndbakar V.Gerszrenkorn A,et al.Coherency calculations in the presence of structural dip[J].Geophysics.1999,64(1):104-111.
[41]汤磊.地质约束相干方法研究及应用[D].北京:中国石油勘探开发研院,2002.
[42]汤祥武,宋中应.三维地震勘探小断层的解释方法及应用[J].煤矿开采,2002,7(3):12-13.
[43]李玲,冯许魁.用地震相干数据体进行断层自动解释[J].石油地球物理勘探.1998,33(1):105-111.
[44]苑书金.地震相干体技术的研究综述[J].勘探地球物理进展,2007,30(1):7-15.
[45]闫德庆,杨飞鹏,郭淑敏,等.应用改进的相干算法提高三维地震资料解释精度[J].物探化探计算技术,2001,23(4):314-317.
[46]张军华,王永刚.相干技术算法改进及其在TJH地区的应用[J].物探与化探,2002,26(1):50-52.
[47]闫德庆,庞留彦,仲其涛,等.C1相干算法及其三维地震勘探解释中的[J].物探化探计算技术,2001,23(3):206-209.
[48]陈凤云,杭远,康建林.相干和方差数据体的算法研究及应用[J].物探与化探,2006,30(3):250-253.
[49]叶增炉,何建军.地震相干体算法的改进及应用[J].物探化探计算方法,2006,28(2):121-123
[50]邓小英.Ricker子波最小二乘支持向量回归机及其在地震勘探资料去噪中的应用研究[D].长春:吉林大学博士论文,2007.
[51]Sheriff R.E.Encyclopedic dictionary of exploration geophy-sics[M].Tulsa:Society of Exploration Geophysicists,1994.
[52]李庆忠.走向精确勘探的道路—多分辨地震勘探系统工程剖析[M].北京:石油工业出版社,1994.
[53]求是科技.Matlab7.0从入门到精通[M].北京:人民邮电出版社,2006.
[54]邓薇.Matlab函数速查手册[M].北京:人民邮电出版社,2006.