CN110068863A - A kind of determination method and device of common depth point CDP - Google Patents

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

A method and device for determining a common depth point CDP

technical field

The present application relates to the technical field of seismic exploration, and in particular, to a method and device for determining a common depth point CDP.

Background technique

At present, seismic data is generally organized in seismic traces and stored in segy file format. The segy format is one of the standard tape data formats proposed by the Society of Exploration Geophysicists (SEG), and it is also one of the most common formats for seismic data in the petroleum exploration industry. A standard segy file generally consists of three parts. The first part is used to store the description information of the seismic data. The second part is used to store the description information of the segy file. The third part is used to store seismic trace information. The seismic trace information includes trace header information and seismic trace data, and the trace header information includes information such as line number, track number, sampling point number, geodetic coordinate, and common depth point (CDP) corresponding to the seismic trace.

Therefore, a method for determining CDP is urgently needed.

SUMMARY OF THE INVENTION

In order to solve the problems in the prior art, the embodiments of the present application provide a method and device for determining CDP, which can determine the CDP of a seismic trace. The technical solution is as follows:

In a first aspect, a method for determining CDP is provided, the method comprising:

Obtain the coordinates of each shot point in the preset work area coordinate system in the target seismic work area, the first coordinate of each detection point in the preset work area coordinate system, and the corresponding relationship between the shot point and the detection point;

According to the coordinates corresponding to each shot point, the first coordinate corresponding to each detection point, and the corresponding relationship between the shot point and the detection point, it is determined that the CMP of each seismic trace in the target seismic work area is in the preset The first coordinate in the work area coordinate system of ;

According to the preset straight line fitting algorithm, perform straight line fitting on the first coordinates corresponding to each CMP to obtain the target straight line;

A first coordinate system is established with the target straight line as the abscissa axis, and the second coordinate of each detection point in the first coordinate system and the second coordinate of each CMP in the first coordinate system are determined ;

According to the second coordinate corresponding to each detection point and the preset bin division rule, determine a plurality of bins;

For each surface element, according to the second coordinate corresponding to the CMP contained in the surface element, determine the first coordinate in the first coordinate system of the CDP of the seismic trace corresponding to the CMP contained in the surface element;

According to the first coordinate corresponding to the CDP of each seismic trace, the second coordinate of the CDP of each seismic trace in the preset work area coordinate system is determined.

Optionally, the straight line fitting algorithm is a least squares method.

Optionally, determining a plurality of surfels according to the second coordinates corresponding to the respective detection points and a preset surfel division rule, including:

Determine the average value of the distances between the abscissas of the second coordinates corresponding to the adjacent detection points;

Taking the ratio of the average value to the preset constant as the bin step size;

A plurality of bins are determined according to the minimum value and the maximum value in the abscissa of the second coordinate corresponding to each detection point and the bin step size.

Optionally, for each surface element, according to the second coordinate corresponding to the CMP included in the surface element, determine the first coordinate system of the CDP of the seismic trace corresponding to the CMP included in the surface element. a coordinate, including:

For each panel, determine the average value of the abscissa and the average value of the ordinate of the second coordinate corresponding to the CMP included in the panel, and use the average value of the abscissa as the CMP corresponding to the panel included in the panel. The abscissa of the first coordinate of the CDP of the seismic trace in the first coordinate system, and the average value of the ordinates is taken as the ordinate of the first coordinate corresponding to the CDP of the seismic trace corresponding to the CMP contained in the bin. coordinate.

Optionally, the method further includes:

Rounding the ratio of the abscissa of the second coordinate corresponding to the first CMP to the step size of the panel to obtain the first integer;

rounding the ratio of the abscissa of the first coordinate corresponding to the first CDP to the step size of the bin to obtain a second integer;

If the first integer is the same as the second integer, it is determined that the CDP corresponding to the seismic trace corresponding to the first CMP is the first CDP.

In a second aspect, a device for determining CDP is provided, the device comprising:

The acquisition module is used to acquire the coordinates of each shot point in the preset work area coordinate system in the target seismic work area, the first coordinate of each detection point in the preset work area coordinate system, and the correspondence between the shot point and the detection point relation;

The first determination module is used to determine each seismic trace in the target seismic work area according to the coordinates corresponding to each shot point, the first coordinate corresponding to each of the detection points, and the correspondence between the shot point and the detection point the first coordinate of the CMP in the preset work area coordinate system;

a fitting module, configured to perform straight line fitting on the first coordinates corresponding to each CMP according to a preset straight line fitting algorithm to obtain a target straight line;

The second determination module is configured to establish a first coordinate system with the target straight line as the abscissa axis, and determine the second coordinates of the detection points in the first coordinate system and the CMPs in the first coordinate system. the second coordinate in the coordinate system;

a third determination module, configured to determine a plurality of surfels according to the second coordinates corresponding to the detection points and a preset surfel division rule;

The fourth determination module is configured to, for each surface element, determine the CDP of the seismic trace corresponding to the CMP included in the surface element in the first coordinate system according to the second coordinate corresponding to the CMP included in the surface element. first coordinate;

The fifth determination module is configured to determine the second coordinate of the CDP of each seismic trace in the preset work area coordinate system according to the first coordinate corresponding to the CDP of each seismic trace.

Optionally, the straight line fitting algorithm is a least squares method.

Optionally, the third determining module is specifically used for:

Determine the average value of the distances between the abscissas of the second coordinates corresponding to the adjacent detection points;

Taking the ratio of the average value to the preset constant as the bin step size;

A plurality of bins are determined according to the minimum value and the maximum value in the abscissa of the second coordinate corresponding to each detection point and the bin step size.

Optionally, the fourth determining module is specifically used for:

For each panel, determine the average value of the abscissa and the average value of the ordinate of the second coordinate corresponding to the CMP included in the panel, and use the average value of the abscissa as the CMP corresponding to the panel included in the panel. The abscissa of the first coordinate of the CDP of the seismic trace in the first coordinate system, and the average value of the ordinates is taken as the ordinate of the first coordinate corresponding to the CDP of the seismic trace corresponding to the CMP contained in the bin. coordinate.

Optionally, the device further includes:

The first rounding module is used for rounding the ratio of the abscissa of the second coordinate corresponding to the first CMP to the step size of the panel to obtain the first integer;

The second rounding module is used for rounding the ratio of the abscissa of the first coordinate corresponding to the first CDP to the step size of the panel to obtain the second integer;

A sixth determining module, configured to determine that the CDP corresponding to the seismic trace corresponding to the first CMP is the first CDP if the first integer is the same as the second integer.

The present application provides a CDP determination method and device. According to the coordinates corresponding to each shot point, the first coordinate corresponding to each detection point, and the corresponding relationship between the shot point and the detection point, it is determined that the CMP of each seismic trace in the target seismic work area is at The first coordinate in the preset work area coordinate system. Then, according to a preset straight line fitting algorithm, perform straight line fitting on the first coordinates corresponding to each CMP to obtain a target straight line, and establish a first coordinate system with the target straight line as the abscissa axis, and determine that each detection point is at the first coordinate The second coordinate in the system and the second coordinate of each CMP in the first coordinate system. Then, according to the second coordinate corresponding to each detection point and the preset bin division rule, a plurality of bins are determined, and for each bin, the bin is determined according to the second coordinate corresponding to the CMP included in the bin. The first coordinate in the first coordinate system of the CDP of the seismic trace corresponding to the CMP contained in the element. Finally, according to the first coordinate corresponding to the CDP of each seismic trace, the second coordinate of the CDP of each seismic trace in the preset work area coordinate system is determined.

Description of drawings

In order to illustrate the technical solutions in the embodiments of the present application more clearly, the following briefly introduces the drawings that are used in the description of the embodiments. Obviously, the drawings in the following description are only some embodiments of the present application. For those of ordinary skill in the art, other drawings can also be obtained from these drawings without creative effort.

1 is a flowchart of a method for determining a CDP according to an embodiment of the present application;

2 is a schematic diagram of a straight line fitting provided by an embodiment of the present application;

3 is a schematic structural diagram of an apparatus for determining a CDP according to an embodiment of the present application;

FIG. 4 is a schematic structural diagram of an apparatus for determining a CDP according to an embodiment of the present application.

Detailed ways

In order to make the objectives, technical solutions and advantages of the present application clearer, the embodiments of the present application will be further described in detail below with reference to the accompanying drawings.

The present application provides a method for determining CDP, the execution subject of the method may be a network device such as a server and a host, and the method can be applied to the technical field of seismic exploration. Based on this method, the CDP corresponding to the seismic trace can be determined. The embodiments of the present application are described by taking the execution subject as the server as an example, and other situations are similar.

A method for determining a CDP provided by the embodiment of the present application will be described in detail below with reference to the specific implementation manner, as shown in FIG. 1 , and the specific steps are as follows:

Step 101: Obtain the coordinates of each shot point in the preset work area coordinate system in the target seismic work area, the first coordinate of each detection point in the preset work area coordinate system, and the corresponding relationship between the shot point and the detection point.

In implementation, a work area coordinate system may be pre-stored in the server, and the work area coordinate system may be set by a technician according to experience. When it is necessary to determine the CDP of each seismic trace in a certain seismic work area (that is, the target seismic work area), the server can obtain the coordinates of each shot point in the target seismic work area in the preset work area coordinate system, and each detection point in the preset work area. The first coordinate in the coordinate system of the work area, and the correspondence between the shot point and the receiver point.

Step 102, according to the coordinates corresponding to each shot point, the first coordinate corresponding to each detection point, and the correspondence between the shot point and the detection point, determine the CMP of each seismic trace in the target seismic work area in the preset work area coordinate system. a coordinate.

In the implementation, the server obtains the coordinates of each shot point in the preset work area coordinate system in the target seismic work area, the first coordinate of each detection point in the preset work area coordinate system, and the correspondence between the shot point and the detection point After the relationship, the common center point (Common Middle Point, CMP) of each seismic trace in the target seismic work area can be further determined according to the coordinates corresponding to each shot point, the first coordinate corresponding to each detection point, and the corresponding relationship between the shot point and the detection point. ) is the first coordinate in the preset work area coordinate system. The server determines, according to the coordinates corresponding to each shot point, the first coordinate corresponding to each detection point, and the corresponding relationship between the shot point and the detection point, the first CMP of each seismic trace in the target seismic work area in the preset work area coordinate system. The processing process of one coordinate and the first coordinate corresponding to the CMP of each seismic trace are determined according to the coordinates corresponding to each shot point, the first coordinate corresponding to each detection point, and the corresponding relationship between the shot point and the detection point disclosed in the existing scheme. The processing process is similar and will not be repeated here.

Step 103, according to a preset straight line fitting algorithm, perform straight line fitting on the first coordinates corresponding to each CMP to obtain a target straight line.

In implementation, the straight line fitting algorithm may be pre-stored in the server. The straight line fitting algorithm can be set by the skilled person based on experience. Optionally, the straight line fitting algorithm may be the least squares method, or may be other straight line fitting algorithms, which are not limited in this application.

After the server obtains the first coordinates of the CMP of each seismic trace in the target seismic work area in the preset work area coordinate system, it can perform line fitting on the first coordinates corresponding to each CMP according to the preset line fitting algorithm to obtain the target straight line.

FIG. 2 is a schematic diagram of a straight line fitting provided by an embodiment of the present application. As shown in FIG. 2 , in the preset coordinate system of the work area, the server performs a straight line on the first coordinates corresponding to each CMP according to a preset straight line fitting algorithm. Fit to get the target straight line.

Step 104, establishing a first coordinate system with the target straight line as the abscissa axis, and determining the second coordinates of each detection point in the first coordinate system and the second coordinates of each CMP in the first coordinate system.

In implementation, after obtaining the target straight line, the server may further establish the first coordinate system with the target straight line as the abscissa axis. The origin of the first coordinate system may be the intersection of the target straight line and the abscissa axis of the work area coordinate system, or may be the intersection of the target straight line and the ordinate axis of the work area coordinate system, or may be other points on the target straight line , the embodiments of the present application are not limited. After the server establishes the first coordinate system, it may further determine the second coordinate of each detection point in the first coordinate system and the second coordinate of each CMP in the first coordinate system.

Step 105: Determine a plurality of bins according to the second coordinates corresponding to each detection point and a preset bin division rule.

In implementation, the surfel division rule may be pre-stored in the server, and the surfel division rule may be set by a technician according to experience. After obtaining the second coordinates of each detection point in the first coordinate system, the server may determine a plurality of bins according to the second coordinates corresponding to each detection point and the preset bin division rule.

Optionally, the server determines, according to the second coordinate corresponding to each detection point and the preset surfel division rule, the processing process for determining multiple surfels as follows:

Step 1: Determine the average value of the distances between the abscissas of the second coordinates corresponding to the adjacent detection points.

In implementation, the server may determine adjacent detection points according to the abscissa of the second coordinate corresponding to each detection point. Then, the server may determine the distance between the abscissas of the second coordinates corresponding to the adjacent detection points. After that, the server may determine the average value of the distances between the abscissas of the second coordinates corresponding to all adjacent detection points.

In step 2, the ratio of the average value to the preset constant is used as the bin step size.

In implementation, preset constants may be pre-stored in the server. The preset constant can be set by the technician according to experience. Optionally, the preset constant may be 2, or may be other constants, which are not limited in this embodiment of the present application. After obtaining the average value of the distances between the abscissas of the second coordinates corresponding to the adjacent detection points, the server may further calculate the ratio of the average value to the preset constant, and use the ratio as the bin step size. For example, the average value of the distance between the abscissas of the second coordinates corresponding to the adjacent detection points is 4, the preset constant is 2, and the bin step size is 2.

Step 3: Determine a plurality of bins according to the minimum value, the maximum value and the bin step size in the abscissa of the second coordinate corresponding to each detection point.

In implementation, after obtaining the bin step size, the server may further use the minimum value in the abscissa of the second coordinate corresponding to each detection point as the starting point, and use the maximum value in the abscissa of the second coordinate corresponding to each detection point as At the end point, the surfels are divided according to the surfel step size, and the range of each surfel is determined. For example, the minimum value in the abscissa of the second coordinate corresponding to each detection point is 2, the maximum value in the abscissa of the second coordinate corresponding to each detection point is 10, and the bin step size is 2, then the first surface The range of the bins is 2 to 4, the range of the second bin is 4 to 6, the range of the third bin is 6 to 8, and the range of the fourth bin is 8 to 10.

Step 106 , for each bin, according to the second coordinate corresponding to the CMP contained in the bin, determine the first coordinate in the first coordinate system of the CDP of the seismic trace corresponding to the CMP contained in the bin.

In implementation, for each surfel, the server may determine the first coordinate in the first coordinate system of the CDP of the seismic trace corresponding to the CMP contained in the surfel according to the second coordinate corresponding to the CMP contained in the surfel .

Optionally, for each surfel, the server determines, according to the second coordinate corresponding to the CMP contained in the surfel, the first coordinate of the CDP of the seismic trace corresponding to the CMP contained in the surfel in the first coordinate system. The processing process is: for each surface element, determine the average value of the abscissa and the average value of the vertical coordinate of the second coordinate corresponding to the CMP included in the surface element, and use the average value of the abscissa as the CMP included in the surface element. The abscissa of the first coordinate of the CDP of the corresponding seismic trace in the first coordinate system, and the average value of the ordinates is taken as the ordinate of the first coordinate corresponding to the CDP of the seismic trace corresponding to the CMP contained in the bin.

In implementation, for each surfel, the server may determine the average value of the abscissa and the average value of the ordinate of the second coordinate corresponding to the CMP included in the surfel. Then, the server may take the average value of the abscissa as the abscissa of the first coordinate in the first coordinate system of the CDP of the seismic trace corresponding to the CMP included in the bin, and take the average value of the ordinate as the The ordinate of the first coordinate corresponding to the CDP of the seismic trace corresponding to the included CMP.

Step 107: Determine the second coordinate of the CDP of each seismic trace in the preset work area coordinate system according to the first coordinate corresponding to the CDP of each seismic trace.

In implementation, after obtaining the first coordinate of the CDP of each seismic trace in the first coordinate system, the server may further determine the second coordinate of the CDP of each seismic trace in the preset work area coordinate system (that is, the CDP).

Optionally, the server can also determine the corresponding relationship between the seismic trace and the CDP according to the second coordinate corresponding to the CMP of each seismic trace and the first coordinate corresponding to each CDP, and the processing process is as follows:

In step 1, the ratio of the abscissa of the second coordinate corresponding to the first CMP to the bin step size is rounded to obtain the first integer.

In implementation, for a certain CMP (that is, the first CMP) in the CMPs of each seismic trace, the server may round the ratio of the abscissa of the second coordinate corresponding to the first CMP to the bin step size to obtain the first CMP. Integer. For example, if the abscissa of the second coordinate corresponding to the first CMP is 10, and the bin step is 3, the first integer is 3.

In step 2, the ratio of the abscissa of the first coordinate corresponding to the first CDP to the bin step size is rounded to obtain a second integer.

In implementation, for a certain CDP (ie, the first CDP), the server may round the ratio of the abscissa of the first coordinate corresponding to the first CDP to the bin step size to obtain the second integer. For example, the abscissa of the first coordinate corresponding to the first CDP is 11, the bin step is 3, and the second integer is 3.

Step 3: If the first integer is the same as the second integer, determine that the CDP corresponding to the seismic trace corresponding to the first CMP is the first CDP.

In an implementation, after obtaining the first integer and the second integer, the server may determine whether the first integer and the second integer are the same. If the first integer is the same as the second integer, it means that the first CMP and the first CDP are in the same bin, and the server can determine that the CDP corresponding to the seismic trace corresponding to the first CMP is the first CDP.

The present application provides a method for determining CDP. According to the coordinates corresponding to each shot point, the first coordinate corresponding to each detection point, and the corresponding relationship between the shot point and the detection point, it is determined that the CMP of each seismic trace in the target seismic work area is at a preset value. The first coordinate in the work area coordinate system of . Then, according to a preset straight line fitting algorithm, perform straight line fitting on the first coordinates corresponding to each CMP to obtain a target straight line, and establish a first coordinate system with the target straight line as the abscissa axis, and determine that each detection point is at the first coordinate The second coordinate in the system and the second coordinate of each CMP in the first coordinate system. Then, according to the second coordinate corresponding to each detection point and the preset bin division rule, a plurality of bins are determined, and for each bin, the bin is determined according to the second coordinate corresponding to the CMP included in the bin. The first coordinate in the first coordinate system of the CDP of the seismic trace corresponding to the CMP contained in the element. Finally, according to the first coordinate corresponding to the CDP of each seismic trace, the second coordinate of the CDP of each seismic trace in the preset work area coordinate system is determined.

Based on the same technical concept, as shown in FIG. 3 , an embodiment of the present application further provides a CDP determination device, which includes:

The acquisition module 310 is used to acquire the coordinates of each shot point in the preset work area coordinate system in the target seismic work area, the first coordinate of each detection point in the preset work area coordinate system, and the corresponding relationship between the shot point and the detection point ;

The first determination module 320 is used to determine that the CMP of each seismic trace in the target seismic work area is in the preset work area according to the coordinates corresponding to each shot point, the first coordinate corresponding to each detection point, and the corresponding relationship between the shot point and the detection point. the first coordinate in the coordinate system;

The fitting module 330 is configured to perform straight line fitting on the first coordinates corresponding to each CMP according to a preset straight line fitting algorithm to obtain a target straight line;

The second determination module 340 is used for establishing the first coordinate system with the target straight line as the abscissa axis, and determining the second coordinate of each detection point in the first coordinate system and the second coordinate of each CMP in the first coordinate system;

The third determining module 350 is configured to determine a plurality of bins according to the second coordinates corresponding to each detection point and a preset bin division rule;

The fourth determination module 360 is configured to, for each surface element, determine, according to the second coordinate corresponding to the CMP contained in the surface element, the first coordinate of the CDP of the seismic trace corresponding to the CMP contained in the surface element in the first coordinate system. a coordinate;

The fifth determination module 370 is configured to determine the second coordinate of the CDP of each seismic trace in the preset work area coordinate system according to the first coordinate corresponding to the CDP of each seismic trace.

Optionally, the straight line fitting algorithm is the least squares method.

Optionally, the third determining module 350 is specifically configured to:

Determine the average value of the distances between the abscissas of the second coordinates corresponding to the adjacent detection points;

Take the ratio of the average value to the preset constant as the bin step size;

According to the minimum value, the maximum value and the bin step size in the abscissa of the second coordinate corresponding to each detection point, a plurality of bins are determined.

Optionally, the fourth determining module 360 is specifically configured to:

For each panel, determine the average value of the abscissa and the average value of the ordinate of the second coordinate corresponding to the CMP contained in the panel, and use the average value of the abscissa as the seismic trace corresponding to the CMP contained in the panel. The abscissa of the first coordinate of the CDP in the first coordinate system, and the average value of the ordinates is taken as the ordinate of the first coordinate corresponding to the CDP of the seismic trace corresponding to the CMP contained in the bin.

Optionally, as shown in Figure 4, the device further includes:

The first rounding module 380 is used to round the ratio of the abscissa of the second coordinate corresponding to the first CMP to the bin step size to obtain the first integer;

The second rounding module 390 is used for rounding the ratio of the abscissa of the first coordinate corresponding to the first CDP to the bin step size to obtain the second integer;

The sixth determining module 3100 is configured to determine the CDP corresponding to the seismic trace corresponding to the first CMP as the first CDP if the first integer is the same as the second integer.

The application provides a CDP determination device, according to the coordinates corresponding to each shot point, the first coordinate corresponding to each detection point, and the corresponding relationship between the shot point and the detection point, it is determined that the CMP of each seismic trace in the target seismic work area is at a preset value The first coordinate in the work area coordinate system of . Then, according to a preset straight line fitting algorithm, perform straight line fitting on the first coordinates corresponding to each CMP to obtain a target straight line, and establish a first coordinate system with the target straight line as the abscissa axis, and determine that each detection point is at the first coordinate The second coordinate in the system and the second coordinate of each CMP in the first coordinate system. Then, according to the second coordinate corresponding to each detection point and the preset bin division rule, a plurality of bins are determined, and for each bin, the bin is determined according to the second coordinate corresponding to the CMP included in the bin. The first coordinate in the first coordinate system of the CDP of the seismic trace corresponding to the CMP contained in the element. Finally, according to the first coordinate corresponding to the CDP of each seismic trace, the second coordinate of the CDP of each seismic trace in the preset work area coordinate system is determined.

It should be noted that the system for determining an investment portfolio provided by the above-mentioned embodiments is only illustrated by the division of the above-mentioned functional modules. The internal structure is divided into different functional modules to complete all or part of the functions described above. In addition, the system for determining an investment portfolio provided by the above embodiment and the method for determining an investment portfolio belong to the same concept, and the specific implementation process is detailed in the method embodiment, which will not be repeated here.

Those of ordinary skill in the art can understand that all or part of the steps of implementing the above embodiments can be completed by hardware, or can be completed by instructing relevant hardware through a program, and the program can be stored in a computer-readable storage medium. The storage medium mentioned may be a read-only memory, a magnetic disk or an optical disk, etc.

The above descriptions are only preferred embodiments of the present application, and are not intended to limit the present application. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present application shall be included in the protection of the present application. within the range.

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.