CN110579799A - Seismic acquisition observation method and system with equal travel time intervals - Google Patents
Seismic acquisition observation method and system with equal travel time intervals Download PDFInfo
- Publication number
- CN110579799A CN110579799A CN201910831794.9A CN201910831794A CN110579799A CN 110579799 A CN110579799 A CN 110579799A CN 201910831794 A CN201910831794 A CN 201910831794A CN 110579799 A CN110579799 A CN 110579799A
- Authority
- CN
- China
- Prior art keywords
- travel time
- coordinate
- seismic
- wave detection
- earth surface
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/18—Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
- G01V1/181—Geophones
- G01V1/184—Multi-component geophones
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
Abstract
The invention relates to a method and a system for seismic acquisition and observation at equal travel time intervals, which are characterized by comprising the following steps: 1) determining a target stratum in the three-dimensional grid-shaped seismic wave velocity model; 2) selecting a certain earth surface seismic source point on the earth surface, and determining x coordinate data and y coordinate data of an earth surface seismic source point corresponding to the earth surface seismic source point; 3) calculating the corresponding relation between the seismic wave travel time reflected by the target stratum and the emergence position of the seismic wave propagation ray on the earth surface by adopting a ray tracing method; 4) determining an x coordinate sequence or a y coordinate sequence of the surface wave detection point by adopting an equal travel time interval method; 5) correspondingly determining a y coordinate sequence or an x coordinate sequence of the earth surface wave detection point by adopting a conventional seismic acquisition method or an equal travel time interval method; 6) the x coordinate sequence and the y coordinate sequence of the earth surface wave detection points are arranged and combined to obtain the coordinate sequence of the earth surface wave detection points.
Description
Technical Field
The invention relates to a seismic acquisition observation method and system with equal travel time intervals, belonging to the technical field of petroleum seismic exploration.
Background
Seismic exploration is one of the main tools for exploration of petroleum and natural gas resources, and the main processes of seismic exploration comprise data acquisition, seismic data processing and data interpretation. The data acquisition comprises three operations of determining a seismic observation system on a land seismic exploration data acquisition field, arranging a seismic source and a geophone in the field according to a determination scheme, and exciting and receiving seismic waves, wherein the determination of the seismic observation system is firstly carried out indoors to determine the optimal arrangement positions of a seismic source point and an earth surface geophone point, and then the seismic source and the geophone are arranged in the field according to the determination scheme. The seismic source point of on-road seismic exploration generally adopts an explosive source, and a plurality of detectors are arranged at equal intervals along a seismic survey line to receive seismic wave signals, and the number of the detectors in modern seismic exploration can reach more than 10000. The seismic source point generates seismic waves after explosion, the seismic waves are reflected back by a rock stratum interface and are received by a detector and transmitted to an instrument vehicle, and the instrument vehicle records signals transmitted by the detector to obtain a seismic record for researching the underground oil and gas burial condition. The seismic data processing is to input the seismic records acquired in the data acquisition stage into a special electronic computer, perform processing operation by adopting a series of programs with different functions according to different requirements, classify and arrange the seismic records, highlight effective seismic records, remove ineffective and interfered seismic records, and finally stack and shift the seismic records after various processing to finally obtain a two-dimensional or three-dimensional seismic data volume file. The data interpretation is a process of changing the processed seismic record into a geological result, and comprises the steps of applying fluctuation theory and geological knowledge, integrating various data such as geology, well drilling, well logging and the like, making structure interpretation, stratum interpretation, lithology and hydrocarbon detection interpretation and comprehensive interpretation, drawing a related result drawing, making oil-gas-containing evaluation on an exploration area, proposing a well drilling position and the like.
The determination of the seismic survey system is the first link of the overall seismic survey and is the basis of the subsequent seismic data processing and data interpretation steps. And (4) determining the earthquake observation system, namely reasonably planning the optimal placement positions of the earthquake source point and the earth surface detection point so as to obtain the earthquake data with the highest quality. In conventional seismic survey system determination, as shown in fig. 1, a seismic source point and a surface detector point are generally arranged in an equally spaced sampling manner. After the space sampling range and the density of the seismic source points and the earth surface wave detection points are given, all the seismic source points and the earth surface wave detection points are sequentially arranged on the ground in an equally-spaced mode. The equal-interval sampling mode can be suitable for solving most earthquake geological conditions and has the advantage of easy field construction, but the optimal sampling mode under all conditions is not always available, and particularly when the sampling density of a surface wave detection point is low, the quality of the earthquake record cannot be effectively improved under the condition of not increasing the earthquake acquisition cost.
Disclosure of Invention
In view of the above problems, an object of the present invention is to provide a seismic acquisition observation method and system with equal travel time intervals, which can effectively improve the quality of seismic records without increasing the seismic acquisition cost.
In order to achieve the purpose, the invention adopts the following technical scheme: a seismic acquisition observation method with equal travel time intervals is characterized by comprising the following steps: 1) determining a target stratum in the three-dimensional grid-shaped seismic wave velocity model, and determining a function of the depth of the target stratum along with the change of a horizontal coordinate; 2) selecting a certain earth surface seismic source point on the earth surface, and determining x coordinate data and y coordinate data of an earth surface seismic source point corresponding to the earth surface seismic source point; 3) calculating the corresponding relation between the seismic wave travel time reflected at the target stratum and the emergence position of the seismic wave propagation ray on the earth surface according to the function of the depth of the target stratum along with the change of the horizontal coordinate by adopting a ray tracing method; 4) determining an x coordinate sequence or a y coordinate sequence of the earth surface wave detection point by adopting an equal travel time interval method according to the obtained corresponding relation and the x coordinate data or the y coordinate data; 5) correspondingly determining a y coordinate sequence or an x coordinate sequence of the earth surface wave detection point by adopting a conventional seismic acquisition method or an equal travel time interval method; 6) arranging and combining the x coordinate sequence and the y coordinate sequence of the earth surface wave detection points to obtain a coordinate sequence of the earth surface wave detection points; 7) and repeating the steps 2) to 6) until coordinate sequences of all the surface wave detection points corresponding to all the surface seismic source points are obtained, and completing seismic acquisition and observation of the target stratum at equal travel time intervals.
Further, the x coordinate data includes an x coordinate range of the surface seismic source point corresponding to the surface detector point and the number of the surface detector points in the x direction, and the y coordinate data includes a y coordinate range of the surface seismic source point corresponding to the surface detector point and the number of the surface detector points in the y direction.
further, the specific process of step 3) is as follows: adopting ray tracing method, according to the function z of depth of said target stratum along with horizontal coordinate changet(x, y), taking the earth surface seismic source point as the starting point of the seismic wave propagation ray, and calculating the seismic wave travel time t reflected at the target stratum and the emergence position x of the seismic wave propagation ray on the earth surface under different seismic wave propagation pathsdThe corresponding relation between them.
Further, the specific process of the step 4) is as follows: 4.1) the x coordinate range (x) of the earth surface wave detection point according to the obtained corresponding relationmin，xmax) Or y coordinate range (y)min，ymax) Determining the travel time range (t) corresponding to the x coordinate range of the surface wave detection pointx，min，tx，max) Or the travel time range (t) corresponding to the y coordinate range of the surface wave detection pointy，min，ty，max) (ii) a 4.2) number N of surface wave detection points in x directionxOr the number N of surface detection points in the y-directionyThe calculated travel time range is divided into travel time intervals delta txOr Δ tyEqually dividing to determine the travel time sequence corresponding to the surface wave detection point in the x directionor travel time sequence corresponding to the earth surface wave detection point in the y directionWherein, tn＝tx，min+(n-1)Δtxand is and Is the Nthxthe travel time corresponding to each earth surface wave detection point; or, tn＝ty，min+(n-1)ΔtyAnd is and Is the NthyThe travel time corresponding to each earth surface wave detection point; 4.3) according to the obtained corresponding relation and the travel time sequence corresponding to the surface wave detection point in the x directionOr travel time sequence corresponding to the earth surface wave detection point in the y directionCalculating the x coordinate sequence { x) of the earth surface wave detection point1，x2，x3，...，xNxSequence of } or y coordinates
further, the specific process of the step 1) is as follows: giving a three-dimensional grid-shaped seismic wave velocity model V (x, y, z), and determining a certain target stratum in the three-dimensional grid-shaped seismic wave velocity model V (x, y, z); determining a function T of the depth of the target formation as a function of horizontal coordinatesz(x，y)。
further, the conventional seismic acquisition method is an equal interval sampling method.
An equal travel time interval seismic acquisition observation system, comprising: the target stratum determining module is used for determining a target stratum in the three-dimensional grid-shaped seismic wave velocity model and determining a function of the depth of the target stratum along with the change of the horizontal coordinate; the coordinate data determination module is used for selecting a certain earth surface seismic source point on the earth surface and determining x coordinate data and y coordinate data of an earth surface detection point corresponding to the earth surface seismic source point; the corresponding relation calculation module is used for calculating the corresponding relation between the seismic wave travel time reflected at the target stratum and the emergence position of the seismic wave propagation ray on the earth surface according to the function of the depth of the target stratum changing along with the horizontal coordinate by adopting a ray tracing method; the first coordinate sequence determination module is used for determining an x coordinate sequence or a y coordinate sequence of the earth surface wave detection point according to the obtained corresponding relation and the x coordinate data or the y coordinate data by adopting an equal travel time interval method; the second coordinate sequence determination module is used for correspondingly determining a y coordinate sequence or an x coordinate sequence of the earth surface wave detection point by adopting a conventional seismic acquisition method or an equal travel time interval method; and the arrangement and combination module is used for carrying out arrangement and combination on the x coordinate sequence and the y coordinate sequence of the surface wave detection points to obtain the coordinate sequence of the surface wave detection points.
Further, the first coordinate sequence determination module includes: a travel time range determining unit for determining the x coordinate range (x) of the earth surface probe point according to the obtained corresponding relationmin，xmax) Or y coordinate range (y)min，ymax) Determining the travel time range (t) corresponding to the x coordinate range of the surface wave detection pointx，min，tx，max) Or the travel time range (t) corresponding to the y coordinate range of the surface wave detection pointy，min，ty，max) (ii) a A travel time sequence determination unit for determining the number N of surface wave detection points in the x directionxOr the number N of surface detection points in the y-directionyThe calculated travel time range is divided into travel time intervals delta txor Δ tyEqually dividing to determine the travel time sequence corresponding to the surface wave detection point in the x directionOr travel time sequence corresponding to the earth surface wave detection point in the y directionWherein, tn＝tx，min+(n-1)Δtxand is and Is the NthxThe travel time corresponding to each earth surface wave detection point; or tn＝ty，min+(n-1)ΔtyAnd is and Is the NthyThe travel time corresponding to each earth surface wave detection point; a coordinate sequence calculation unit for calculating the travel time sequence corresponding to the surface wave detection point in the x direction according to the obtained corresponding relationor travel time sequence corresponding to the earth surface wave detection point in the y directionCalculating the x coordinate sequence { x) of the earth surface wave detection point1，x2，x3，...，xNxSequence of } or y coordinates
Due to the adoption of the technical scheme, the invention has the following advantages: compared with the traditional seismic acquisition mode with equal geophone point intervals, the method adopts the equal travel time interval method, can effectively strengthen the pertinence of the existing seismic acquisition observation to the target stratum under the condition of not increasing the seismic acquisition cost, effectively improves the quality of seismic data, improves the seismic migration imaging effect, has important significance for improving the success rate of oil and gas exploration, and can be widely applied to the technical field of oil and seismic exploration.
Drawings
FIG. 1 is a schematic diagram of a conventional equal-interval sampling method in the prior art;
FIG. 2 is a schematic illustration of equal travel time interval sampling in the method of the present invention;
FIG. 3 is a schematic diagram of a seismic acquisition observation system with equal travel time interval in a complex medium, wherein the abscissa is position (km) and the ordinate is depth (km), obtained by the method of the present invention.
Detailed Description
the present invention is described in detail below with reference to the attached drawings. It is to be understood, however, that the drawings are provided solely for the purposes of promoting an understanding of the invention and that they are not to be construed as limiting the invention.
The invention provides a seismic acquisition observation method with equal travel time intervals, which comprises the following steps:
1) given a three-dimensional grid-shaped seismic wave velocity model V (x, y, z), determining a certain target stratum in the three-dimensional grid-shaped seismic wave velocity model V (x, y, z), and determining a function T of the depth of the target stratum along with the change of horizontal coordinatesz(x，y)。
2) selecting a certain earth surface seismic source point on the earth surface, and determining the x coordinate range (x) of the earth surface seismic source point corresponding to the earth surface detection pointmin，xmax) Y coordinate range (y)min，ymax) The number N of surface wave detection points in the x directionxAnd the number N of surface detection points in the y directionyThe x coordinate range, the y coordinate range and the number of surface wave detection points can be obtained by a traditional seismic acquisition and evaluation method, and then are directly used as input parameters of the method, and the specific analysis process is not repeated here. As shown in FIG. 2, the minimum value of the x-coordinate range is 0, and the number N of surface detection points in the x-directionx＝4。
3) Adopting ray tracing method, according to the function z of depth of said target stratum along with horizontal coordinate changet(x, y), taking the earth surface seismic source point as the starting point of the seismic wave propagation ray, and calculating the seismic wave travel time t reflected at the target stratum and the emergence position x of the seismic wave propagation ray on the earth surface under different seismic wave propagation pathsdThe ray tracing method is a method disclosed in the prior art, and the specific process is not described herein. FIG. 3 illustrates an example of ray tracing under a complex velocity model, wherein rays of seismic waves are seismically initiated from the earth's surfacethe emission is carried out on the target stratum from the source point along the curved propagation path to reach the surface emergence position xd。
4) Determining an x coordinate sequence and a y coordinate sequence of a surface wave detection point by adopting an equal travel time interval method, which specifically comprises the following steps:
4.1) according to the obtained seismic wave travel time t and the emergence position x of the seismic wave propagation ray on the ground surfacedcorresponding relation between them, and x coordinate range (x) of surface wave detection pointmin，xmax) And y coordinate range (y)min，ymax) Respectively calculating the x coordinate range (x) of the earth surface wave detection pointmin，xmax) Corresponding travel time range (t)x，min，tx，max) And the y coordinate range (y) of the earth's surface probe pointmin，ymax) Corresponding travel time range (t)y，min，ty，max)。
4.2) number N of surface wave detection points in x directionxAnd the number N of surface detection points in the y directionyrange of travel time (t) to be calculatedx，min，tx，max) And (t)y，min，ty，max) According to the corresponding travel time interval Δ txAnd Δ tyEqually dividing to determine the travel time sequence corresponding to the surface wave detection point in the x directionTravel time sequence corresponding to surface wave detection point in y directionWherein, tn＝tx，min+(n-1)ΔtxAnd is and Is the NthxThe travel time corresponding to each earth surface wave detection point; or tn＝ty，min+(n-1)ΔtyAnd is and Is the NthyAnd the travel time corresponding to each earth surface wave detection point.
4.3) according to the seismic wave travel time t and the emergence position x of the seismic wave propagation ray on the earth surfacedCorresponding relation between them, travel time series corresponding to surface wave detection point in x directiontravel time sequence corresponding to surface wave detection point in y directionCorrespondingly calculating an x coordinate sequence { x) of the earth surface detection point1，x2，x3，...，xNxAnd y coordinate sequenceAs shown in fig. 3, the isoreflection angle rays from the underground target point reach the coordinate position of the earth surface, that is, the coordinates of the earth surface detection point.
5) Arranging and combining the x coordinate sequence and the y coordinate sequence of the surface wave detection points to obtain the coordinate sequence of the surface wave detection points:
6) And repeating the steps 2) to 5) until coordinate sequences of all the surface wave detection points corresponding to all the surface seismic source points are obtained, and completing seismic acquisition and observation at equal travel time intervals for a certain target stratum.
In the above steps 4) to 6), the coordinate sequence of the surface wave detection points at equal travel time intervals may be observed only in one direction (for example, the x direction), and the coordinate sequence of the surface wave detection points in the other direction may be obtained by using a conventional seismic acquisition method, i.e., an equal-interval sampling method.
based on the method for acquiring and observing the earthquake with the equal travel time interval, the invention also provides an earthquake acquisition and observation system with the equal travel time interval, which comprises the following steps:
the target stratum determining module is used for determining a target stratum in the three-dimensional grid-shaped seismic wave velocity model and determining a function of the depth of the target stratum along with the change of the horizontal coordinate; the coordinate data determination module is used for selecting a certain earth surface seismic source point on the earth surface and determining x coordinate data and y coordinate data of an earth surface detection point corresponding to the earth surface seismic source point; the corresponding relation calculation module is used for calculating the corresponding relation between the seismic wave travel time reflected at the target stratum and the emergence position of the seismic wave propagation ray on the earth surface according to the function of the depth of the target stratum changing along with the horizontal coordinate by adopting a ray tracing method; the first coordinate sequence determination module is used for determining an x coordinate sequence or a y coordinate sequence of the earth surface wave detection point according to the obtained corresponding relation and the x coordinate data or the y coordinate data by adopting an equal travel time interval method; the second coordinate sequence determination module is used for correspondingly determining a y coordinate sequence or an x coordinate sequence of the earth surface wave detection point by adopting a conventional seismic acquisition method or an equal travel time interval method; and the arrangement and combination module is used for carrying out arrangement and combination on the x coordinate sequence and the y coordinate sequence of the surface wave detection points to obtain the coordinate sequence of the surface wave detection points.
in a preferred embodiment, the first coordinate sequence determination module includes:
A travel time range determining unit for determining the x coordinate range (x) of the earth surface probe point according to the obtained corresponding relationmin，xmax) Or y coordinate range (y)min，ymax) Determining the travel time range (t) corresponding to the x coordinate range of the surface wave detection pointx，min，tx，max) Or the travel time range (t) corresponding to the y coordinate range of the surface wave detection pointy，min，ty，max)；
A travel time sequence determination unit for determining the number N of surface wave detection points in the x directionxor the number N of surface detection points in the y-directionyThe calculated travel time range is according to the travelTime interval Δ txOr Δ tyEqually dividing to determine the travel time sequence corresponding to the surface wave detection point in the x directionOr travel time sequence corresponding to the earth surface wave detection point in the y directionwherein, tn＝tx，min+(n-1)ΔtxAnd is and Is the NthxThe travel time corresponding to each earth surface wave detection point; or tn＝ty，min+(n-1)ΔtyAnd is and is the NthyThe travel time corresponding to each earth surface wave detection point;
A coordinate sequence calculation unit for calculating the travel time sequence corresponding to the surface wave detection point in the x direction according to the obtained corresponding relationOr travel time sequence corresponding to the earth surface wave detection point in the y directionCalculating the x coordinate sequence { x) of the earth surface wave detection point1，x2，x3，...，xNxSequence of } or y coordinates
The above embodiments are only used for illustrating the present invention, and the structure, connection mode, manufacturing process, etc. of the components may be changed, and all equivalent changes and modifications performed on the basis of the technical solution of the present invention should not be excluded from the protection scope of the present invention.
Claims (8)
1. A seismic acquisition observation method with equal travel time intervals is characterized by comprising the following steps:
1) Determining a target stratum in the three-dimensional grid-shaped seismic wave velocity model, and determining a function of the depth of the target stratum along with the change of a horizontal coordinate;
2) Selecting a certain earth surface seismic source point on the earth surface, and determining x coordinate data and y coordinate data of an earth surface seismic source point corresponding to the earth surface seismic source point;
3) Calculating the corresponding relation between the seismic wave travel time reflected at the target stratum and the emergence position of the seismic wave propagation ray on the earth surface according to the function of the depth of the target stratum along with the change of the horizontal coordinate by adopting a ray tracing method;
4) Determining an x coordinate sequence or a y coordinate sequence of the earth surface wave detection point by adopting an equal travel time interval method according to the obtained corresponding relation and the x coordinate data or the y coordinate data;
5) correspondingly determining a y coordinate sequence or an x coordinate sequence of the earth surface wave detection point by adopting a conventional seismic acquisition method or an equal travel time interval method;
6) arranging and combining the x coordinate sequence and the y coordinate sequence of the earth surface wave detection points to obtain a coordinate sequence of the earth surface wave detection points;
7) And repeating the steps 2) to 6) until coordinate sequences of all the surface wave detection points corresponding to all the surface seismic source points are obtained, and completing seismic acquisition and observation of the target stratum at equal travel time intervals.
2. The method of claim, wherein the x-coordinate data comprises an x-coordinate range of the seismic source point corresponding to the surface detector points and a number of surface detector points in the x-direction, and the y-coordinate data comprises a y-coordinate range of the seismic source point corresponding to the surface detector points and a number of surface detector points in the y-direction.
3. The method for seismic acquisition and observation with equal travel time intervals as claimed in claim 1, wherein the specific process of the step 3) is as follows:
Adopting ray tracing method, according to the function z of depth of said target stratum along with horizontal coordinate changet(x, y), taking the earth surface seismic source point as the starting point of the seismic wave propagation ray, and calculating the seismic wave travel time t reflected at the target stratum and the emergence position x of the seismic wave propagation ray on the earth surface under different seismic wave propagation pathsdThe corresponding relation between them.
4. The method for seismic acquisition and observation with equal travel time intervals as claimed in claim 2, wherein the specific process of the step 4) is as follows:
4.1) the x coordinate range (x) of the earth surface wave detection point according to the obtained corresponding relationmin，xmax) Or y coordinate range (y)min，ymax) Determining the travel time range (t) corresponding to the x coordinate range of the surface wave detection pointx，min，tx，max) Or the travel time range (t) corresponding to the y coordinate range of the surface wave detection pointy，min，ty，max)；
4.2) number N of surface wave detection points in x directionxOr the number N of surface detection points in the y-directionyThe calculated travel time range is divided into travel time intervals delta txor Δ tyequally dividing to determine the travel time sequence corresponding to the surface wave detection point in the x directionOr travel time sequence corresponding to the earth surface wave detection point in the y directionWherein, tn＝tx，min+(n-1)ΔtxAnd is and Is the NthxThe travel time corresponding to each earth surface wave detection point; or, tn＝ty，min+(n-1)ΔtyAnd is and is the Nthythe travel time corresponding to each earth surface wave detection point;
4.3) according to the obtained corresponding relation and the travel time sequence corresponding to the surface wave detection point in the x directionOr travel time sequence corresponding to the earth surface wave detection point in the y directionCalculating the x coordinate sequence { x) of the earth surface wave detection point1，x2，x3，...，xNxSequence of } or y coordinates
5. The method for seismic acquisition and observation with equal travel time intervals as claimed in claim 1, wherein the specific process of the step 1) is as follows:
Giving a three-dimensional grid-shaped seismic wave velocity model V (x, y, z), and determining a certain target stratum in the three-dimensional grid-shaped seismic wave velocity model V (x, y, z);
Determining a function T of the depth of the target formation as a function of horizontal coordinatesz(x，y)。
6. The method of claim 1, wherein the conventional seismic acquisition method is an equal interval sampling method.
7. An equal travel time interval seismic acquisition observation system, comprising:
The target stratum determining module is used for determining a target stratum in the three-dimensional grid-shaped seismic wave velocity model and determining a function of the depth of the target stratum along with the change of the horizontal coordinate;
the coordinate data determination module is used for selecting a certain earth surface seismic source point on the earth surface and determining x coordinate data and y coordinate data of an earth surface detection point corresponding to the earth surface seismic source point;
The corresponding relation calculation module is used for calculating the corresponding relation between the seismic wave travel time reflected at the target stratum and the emergence position of the seismic wave propagation ray on the earth surface according to the function of the depth of the target stratum changing along with the horizontal coordinate by adopting a ray tracing method;
the first coordinate sequence determination module is used for determining an x coordinate sequence or a y coordinate sequence of the earth surface wave detection point according to the obtained corresponding relation and the x coordinate data or the y coordinate data by adopting an equal travel time interval method;
The second coordinate sequence determination module is used for correspondingly determining a y coordinate sequence or an x coordinate sequence of the earth surface wave detection point by adopting a conventional seismic acquisition method or an equal travel time interval method;
And the arrangement and combination module is used for carrying out arrangement and combination on the x coordinate sequence and the y coordinate sequence of the surface wave detection points to obtain the coordinate sequence of the surface wave detection points.
8. The seismic acquisition observation system for equal travel time intervals of claim 7, wherein the first coordinate sequence determination module comprises:
A travel time range determining unit for determining the x coordinate range (x) of the earth surface probe point according to the obtained corresponding relationmin，xmax) Or a y coordinate rangeEnclose (y)min，ymax) Determining the travel time range (t) corresponding to the x coordinate range of the surface wave detection pointx，min，tx，max) Or the travel time range (t) corresponding to the y coordinate range of the surface wave detection pointy，min，ty，max)；
A travel time sequence determination unit for determining the number N of surface wave detection points in the x directionxOr the number N of surface detection points in the y-directionyThe calculated travel time range is divided into travel time intervals delta txor Δ tyEqually dividing to determine the travel time sequence corresponding to the surface wave detection point in the x directionOr travel time sequence corresponding to the earth surface wave detection point in the y directionWherein, tn＝tx，min+(n-1)ΔtxAnd is and is the NthxThe travel time corresponding to each earth surface wave detection point; or tn＝ty，min+(n-1)ΔtyAnd is and Is the NthyThe travel time corresponding to each earth surface wave detection point;
A coordinate sequence calculation unit for calculating the travel time sequence corresponding to the surface wave detection point in the x direction according to the obtained corresponding relationOr travel time sequence corresponding to the earth surface wave detection point in the y directioncalculating the x coordinate sequence { x) of the earth surface wave detection point1，x2，x3，...，xNxsequence of } or y coordinates
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910831794.9A CN110579799B (en) | 2019-09-04 | 2019-09-04 | Seismic acquisition observation method and system with equal travel time intervals |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910831794.9A CN110579799B (en) | 2019-09-04 | 2019-09-04 | Seismic acquisition observation method and system with equal travel time intervals |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110579799A true CN110579799A (en) | 2019-12-17 |
| CN110579799B CN110579799B (en) | 2020-10-27 |
Family
ID=68812500
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910831794.9A Active CN110579799B (en) | 2019-09-04 | 2019-09-04 | Seismic acquisition observation method and system with equal travel time intervals |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110579799B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111856581A (en) * | 2020-07-27 | 2020-10-30 | 广州海洋地质调查局 | OBS clock drift correction method and processing terminal |
| CN112505749A (en) * | 2020-10-19 | 2021-03-16 | 中国地质调查局南京地质调查中心(华东地质科技创新中心) | Micro-motion data acquisition method based on multiple coverage of linear array |
| CN112505749B (en) * | 2020-10-19 | 2024-04-26 | 中国地质调查局南京地质调查中心(华东地质科技创新中心) | Micro-motion data acquisition method based on linear array multiple coverage |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030208321A1 (en) * | 2001-07-31 | 2003-11-06 | Ruben Martinez | Relative true amplitude migration |
| CN101173989A (en) * | 2006-11-03 | 2008-05-07 | 中国石油集团东方地球物理勘探有限责任公司 | Quantitative analysis method for earthquake-capturing three-dimensional observation system |
| CN102636809A (en) * | 2012-03-27 | 2012-08-15 | 中国科学院地质与地球物理研究所 | Method for generating spreading angle domain common image point gathers |
| CN104090297A (en) * | 2014-06-14 | 2014-10-08 | 吉林大学 | Reverse illumination method for optimizing earthquake collection observing system |
| WO2014164354A1 (en) * | 2013-04-02 | 2014-10-09 | Halliburton Energy Services, Inc. | Anisotropy analysis using direct and reflected arrivals in seismic survey data |
| CN104267434A (en) * | 2014-09-12 | 2015-01-07 | 中国科学院地质与地球物理研究所 | Three-dimensional multi-component earthquake observation system geophone offset distribution obtaining method and device |
| CN104297790A (en) * | 2014-09-19 | 2015-01-21 | 中国海洋石油总公司 | Quantitative evaluation method for offset distribution in fluctuating submarine earthquake observation system |
| CN104360388A (en) * | 2014-10-31 | 2015-02-18 | 中国石油化工股份有限公司 | Method for evaluating three-dimensional seismic observation systems |
| CN104570069A (en) * | 2013-10-12 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | Seismic imaging method and device for direct imaging from undulating surface |
| CN105388513A (en) * | 2015-10-15 | 2016-03-09 | 中国石油天然气集团公司 | Method and device for establishing earthquake acquisition and observation system |
| CN105974463A (en) * | 2016-06-21 | 2016-09-28 | 吉林大学 | Seismic exploration method based on source inspection bi-orienting |
| CN106324702A (en) * | 2015-06-16 | 2017-01-11 | 核工业北京地质研究院 | Quantitative evaluation method of seismic interferometry imaging observation system design |
| CN106932821A (en) * | 2015-12-31 | 2017-07-07 | 上海青凤致远地球物理地质勘探科技有限公司 | A kind of direction ray tracer technique in seismic tomography inverting |
-
2019
- 2019-09-04 CN CN201910831794.9A patent/CN110579799B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030208321A1 (en) * | 2001-07-31 | 2003-11-06 | Ruben Martinez | Relative true amplitude migration |
| CN101173989A (en) * | 2006-11-03 | 2008-05-07 | 中国石油集团东方地球物理勘探有限责任公司 | Quantitative analysis method for earthquake-capturing three-dimensional observation system |
| CN102636809A (en) * | 2012-03-27 | 2012-08-15 | 中国科学院地质与地球物理研究所 | Method for generating spreading angle domain common image point gathers |
| WO2014164354A1 (en) * | 2013-04-02 | 2014-10-09 | Halliburton Energy Services, Inc. | Anisotropy analysis using direct and reflected arrivals in seismic survey data |
| CN104570069A (en) * | 2013-10-12 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | Seismic imaging method and device for direct imaging from undulating surface |
| CN104090297A (en) * | 2014-06-14 | 2014-10-08 | 吉林大学 | Reverse illumination method for optimizing earthquake collection observing system |
| CN104267434A (en) * | 2014-09-12 | 2015-01-07 | 中国科学院地质与地球物理研究所 | Three-dimensional multi-component earthquake observation system geophone offset distribution obtaining method and device |
| CN104297790A (en) * | 2014-09-19 | 2015-01-21 | 中国海洋石油总公司 | Quantitative evaluation method for offset distribution in fluctuating submarine earthquake observation system |
| CN104360388A (en) * | 2014-10-31 | 2015-02-18 | 中国石油化工股份有限公司 | Method for evaluating three-dimensional seismic observation systems |
| CN106324702A (en) * | 2015-06-16 | 2017-01-11 | 核工业北京地质研究院 | Quantitative evaluation method of seismic interferometry imaging observation system design |
| CN105388513A (en) * | 2015-10-15 | 2016-03-09 | 中国石油天然气集团公司 | Method and device for establishing earthquake acquisition and observation system |
| CN106932821A (en) * | 2015-12-31 | 2017-07-07 | 上海青凤致远地球物理地质勘探科技有限公司 | A kind of direction ray tracer technique in seismic tomography inverting |
| CN105974463A (en) * | 2016-06-21 | 2016-09-28 | 吉林大学 | Seismic exploration method based on source inspection bi-orienting |
Non-Patent Citations (2)
| Title |
|---|
| 刘晨: ""基于地质目标的观测系统优化设计方法研究"", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 蒋先艺: ""基于二维和三维复杂结构模型正演的地震数据采集设计方法研究"", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111856581A (en) * | 2020-07-27 | 2020-10-30 | 广州海洋地质调查局 | OBS clock drift correction method and processing terminal |
| CN111856581B (en) * | 2020-07-27 | 2022-02-22 | 广州海洋地质调查局 | OBS clock drift correction method and processing terminal |
| CN112505749A (en) * | 2020-10-19 | 2021-03-16 | 中国地质调查局南京地质调查中心(华东地质科技创新中心) | Micro-motion data acquisition method based on multiple coverage of linear array |
| CN112505749B (en) * | 2020-10-19 | 2024-04-26 | 中国地质调查局南京地质调查中心(华东地质科技创新中心) | Micro-motion data acquisition method based on linear array multiple coverage |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110579799B (en) | 2020-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10577926B2 (en) | Detecting sub-terranean structures | |
| Verdon et al. | Microseismic monitoring using a fiber-optic distributed acoustic sensor array | |
| CN102395902B (en) | Seismic imaging systems and methods employing a fast target-oriented illumination calculation | |
| US7952960B2 (en) | Seismic imaging with natural Green's functions derived from VSP data | |
| CN102282481B (en) | Data acquisition and prestack migration based on seismic visibility analysis | |
| US9869783B2 (en) | Structure tensor constrained tomographic velocity analysis | |
| AU742887B2 (en) | Hydrocarbon edge detection using seismic amplitude | |
| Liu et al. | A new 3D observation system designed for a seismic ahead prospecting method in tunneling | |
| CN110529087B (en) | Method and device for evaluating hydraulic fracturing effect of stratum | |
| AU2015253627B2 (en) | Method for using semblance of corrected amplitudes due to source mechanisms for microseismic event detection and location | |
| US20190094397A1 (en) | Surface detection and location of microseismic events and earthquakes without the use of a velocity model | |
| Zhang et al. | Microseismic hydraulic fracture imaging in the Marcellus Shale using head waves | |
| Zhang et al. | Simultaneous inversion for microseismic event location and velocity model in Vaca Muerta Formation | |
| US9658354B2 (en) | Seismic imaging systems and methods employing correlation-based stacking | |
| US6249746B1 (en) | Automated seismic isochron analysis | |
| CN113805237B (en) | Method and system for offset land crossrange seismic using compressed sensing models | |
| Foulger et al. | Earthquakes and errors: Methods for industrial applications | |
| Xuan et al. | Probabilistic microearthquake location for reservoir monitoring | |
| CN110579799B (en) | Seismic acquisition observation method and system with equal travel time intervals | |
| US20120269035A1 (en) | Evaluating Prospects from P-Wave Seismic Data Using S-Wave Vertical Shear Profile Data | |
| CN110579798B (en) | Seismic acquisition observation method and system with equal reflection angle intervals | |
| CN108375794B (en) | VSP (vertical seismic profiling) slit-hole diffraction imaging technical method based on symmetrical observation | |
| EP3669210B1 (en) | Method for obtaining estimates of a model parameter so as to characterise the evolution of a subsurface volume over a time period using time-lapse seismic | |
| CN116736383A (en) | Updating method and device of seismic wave velocity model, electronic equipment and storage medium | |
| CN115793044A (en) | Seismic source positioning method, processor and device based on seismic wave dynamic parameters |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |