CN111007580B - Method and system for evaluating relative consistency of surface elements of three-dimensional seismic observation system - Google Patents
Method and system for evaluating relative consistency of surface elements of three-dimensional seismic observation system Download PDFInfo
- Publication number
- CN111007580B CN111007580B CN201811168459.7A CN201811168459A CN111007580B CN 111007580 B CN111007580 B CN 111007580B CN 201811168459 A CN201811168459 A CN 201811168459A CN 111007580 B CN111007580 B CN 111007580B
- Authority
- CN
- China
- Prior art keywords
- dimensional space
- surface element
- dimensional
- bin
- observation system
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
Abstract
A method and a system for evaluating the relative consistency of three-dimensional seismic observation system bins are disclosed. The method comprises the following steps: 1) inputting an azimuth angle-offset data pair of each surface element of the three-dimensional earthquake observation system; 2) mapping each data pair in respective surface element two-dimensional space; 3) taking one surface element as a center, mapping data pairs of the surrounding 9 surface elements including the surface element into a two-dimensional space of a contrast surface element; 4) calculating the correlation degree by using the data of the two-dimensional space of the surface element serving as the center in the step 3) and the two-dimensional space of the contrast surface element; 5) and repeating the steps 3) and 4), traversing all the surface elements of the three-dimensional observation system, and calculating the correlation of each surface element. The invention quantitatively calculates the correlation degree of the two-dimensional space distribution of each surface element azimuth angle and shot-geophone distance with a compared surface element azimuth angle and shot-geophone distance, quantitatively evaluates the consistency of the distribution of the azimuth angle and shot-geophone distance of the surface element of the observation system, and can give out an evaluation result through an intuitive data graph.
Description
Technical Field
The invention relates to the field of three-dimensional seismic exploration and observation systems, in particular to a method and a system for evaluating the relative bin consistency of a three-dimensional seismic observation system.
Background
The three-dimensional seismic exploration technology is the mainstream technology of petroleum exploration at present. The quality of the three-dimensional seismic observation system directly influences the quality of the acquired seismic data, and further influences the seismic exploration effect. Therefore, evaluation of the three-dimensional seismic observation system is important. The bin property of the three-dimensional seismic observation system is one of important means for evaluating the three-dimensional observation system, wherein the distribution of bin azimuth angles and shot-geophone distances is the most important bin property.
From the view of collection and imaging, an ideal three-dimensional observation system requires all-directional bin azimuth angles and uniform distribution of shot-geophone distances. The observation system with wide azimuth and offset evenly distributed can effectively attenuate and suppress multiple waves, ground roll waves, offset noise and other various random interferences and noises.
At present, evaluation of the azimuth angle and the offset of a surface element of a three-dimensional observation system is mainly represented by an azimuth spider graph, an offset line graph and the like of a single surface element. On one hand, the representation method cannot quantitatively evaluate the difference between different surface elements and basically depends on the observation and experience of designers, and on the other hand, the comprehensive evaluation of offset and azimuth is not carried out. In recent years, quantitative evaluation of bin offset and azimuth distributions by one-dimensional cross-correlation has also been attempted. However, these methods either only consider the distribution of offsets, or even if the azimuth problem is considered, consider the offsets first, and then evaluate the distribution of offsets at different azimuths one by one, without comprehensively considering the azimuths of the offsets as a whole. Therefore, it is expected that a method and a system for evaluating the relative consistency of surface elements can be designed for a seismic exploration observation system.
The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides an evaluation method and system based on the consistency of azimuth angles and offset distances, and the optimization and selection of a seismic exploration observation system are solved.
According to one aspect of the invention, the method for evaluating the relative bin consistency of the three-dimensional seismic observation system comprises the following steps:
1) inputting an azimuth angle-offset data pair of each surface element of the three-dimensional earthquake observation system;
2) mapping the azimuth angle-offset data pairs of each surface element into respective surface element two-dimensional space;
3) taking one surface element, taking the surface element as a center, and mapping azimuth angle-offset data pairs of the surrounding 9 surface elements including the surface element into a two-dimensional space of a contrast surface element;
4) calculating the correlation degree by using the data of the two-dimensional space of the surface element serving as the center in the step 3) and the two-dimensional space of the contrast surface element;
5) and repeating the steps 3) and 4), traversing all surface elements of the three-dimensional observation system, and calculating the correlation of each surface element.
Preferably, the step 2) includes:
2.1) constructing an M N two-dimensional space G for one of the binsavoAnd is GavoSetting a counter for each element, where M is 360/d theta, d theta is azimuth step, and N is Xmax/dx，XmaxThe maximum offset of the observation system is adopted, and dx is offset step length;
2.2) traverse all azimuth-offset data pairs of the bins, mapping each data pair to Gavo[idx,idy]Each time the mapping is successful, adding 1 to a counter of the corresponding element, wherein idx is θ/d θ, and idy is x/dx;
2.3) traversing all the bins, repeating the steps 2.1) -2.2), and mapping the azimuth angle-offset data pairs of each bin in the respective two-dimensional space.
Preferably, the mapping of the azimuth-offset data pairs of the surrounding 9 bins including the bin into the two-dimensional space of the contrast bin in the step 3) comprises:
3.1) construction of an MxN two-dimensional space G of contrast binsstaAnd is GstaSetting a counter for each element;
3.2) mapping the azimuth-offset data pair of 9 bins to the contrast bin two-dimensional space Gsta，GstaEach element in (1) is mapped once successfully, then pairAdding 1 to a corresponding calculator;
3.3) mixing GstaThe counter result for each element of (a) is divided by 9 to take the average.
Preferably, the step 4) includes:
4.1) calculating the mean value of each element in the two-dimensional space of the contrast surface element based on the formula (1)
4.2) calculating the mean value of each element in the bin two-dimensional space based on the formula (2)
4.3) calculating the correlation R of the two-dimensional space of the contrast surface element and the two-dimensional space of the surface element taken in the step 3) based on the formula (3):
preferably, the method further comprises a step 6), which specifically comprises the following steps:
6.1) selecting the maximum value R of the correlation R of the two-dimensional space of all the surface elements and the two-dimensional space of the contrast surface elementmax；
6.2) mixing of [0, Rmax]Dividing the range into L sections, traversing the correlation degree R of all surface elements, and counting the number R of the correlation degree belonging to each sectionnum[i]Wherein i is 0,1,2, …, L.
According to another aspect of the present invention, a three-dimensional seismic observation system bin relative consistency evaluation system is provided, which is stored in a computer program, and the program is executed by a processor to perform the following steps:
1) inputting an azimuth angle-offset data pair of each surface element of the three-dimensional earthquake observation system;
2) mapping the azimuth angle-offset data pairs of each surface element into respective surface element two-dimensional space;
3) taking one surface element as a center, and mapping azimuth angle-offset data pairs of the surrounding 9 surface elements including the surface element into a two-dimensional space of a contrast surface element;
4) calculating the correlation degree by using the data of the two-dimensional space of the surface element serving as the center in the step 3) and the two-dimensional space of the contrast surface element;
5) and repeating the steps 3) and 4), traversing all surface elements of the three-dimensional observation system, and calculating the correlation of each surface element.
Preferably, the step 2) includes:
2.1) constructing an M N two-dimensional space G for one of the binsavoAnd is GavoSetting a counter for each element, where M is 360/d theta, d theta is azimuth step, and N is Xmax/dx，XmaxThe maximum offset of the observation system is adopted, and dx is offset step length;
2.2) traverse all azimuth-offset data pairs of the bins, mapping each data pair to Gavo[idx,idy]Each time the mapping is successful, adding 1 to a counter of the corresponding element, wherein idx is θ/d θ, and idy is x/dx;
2.3) traversing all the bins, repeating the steps 2.1) -2.2), and mapping the azimuth angle-offset data pairs of each bin in the respective two-dimensional space.
Preferably, the mapping of the azimuth-offset data pairs of the surrounding 9 bins including the bin into the two-dimensional space of the contrast bin in the step 3) comprises:
3.1) construction of an MxN two-dimensional space G of contrast binsstaAnd is GstaSetting a counter for each element;
3.2) mapping the azimuth-offset data pair of 9 bins to the contrast bin two-dimensional space Gsta，GstaEach element in (1) is mapped once successfully, and the corresponding calculator is added with 1;
3.3) mixing GstaThe counter result for each element of (a) is divided by 9 to take the average.
Preferably, the step 4) includes:
4.1) calculating the mean value of each element in the two-dimensional space of the contrast surface element based on the formula (1)
4.2) calculating the mean value of each element in the bin two-dimensional space based on the formula (2)
4.3) calculating the correlation R of the two-dimensional space of the contrast surface element and the two-dimensional space of the surface element taken in the step 3) based on the formula (3):
preferably, the processor executes the program to further implement the following steps:
6.1) selecting the maximum value R of the correlation R of the two-dimensional space of all the surface elements and the two-dimensional space of the contrast surface elementmax；
6.2) mixing of [0, Rmax]Dividing the range into L sections, traversing the correlation degree R of all surface elements, and counting the number R of the correlation degree belonging to each sectionnum[i]Wherein i is 0,1,2, …, L.
According to the invention, the bin azimuth angle and the offset distance data are mapped in a two-dimensional space, the correlation degree of the two-dimensional space of each bin azimuth angle and offset distance with the two-dimensional space distribution of a compared bin azimuth angle and offset distance is calculated quantitatively by a two-dimensional space cross-correlation method, the consistency of the azimuth angle and offset distance distribution of the observation system bin is evaluated quantitatively, and an evaluation result can be given by an intuitive data graph. The designer can adjust the designed observation system based on the result, or in many alternatives, prefer the best. The method and the system are feasible and can solve the problems of evaluation and optimization of the binning consistency of the observation system in practical application.
The method and system of the present invention have other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts.
FIG. 1 illustrates a flow chart of a three-dimensional seismic observation system bin relative consistency evaluation method according to an exemplary embodiment of the present invention;
FIG. 2 shows a schematic of offset data in one bin mapped to two-dimensional space;
3a-3c show the correlation output of the 34L8S240T orthogonal observation system 1;
4a-4c show the correlation output of the 34L8S240T orthogonal observation system 2;
fig. 5a-5c show the correlation output results for the 34L8S240T orthogonal observation system 3.
Detailed Description
The invention will be described in more detail below with reference to the accompanying drawings. While the preferred embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
FIG. 1 shows a flow diagram of a three-dimensional seismic observation system bin relative consistency evaluation method according to an exemplary embodiment of the invention. As shown in FIG. 1, the method includes steps S1-S5.
In step S1, azimuth-offset data pairs for each bin of the three-dimensional seismic survey system are input.
The three-dimensional seismic observation system is an observation mode for describing the spatial position relationship between excitation points and receiving points distributed on the area during three-dimensional seismic data acquisition. The bin is a basic unit (like a pixel) of three-dimensional imaging, each shot point and wave pair are corresponding to different bins, the shot-receiver pair combination of the same shot but different channels is corresponding to different bins, and the shot-receiver pair combination of the same shot but different shots is also corresponding to different bins.
In step S2, the azimuth-offset data pairs for each bin are mapped into the respective bin two-dimensional space. Referring to fig. 2, Xmax is 3500 m, dx is 100 m, and d θ is 15 °.
Specifically, the mapping of azimuth-offset data to bin two-dimensional space can be realized by:
firstly, an azimuth angle step length d theta is specified, and the azimuth angle is equally divided into M equal parts, wherein M is 360/d theta; and designating a offset step dx, and dividing the offset into N equal parts, wherein N is Xmax/dx, and Xmax is the maximum offset of the observation system.
One surface element is specified, and an M multiplied by N two-dimensional space G is constructedavoAnd is GavoEach element is provided with a counter.
Traversing all azimuth-offset data pairs of the bin, mapping each avo data pair to G using the following formulaavoCorresponding element Gavo[idx,idy]In each mapping success, make Gavo[idx,idy]Counter + 1:
and traversing each surface element, repeating the steps in the two sections, constructing a two-dimensional space for each surface element, and mapping the azimuth angle-offset data to the space.
In step S3, the azimuth-offset data pairs of the surrounding 9 bins including one of the bins are mapped into a two-dimensional space of contrast bins centered on the one of the bins.
Specifically, similarly to the content of step S2, a contrast bin two-dimensional space G is first createdstaThe same size is M × N and GstaEach element is provided with a counter.
For the selected bin, the azimuth-offset data of 9 bins, which are the bin and 8 bins surrounding the bin, are mapped to the contrast bin space G, using the bin as the center, in a manner similar to that in step S2sta。
G is to bestaIs averaged with the counter result of each element of (1), i.e. Gsta[i][j]And/9, wherein i is 1,2, … M, and j is 1,2, …, N.
In step S4, a correlation is calculated using the data of the bin two-dimensional space and the contrast bin two-dimensional space, which are the centers in step S3.
Specifically, the correlation degree may be calculated by the following method;
for the selected surface element, calculating the mean value of the contrast surface element according to the formula (1)
For the selected surface element, calculating the mean value of the surface element according to the formula (2)
Calculating the correlation value R of the selected surface element according to the formula (3):
in step S5, steps S3 and S4 are repeated, and all bins of the three-dimensional observation system are traversed to calculate the correlation of each bin.
In one example, the method further includes step S6, which specifically includes:
6.1) selecting the maximum value R of the correlation R of the two-dimensional space of all the surface elements and the two-dimensional space of the contrast surface elementmax；
6.2) mixing of [0, Rmax]Dividing the range into L sections, traversing the correlation degree R of all surface elements, and counting the number R of the correlation degree belonging to each sectionnum[i]Wherein i is 0,1,2, …, L.
The statistical results can be visually displayed in the following graphical manner:
1. and drawing a color plane graph by taking the value of the correlation R of the surface element as a color scale and the central position of the surface element as a coordinate, wherein the color depth of the plane graph indicates the height of the value of the correlation R of the surface element. Those skilled in the art will understand that the degree of correlation R can be represented by gray scale;
2. r counted in step 6.2)numDrawing a bar graph as a result, wherein each bar in the bar graph represents the number of the surface elements occupied by the segment;
3. r counted in step 6.2)numThe result is a sector plot, where each sector in the sector plot represents the percentage of the segmented bin number to the total bin number.
The invention also provides a three-dimensional seismic observation system bin relative consistency evaluation system, which is stored in a computer program, wherein the program is executed by a processor and comprises the following steps:
1) inputting an azimuth angle-offset data pair of each surface element of the three-dimensional earthquake observation system;
2) mapping the azimuth angle-offset data pairs of each surface element into respective surface element two-dimensional space;
3) taking one surface element as a center, and mapping azimuth angle-offset data pairs of the surrounding 9 surface elements including the surface element into a two-dimensional space of a contrast surface element;
4) calculating the correlation degree by using the data of the two-dimensional space of the surface element serving as the center in the step 3) and the two-dimensional space of the contrast surface element;
5) and repeating the steps 3) and 4), traversing all surface elements of the three-dimensional observation system, and calculating the correlation of each surface element.
In one example, the step 2) includes:
2.1) constructing an M N two-dimensional space G for one of the binsavoAnd is GavoSetting a counter for each element, where M is 360/d theta, d theta is azimuth step, and N is Xmax/dx，XmaxThe maximum offset of the observation system is adopted, and dx is offset step length;
2.2) traverse all azimuth-offset data pairs of the bins, mapping each data pair to Gavo[idx,idy]Each time the mapping is successful, adding 1 to a counter of the corresponding element, wherein idx is θ/d θ, and idy is x/dx;
2.3) traversing all the bins, repeating the steps 2.1) -2.2), and mapping the azimuth angle-offset data pairs of each bin in the respective two-dimensional space.
In one example, the mapping of the azimuth-offset data pair of the surrounding 9 bins including the bin into the contrast bin two-dimensional space in step 3) comprises:
3.1) construction of an MxN two-dimensional space G of contrast binsstaAnd is GstaSetting a counter for each element;
3.2) the azimuth-offset data pairs of the 9 binsMapping to the contrast bin two-dimensional space Gsta，GstaEach element in (1) is mapped once successfully, and the corresponding calculator is added with 1;
3.3) mixing GstaThe counter result for each element of (a) is divided by 9 to take the average.
In one example, the step 4) includes:
4.1) calculating the mean value of each element in the two-dimensional space of the contrast surface element based on the formula (1)
4.2) calculating the mean value of each element in the bin two-dimensional space based on the formula (2)
4.3) calculating the correlation R of the two-dimensional space of the contrast surface element and the two-dimensional space of the surface element taken in the step 3) based on the formula (3):
in one example, the program is further executable by the processor to:
6.1) selecting the maximum value R of the correlation R of the two-dimensional space of all the surface elements and the two-dimensional space of the contrast surface elementmax；
6.2) mixing of [0, Rmax]Dividing the range into L sections, traversing the correlation degree R of all surface elements, and counting the number R of the correlation degree belonging to each sectionnum[i]Wherein i is 0,1,2, …, L.
Application example
The three-dimensional seismic observation system bin relative consistency evaluation method is applied to the 34L8S240T observation systems 1,2 and 3 respectively. The parameters of the observation system are shown in the following table, and the parameters of the three systems are the same except for the distance between the gun lines.
| |
|
|
|
| Road pitch (Rice) | 50 | 50 | 50 |
| Receiving line distance (rice) | 400 | 400 | 400 |
| Distance of fire point (rice) | 50 | 50 | 50 |
| Distance between gun lines (rice) | 300 | 400 | 500 |
| Number of received lines | 34 | 34 | 34 |
| Number of single-row tracks | 240 | 240 | 240 |
| Number of times of coverage | 304 | 255 | 204 |
In order to evaluate the bin consistency of the three observation systems, the bin correlation degrees of the three observation systems are respectively calculated by using the three-dimensional seismic observation system bin relative consistency evaluation method, so that the output results shown in fig. 3a-3c, fig. 4a-4c and fig. 5a-5c are obtained.
Comparing the three systems, from system 1 to system 3, the phenomenon of collecting footprints in fig. 3a, 4a and 5a is gradually enhanced, which shows that the correlation degree is gradually reduced. Similarly, from the bar graphs of fig. 3b, 4b, and 5b and the sector graphs of fig. 3c, 4c, and 5c, it can be seen that the proportion of the number of bins with high correlation of the system 1 is also significantly higher than that of the systems 2 and 3, and comparing the systems 2 and 3, the correlation of the system 2 is higher than that of the system 3. Therefore, the observation system 1 is superior to the systems 2 and 3 in bin consistency in terms of the correlation index in this example.
In conclusion, the application example proves that the method is feasible, and the problems of evaluation and optimization of the observation system binning consistency can be solved in practical application.
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention is intended only to illustrate the benefits of embodiments of the invention and is not intended to limit embodiments of the invention to any examples given.
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The terminology used herein is chosen in order to best explain the principles of the embodiments, the practical application, or improvements made to the technology in the marketplace, or to enable others of ordinary skill in the art to understand the embodiments disclosed herein.
Claims (10)
1. A three-dimensional seismic observation system surface element relative consistency evaluation method is characterized by comprising the following steps:
1) inputting an azimuth angle-offset data pair of each surface element of the three-dimensional earthquake observation system;
2) mapping the azimuth angle-offset data pairs of each surface element into respective surface element two-dimensional space;
3) taking one surface element as a center, and mapping azimuth angle-offset data pairs of the surrounding 9 surface elements including the surface element into a two-dimensional space of a contrast surface element;
4) calculating the correlation degree by using the data of the two-dimensional space of the surface element serving as the center in the step 3) and the two-dimensional space of the contrast surface element;
5) and repeating the steps 3) and 4), traversing all surface elements of the three-dimensional observation system, and calculating the correlation of each surface element.
2. The three-dimensional seismic observation system bin relative consistency evaluation method according to claim 1, wherein the step 2) comprises:
2.1) constructing an M N two-dimensional space G for one of the binsavoAnd is GavoSetting a counter for each element, where M is 360/d theta, d theta is azimuth step, and N is Xmax/dx，XmaxFor observation systemThe unified maximum offset, dx is offset step length;
2.2) traverse all azimuth-offset data pairs of the bins, mapping each data pair to Gavo[idx,idy]Each time the mapping is successful, adding 1 to a counter of the corresponding element, wherein idx is θ/d θ, and idy is x/dx;
2.3) traversing all the bins, repeating the steps 2.1) -2.2), and mapping the azimuth angle-offset data pairs of each bin in the respective two-dimensional space.
3. The method of claim 2, wherein the step 3) of mapping the azimuth-offset pairs of the surrounding 9 bins including the bin into a two-dimensional space of a contrast bin comprises:
3.1) construction of an MxN two-dimensional space G of contrast binsstaAnd is GstaSetting a counter for each element;
3.2) mapping the azimuth-offset data pair of 9 bins to the contrast bin two-dimensional space Gsta，GstaEach element in (1) is mapped once successfully, and the corresponding calculator is added with 1;
3.3) mixing GstaThe counter result for each element of (a) is divided by 9 to take the average.
4. The three-dimensional seismic observation system bin relative consistency evaluation method according to claim 3, wherein the step 4) comprises:
4.1) calculating the mean value of each element in the two-dimensional space of the contrast surface element based on the formula (1)
4.2) calculating the surface based on equation (2)Mean of individual elements in a binary space
4.3) calculating the correlation R of the two-dimensional space of the contrast surface element and the two-dimensional space of the surface element taken in the step 3) based on the formula (3):
5. the three-dimensional seismic observation system bin relative consistency evaluation method according to claim 1, further comprising a step 6), specifically comprising:
6.1) selecting the maximum value R of the correlation R of the two-dimensional space of all the surface elements and the two-dimensional space of the contrast surface elementmax；
6.2) mixing of [0, Rmax]Dividing the range into L sections, traversing the correlation degree R of all surface elements, and counting the number R of the correlation degree belonging to each sectionnum[i]Wherein i is 0,1,2, …, L.
6. A three-dimensional seismic observation system bin relative consistency evaluation system, the system comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the following steps when executing the program:
1) inputting an azimuth angle-offset data pair of each surface element of the three-dimensional earthquake observation system;
2) mapping the azimuth angle-offset data pairs of each surface element into respective surface element two-dimensional space;
3) taking one surface element as a center, and mapping azimuth angle-offset data pairs of the surrounding 9 surface elements including the surface element into a two-dimensional space of a contrast surface element;
4) calculating the correlation degree by using the data of the two-dimensional space of the surface element serving as the center in the step 3) and the two-dimensional space of the contrast surface element;
5) and repeating the steps 3) and 4), traversing all surface elements of the three-dimensional observation system, and calculating the correlation of each surface element.
7. The three-dimensional seismic observation system bin relative consistency evaluation system according to claim 6, wherein the step 2) comprises:
2.1) constructing an M N two-dimensional space G for one of the binsavoAnd is GavoSetting a counter for each element, where M is 360/d theta, d theta is azimuth step, and N is Xmax/dx，XmaxThe maximum offset of the observation system is adopted, and dx is offset step length;
2.2) traverse all azimuth-offset data pairs of the bins, mapping each data pair to Gavo[idx,idy]Each time the mapping is successful, adding 1 to a counter of the corresponding element, wherein idx is θ/d θ, and idy is x/dx;
2.3) traversing all the bins, repeating the steps 2.1) -2.2), and mapping the azimuth angle-offset data pairs of each bin in the respective two-dimensional space.
8. A three-dimensional seismic observation system bin relative consistency evaluation system according to claim 7, wherein the mapping of azimuth-offset data pairs of the surrounding 9 bins including the bin into a contrast bin two-dimensional space in step 3) comprises:
3.1) construction of an MxN two-dimensional space G of contrast binsstaAnd is GstaSetting a counter for each element;
3.2) mapping the azimuth-offset data pair of 9 bins to the contrast bin two-dimensional space Gsta，GstaEach element in (1) is mapped once successfully, and the corresponding calculator is added with 1;
3.3) mixing GstaThe counter result for each element of (a) is divided by 9 to take the average.
9. The three-dimensional seismic observation system bin relative consistency evaluation system according to claim 8, wherein the step 4) comprises:
4.1) calculating the mean value of each element in the two-dimensional space of the contrast surface element based on the formula (1)
4.2) calculating the mean value of each element in the bin two-dimensional space based on the formula (2)
4.3) calculating the correlation R of the two-dimensional space of the contrast surface element and the two-dimensional space of the surface element taken in the step 3) based on the formula (3):
10. the three-dimensional seismic observation system bin relative consistency evaluation system according to claim 6, wherein the processor, when executing the program, further performs the steps of:
6.1) selecting the maximum value R of the correlation R of the two-dimensional space of all the surface elements and the two-dimensional space of the contrast surface elementmax；
6.2) mixing of [0, Rmax]The range is divided into L sections and traversedThe correlation degree R of the surface element is counted, and the number R of the correlation degrees belonging to each equal division is countednum[i]Wherein i is 0,1,2, …, L.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811168459.7A CN111007580B (en) | 2018-10-08 | 2018-10-08 | Method and system for evaluating relative consistency of surface elements of three-dimensional seismic observation system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811168459.7A CN111007580B (en) | 2018-10-08 | 2018-10-08 | Method and system for evaluating relative consistency of surface elements of three-dimensional seismic observation system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111007580A CN111007580A (en) | 2020-04-14 |
| CN111007580B true CN111007580B (en) | 2021-11-05 |
Family
ID=70111094
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811168459.7A Active CN111007580B (en) | 2018-10-08 | 2018-10-08 | Method and system for evaluating relative consistency of surface elements of three-dimensional seismic observation system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111007580B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2308993A1 (en) * | 1999-05-25 | 2000-11-25 | Baker Hughes Incorporated | Method of determining merit of offset and azimuth distributions in seismic data acquisition |
| CN104142514A (en) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | Three-dimensional earthquake observing system quantitative designing method |
| CN105044799A (en) * | 2015-07-10 | 2015-11-11 | 中国石油大学(北京) | Method for determining bin attribute uniformity degree and uniformization of three-dimensional earthquake observation system |
| CN108169818A (en) * | 2017-12-15 | 2018-06-15 | 西南石油大学 | A kind of novel observation system quality evaluating method |
-
2018
- 2018-10-08 CN CN201811168459.7A patent/CN111007580B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2308993A1 (en) * | 1999-05-25 | 2000-11-25 | Baker Hughes Incorporated | Method of determining merit of offset and azimuth distributions in seismic data acquisition |
| CN104142514A (en) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | Three-dimensional earthquake observing system quantitative designing method |
| CN105044799A (en) * | 2015-07-10 | 2015-11-11 | 中国石油大学(北京) | Method for determining bin attribute uniformity degree and uniformization of three-dimensional earthquake observation system |
| CN108169818A (en) * | 2017-12-15 | 2018-06-15 | 西南石油大学 | A kind of novel observation system quality evaluating method |
Non-Patent Citations (2)
| Title |
|---|
| 三维观测系统采样均匀性分析方法的研究;碗学俭，葛雪钦，张松平等;《中国石油勘探》;20091215;第14卷(第6期);1、42-48 * |
| 基于面元属性分析的三维地震观测系统设计;刘二鹏，李瑞强，李宇;《山西煤炭》;20110531;第31卷(第5期);54-57 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111007580A (en) | 2020-04-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109752715B (en) | SAR data total-dispersion body detection method and device | |
| US9366771B2 (en) | System and method for migration velocity modeling | |
| CN108318921A (en) | A kind of quick earthquake stochastic inversion methods based on lateral confinement | |
| WO2012177335A1 (en) | Improved dispersion estimation by nonlinear optimization of beam-formed fields | |
| CN104657955A (en) | Displacement field iteration smoothing method of kernel function based digital image correlation method | |
| CN111597753B (en) | Data depth change characteristic self-adaptive two-dimensional resistivity inversion method and system | |
| CN111007580B (en) | Method and system for evaluating relative consistency of surface elements of three-dimensional seismic observation system | |
| CN104316961B (en) | Method for obtaining geological parameters of weathered layer | |
| CN111007581B (en) | Method and system for evaluating bin uniformity of three-dimensional seismic observation system | |
| CN105277974A (en) | Stratum data interpolation method | |
| CN109507750B (en) | Method and system for evaluating bin consistency of three-dimensional seismic observation system | |
| CN116449305A (en) | Dense time-varying array construction method and system based on controllable variation self-encoder | |
| CN104656126B (en) | A kind of method generating seismic observation system | |
| CN104570062B (en) | A kind of VSP observation system design methods centered on exciting | |
| JP6960834B2 (en) | Ground survey position determination method, determination device, ground estimation method and ground estimation device | |
| Behrens et al. | Incorporating seismic data of intermediate vertical resolution into 3D reservoir models | |
| CN112630831A (en) | Method and system for calculating longitudinal dimension of carbonate karst cave | |
| Huang* et al. | Non-local means filter for trim statics | |
| Behrens et al. | Incorporating seismic data of intermediate vertical resolution into three-dimensional reservoir models: a new method | |
| Demirhan et al. | Performance evaluation of spatial interpolation methods in the presence of noise | |
| CN105844086A (en) | Arrangement structure of pressure sensor | |
| CN111538078A (en) | Observation mode determination method and device of two-dimensional wide-line seismic observation system | |
| CN113805233A (en) | Calculation method of point spread function | |
| CN112835122B (en) | Discontinuous three-dimensional joint inversion method based on smooth focusing regularization | |
| CN110308428A (en) | A method of establishing low latitude clutter distributed simulation model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |