[go: up one dir, main page]

CN101339252B - Directional lighting seismic exploration method based on single earthquake focus - Google Patents

Directional lighting seismic exploration method based on single earthquake focus Download PDF

Info

Publication number
CN101339252B
CN101339252B CN200810051062XA CN200810051062A CN101339252B CN 101339252 B CN101339252 B CN 101339252B CN 200810051062X A CN200810051062X A CN 200810051062XA CN 200810051062 A CN200810051062 A CN 200810051062A CN 101339252 B CN101339252 B CN 101339252B
Authority
CN
China
Prior art keywords
shot
seismic
distance
point
points
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN200810051062XA
Other languages
Chinese (zh)
Other versions
CN101339252A (en
Inventor
姜弢
林君
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jilin University
Original Assignee
Jilin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jilin University filed Critical Jilin University
Priority to CN200810051062XA priority Critical patent/CN101339252B/en
Publication of CN101339252A publication Critical patent/CN101339252A/en
Application granted granted Critical
Publication of CN101339252B publication Critical patent/CN101339252B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention relates to a seismic prospecting method based on a single seismic focus directional lighting. An offset distance and a track pitch are determined according to the traditional geologic prospecting method; a wave-detector array is arranged in a collinear way; according to the different offset distance of the shot-points in the same group and the same offset distance of the m shot-points in the different groups with the serial number differing in n, seismic data in a single shot on an entire cross-section is gathered; a uniformly-spaced time delay is made for the seismic data in the single-shot in the same group according to the sequence of the serial number; a linear superposition is made for the delayed seismic data in the single-shot so as to compose a directional lighting seismic signal with directivity, and the directional lighting seismic signal is corresponding to a directional seismic wave field; by adopting a routine seismic data processing method, a directional lighting seismic cross-section based on the single seismic focus can be obtained. After a test, the direction of the directional lighting seismic wave field that is composed can be controlled manually according to a target occurrence and has better prospecting effect to the target body of a steep dip, as compared with the common seismic prospecting method based on the single seismic focus; the seismic prospecting method of the invention also has the advantages of low operation cost and simple operation as compared with a phased seismic prospecting method.

Description

基于单震源的定向照明地震勘探方法 Directional illumination seismic exploration method based on single source

技术领域:Technical field:

本发明涉及一种地震勘探方法,尤其适于背景噪声较强,信噪比较低区域的单震源定向照明地震勘探方法。The invention relates to a seismic prospecting method, which is especially suitable for the seismic prospecting method of single seismic source directional illumination in areas with strong background noise and low signal-to-noise ratio.

背景技术:Background technique:

现有的可控震源地震勘探方法按参与工作的震源数目可分为两类。一类是单可控震源地震勘探;另一类是阵列可控震源地震勘探。单可控震源在地下半空间激发的地震波场可近似认为球面波场,无方向性;阵列可控震源激发的地震波场具有明显的方向性,其中组合震源激发的地震波场主波束方向垂直地面向下,相控震源激发的地震波场也存在方向性,且方向可控;一般来说,合理控制地震波方向,采用相控震源可获得更高信噪比的地震资料。但是,相控震源地震勘探方法具有成本高,控制复杂,并且存在多个震源的不一致性问题,上述问题限制了相控震源地震勘探方法的推广。为此,考虑仍采用单震源作为激发源,以基于单震源的单炮地震信号为数据源,采用定向照明地震信号合成方法合成定向照明单道地震信号,类似地再合成定向照明单炮地震数据,该数据对应于定向地震波场。与相控地震波场相似之处在于,两种波场都是定向波场,因此对提高地震资料信噪比都是有效的;与相控地震波场不同之处在于,该波场为人工合成的波场,而非物理波场。因此定向照明地震勘探方法,在保证波场方向性特点的基础上,克服了阵列震源的不一致性和高成本的缺点,势必为高分辨率地震勘探提供一种有效方法。The existing vibrator seismic exploration methods can be divided into two categories according to the number of seismic sources involved in the work. One is single vibrator seismic exploration; the other is array vibroseis seismic exploration. The seismic wavefield excited by a single vibroseis in the underground half-space can be approximated as a spherical wavefield with no directionality; the seismic wavefield excited by an array vibroseis has obvious directionality, and the main beam direction of the seismic wavefield excited by the combined vibrator is perpendicular to the ground. Under normal circumstances, the seismic wave field excited by the phase-controlled vibrator also has directionality, and the direction is controllable; generally speaking, the seismic wave direction can be reasonably controlled, and the seismic data with a higher signal-to-noise ratio can be obtained by using the phase-controlled vibrator. However, the phase-controlled seismic source seismic exploration method has high cost, complex control, and the inconsistency of multiple sources, which limit the promotion of the phase-controlled seismic source seismic exploration method. For this reason, it is still considered to use a single seismic source as the excitation source, and use the single-shot seismic signal based on the single seismic source as the data source, and use the directional illumination seismic signal synthesis method to synthesize the directional illumination single-channel seismic signal, and then synthesize the directional illumination single-shot seismic data similarly , the data correspond to a directional seismic wavefield. The similarity with the phased seismic wavefield is that both wavefields are directional wavefields, so both are effective for improving the signal-to-noise ratio of seismic data; the difference with the phased seismic wavefield is that the wavefield is artificially synthesized Wavefields, not physical wavefields. Therefore, the directional illumination seismic exploration method, on the basis of ensuring the directional characteristics of the wave field, overcomes the shortcomings of the inconsistency and high cost of the array source, and is bound to provide an effective method for high-resolution seismic exploration.

发明内容:Invention content:

本发明的目的是针对普通单震源和阵列震源地震勘探方法的缺点和不足,提供一种通过合成定向地震波场,改善反射波信号信噪比的基于单震源的定向照明地震勘探方法。The purpose of the present invention is to provide a single-source directional lighting seismic exploration method that improves the signal-to-noise ratio of reflected wave signals by synthesizing directional seismic wavefields and aiming at the shortcomings and deficiencies of common single-source and array-source seismic prospecting methods.

本发明的目的是通过以下方式实现的:The purpose of the present invention is achieved in the following manner:

基于单震源的定向照明地震勘探方法包括下列顺序和步骤:The directional illumination seismic exploration method based on single source includes the following sequence and steps:

a、野外实施勘探时,依据地质资料按照传统地震勘探方法确定偏移距p、炮点距d和道间距及每一炮的炮点位置;这里p是指第1个炮点距第1个检波器的距离,即偏移距p;a. When carrying out exploration in the field, according to the geological data, determine the offset distance p, shot point distance d, trace spacing and the shot point position of each shot according to the traditional seismic prospecting method; here p refers to the distance between the first shot point and the first shot point The distance of the detector, that is, the offset p;

b、依据地质背景和勘探要求,确定定向照明组数m和合成一个定向照明炮点所需的单震源炮点数n,因此定向照明地震勘探包含的总炮数N=n×m,即将定向照明地震勘探总炮数分成m组,每组包含n炮单震源地震数据;b. According to the geological background and exploration requirements, determine the number of directional lighting groups m and the number of single-source shot points n required to synthesize a directional lighting shot point. Therefore, the total number of shots included in the directional lighting seismic survey is N=n×m, that is, the directional lighting The total number of seismic exploration shots is divided into m groups, and each group contains n shots single-source seismic data;

c、按照步骤a中设定的的偏移距p,炮点距d和道间距,在测线上布置检波器阵列和炮点;c. According to the offset p set in step a, the shot distance d and the track spacing, arrange the detector array and the shot point on the survey line;

d、若选择n=3个单震源炮点数据合成单炮定向照明炮点数据,首先布置第1至第n个炮点,这n个炮点可视为定向照明勘探的第1组炮点,按照传统地质勘探方法布置检波器阵列,并根据指定偏移距p,设置第1个炮点,该炮点距离第一个检波器的距离为:p1=p,据此取得可控震源第1炮地震数据,保持检波器阵列位置不变,部置第2个炮点,该炮点距离第一个检波器的距离为p2=p1-d,取得可控震源取得第2炮地震数据;保持检波器阵列位置不变,部置第3个炮点,该炮点距离第一个检波器的距离为p3=p2-d=p1-2d,取得可控震源第3炮地震数据,此时完成了第1组n个炮点的地震数据采集;d. If n=3 single-source shot point data are selected to synthesize single-shot directional lighting shot point data, the first to nth shot points are arranged first, and these n shot points can be regarded as the first group of shot points for directional lighting exploration , arrange the geophone array according to the traditional geological exploration method, and set the first shot point according to the specified offset p, the distance between the shot point and the first geophone is: p 1 = p, and the vibroseis For the seismic data of the first shot, keep the position of the geophone array unchanged, install the second shot point, the distance between the shot point and the first geophone is p 2 =p 1 -d, obtain the vibrator and get the second shot Seismic data: keep the location of the geophone array unchanged, install the third shot point, the distance between the shot point and the first geophone is p 3 =p 2 -d=p 1 -2d, and obtain the third vibroseis Shot seismic data. At this time, the seismic data acquisition of the first group of n shot points has been completed;

e、按照传统地震勘探方法,滚动检波器阵列,接着布置第n+1至第2n个炮点,这n个炮点可视为定向照明勘探的第2组炮点,在同组各炮点的数据采集过程中,检波器阵列道间距始终保持不变,第n+1个炮点设置在距离第1个检波器的距离为:p1=p的位置,即与第1组的第1个炮点偏移距相同,取得可控震源第4炮地震数据;第n+2个炮点设置在距距离第1个检波器的距离为:p2=p1-d的位置,即与第1组的第2个炮点偏移距相同,取得可控震源第5炮地震数据;第n+3个炮点设置在距距离第1个检波器的距离为:p3=p2-d=p1-2d的位置,取得可控震源第6炮地震数据;此时完成了第2组n个炮点的地震数据采集;即同一组内各炮点的偏移距不同,不同组序号相差n的m个炮点的偏移距相同;e. According to the traditional seismic exploration method, roll the geophone array, and then arrange the n+1th to 2nth shot points. These n shot points can be regarded as the second group of shot points for directional lighting exploration. During the data acquisition process of , the track spacing of the geophone array remains unchanged, and the n+1th shot point is set at a distance from the 1st geophone: p 1 =p, that is, the distance from the 1st geophone of the 1st group The offsets of the two shot points are the same, and the seismic data of the fourth shot of the vibrator are obtained; the n+2th shot point is set at a distance from the first geophone: p 2 =p 1 -d, that is, the same as The offset distance of the second shot point in the first group is the same, and the seismic data of the fifth shot of the vibrator is obtained; the n+3th shot point is set at a distance from the first geophone: p 3 =p 2 - At the position of d=p 1 -2d, the seismic data of the sixth shot of the vibrator is obtained; at this time, the seismic data acquisition of n shot points in the second group is completed; that is, the offset distances of each shot point in the same group are different, and the The offset distances of the m shot points whose sequence numbers differ by n are the same;

f、依次类推,继续滚动检波器阵列;f, and so on, continue to scroll the detector array;

g、布置下一组的n个炮点,设该组的编号为i,对应第i组的n个炮点在总炮点中的排列序号即为第(i-1)×n+1至第(i-1)×n+n个炮点,布置第(i-1)×n+1个炮点,该炮点距离第一个检波器的距离为p1=p;再布置第(i-1)×n+2个炮点,该炮点距第一个检波器的距离为p2=p1-d;……,再布置第(i-1)×n+n个炮点,该炮点距第一个检波器的距离为:g, arranging the n shot points of the next group, assuming that the group number is i, the sequence number of the n shot points corresponding to the i group in the total shot points is (i-1)×n+1 to For the (i-1)×n+n shot point, arrange the (i-1)×n+1 shot point, the distance between the shot point and the first geophone is p 1 =p; then arrange the ( i-1)×n+2 shot points, the distance between the shot point and the first geophone is p 2 =p 1 -d; ..., then arrange the (i-1)×n+nth shot point , the distance between the shot point and the first geophone is:

pn=pn-1-d=p1-(n-1)×d;在以上n个炮点取得第(i-1)×n+1至第(i-1)×n+n个炮点的地震数据;p n =p n-1 -d=p 1 -(n-1)×d; obtain the (i-1)×n+1th to (i-1)×n+nth shots at the above n shot points Seismic data of the shot point;

h、重复步骤f、g直至取得所有N炮地震数据h. Repeat steps f and g until all N shot seismic data are obtained

i、将取得的单炮地震数据分成m组,对属于同一组的单炮地震数据按编号顺序依次作等间隔延时,将延时后的n个单炮地震数据做线性叠加,合成具有方向性的定向照明地震信号,该信号对应定向地震波场,利用常规地震数据处理方法就能得到基于单震源的定向照明地震剖面。i. Divide the obtained single-shot seismic data into m groups, perform equal-interval time delay on the single-shot seismic data belonging to the same group according to the order of numbering, do linear superposition of the n single-shot seismic data after the delay, and synthesize a direction The directional illumination seismic signal corresponds to the directional seismic wave field, and the directional illumination seismic section based on a single source can be obtained by using the conventional seismic data processing method.

本发明的目的还可以通过以下方式实现:The purpose of the present invention can also be achieved in the following ways:

炮点距d等于或不等于道间距;步骤b中的m为>1的奇数。The shot pitch d is equal to or not equal to the track pitch; m in step b is an odd number >1.

有益效果:经试验,与普通单震源地震勘探方法相比,定向照明地震勘探方法可使来自目标体的反射波信号的信噪比得到改善;由于合成定向照明地震波场时方向可根据目标产状人为控制,对于陡倾角目标体勘探效果更好;与相控地震方法相比,具有工作成本低,控制简单的优点。Beneficial effects: After testing, compared with the ordinary single source seismic exploration method, the directional illumination seismic exploration method can improve the signal-to-noise ratio of the reflected wave signal from the target body; because the direction of the synthetic directional illumination seismic wave field can be adjusted according to the occurrence of the target Artificial control is better for the exploration of objects with steep dip angles; compared with the phase-controlled seismic method, it has the advantages of low work cost and simple control.

附图说明Description of drawings

图1:定向照明地震勘探不同组炮点和检波器布置图Figure 1: Layout of different groups of shot points and geophones for directional illumination seismic exploration

图2:定向照明同一组炮点和检波器布置图Figure 2: Layout of the same set of shot points and detectors for directional lighting

m:定向照明组数,即定向照明地震方法中对总炮数的分组数m: the number of directional lighting groups, that is, the grouping number of the total shot number in the directional lighting seismic method

n:合成一个定向照明炮点所需的单震源炮点数,即总炮数分组中每定向照明合成组的炮点数n: the number of single-source shot points required to synthesize a directional lighting shot point, that is, the number of shot points per directional lighting synthesis group in the total shot number group

N:测线上总炮点数N: total number of shot points on the survey line

d:炮点距d: shot point distance

P:偏移距P: Offset

具体实施方式:Detailed ways:

下面结合附图和实施例作进一步详细说明:Below in conjunction with accompanying drawing and embodiment describe in further detail:

基于单震源的定向照明地震勘探方法包括下列顺序和步骤:The directional illumination seismic exploration method based on single source includes the following sequence and steps:

a、野外实施勘探时,依据地质资料按照传统地震勘探方法确定偏移距p、炮点距d和道间距及每一炮的炮点位置;这里p是指第1个炮点距第1个检波器的距离,即偏移距p;a. When carrying out exploration in the field, according to the geological data, determine the offset distance p, shot point distance d, trace spacing and the shot point position of each shot according to the traditional seismic prospecting method; here p refers to the distance between the first shot point and the first shot point The distance of the detector, that is, the offset p;

b、依据地质背景和勘探要求,确定定向照明组数m和合成一个定向照明炮点所需的单震源炮点数n,因此定向照明地震勘探包含的总炮数N=n×m,即将定向照明地震勘探总炮数分成m组,每组包含n炮单震源地震数据;b. According to the geological background and exploration requirements, determine the number of directional lighting groups m and the number of single-source shot points n required to synthesize a directional lighting shot point. Therefore, the total number of shots included in the directional lighting seismic survey is N=n×m, that is, the directional lighting The total number of seismic exploration shots is divided into m groups, and each group contains n shots single-source seismic data;

c、按照步骤a中设定的的偏移距p,炮点距d和道间距,在测线上布置检波器阵列和炮点;c. According to the offset p set in step a, the shot distance d and the track spacing, arrange the detector array and the shot point on the survey line;

d、若选择n=3个单震源炮点数据合成单炮定向照明炮点数据,首先布置第1至第n个炮点,这n个炮点可视为定向照明勘探的第1组炮点,按照传统地质勘探方法布置检波器阵列,并根据指定偏移距p,设置第1个炮点,该炮点距离第一个检波器的距离为:p1=p,据此取得可控震源第1炮地震数据,保持检波器阵列位置不变,部置第2个炮点,该炮点距离第一个检波器的距离为p2=p1-d,取得可控震源取得第2炮地震数据;保持检波器阵列位置不变,部置第3个炮点,该炮点距离第一个检波器的距离为p3=p2-d=p1-2d,取得可控震源第3炮地震数据,此时完成了第1组n个炮点的地震数据采集;d. If n=3 single-source shot point data are selected to synthesize single-shot directional lighting shot point data, the first to nth shot points are arranged first, and these n shot points can be regarded as the first group of shot points for directional lighting exploration , arrange the geophone array according to the traditional geological exploration method, and set the first shot point according to the specified offset p, the distance between the shot point and the first geophone is: p 1 = p, and the vibroseis For the seismic data of the first shot, the position of the geophone array remains unchanged, and the second shot point is placed. The distance between the shot point and the first geophone is p 2 =p 1 -d, and the vibroseis is obtained to obtain the second shot Seismic data: keep the location of the geophone array unchanged, install the third shot point, the distance between the shot point and the first geophone is p 3 =p 2 -d=p 1 -2d, and obtain the third shot point of the vibroseis Shot seismic data. At this time, the seismic data acquisition of the first group of n shot points has been completed;

e、按照传统地震勘探方法,滚动检波器阵列,接着布置第n+1至第2n个炮点,这n个炮点可视为定向照明勘探的第2组炮点,在同组各炮点的数据采集过程中,检波器阵列道间距始终保持不变,第n+1个炮点设置在距离第1个检波器的距离为:p1=p的位置,即与第1组的第1个炮点偏移距相同,取得可控震源第4炮地震数据;第n+2个炮点设置在距距离第1个检波器的距离为:p2=p1-d的位置,即与第1组的第2个炮点偏移距相同,取得可控震源第5炮地震数据;第n+3个炮点设置在距距离第1个检波器的距离为:p3=p2-d=p1-2d的位置,取得可控震源第6炮地震数据;此时完成了第2组n个炮点的地震数据采集;即同一合成组内各炮点的偏移距不同,不同合成组序号相差n的m个的炮点偏移距相同;e. According to the traditional seismic exploration method, roll the geophone array, and then arrange the n+1th to 2nth shot points. These n shot points can be regarded as the second group of shot points for directional lighting exploration. During the data acquisition process of , the track spacing of the geophone array remains unchanged, and the n+1th shot point is set at a distance from the 1st geophone: p 1 =p, that is, it is the same as the 1st The offsets of the two shot points are the same, and the seismic data of the fourth shot of the vibrator is obtained; the n+2th shot point is set at a distance from the first geophone: p 2 =p 1 -d, that is, the same as The offset distance of the second shot point in the first group is the same, and the seismic data of the fifth shot of the vibrator is obtained; the n+3th shot point is set at a distance from the first geophone: p 3 =p 2 - At the position of d=p 1 -2d, the seismic data of the sixth shot of the vibrator is obtained; at this time, the seismic data acquisition of n shot points in the second group is completed; that is, the offsets of each shot point in the same synthetic group are different, The offset distances of the m shot points with a difference of n in the synthesized group numbers are the same;

f、依次类推,继续滚动检波器阵列;f, and so on, continue to scroll the detector array;

g、布置下一组的n个炮点,设该组的编号为i,对应第i组的n个炮点在总炮点中的排列序号即为第(i-1)×n+1至第(i-1)×n+n个炮点,布置第(i-1)×n+1个炮点,该炮点距离第一个检波器的距离为p1=p;再布置第(i-1)×n+2个炮点,该炮点距第一个检波器的距离为p2=p1-d;……,再布置第(i-1)×n+n个炮点,该炮点距第一个检波器的距离为:g, arranging the n shot points of the next group, assuming that the group number is i, the sequence number of the n shot points corresponding to the i group in the total shot points is (i-1)×n+1 to For the (i-1)×n+n shot point, arrange the (i-1)×n+1 shot point, the distance between the shot point and the first geophone is p 1 =p; then arrange the ( i-1)×n+2 shot points, the distance between the shot point and the first geophone is p 2 =p 1 -d; ..., then arrange the (i-1)×n+nth shot point , the distance between the shot point and the first geophone is:

pn=pn-1-d=p1-(n-1)×d;在以上n个炮点取得第(i-1)×n+1至第(i-1)×n+n个炮点的地震数据;p n =p n-1 -d=p 1 -(n-1)×d; obtain the (i-1)×n+1th to (i-1)×n+nth shots at the above n shot points Seismic data of the shot point;

h、重复步骤f、g直至取得所有N炮地震数据h. Repeat steps f and g until all N shot seismic data are obtained

i、将取得的单炮地震数据分成m组,对属于同一组的单炮地震数据按编号顺序依次作等间隔延时,将延时后的n个单炮地震数据做线性叠加,合成具有方向性的定向照明地震信号,该信号对应定向地震波场,利用常规地震数据处理方法就能得到基于单震源的定向照明地震剖面。i. Divide the obtained single-shot seismic data into m groups, perform equal-interval time delay on the single-shot seismic data belonging to the same group according to the order of numbering, do linear superposition of the n single-shot seismic data after the delay, and synthesize a direction The directional illumination seismic signal corresponds to the directional seismic wave field, and the directional illumination seismic section based on a single source can be obtained by using the conventional seismic data processing method.

m:定向照明组数,即定向照明地震方法中对总炮数的分组数m: the number of directional lighting groups, that is, the grouping number of the total shot number in the directional lighting seismic method

n:合成一个定向照明炮点所需的单震源炮点数,即总炮数分组中每定向照明合成组的炮点数n: the number of single-source shot points required to synthesize a directional lighting shot point, that is, the number of shot points per directional lighting synthesis group in the total shot number group

N:测线上总炮点数N: total number of shot points on the survey line

d:炮点距d: shot point distance

P:偏移距P: Offset

本发明的目的还可以通过以下方式实现:The purpose of the present invention can also be achieved in the following ways:

炮点距d等于或不等于道间距;步骤b中的m为>1的奇数。The shot pitch d is equal to or not equal to the track pitch; m in step b is an odd number >1.

同一组内各炮点的偏移距是不同的,第一炮的偏移距为p,第二炮的偏移距为p2=p-d,第三炮的偏移距为p3=p-2d,第四炮的偏移距为p4=p-3d。不同组各炮点的偏移距是相同的,第一组第一炮偏移距为p,第二组第一炮偏移距也为p,第一组第二炮的偏移距为p2=p-d,第二组第二炮的偏移距也为p2=p-d,第一组第三炮的偏移距为p3=p-2d,第二组第三炮的偏移距为p3=p-2d.......。The offset distance of each shot point in the same group is different, the offset distance of the first shot is p, the offset distance of the second shot is p2=p-d, the offset distance of the third shot is p3=p-2d, The offset of the fourth shot is p4=p-3d. The offset distance of each shot point in different groups is the same, the offset distance of the first shot in the first group is p, the offset distance of the first shot in the second group is also p, and the offset distance of the second shot in the first group is p2 =p-d, the offset distance of the second group of second guns is also p2=p-d, the offset distance of the first group of third guns is p3=p-2d, the offset distance of the second group of third guns is p3=p -2d.......

根据常规反射地震勘探方法设计炮检距和道间距;参考测区地质资料大致估计目标深度和形态,粗略确定定向照明延时控制参数τ的取值范围;选取定向照明叠加控制参数n,一般地n的取值与测区噪声强度有关,噪声越强n取值越大,而伴随n取值增大,需要获取的原始单炮记录和合成定向照明地震信号的算法运算量随之增大,可见n的取值应根据地震勘探分辨率的要求与工作量及成本两方面折衷考虑。Design the offset and trace spacing according to the conventional reflection seismic exploration method; roughly estimate the depth and shape of the target with reference to the geological data of the survey area, and roughly determine the value range of the directional lighting delay control parameter τ; select the directional lighting superposition control parameter n, generally The value of n is related to the noise intensity of the survey area. The stronger the noise is, the larger the value of n is, and with the increase of the value of n, the computational complexity of the original single-shot records and synthetic directional illumination seismic signals that need to be acquired increases accordingly. It can be seen that the value of n should be considered as a compromise between the requirements of seismic exploration resolution and workload and cost.

在测线上设置m组炮点,每组包含n炮,即总共包含N=n×m个炮点,各炮点依次标号为1,2,…,N。检波器共线均匀排列,检波器阵列与炮点共线。同一合成组检波器阵列位置保持不变,炮点距第一个检波器的距离是变的,同一合成组内各炮点的偏移距不同,即同一组内的n个炮点偏移距依次为p,p-d,p-2d,…,p-(n-1)d。不同合成组的炮点和检波器阵列位置按设定的炮点距d和道间距滚动。不同组序号相差n的m个炮点的偏移距相同,即第1,n+1,2n+1,…,(m-1)×n+1这m个炮点的偏移距为p;第2,n+2,2n+2,…,(m-1)×n+2这m个炮点的偏移距为p-d;第3,n+3,2n+3,…,(m-1)×n+3这m个炮点的偏移距为p-2d;…,以此类推,第n,2n,3n,…,m×n这n个炮点的偏移距为p-(n-1)d。Set up m groups of shot points on the survey line, each group contains n shots, that is, it contains N=n×m shot points in total, and each shot point is numbered 1, 2, ..., N in turn. The geophones are arranged collinearly and evenly, and the geophone array and the shot point are collinear. The position of the geophone array in the same synthesis group remains unchanged, the distance between the shot point and the first geophone is variable, and the offset distances of each shot point in the same synthesis group are different, that is, the offset distance of n shot points in the same group In turn, p, p-d, p-2d, ..., p-(n-1)d. The positions of shot points and geophone arrays of different synthesis groups are scrolled according to the set shot point distance d and track spacing. The offset distances of the m shot points with a difference of n in different groups of serial numbers are the same, that is, the offset distance of the m shot points of the 1st, n+1, 2n+1, ..., (m-1)×n+1 is p ; The 2nd, n+2, 2n+2, ..., (m-1) × n+2 the offset distance of the m shot points is p-d; the 3rd, n+3, 2n+3, ..., (m -1)×n+3, the offset distance of the m shot points is p-2d; ..., and so on, the offset distance of the nth, 2n, 3n, ..., m×n shot points is p -(n-1)d.

具体地,使用单可控震源按以上排列方法,分别获取m组总共N炮单可控震源单炮地震记录。采用定向照明地震信号合成方法,利用第1组标号为1至n的n个单震源单炮记录合成第1炮定向照明单炮记录,该记录对应的等效定向照明炮点位置为第(n+1)/2个炮点所在位置,即第1组中间炮点的位置;利用第2组标号为n+1至2n的n个单震源单炮记录合成第2炮定向照明单炮记录,该记录对应的等效定向照明炮点位置为第n+(n+1)/2个炮点所在位置,即第2组中间炮点的位置;…;类似地,利用第m组标号为(m-1)n+1至m×n的n个单震源单炮记录合成第m炮定向照明单炮记录,该记录对应的等效定向照明炮点位置为第(m-1)n+(n+1)/2个炮点所在位置,即第m组中间炮点的位置;按照上述方法可合成m组定向照明单炮记录。其合成方法已被CN200810050992.3公开。Specifically, single-vibrator single-shot seismic records with a total of N shots in m groups were obtained using the single-vibrator according to the above arrangement method. Using the directional illumination seismic signal synthesis method, using the first group of n single-source single-shot records labeled 1 to n to synthesize the first shot directional illumination single-shot record, the corresponding equivalent directional illumination shot point position of the record is (n +1)/2 The position of the shot point is the position of the middle shot point of the first group; use the n single-source single-shot records labeled n+1 to 2n in the second group to synthesize the second shot directional illumination single-shot record, The position of the equivalent directional lighting shot point corresponding to this record is the position of the n+(n+1)/2 shot point, that is, the position of the middle shot point of the second group; ...; similarly, using the m group label as (m -1) n single-source single-shot records from n+1 to m×n are combined to synthesize the m-th shot directional lighting single-shot record, and the corresponding equivalent directional lighting shot point position of the record is (m-1)n+(n+ 1)/2 The position of the shot point, that is, the position of the middle shot point of the m group; according to the above method, the m group of directional illumination single shot records can be synthesized. Its synthesis method has been disclosed by CN200810050992.3.

通过对定向照明地震信号实施互相关检测,再利用常规地震数据处理方法及软件就可以得到地震剖面。Seismic sections can be obtained by performing cross-correlation detection on seismic signals with directional illumination and using conventional seismic data processing methods and software.

定向照明地震波场方向的计算方法:Directional Lighting Seismic Wavefield Direction Calculation method:

根据测区地质资料或单震源地震数据,可估计定向照明方向。对于近似水平的目标层,设炮检距为x,目标体深度为h,此时定向照明方向可按下式估计According to the geological data of the survey area or the seismic data of a single source, the direction of directional lighting can be estimated. For an approximately horizontal target layer, set the offset as x and the depth of the target as h, and at this time the directional lighting direction can be estimated by

Figure G200810051062XD00053
(公式1),
Figure G200810051062XD00053
(Formula 1),

定向照明地震延时控制参数τ的估算方法:Estimation method of directional lighting seismic delay control parameter τ:

设炮点距为d,对于近似水平的目标层,目标层的上覆介质平均速度为v时,延时参数可按下式估算。Assuming that the shot distance is d, for an approximately horizontal target layer, when the average velocity of the overlying medium of the target layer is v, the delay parameter can be estimated by the following formula.

Figure G200810051062XD00054
(公式2)
Figure G200810051062XD00054
(Formula 2)

定向照明单炮地震记录合成方法:Synthesis method of single-shot seismic records with directional illumination:

采用n炮单震源数据合成1炮定向照明地震数据,合成方法已被CN200810050992.3公开。Synthesize 1-shot directional illumination seismic data by using n-shot single-source data, and the synthesis method has been disclosed by CN200810050992.3.

首先用n炮单震源数据合成1炮定向照明地震数据,分析定向照明地震信号质量。若该信号满足预期的信噪比要求,则确定n炮单震源数据为一个合成组,合成1炮定向照明地震数据;否则,调整参数,重新确定一个合成组所需定向照明的炮数。最后,利用合成的m组炮定向照明地震数据,采用常规地震数据处理方法合成地震剖面。Firstly, single source data of n shots are used to synthesize 1 shot directional illumination seismic data, and the quality of directional illumination seismic signal is analyzed. If the signal meets the expected signal-to-noise ratio requirements, determine n-shot single source data as a composite group, and synthesize 1-shot directional illumination seismic data; otherwise, adjust the parameters to re-determine the number of directional illumination shots required for a composite group. Finally, using the synthesized m-shot directional illumination seismic data, the conventional seismic data processing method is used to synthesize the seismic section.

Claims (3)

1.一种基于单震源的定向照明地震勘探方法,其特征在于,包括下列顺序和步骤:1. A method for directional illumination seismic exploration based on a single source, characterized in that it comprises the following sequence and steps: a、野外实施勘探时,依据地质资料按照传统地震勘探方法确定偏移距p、炮点距d和道间距及每一炮的炮点位置;这里p是指第1个炮点距第1个检波器的距离,即偏移距p;a. When carrying out exploration in the field, according to the geological data, determine the offset distance p, shot point distance d, trace spacing and the shot point position of each shot according to the traditional seismic prospecting method; here p refers to the distance between the first shot point and the first shot point The distance of the detector, that is, the offset p; b、依据地质背景和勘探要求,确定定向照明组数m和合成一个定向照明炮点所需的单震源炮点数n,因此定向照明地震勘探包含的总炮数N=n×m,即将定向照明地震勘探总炮数分成m组,每组包含n炮单震源地震数据;b. According to the geological background and exploration requirements, determine the number of directional lighting groups m and the number of single-source shot points n required to synthesize a directional lighting shot point. Therefore, the total number of shots included in the directional lighting seismic survey is N=n×m, that is, the directional lighting The total number of seismic exploration shots is divided into m groups, and each group contains n shots single-source seismic data; c、按照步骤a中设定的的偏移距p,炮点距d和道间距,在测线上布置检波器阵列和炮点;c. According to the offset p set in step a, the shot distance d and the track spacing, arrange the detector array and the shot point on the survey line; d、利用n=3个单震源炮点数据合成单炮定向照明炮点数据,首先布置第1至第n个炮点,这n个炮点可视为定向照明勘探的第1组炮点,按照传统地质勘探方法布置检波器阵列,并根据指定偏移距p,设置第1个炮点,该炮点距离第一个检波器的距离为:p1=p,据此取得可控震源第1炮地震数据,保持检波器阵列位置不变,布置第2个炮点,该炮点距离第一个检波器的距离为p2=p1-d,取得可控震源第2炮地震数据;保持检波器阵列位置不变,布置第3个炮点,该炮点距离第一个检波器的距离为p3=p2-d=p1-2d,取得可控震源第3炮地震数据,此时完成了第1组n个炮点的地震数据采集;d. Synthesize single-shot directional lighting shot point data by using n=3 single-source shot point data, and first arrange the 1st to nth shot points, and these n shot points can be regarded as the first group of shot points for directional lighting exploration, The geophone array is arranged according to the traditional geological prospecting method, and the first shot point is set according to the specified offset p. The distance between the shot point and the first geophone is: p 1 = p. For the seismic data of the first shot, keep the position of the geophone array unchanged, arrange the second shot point, the distance between the shot point and the first geophone is p 2 =p 1 -d, and obtain the seismic data of the second shot of the vibrator; Keep the location of the geophone array unchanged, arrange the third shot point, the distance between the shot point and the first geophone is p 3 =p 2 -d=p 1 -2d, and obtain the seismic data of the third shot of the vibrator, At this point, the seismic data acquisition of the first group of n shot points is completed; e、按照传统地震勘探方法,滚动检波器阵列,接着布置第n+1至第2n个炮点,这n个炮点可视为定向照明勘探的第2组炮点,在同组各炮点的数据采集过程中,检波器阵列道间距始终保持不变,第n+1个炮点设置在距离第1个检波器的距离为:p1=p的位置,即与第1组的第1个炮点偏移距相同,取得可控震源第4炮地震数据;第n+2个炮点设置在距离第1个检波器的距离为:p2=p1-d的位置,即与第1组的第2个炮点偏移距相同,取得可控震源第5炮地震数据;第n+3个炮点设置在距离第1个检波器的距离为:p3=p2-d=p1-2d的位置,取得可控震源第6炮地震数据;此时完成了第2组n个炮点的地震数据采集,即同一组内各炮点的偏移距不同,不同合成组序号相差n的m个炮点的偏移距相同;e. According to the traditional seismic exploration method, roll the geophone array, and then arrange the n+1th to 2nth shot points. These n shot points can be regarded as the second group of shot points for directional lighting exploration. During the data acquisition process of , the track spacing of the geophone array remains unchanged, and the n+1th shot point is set at a distance from the 1st geophone: p 1 =p, that is, the distance from the 1st geophone of the 1st group The offset distances of the two shot points are the same, and the seismic data of the fourth shot of the vibrator are obtained ; the n+2th shot point is set at a distance from the first geophone: The offset distance of the second shot point in group 1 is the same, and the seismic data of the fifth shot of the vibrator is obtained; the n+3th shot point is set at a distance from the first geophone: p 3 =p 2 -d= At the position of p 1 -2d, the seismic data of the sixth shot of the vibrator is obtained; at this time, the seismic data acquisition of the second group of n shot points is completed, that is, the offset distances of each shot point in the same group are different, and the sequence numbers of different composite groups The offsets of the m shot points with a difference of n are the same; f、依次类推,继续滚动检波器阵列;f, and so on, continue to scroll the detector array; g、布置下一组的n个炮点,设该组的编号为i,对应第i组的n个炮点在总炮点中的排列序号即为第(i-1)×n+1至第(i-1)×n+n个炮点,布置第(i-1)×n+1个炮点,该炮点距离第一个检波器的距离为p1=p;再布置第(i-1)×n+2个炮点,该炮点距第一个检波器的距离为p2=p1-d;……,再布置第(i-1)×n+n个炮点,该炮点距第一个检波器的距离为:g, arranging the n shot points of the next group, assuming that the group number is i, the sequence number of the n shot points corresponding to the i group in the total shot points is (i-1)×n+1 to For the (i-1)×n+n shot point, arrange the (i-1)×n+1 shot point, the distance between the shot point and the first geophone is p 1 =p; then arrange the ( i-1)×n+2 shot points, the distance between the shot point and the first geophone is p 2 =p 1 -d; ..., then arrange the (i-1)×n+nth shot point , the distance between the shot point and the first geophone is: pn=pn-1-d=p1-(n-1)×d;在以上n个炮点取得第(i-1)×n+1至第(i-1)×n+n个炮点的地震数据;p n =p n-1 -d=p 1 -(n-1)×d; obtain the (i-1)×n+1th to (i-1)×n+nth shots at the above n shot points Seismic data of the shot point; h、重复步骤f、g直至取得所有N炮地震数据;h, repeating steps f and g until all N shot seismic data are obtained; i、将取得的单炮地震数据分成m组,对属于同一组的单炮地震数据按编号顺序依次作等间隔延时,将延时后的n个单炮地震数据做线性叠加,合成具有方向性的定向照明地震信号,该信号对应定向地震波场,利用常规地震数据处理方法就能得到基于单震源的定向照明地震剖面。i. Divide the obtained single-shot seismic data into m groups, perform equal-interval time delay on the single-shot seismic data belonging to the same group according to the order of numbering, do linear superposition of the n single-shot seismic data after the delay, and synthesize a direction The directional illumination seismic signal corresponds to the directional seismic wave field, and the directional illumination seismic section based on a single source can be obtained by using the conventional seismic data processing method. 2.按照权利要求1所述的基于单震源的定向照明地震勘探方法,其特征在于,炮点距d等于或不等于道间距。2. The single-source-based directional illumination seismic exploration method according to claim 1, characterized in that the shot distance d is equal to or not equal to the trace distance. 3.按照权利要求1所述的基于单震源的定向照明地震勘探方法,其特征在于,步骤b中的m为>1的奇数。3. The single-source-based directional illumination seismic exploration method according to claim 1, characterized in that m in step b is an odd number >1.
CN200810051062XA 2008-08-08 2008-08-08 Directional lighting seismic exploration method based on single earthquake focus Expired - Fee Related CN101339252B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200810051062XA CN101339252B (en) 2008-08-08 2008-08-08 Directional lighting seismic exploration method based on single earthquake focus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200810051062XA CN101339252B (en) 2008-08-08 2008-08-08 Directional lighting seismic exploration method based on single earthquake focus

Publications (2)

Publication Number Publication Date
CN101339252A CN101339252A (en) 2009-01-07
CN101339252B true CN101339252B (en) 2010-08-18

Family

ID=40213366

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200810051062XA Expired - Fee Related CN101339252B (en) 2008-08-08 2008-08-08 Directional lighting seismic exploration method based on single earthquake focus

Country Status (1)

Country Link
CN (1) CN101339252B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102590861B (en) * 2012-01-17 2014-11-19 吉林大学 Seismic beamforming method based on receiving array
CN103336302B (en) * 2013-06-27 2015-07-08 吉林大学 Seismic beamforming method based on high-order cosine amplitude weighting
CN103984007B (en) * 2014-06-09 2016-07-13 吉林大学 Optimal design method of directional seismic wave time-delay parameters
CN104570121B (en) * 2014-12-26 2017-01-04 吉林大学 Directionally seismic wave distorted signal removing method
CN104793243B (en) * 2015-04-23 2017-05-10 吉林大学 Directional seismic data processing method based on Nth root stacking
CN105629300B (en) * 2015-12-30 2018-05-15 中国石油集团东方地球物理勘探有限责任公司 The method for improving complicated structure offset data signal-to-noise ratio
CN111965729B (en) * 2019-05-20 2023-04-07 中国石油天然气集团有限公司 Real-time monitoring method, system and device for vibroseis combination center
CN111221040B (en) * 2020-02-28 2021-06-22 吉林大学 A kind of formation dip detection method and system
CN112505749B (en) * 2020-10-19 2024-04-26 中国地质调查局南京地质调查中心(华东地质科技创新中心) Micro-motion data acquisition method based on linear array multiple coverage
CN112965193A (en) * 2021-02-04 2021-06-15 广州海洋地质调查局 Ocean seismic streamer based on DAS optical fiber

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1185063A (en) * 1966-03-21 1970-03-18 Texas Instruments Inc Seismic Exploration
CN1259208A (en) * 1997-06-13 2000-07-05 菲利浦石油公司 Method for sorting seismic data
CN1560651A (en) * 2004-03-03 2005-01-05 吉林大学 Phased Array Seismic Exploration Method
CN1560650A (en) * 2004-03-03 2005-01-05 吉林大学 Phase Control Vibroseis System
CN1797031A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for shifting depth before superposition in combined earthquake data from multiple times of shooting

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1185063A (en) * 1966-03-21 1970-03-18 Texas Instruments Inc Seismic Exploration
CN1259208A (en) * 1997-06-13 2000-07-05 菲利浦石油公司 Method for sorting seismic data
CN1560651A (en) * 2004-03-03 2005-01-05 吉林大学 Phased Array Seismic Exploration Method
CN1560650A (en) * 2004-03-03 2005-01-05 吉林大学 Phase Control Vibroseis System
CN1797031A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for shifting depth before superposition in combined earthquake data from multiple times of shooting

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
姜弢.基于相控震源的地震波定向方法研究.《中国优秀博硕士学位论文全文数据库(博士)基础科学辑》.2006,(第10期),1-104. *
杨冬,姜弢.可控震源定向照明方法的仿真研究.《地球物理学进展》.2008,第23卷(第3期),822-825. *

Also Published As

Publication number Publication date
CN101339252A (en) 2009-01-07

Similar Documents

Publication Publication Date Title
CN101339252B (en) Directional lighting seismic exploration method based on single earthquake focus
US7706211B2 (en) Method of determining a seismic velocity profile
AU2004304332B2 (en) Methods for acquiring and processing seismic data from quasi-simultaneously activated translating energy sources
US6574564B2 (en) 3D prestack seismic data migration method
US9052410B2 (en) Multiple seismic signal inversion
US9310503B2 (en) Methods to process seismic data contaminated by coherent energy radiated from more than one source
US8665667B2 (en) Vertical seismic profiling velocity estimation method
US7929373B2 (en) Method of wavefield extrapolation for single-station, dual-sensor towed streamer signals
Vinje et al. Shooting over the seismic spread
US20080130411A1 (en) Seismic imaging with natural green's functions derived from vsp data
US20090147620A1 (en) Method for processing marine towed streamer seismic data from regular multi-azimuth surveys
US20030067842A1 (en) Helix streamer acquisition of seismic data
US7969818B2 (en) Method for regularizing offset distribution in towed seismic streamer data
CN103984019B (en) Local relevant weighted earthquake beam synthesis method
CN101320095B (en) Seismic Signal Synthesis Method for Single Source Directional Illumination
Hegna et al. Making the transition from discrete shot records to continuous seismic records and source wavefields, and its potential impact on survey efficiency and environmental footprint
CN1220884C (en) Phased array seismic prospecting method
Alfaro et al. Reducing exploration risk
GB2368911A (en) Computing a stacked seismic line by interpolation between known stacks
Nasyrov et al. Multiple migration of VSP data for velocity analysis
Poole et al. Advanced imaging solutions for tailored multi-source and multi-vessel surveys
Pecholcs A history of seismic acquisition design in Saudi Arabia-lessons learned
CN117761758A (en) Wire harness seismic exploration method and system
JSC et al. Advantages and Disadvantages of Surface and Downhole Seismic Illustrated by Processing Results of 3D VSP and 3D+ VSP
Landa et al. Multifocusing stack-a new time imaging method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20100818

Termination date: 20130808