CN110068863A - A kind of determination method and device of common depth point CDP - Google Patents
A kind of determination method and device of common depth point CDP Download PDFInfo
- Publication number
- CN110068863A CN110068863A CN201910271485.0A CN201910271485A CN110068863A CN 110068863 A CN110068863 A CN 110068863A CN 201910271485 A CN201910271485 A CN 201910271485A CN 110068863 A CN110068863 A CN 110068863A
- Authority
- CN
- China
- Prior art keywords
- coordinate
- cmp
- cdp
- face element
- geophone station
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000007689 inspection Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- 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. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
This application discloses the determination method and devices of CDP a kind of, are related to seismic exploration technique field.This method is included in the coordinate system of work area, according to the corresponding relationship of the coordinate of shot point each in seismic target earthquakes work area, the first coordinate of each geophone station and shot point and geophone station, determines the first coordinate of the CMP of each seismic channel；The first coordinate corresponding to each CMP carries out straight line fitting, obtains target line；The first coordinate system is established by axis of abscissas of target line, determines the second coordinate of each geophone station and the second coordinate of each CMP；According to corresponding second coordinate of each geophone station and preset binning rule, multiple face elements are determined；The first coordinate of the CDP of the corresponding seismic channel of the CMP for including in the face element is determined according to corresponding second coordinate of the CMP for including in the face element for each face element；According to corresponding first coordinate of the CDP of each seismic channel, second coordinate of the CDP of each seismic channel in preset work area coordinate system is determined.The CDP of seismic channel can be determined using the application.
Description
Technical field
This application involves seismic exploration technique field, in particular to a kind of determination method and device of common depth point CDP.
Background technique
Currently, seismic data carries out tissue generally as unit of seismic channel, stored using segy file format.segy
Format is proposed by Exploration Geophysicists Association (Society of Exploration Geophysicists, SEG)
One of the format the most universal of one of standard reference tape data format and oil exploration industry seismic data.Standard segy text
Part generally comprises three parts.First part is used to store the description information of seismic data.Second part is for storing segy file
Description information.Part III is for storing seismic channel information.Wherein, seismic channel information includes trace header information and seismic channel number
According to trace header information includes the corresponding wire size of seismic channel, Taoist monastic name, sampling number, geodetic coordinates, common depth point (Common Depth
Point, CDP) etc. information.
It would therefore be highly desirable to a kind of method of determining CDP.
Summary of the invention
In order to solve problems in the prior art, the embodiment of the present application provides the determination method and device of CDP a kind of, can be with
Determine the CDP of seismic channel.The technical solution is as follows:
In a first aspect, providing the determination method of CDP a kind of, which comprises
Coordinate, each geophone station of each shot point in preset work area coordinate system in seismic target earthquakes work area is obtained to preset described
Work area coordinate system in the first coordinate and shot point and geophone station corresponding relationship；
According to the corresponding coordinate of each shot point, each corresponding first coordinate of geophone station and the shot point and inspection
The corresponding relationship of wave point determines the CMP of each seismic channel in the seismic target earthquakes work area in the preset work area coordinate system
First coordinate；
According to preset Algorithm of fitting a straight line, the first coordinate corresponding to each CMP carries out straight line fitting, and it is straight to obtain target
Line；
The first coordinate system is established by axis of abscissas of the target line, determines each geophone station in first coordinate
The second coordinate of the second coordinate and each CMP in first coordinate system in system；
According to each corresponding second coordinate of geophone station and preset binning rule, multiple face elements are determined；
For each face element, according to corresponding second coordinate of the CMP for including in the face element, determines in the face element and include
First coordinate of the CDP of the corresponding seismic channel of CMP in first coordinate system；
According to corresponding first coordinate of the CDP of each seismic channel, determine the CDP of each seismic channel described preset
The second coordinate in the coordinate system of work area.
Optionally, the Algorithm of fitting a straight line is least square method.
Optionally, described according to each corresponding second coordinate of geophone station and preset binning rule, determination is more
A face element, comprising:
Determine the average value of the distance between abscissa of corresponding second coordinate of adjacent geophone station；
Using the average value and the ratio of preset constant as face element step-length；
According to minimum value, maximum value and the face element step in the abscissa of corresponding second coordinate of each geophone station
It is long, determine multiple face elements.
Optionally, described that the face is determined according to corresponding second coordinate of the CMP for including in the face element for each face element
First coordinate of the CDP of the corresponding seismic channel of the CMP for including in member in first coordinate system, comprising:
For each face element, determines the average value of the abscissa of corresponding second coordinate of the CMP for including in the face element and indulge
The average value of coordinate, using the average value of the abscissa as the CDP for the corresponding seismic channel of CMP for including in the face element described
The abscissa of the first coordinate in first coordinate system, and using the average value of the ordinate as CMP pairs for including in the face element
The ordinate of corresponding first coordinate of the CDP for the seismic channel answered.
Optionally, the method also includes:
The ratio of abscissa and the face element step-length to corresponding second coordinate of the first CMP is rounded, and obtains the first integer；
The ratio of abscissa and the face element step-length to corresponding first coordinate of the first CDP is rounded, and obtains the second integer；
If first integer is identical as second integer, it is determined that the corresponding seismic channel of the first CMP is corresponding
CDP be the first CDP.
Second aspect, provides the determining device of CDP a kind of, and described device includes:
Module is obtained, for obtaining coordinate in preset work area coordinate system of each shot point in seismic target earthquakes work area, each inspection
The corresponding relationship of first coordinate and shot point and geophone station of the wave point in the preset work area coordinate system；
First determining module, for according to the corresponding coordinate of each shot point, corresponding first coordinate of each geophone station,
And the corresponding relationship of the shot point and geophone station, determine the CMP of each seismic channel in the seismic target earthquakes work area described default
Work area coordinate system in the first coordinate；
Fitting module, for according to preset Algorithm of fitting a straight line, it is quasi- that the first coordinate corresponding to each CMP carries out straight line
It closes, obtains target line；
Second determining module determines each inspection for establishing the first coordinate system by axis of abscissas of the target line
Wave point is in the second coordinate of the second coordinate and each CMP in first coordinate system in first coordinate system；
Third determining module, for regular according to each corresponding second coordinate of geophone station and preset binning,
Determine multiple face elements；
4th determining module, for being directed to each face element, according to corresponding second coordinate of the CMP for including in the face element, really
First coordinate of the CDP for the corresponding seismic channel of CMP for including in the fixed face element in first coordinate system；
5th determining module determines each seismic channel for corresponding first coordinate of CDP according to each seismic channel
Second coordinate of the CDP in the preset work area coordinate system.
Optionally, the Algorithm of fitting a straight line is least square method.
Optionally, the third determining module, is specifically used for:
Determine the average value of the distance between abscissa of corresponding second coordinate of adjacent geophone station；
Using the average value and the ratio of preset constant as face element step-length；
According to minimum value, maximum value and the face element step in the abscissa of corresponding second coordinate of each geophone station
It is long, determine multiple face elements.
Optionally, the 4th determining module, is specifically used for:
For each face element, determines the average value of the abscissa of corresponding second coordinate of the CMP for including in the face element and indulge
The average value of coordinate, using the average value of the abscissa as the CDP for the corresponding seismic channel of CMP for including in the face element described
The abscissa of the first coordinate in first coordinate system, and using the average value of the ordinate as CMP pairs for including in the face element
The ordinate of corresponding first coordinate of the CDP for the seismic channel answered.
Optionally, described device further include:
First floor module, the ratio for abscissa and the face element step-length to corresponding second coordinate of the first CMP
It is rounded, obtains the first integer；
Second floor module, the ratio for abscissa and the face element step-length to corresponding first coordinate of the first CDP
It is rounded, obtains the second integer；
6th determining module, if identical as second integer for first integer, it is determined that the first CMP
The corresponding CDP of corresponding seismic channel is the first CDP.
The application provides the determination method and device of CDP a kind of, corresponding according to the corresponding coordinate of each shot point, each geophone station
The corresponding relationship of first coordinate and shot point and geophone station determines the CMP of each seismic channel in seismic target earthquakes work area in preset work
The first coordinate in area's coordinate system.Then, according to preset Algorithm of fitting a straight line, the first coordinate corresponding to each CMP carries out straight
Line fitting, obtains target line, and establish the first coordinate system by axis of abscissas of target line, determines that each geophone station is sat first
The second coordinate of the second coordinate and each CMP in the first coordinate system in mark system.Later, it is sat according to each geophone station corresponding second
Mark and preset binning rule, determine multiple face elements, and be directed to each face element, corresponding according to the CMP for including in the face element
The second coordinate, determine first coordinate of the CDP of the corresponding seismic channel of the CMP for including in the face element in the first coordinate system.Most
Afterwards, according to corresponding first coordinate of the CDP of each seismic channel, of the CDP of each seismic channel in preset work area coordinate system is determined
Two coordinates.
Detailed description of the invention
In order to more clearly explain the technical solutions in the embodiments of the present application, make required in being described below to embodiment
Attached drawing is briefly described, it should be apparent that, the drawings in the following description are only some examples of the present application, for
For those of ordinary skill in the art, without creative efforts, it can also be obtained according to these attached drawings other
Attached drawing.
Fig. 1 is a kind of flow chart of the determination method of CDP provided by the embodiments of the present application；
Fig. 2 is the schematic diagram of straight line fitting provided by the embodiments of the present application；
Fig. 3 is a kind of structural schematic diagram of the determining device of CDP provided by the embodiments of the present application；
Fig. 4 is a kind of structural schematic diagram of the determining device of CDP provided by the embodiments of the present application.
Specific embodiment
To keep the purposes, technical schemes and advantages of the application clearer, below in conjunction with attached drawing to the application embodiment party
Formula is described in further detail.
This application provides the determination method of CDP a kind of, the executing subject of this method can be the networks such as server, host
Equipment, this method can be applied to seismic exploration technique field.The corresponding CDP of seismic channel can be determined based on this method.This Shen
Please embodiment by executing subject be server for be introduced, other situations are similar therewith.
It is detailed to a kind of determination method progress of CDP provided by the embodiments of the present application below in conjunction with specific embodiment
Explanation, as shown in Figure 1, the specific steps are as follows:
Step 101, coordinate, each geophone station of each shot point in preset work area coordinate system in seismic target earthquakes work area is obtained to exist
The corresponding relationship of the first coordinate and shot point and geophone station in preset work area coordinate system.
In an implementation, work area coordinate system can be previously stored in server, which can be by technical staff
Rule of thumb it is configured.When it needs to be determined that when the CDP of each seismic channel in a certain earthquake work area (namely seismic target earthquakes work area), clothes
Coordinate of each shot point in preset work area coordinate system, each geophone station be engaged in the available seismic target earthquakes work area of device preset
The corresponding relationship of the first coordinate and shot point and geophone station in the coordinate system of work area.
Step 102, according to the corresponding coordinate of each shot point, corresponding first coordinate of each geophone station and shot point and geophone station
Corresponding relationship, determine first coordinate of the CMP of each seismic channel in seismic target earthquakes work area in preset work area coordinate system.
In an implementation, server gets seat of each shot point in preset work area coordinate system in the seismic target earthquakes work area
Mark, each geophone station are after the corresponding relationship of the first coordinate and shot point and geophone station in preset work area coordinate system, Ke Yijin
One step is according to the corresponding relationship of the corresponding coordinate of each shot point, corresponding first coordinate of each geophone station and shot point and geophone station, really
Set the goal each seismic channel in earthquake work area common point (Common Middle Point, CMP) in preset work area coordinate system
In the first coordinate.Wherein, server is according to the corresponding coordinate of each shot point, corresponding first coordinate of each geophone station and shot point
With the corresponding relationship of geophone station, first of the CMP of each seismic channel in seismic target earthquakes work area in preset work area coordinate system is determined
The treatment process of coordinate with disclosed in existing scheme according to the corresponding coordinate of each shot point, corresponding first coordinate of each geophone station,
And the corresponding relationship of shot point and geophone station, determine that the treatment process of corresponding first coordinate of the CMP of each seismic channel is similar, herein
It repeats no more.
Step 103, according to preset Algorithm of fitting a straight line, the first coordinate corresponding to each CMP carries out straight line fitting, obtains
Target line.
In an implementation, Algorithm of fitting a straight line can be previously stored in server.The Algorithm of fitting a straight line can be by technology
Personnel are rule of thumb configured.Optionally, which can be least square method, can also be quasi- for other straight lines
Hop algorithm, the application are not construed as limiting.
Server obtains first coordinate of the CMP of each seismic channel in seismic target earthquakes work area in preset work area coordinate system
Afterwards, straight line fitting can be carried out, target line is obtained according to preset Algorithm of fitting a straight line, the first coordinate corresponding to each CMP.
Fig. 2 is the schematic diagram of straight line fitting provided by the embodiments of the present application, as shown in Fig. 2, in preset work area coordinate system
In, server carries out straight line fitting, it is straight to obtain target according to preset Algorithm of fitting a straight line, the first coordinate corresponding to each CMP
Line.
Step 104, the first coordinate system is established by axis of abscissas of target line, determines each geophone station in the first coordinate system
The second coordinate in the first coordinate system of the second coordinate and each CMP.
In an implementation, after server obtains target line, first further can be established by axis of abscissas of target line
Coordinate system.Wherein, the origin of first coordinate system can be the intersection point of target line and the axis of abscissas of work area coordinate system, can also
Think the intersection point of the axis of ordinates of target line and work area coordinate system, can also be other points in the target line, the application
Embodiment is not construed as limiting.After server establishes the first coordinate system, it may further determine that each geophone station in the first coordinate system
The second coordinate of second coordinate and each CMP in the first coordinate system.
Step 105, according to corresponding second coordinate of each geophone station and preset binning rule, multiple face elements are determined.
In an implementation, binning rule can be previously stored in server, which can be by technology
Personnel are rule of thumb configured.It, can be according to each after server obtains the second coordinate of each geophone station in the first coordinate system
Corresponding second coordinate of geophone station and preset binning rule, determine multiple face elements.
Optionally, server determines multiple according to corresponding second coordinate of each geophone station and preset binning rule
The treatment process of face element is as follows:
Step 1 determines the average value of the distance between abscissa of corresponding second coordinate of adjacent geophone station.
In an implementation, server can determine adjacent detection according to the abscissa of corresponding second coordinate of each geophone station
Point.Then, server can determine the distance between the abscissa of corresponding second coordinate of adjacent geophone station.Later, it services
Device can determine the average value of the distance between abscissa of corresponding second coordinate of all adjacent geophone stations.
Step 2, using average value and the ratio of preset constant as face element step-length.
In an implementation, preset constant can be previously stored in server.The preset constant can by technical staff according to
Experience is configured.Optionally, which can be 2, or other constants, the embodiment of the present application are not construed as limiting.
After server obtains the average value of the distance between abscissa of corresponding second coordinate of adjacent geophone station, can further it count
The ratio of the average value and preset constant is calculated, and using the ratio as face element step-length.For example, adjacent geophone station corresponding second
The average value of the distance between the abscissa of coordinate is 4, preset constant 2, then face element step-length is 2.
Step 3, according to minimum value, maximum value and the face element step-length in the abscissa of corresponding second coordinate of each geophone station,
Determine multiple face elements.
It in an implementation, can be further with the cross of corresponding second coordinate of each geophone station after server obtains face element step-length
Minimum value in coordinate is starting point, using the maximum value in the abscissa of corresponding second coordinate of each geophone station as terminal, according to face
First step-length divides face element, determines the range of each face element.For example, the minimum value in the abscissa of corresponding second coordinate of each geophone station
It is 2, with the maximum value in the abscissa of corresponding second coordinate of each geophone station for 10, face element step-length is 2, then the first binning range
It is 2~4, the second binning range is 4~6, and third binning range is 6~8, and fourth face member range is 8~10.
Step 106, it is determined in the face element for each face element according to corresponding second coordinate of the CMP for including in the face element
First coordinate of the CDP of the corresponding seismic channel of the CMP for including in the first coordinate system.
In an implementation, for each face element, server can according to corresponding second coordinate of the CMP for including in the face element,
Determine first coordinate of the CDP of the corresponding seismic channel of the CMP for including in the face element in the first coordinate system.
Optionally, for each face element, for server according to corresponding second coordinate of the CMP for including in the face element, determining should
The treatment process of first coordinate of the CDP of the corresponding seismic channel of the CMP for including in face element in the first coordinate system are as follows: for each
Face element determines the average value of the abscissa of corresponding second coordinate of the CMP for including in the face element and the average value of ordinate, will be horizontal
CDP first coordinate in first coordinate system of the average value of coordinate as the corresponding seismic channel of CMP for including in the face element
Abscissa, and using the average value of ordinate as corresponding first coordinate of the CDP for the corresponding seismic channel of CMP for including in the face element
Ordinate.
In an implementation, for each face element, server can determine corresponding second coordinate of the CMP for including in the face element
The average value of abscissa and the average value of ordinate.Then, server can be using the average value of abscissa as wrapping in the face element
The abscissa of first coordinate of the CDP of the corresponding seismic channel of the CMP contained in the first coordinate system, and the average value of ordinate is made
For the ordinate of corresponding first coordinate of CDP for the corresponding seismic channel of CMP for including in the face element.
Step 107, according to corresponding first coordinate of the CDP of each seismic channel, determine the CDP of each seismic channel in preset work area
The second coordinate in coordinate system.
It in an implementation, can be further after server obtains the first coordinate of the CDP of each seismic channel in the first coordinate system
Determine second coordinate (namely CDP of each seismic channel) of the CDP of each seismic channel in preset work area coordinate system.
Optionally, server can also be sat according to corresponding second coordinate of CMP of each seismic channel and each CDP corresponding first
The corresponding relationship for determining seismic channel and CDP is marked, treatment process is as follows:
The ratio of step 1, abscissa and face element step-length to corresponding second coordinate of the first CMP is rounded, and it is whole to obtain first
Number.
In an implementation, for a certain CMP (namely the first CMP) in the CMP of each seismic channel, server can be to first
The abscissa of corresponding second coordinate of CMP and the ratio of face element step-length are rounded, and obtain the first integer.For example, the first CMP is corresponding
The abscissa of second coordinate is 10, and face element step-length is 3, then the first integer is 3.
The ratio of step 2, abscissa and face element step-length to corresponding first coordinate of the first CDP is rounded, and it is whole to obtain second
Number.
In an implementation, for a certain CDP (namely the first CDP), server can be to corresponding first coordinate of the first CDP
Abscissa and the ratio of face element step-length are rounded, and obtain the second integer.For example, the abscissa of corresponding first coordinate of the first CDP is
11, face element step-length is 3, then the second integer is 3.
Step 3, if the first integer is identical as the second integer, it is determined that the corresponding CDP of the corresponding seismic channel of the first CMP
For the first CDP.
In an implementation, after server obtains the first integer and the second integer, it can be determined that the first integer is with the second integer
It is no identical.If the first integer is identical as the second integer, illustrate that the first CMP and the first CDP in the same face element, are serviced
Device can determine that the corresponding CDP of the corresponding seismic channel of the first CMP is the first CDP.
The application provides the determination method of CDP a kind of, is sat according to the corresponding coordinate of each shot point, each geophone station corresponding first
The corresponding relationship of mark and shot point and geophone station, determines the CMP of each seismic channel in seismic target earthquakes work area in preset work area coordinate
The first coordinate in system.Then, according to preset Algorithm of fitting a straight line, the first coordinate corresponding to each CMP carries out straight line fitting,
Target line is obtained, and establishes the first coordinate system by axis of abscissas of target line, determines each geophone station in the first coordinate system
The second coordinate in the first coordinate system of the second coordinate and each CMP.Later, according to corresponding second coordinate of each geophone station and in advance
If binning rule, multiple face elements are determined, and be directed to each face element, according to the CMP for including in the face element corresponding second
Coordinate determines first coordinate of the CDP of the corresponding seismic channel of the CMP for including in the face element in the first coordinate system.Finally, according to
Corresponding first coordinate of the CDP of each seismic channel determines second coordinate of the CDP of each seismic channel in preset work area coordinate system.
It based on the same technical idea, should as shown in figure 3, the embodiment of the present application also provides the determining device of CDP a kind of
Device includes:
Module 310 is obtained, for obtaining coordinate in preset work area coordinate system of each shot point in seismic target earthquakes work area, each
The corresponding relationship of first coordinate and shot point and geophone station of the geophone station in preset work area coordinate system；
First determining module 320, for according to the corresponding coordinate of each shot point, corresponding first coordinate of each geophone station and
The corresponding relationship of shot point and geophone station determines the CMP of each seismic channel in seismic target earthquakes work area in preset work area coordinate system
First coordinate；
Fitting module 330, for according to preset Algorithm of fitting a straight line, the first coordinate corresponding to each CMP to carry out straight line
Fitting, obtains target line；
Second determining module 340 determines that each geophone station exists for establishing the first coordinate system by axis of abscissas of target line
The second coordinate of the second coordinate and each CMP in the first coordinate system in first coordinate system；
Third determining module 350, for regular according to corresponding second coordinate of each geophone station and preset binning, really
Fixed multiple face elements；
4th determining module 360, for being directed to each face element, according to corresponding second coordinate of the CMP for including in the face element,
Determine first coordinate of the CDP of the corresponding seismic channel of the CMP for including in the face element in the first coordinate system；
5th determining module 370 determines the CDP of each seismic channel for corresponding first coordinate of CDP according to each seismic channel
The second coordinate in preset work area coordinate system.
Optionally, Algorithm of fitting a straight line is least square method.
Optionally, third determining module 350, is specifically used for:
Determine the average value of the distance between abscissa of corresponding second coordinate of adjacent geophone station；
Using average value and the ratio of preset constant as face element step-length；
According to minimum value, maximum value and the face element step-length in the abscissa of corresponding second coordinate of each geophone station, determination is more
A face element.
Optionally, the 4th determining module 360 is specifically used for:
For each face element, determines the average value of the abscissa of corresponding second coordinate of the CMP for including in the face element and indulge
The average value of coordinate, using the average value of abscissa as the CDP for the corresponding seismic channel of CMP for including in the face element in the first coordinate
The abscissa of the first coordinate in system, and using the average value of ordinate as the corresponding seismic channel of CMP for including in the face element
The ordinate of corresponding first coordinate of CDP.
Optionally, as shown in figure 4, the device further include:
First floor module 380, the ratio for abscissa and face element step-length to corresponding second coordinate of the first CMP take
It is whole, obtain the first integer；
Second floor module 390, the ratio for abscissa and face element step-length to corresponding first coordinate of the first CDP take
It is whole, obtain the second integer；
6th determining module 3100, if identical as the second integer for the first integer, it is determined that the first CMP is correspondingly
Shaking the corresponding CDP in road is the first CDP.
The application provides the determining device of CDP a kind of, is sat according to the corresponding coordinate of each shot point, each geophone station corresponding first
The corresponding relationship of mark and shot point and geophone station, determines the CMP of each seismic channel in seismic target earthquakes work area in preset work area coordinate
The first coordinate in system.Then, according to preset Algorithm of fitting a straight line, the first coordinate corresponding to each CMP carries out straight line fitting,
Target line is obtained, and establishes the first coordinate system by axis of abscissas of target line, determines each geophone station in the first coordinate system
The second coordinate in the first coordinate system of the second coordinate and each CMP.Later, according to corresponding second coordinate of each geophone station and in advance
If binning rule, multiple face elements are determined, and be directed to each face element, according to the CMP for including in the face element corresponding second
Coordinate determines first coordinate of the CDP of the corresponding seismic channel of the CMP for including in the face element in the first coordinate system.Finally, according to
Corresponding first coordinate of the CDP of each seismic channel determines second coordinate of the CDP of each seismic channel in preset work area coordinate system.
It should be understood that the system of determining investment combination provided by the above embodiment is only with above-mentioned each functional module
It divides and carries out for example, can according to need in practical application and be completed by different functional modules above-mentioned function distribution,
The internal structure of system is divided into different functional modules, to complete all or part of the functions described above.In addition,
The system of determining investment combination provided by the above embodiment and the embodiment of the method for determining investment combination belong to same design, have
Body realizes that process is detailed in embodiment of the method, and which is not described herein again.
Those of ordinary skill in the art will appreciate that realizing that all or part of the steps of above-described embodiment can pass through hardware
It completes, relevant hardware can also be instructed to complete by program, the program can store in a kind of computer-readable
In storage medium, storage medium mentioned above can be read-only memory, disk or CD etc..
The foregoing is merely the preferred embodiments of the application, not to limit the application, it is all in spirit herein and
Within principle, any modification, equivalent replacement, improvement and so on be should be included within the scope of protection of this application.
Claims (10)
1. a kind of determination method of common depth point CDP, which is characterized in that the described method includes:
Coordinate of each shot point in preset work area coordinate system, each geophone station are obtained in seismic target earthquakes work area in the preset work
The corresponding relationship of the first coordinate and shot point and geophone station in area's coordinate system；
According to the corresponding coordinate of each shot point, each corresponding first coordinate of geophone station and the shot point and geophone station
Corresponding relationship, determine the common point CMP of each seismic channel in the seismic target earthquakes work area in the preset work area coordinate system
In the first coordinate；
According to preset Algorithm of fitting a straight line, the first coordinate corresponding to each CMP carries out straight line fitting, obtains target line；
The first coordinate system is established by axis of abscissas of the target line, determines each geophone station in first coordinate system
The second coordinate in first coordinate system of the second coordinate and each CMP；
According to each corresponding second coordinate of geophone station and preset binning rule, multiple face elements are determined；
CMP pairs for including in the face element is determined according to corresponding second coordinate of the CMP for including in the face element for each face element
First coordinate of the CDP for the seismic channel answered in first coordinate system；
According to corresponding first coordinate of the CDP of each seismic channel, determine the CDP of each seismic channel in the preset work area
The second coordinate in coordinate system.
2. the method according to claim 1, wherein the Algorithm of fitting a straight line is least square method.
3. the method according to claim 1, wherein it is described according to corresponding second coordinate of each geophone station and
Preset binning rule, determines multiple face elements, comprising:
Determine the average value of the distance between abscissa of corresponding second coordinate of adjacent geophone station；
Using the average value and the ratio of preset constant as face element step-length；
According to minimum value, maximum value and the face element step-length in the abscissa of corresponding second coordinate of each geophone station, really
Fixed multiple face elements.
4. including according in the face element the method according to claim 1, wherein described be directed to each face element
Corresponding second coordinate of CMP determines the CDP of the corresponding seismic channel of the CMP for including in the face element in first coordinate system
First coordinate, comprising:
For each face element, the average value and ordinate of the abscissa of corresponding second coordinate of the CMP for including in the face element are determined
Average value, using the average value of the abscissa as the CDP for the corresponding seismic channel of CMP for including in the face element described first
The abscissa of the first coordinate in coordinate system, and the average value of the ordinate is corresponding as the CMP for including in the face element
The ordinate of corresponding first coordinate of the CDP of seismic channel.
5. according to the method described in claim 3, it is characterized in that, the method also includes:
The ratio of abscissa and the face element step-length to corresponding second coordinate of the first CMP is rounded, and obtains the first integer；
The ratio of abscissa and the face element step-length to corresponding first coordinate of the first CDP is rounded, and obtains the second integer；
If first integer is identical as second integer, it is determined that the corresponding CDP of the corresponding seismic channel of the first CMP
For the first CDP.
6. a kind of determining device of common depth point CDP, which is characterized in that described device includes:
Module is obtained, for obtaining coordinate, each geophone station of each shot point in preset work area coordinate system in seismic target earthquakes work area
The corresponding relationship of the first coordinate and shot point and geophone station in the preset work area coordinate system；
First determining module, for according to the corresponding coordinate of each shot point, corresponding first coordinate of each geophone station and
The corresponding relationship of the shot point and geophone station determines the common point CMP of each seismic channel in the seismic target earthquakes work area described
The first coordinate in preset work area coordinate system；
Fitting module, for according to preset Algorithm of fitting a straight line, the first coordinate corresponding to each CMP to carry out straight line fitting, obtains
To target line；
Second determining module determines each geophone station for establishing the first coordinate system by axis of abscissas of the target line
In the second coordinate of the second coordinate and each CMP in first coordinate system in first coordinate system；
Third determining module, for determining according to each corresponding second coordinate of geophone station and preset binning rule
Multiple face elements；
4th determining module, for being directed to each face element, according to corresponding second coordinate of the CMP for including in the face element, determining should
First coordinate of the CDP of the corresponding seismic channel of the CMP for including in face element in first coordinate system；
5th determining module determines each seismic channel for corresponding first coordinate of CDP according to each seismic channel
Second coordinate of the CDP in the preset work area coordinate system.
7. device according to claim 6, which is characterized in that the Algorithm of fitting a straight line is least square method.
8. device according to claim 6, which is characterized in that the third determining module is specifically used for:
Determine the average value of the distance between abscissa of corresponding second coordinate of adjacent geophone station；
Using the average value and the ratio of preset constant as face element step-length；
According to minimum value, maximum value and the face element step-length in the abscissa of corresponding second coordinate of each geophone station, really
Fixed multiple face elements.
9. device according to claim 6, which is characterized in that the 4th determining module is specifically used for:
For each face element, the average value and ordinate of the abscissa of corresponding second coordinate of the CMP for including in the face element are determined
Average value, using the average value of the abscissa as the CDP for the corresponding seismic channel of CMP for including in the face element described first
The abscissa of the first coordinate in coordinate system, and the average value of the ordinate is corresponding as the CMP for including in the face element
The ordinate of corresponding first coordinate of the CDP of seismic channel.
10. device according to claim 8, which is characterized in that described device further include:
First floor module, the ratio for abscissa and the face element step-length to corresponding second coordinate of the first CMP are rounded,
Obtain the first integer；
Second floor module, the ratio for abscissa and the face element step-length to corresponding first coordinate of the first CDP are rounded,
Obtain the second integer；
6th determining module, if identical as second integer for first integer, it is determined that the first CMP is corresponding
The corresponding CDP of seismic channel be the first CDP.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910271485.0A CN110068863A (en) | 2019-04-04 | 2019-04-04 | A kind of determination method and device of common depth point CDP |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910271485.0A CN110068863A (en) | 2019-04-04 | 2019-04-04 | A kind of determination method and device of common depth point CDP |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110068863A true CN110068863A (en) | 2019-07-30 |
Family
ID=67367159
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910271485.0A Pending CN110068863A (en) | 2019-04-04 | 2019-04-04 | A kind of determination method and device of common depth point CDP |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110068863A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112444847A (en) * | 2019-08-28 | 2021-03-05 | 中国石油天然气集团有限公司 | Method and device for determining outer boundary of seismic data of three-dimensional work area |
| CN114494571A (en) * | 2020-10-23 | 2022-05-13 | 中国石油化工股份有限公司 | Method, device, equipment and medium for extracting arbitrary section of three-dimensional seismic data volume |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5515335A (en) * | 1993-08-16 | 1996-05-07 | Atlantic Richfield Company | Seismic trace overburden correction method |
| CN102707316A (en) * | 2012-06-06 | 2012-10-03 | 中国石油天然气股份有限公司 | 2.5-time three-dimensional seismic exploration method for complex geological target |
| CN103091713A (en) * | 2011-10-28 | 2013-05-08 | 中国石油化工股份有限公司 | Automatic optimization method for undulating surface observation system |
| CN104502967A (en) * | 2014-12-23 | 2015-04-08 | 中国石油天然气集团公司 | Method and device for quickly obtaining surface element information of seismic exploration and observation system |
| CN106383361A (en) * | 2016-08-31 | 2017-02-08 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Speed data grid updating method |
| CN107300717A (en) * | 2016-04-15 | 2017-10-27 | 中国石油化工股份有限公司 | A kind of azimuth deviation of three dimensions is away from the homogenized method of trace gather |
-
2019
- 2019-04-04 CN CN201910271485.0A patent/CN110068863A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5515335A (en) * | 1993-08-16 | 1996-05-07 | Atlantic Richfield Company | Seismic trace overburden correction method |
| CN103091713A (en) * | 2011-10-28 | 2013-05-08 | 中国石油化工股份有限公司 | Automatic optimization method for undulating surface observation system |
| CN102707316A (en) * | 2012-06-06 | 2012-10-03 | 中国石油天然气股份有限公司 | 2.5-time three-dimensional seismic exploration method for complex geological target |
| CN104502967A (en) * | 2014-12-23 | 2015-04-08 | 中国石油天然气集团公司 | Method and device for quickly obtaining surface element information of seismic exploration and observation system |
| CN107300717A (en) * | 2016-04-15 | 2017-10-27 | 中国石油化工股份有限公司 | A kind of azimuth deviation of three dimensions is away from the homogenized method of trace gather |
| CN106383361A (en) * | 2016-08-31 | 2017-02-08 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Speed data grid updating method |
Non-Patent Citations (1)
| Title |
|---|
| 柳旭茏: "束线状三维正交观测系统满覆盖次数点坐标确定方法", 《石油天然气学报（江汉石油学院学报）》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112444847A (en) * | 2019-08-28 | 2021-03-05 | 中国石油天然气集团有限公司 | Method and device for determining outer boundary of seismic data of three-dimensional work area |
| CN114494571A (en) * | 2020-10-23 | 2022-05-13 | 中国石油化工股份有限公司 | Method, device, equipment and medium for extracting arbitrary section of three-dimensional seismic data volume |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Martin et al. | Seismic monitoring leveraging existing telecom infrastructure at the SDASA: Active, passive, and ambient-noise analysis | |
| CN110068863A (en) | A kind of determination method and device of common depth point CDP | |
| CN102970771B (en) | A kind of acquisition of signal acquisition method based on wireless sensor network and system | |
| CN106020757A (en) | Screen splicing method and electronic equipment | |
| Forghani-Arani et al. | An effective noise-suppression technique for surface microseismic data | |
| CN107505653B (en) | A kind of method and apparatus of determining migration before stack time result | |
| CN105204065B (en) | A kind of method and apparatus of first break picking ripple | |
| US20190383958A1 (en) | Method, seismic sensor and system for wireless seismic networking | |
| CN105785475A (en) | Hydraulic fracturing seismo-electric combined detection system, detection method and field work method | |
| CN107948953A (en) | Location determining method and device, the terminal of a kind of wearable device | |
| CN103870143A (en) | Method for data transmission control and electronic device | |
| Satoh | Empirical frequency-dependent radiation pattern of the 1998 Miyagiken-Nanbu earthquake in Japan | |
| CN105974473B (en) | A kind of method and apparatus for the seismic properties determining thin reservoir | |
| CN104299434A (en) | Road condition obtaining-presenting method and device | |
| CN107918623B (en) | Pipeline digital map display method and device | |
| CN106503071A (en) | The processing method and processing device of POI | |
| CN108921838A (en) | Method and apparatus for detection image | |
| US20180213081A1 (en) | Method and device for detecting a broken binder | |
| CN105444751A (en) | LED-lamp-based indoor positioning method, device and system | |
| CN109345610A (en) | The method and system of operational commanding are shown on map | |
| CN112950783B (en) | Three-dimensional dynamic loading method and device for raster data | |
| CN104459772A (en) | Frequency extracting method and device for earthquake digital signals | |
| CN111751880B (en) | Well earthquake data combined display method and device | |
| Stirling et al. | Assessment of the site conditions of precariously balanced rocks in the Mojave Desert, southern California | |
| CN104597498B (en) | A kind of geological data drawing methods and device |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190730 |
|
| WD01 | Invention patent application deemed withdrawn after publication |