CN112462422B - Automatic design method and device for three-dimensional seismic exploration physical point positions in desert area - Google Patents
Automatic design method and device for three-dimensional seismic exploration physical point positions in desert area Download PDFInfo
- Publication number
- CN112462422B CN112462422B CN202011202469.5A CN202011202469A CN112462422B CN 112462422 B CN112462422 B CN 112462422B CN 202011202469 A CN202011202469 A CN 202011202469A CN 112462422 B CN112462422 B CN 112462422B
- Authority
- CN
- China
- Prior art keywords
- point
- slope
- determining
- points
- sample
- 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
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000013461 design Methods 0.000 title claims description 16
- 239000004576 sand Substances 0.000 claims abstract description 65
- 238000010276 construction Methods 0.000 claims abstract description 30
- 238000004590 computer program Methods 0.000 claims description 17
- 238000001514 detection method Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 11
- 238000012876 topography Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 5
- 230000002349 favourable effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000010187 selection method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000005056 compaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000007664 blowing Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008439 repair process Effects 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
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
The invention discloses a method and a device for automatically designing physical point positions of three-dimensional seismic exploration in a desert area, wherein the method comprises the following steps: acquiring elevation data of a construction area and coordinates of preset theoretical points; for each theoretical point, determining a three-dimensional seismic exploration physical point position corresponding to the theoretical point according to the following method: determining the coordinates of all sample points in the selected point range of the theoretical point according to the coordinates of the theoretical point; determining the slope direction, the slope data value, the slope type and the elevation of the sand dunes at each sample point according to the elevation data; calculating the contribution rate of each sample point according to the coordinates and the elevation of each sample point; and selecting the sample points with contribution rates, slope directions and slope types meeting preset conditions from all the sample points as three-dimensional seismic exploration physical points. The invention can reasonably and efficiently select the physical point location of the three-dimensional seismic exploration.
Description
Technical Field
The invention relates to the technical field of land petroleum seismic exploration, in particular to an automatic design method and device for three-dimensional seismic exploration physical points in a desert area.
Background
This section is intended to provide a background or context to the embodiments of the invention that are recited in the claims. The description herein is not admitted to be prior art by inclusion in this section.
In the process of acquiring three-dimensional seismic data in a desert area, factors objectively existing such as sand dune thickness, earth surface compaction degree, topography fluctuation degree and the like directly influence the seismic data acquisition quality and construction difficulty, so that the physical point position which is theoretically is properly optimized within the technical allowable range in the measurement procedure.
At present, the three-dimensional seismic data acquisition physical point location of the mobile desert area preferably mainly adopts a manual point selection method. The manual point selection method is a method for selecting the position of a physical point through subjective judgment according to construction experience after observing topography and topography near a point of theoretical design. Because the motor equipment is constructed in the mobile desert area and needs the bulldozer to repair the road, the manual point selection implementation process is roughly divided into two steps: firstly, roughly, namely, measuring a preferable travelling route of mobile equipment by a worker, observing topography and topography in a visual field range after the technician responsible for selecting points reaches a theoretical point, limiting a route which is easy to pass by the mobile equipment in a technical allowable range, and carrying out road pushing work along a route appointed by the technician by a bulldozer operator; secondly, fine placement is performed, and based on rough placement, a three-dimensional seismic exploration physical point position which is beneficial to reducing construction difficulty and obtaining high-quality seismic data is further optimized in a technical allowable offset range aiming at each theoretical point.
The method mainly has two problems: first, the subjective factors of the person are greatly affected. The construction route and the physical point position quality are determined by subjective judgment of a technical worker, and the point selection effect is closely related to the construction experience and responsibility of the technical worker; second, work efficiency is low. In the vast sand sea, due to the limitation of the visual field range, the judgment and evaluation of the technical workers on the whole topography are affected, and the situation that a rough route and a fine point are reselected often occurs.
In seismic exploration in a mobile desert area, the data quality is closely related to the earth surface condition, and the selection of physical point positions needs to consider the needs of reducing the construction difficulty, improving the seismic data quality, reducing the offset, improving the construction efficiency and the like, so that the current point selection method cannot completely meet the production needs, and a reasonable and efficient physical point position optimization method is necessary to be provided.
Disclosure of Invention
The embodiment of the invention provides an automatic design method for three-dimensional seismic exploration physical points in a desert area, which is used for reasonably and efficiently selecting the three-dimensional seismic exploration physical points, and comprises the following steps:
acquiring elevation data of a construction area and coordinates of preset theoretical points;
for each theoretical point, determining a three-dimensional seismic exploration physical point position corresponding to the theoretical point according to the following method:
Determining the coordinates of all sample points in the selected point range of the theoretical point according to the coordinates of the theoretical point; the selection range of the shot points is set to be a rectangular range with the position of the theoretical shot point as the center InLine m in the direction of 25m and the CrossLine direction of 120 m; the selection range of the detection point is set to be a rectangular range with the detection point theoretical point position as a center and the direction of InLine m and the direction of cross line of 25 m;
Determining the slope direction, the slope data value, the slope type and the elevation of the sand dunes at each sample point according to the elevation data;
calculating the contribution rate of each sample point according to the coordinates and the elevation of each sample point;
Selecting sample points with contribution rate, slope direction and slope position type meeting preset conditions from all the sample points as three-dimensional seismic exploration physical points;
After selecting the sample points with contribution rates, slope directions and slope types meeting preset conditions from all the sample points as three-dimensional seismic exploration physical points, the method further comprises the following steps: inserting virtual shots between every two adjacent three-dimensional seismic exploration physical points by using an interpolation method; taking the virtual shot points as theoretical points, and determining virtual three-dimensional seismic exploration physical points corresponding to each virtual shot point; performing curve fitting by utilizing the three-dimensional seismic exploration physical points and the virtual three-dimensional seismic exploration physical points, and taking the fitted curve as a construction route of the controllable seismic source;
The maximum offset distance of the virtual shot point in the InLine direction is allowed to be 2.5m, the maximum offset distance in the CrossLine direction is allowed to be 30m, and the selection point range of the virtual shot point is set to be a rectangular range of 5m in the direction InLine and 60m in the CrossLine direction with the virtual shot point position as the center.
The embodiment of the invention also provides an automatic design device for the three-dimensional seismic exploration physical point positions of the desert area, which is used for reasonably and efficiently selecting the three-dimensional seismic exploration physical point positions and comprises the following steps:
the acquisition module is used for acquiring elevation data of the construction area and coordinates of preset theoretical points;
For each theoretical point, the point position determining module is used for determining a three-dimensional seismic exploration physical point position corresponding to the theoretical point according to the following method:
Determining the coordinates of all sample points in the selected point range of the theoretical point according to the coordinates of the theoretical point; the selection range of the shot points is set to be a rectangular range with the position of the theoretical shot point as the center InLine m in the direction of 25m and the CrossLine direction of 120 m; the selection range of the detection point is set to be a rectangular range with the detection point theoretical point position as a center and the direction of InLine m and the direction of cross line of 25 m;
Determining the slope direction, the slope data value, the slope type and the elevation of the sand dunes at each sample point according to the elevation data;
calculating the contribution rate of each sample point according to the coordinates and the elevation of each sample point;
Selecting sample points with contribution rate, slope direction and slope position type meeting preset conditions from all the sample points as three-dimensional seismic exploration physical points;
the interpolation module is used for inserting virtual shot points between every two adjacent three-dimensional seismic exploration physical points by using an interpolation method;
The point position determining module is also used for determining virtual three-dimensional seismic exploration physical points corresponding to each virtual shot point by taking the virtual shot point as a theoretical point;
The curve determining module is used for performing curve fitting by utilizing the three-dimensional seismic exploration physical points and the virtual three-dimensional seismic exploration physical points, and taking the fitting curve as a construction route of the controllable seismic source;
The maximum offset distance of the virtual shot point in the InLine direction is allowed to be 2.5m, the maximum offset distance in the CrossLine direction is allowed to be 30m, and the selection point range of the virtual shot point is set to be a rectangular range of 5m in the direction InLine and 60m in the CrossLine direction with the virtual shot point position as the center.
The embodiment of the invention also provides computer equipment, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the computer program is executed by the processor to realize the method for automatically designing the physical point location of the three-dimensional seismic exploration in the desert area.
The embodiment of the invention also provides a computer readable storage medium which stores a computer program, and the computer program realizes the automatic design method of the three-dimensional seismic exploration physical point location of the desert area when being executed by a processor.
In the embodiment of the invention, the slope direction, the slope data value, the slope type and the elevation of the sand dune at each sample point are obtained by processing the elevation data; after the data are further processed, digitized data, such as contribution rate, slope direction, slope position type and the like, which accord with the characteristics of the mobile desert area are obtained, the physical point positions of the three-dimensional seismic exploration of the area are selected by utilizing the digitized data, and the physical point positions are selected at positions which are most favorable for guaranteeing the quality of seismic data and reducing the construction difficulty, so that the requirements of reasonable, accurate and efficient optimization of the physical point positions in the three-dimensional seismic exploration of the mobile desert area are met.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art. In the drawings:
FIG. 1 is a flow chart of an automated design method for physical points of three-dimensional seismic exploration in a desert area according to an embodiment of the invention;
FIG. 2 (a) is a plan view of altitude data obtained by aerial photography in an embodiment of the present invention;
FIG. 2 (b) is a partial 3D elevation view obtained by aerial photography in an embodiment of the present invention;
FIG. 3 is a schematic view of a theoretical point selection range according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of a slope data body according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of a selection point within an offset range according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of a device for automatically designing physical points of three-dimensional seismic exploration in a desert area according to an embodiment of the invention;
fig. 7 is a schematic structural diagram of a computer device according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention will be described in further detail with reference to the accompanying drawings. The exemplary embodiments of the present invention and their descriptions herein are for the purpose of explaining the present invention, but are not to be construed as limiting the invention.
The purpose of the three-dimensional seismic exploration physical point optimization work in the flow desert area comprises the following three aspects: 1. the quality of the seismic data is improved; 2. the construction difficulty is reduced; 3. and the point selection efficiency is improved.
In order to achieve the above object, the following point selection principle is generally adopted:
1. Preferred windward side
The windward side and the leeward side are special micro-topography of the flowing desert area. The windward side is influenced by the continuous blowing of the monsoon, and the relative compaction degree of the ground surface is relatively high; the sand on the surface of the leeward surface is derived from the fine sand particles blown by the monsoon on the top of the dune, the sand in the area is soft, and the slope of the terrain is larger due to the accumulation effect; past seismic data confirm that the windward side is the most preferred point area in the mobile desert area.
2. Selecting a gentle slope region
In consideration of the need to reduce construction difficulty, the physical points, particularly the excitation points, should be selected in a gentle region where the mobile device is easily in place.
3. Preferably low altitude areas
In the past, the exploration data prove that the thickness of a sand dune is one of key factors influencing the quality of seismic data in a desert area, and therefore, physical points are selected at positions with relatively low altitudes based on the requirements for improving the quality of the seismic data.
4. Ensuring uniformity of three-dimensional observation system attributes
When the theoretical design point position (theoretical point) is located on the windward side, the altitude of the physical point after the offset and the theoretical point is Cheng Chaxiao m (threshold value) and the theoretical point position is kept unchanged.
Based on the selection principle, the embodiment of the invention provides an automatic design method for three-dimensional seismic exploration physical points in a desert area, as shown in fig. 1, comprising the following steps:
And 101, acquiring elevation data of a construction area and coordinates of preset theoretical points.
Because the change of the topography and the landform of the mobile desert area is large, aerial photographing data acquisition is carried out before physical point construction is carried out, and the acquired aerial photographing data comprise elevation data, image data and the like.
For example, in a large desert area of a Tarim basin, aerial photographing data are acquired by aerial photographing, and the scale of the acquired aerial photographing data is 1:2000, the plane resolution is 0.2m, the elevation plane sampling interval is 2m, and the imaging control point density is 1km multiplied by 1km to 1.5km multiplied by 1.5km. Referring to fig. 2 (a) and 2 (b), an aerial acquired elevation data plan and a local 3D elevation map are respectively.
Step 102, determining three-dimensional seismic exploration physical points corresponding to the theoretical points according to the methods in the steps 103 to 106 for each theoretical point.
And 103, determining the coordinates of all the sample points in the selected point range of the theoretical point according to the coordinates of the theoretical point.
The coordinate system used in the embodiment of the invention is a self-built coordinate system after the user selects the original point, the coordinate axis direction and the like, and the coordinates of all the sample points in the range of the selected point of the theoretical point are determined under the coordinate system.
In engineering design, the maximum offset distance of the shot point in the transverse line (InLine) direction is allowed to be 12.5m, and the maximum offset distance in the longitudinal line (CrossLine) direction is allowed to be 60m, so that the shot point selection range of the shot point is set to be a rectangular range of 25m in the InLine direction and 120m in the CrossLine direction by taking the shot point position as the center; the maximum offset distance of the allowable detector point in InLine direction is 2.5m, the maximum offset distance in the CrossLine direction is 12.5m, and the selection point range is set as a rectangular range of 5m in InLine direction and 25m in CrossLine direction with the theoretical position of the detector point as the center. As shown in fig. 3, circles on the horizontal measurement lines in fig. 3 are theoretical points of the detectors, dots on the vertical measurement lines are theoretical points of the detectors, one theoretical point of the detector and one theoretical point of the detector are taken as examples, a range of the detector offset is shown in a dotted rectangle in fig. 3, and a better detector and detector are selected from the range of the detector offset and the range of the detector offset.
And 104, determining the slope direction, the gradient data value, the slope position type and the elevation of the sand dune at each sample point according to the elevation data.
Determining the downhill direction and the sand hill trend of the sand hill by using the elevation data, and determining the included angle between the direction with the largest elevation change rate from each pixel to the adjacent pixels in each downhill direction of the sand hill and the designated direction as slope data; determining a reference angle range of each slope according to the sand dune trend, wherein the slope comprises a windward side and a leeward side; determining a windward slope and a leeward slope of the sand dune according to the reference angle range and the slope direction data; and determining the slope direction of the sand dune at each sample point according to the coordinates of each sample point.
The value of each pixel in the grid in the elevation data may indicate the compass direction of the orientation of the surface at the location of each pixel.
In the embodiment of the invention, the downhill direction of the sand hill is radiated from the top of the sand hill to the periphery, so that the included angle between the downhill direction of the sand hill and the designated direction (such as northbound direction, southward direction and the like) is also between 0 and 360 degrees. The north direction is taken as the appointed direction, the measurement is carried out according to the clockwise or anticlockwise direction, the included angle between the initial direction and the north direction is 0 degree, the measurement returns to the north direction after one circle, and the angle is 360 degrees, so that the circle is a complete circle. And taking the measured included angle as slope data.
Each slope has a certain angle range, taking a desert area in a tower as an example, the desert area in the tower is influenced by a monsoon, the sand dunes are relatively strong in regularity, the reference angle range of the leeward surface is distributed between 110 degrees and 270 degrees, the rest areas are windward surface ranges, different values are assigned to different slopes according to the angle ranges, for example, the leeward surface is assigned to 1, the windward surface is assigned to 2, the flat area without downhill direction is assigned to-1, slope data are digitized, and a slope data body shown in fig. 4 can be obtained.
In one implementation, determining the slope type of the dune at each sample point according to the elevation data may be performed as follows: generating a terrain location index (Topographic Position Index, TPI) data volume using the elevation data; generating a terrain slope classification (Slope Position Classification, SPC) data body of the construction area according to the elevation data and the TPI data body; according to the SPC data body and the coordinates of each sample point, determining the slope type of the sand dune at each sample point; the slope types of the sand dunes comprise sand valleys, slope feet, flat areas, medium slope areas, steep slope areas and sand ridges.
The TPI data body refers to that under the existing DEM data, the TPI raster data body is calculated and generated by utilizing software with a TPI calculation function according to specified grid intervals (for example, 5×5) in a gridding way. The SPC data body is characterized in that the software with the terrain slope classification calculation function is utilized, the input DEM data body and the TPI data body are used as the basis, the terrain is divided into Sha Gu types, slope toe types, flat areas, medium slope areas, steep slope areas and sand ridges 6 types, the types are respectively assigned to be 1,2, 3, 4, 5 and 6, and the SPC grid data body is generated.
And 105, calculating the contribution rate of each sample point according to the coordinates and the elevation of each sample point.
Specifically, determining the height difference between each sample point and the theoretical point according to the elevation; determining the plane distance between each sample point and the theoretical point according to the coordinates of each sample point; and determining the ratio of the height difference to the plane distance as the contribution rate of the sample point.
And 106, selecting the sample points with contribution rates, slope directions and slope types meeting preset conditions from all the sample points as three-dimensional seismic exploration physical points.
Specifically, selecting samples with contribution rate greater than 10% from all samples, wherein the slope direction is windward, and the slope type is not the samples of the sand valley, the slope toe and the steep slope area, and the samples are used as initial selected samples; and determining the sample point with the highest contribution rate in the initially selected sample points as a physical point of the three-dimensional seismic exploration. For example, referring to fig. 5, the theoretical point is located on the leeward side, and is not suitable for being used as an actual exploration point, the method provided by the embodiment of the invention is used for selecting points in the offset range of the theoretical point, the selected points are used as three-dimensional seismic exploration physical points, the selected points are located on the windward side, and the condition is suitable for actual exploration.
For each theoretical point, the steps 102 to 106 are repeatedly executed, and the three-dimensional seismic exploration physical point position corresponding to the theoretical point can be determined in the offset range of all the theoretical points.
Because the distance between the three-dimensional seismic exploration physical points is longer, and the route from one physical point to another physical point is more, in the embodiment of the invention, after the sample points with contribution rates, slope directions and slope types meeting the preset conditions are selected from all the sample points to serve as the three-dimensional seismic exploration physical points, a virtual shot point can be inserted between every two adjacent three-dimensional seismic exploration physical points by using an interpolation method; taking the virtual shot points as theoretical points, and determining virtual three-dimensional seismic exploration physical points corresponding to each virtual shot point; and performing curve fitting by utilizing the three-dimensional seismic exploration physical points and the virtual three-dimensional seismic exploration physical points, and taking the fitted curve as a construction route of the controllable seismic source.
It should be noted that, the distance between the virtual shot point and the virtual shot point is short, for example, the shot point distance is set to 5 meters. Because the shot size is small, the dot selection range of the corresponding virtual shot is also set to a small value, such as allowing the maximum offset distance of the virtual shot in the InLine direction to be 2.5m and the maximum offset distance of the virtual shot in the CrossLine direction to be 30m, and thus the dot selection range of the virtual shot is set to a rectangular range of 5m in the InLine direction and 60m in the CrossLine direction with the virtual shot position as the center.
In the embodiment of the invention, the slope direction, the slope data value, the slope type and the elevation of the sand dune at each sample point are obtained by processing the elevation data; after the data are further processed, digitized data, such as contribution rate, slope direction, slope position type and the like, which accord with the characteristics of the mobile desert area are obtained, the physical point positions of the three-dimensional seismic exploration of the area are selected by utilizing the digitized data, and the physical point positions are selected at positions which are most favorable for guaranteeing the quality of seismic data and reducing the construction difficulty, so that the requirements of reasonable, accurate and efficient optimization of the physical point positions in the three-dimensional seismic exploration of the mobile desert area are met.
The embodiment of the invention also provides an automatic design device for the physical point location of the three-dimensional seismic exploration in the desert area, as described in the following embodiment. Because the principle of the device for solving the problems is similar to that of the automatic design method of the three-dimensional seismic exploration physical point positions in the desert area, the implementation of the device can be referred to the implementation of the automatic design method of the three-dimensional seismic exploration physical point positions in the desert area, and repeated parts are not repeated.
As shown in fig. 6, the apparatus 600 includes an acquisition module 601 and a point location determination module 602.
The acquiring module 601 is configured to acquire elevation data of a construction area and coordinates of a preset theoretical point;
for each theoretical point, the point location determining module 602 is configured to determine a three-dimensional seismic exploration physical point location corresponding to the theoretical point according to the following method:
determining the coordinates of all sample points in the selected point range of the theoretical point according to the coordinates of the theoretical point;
Determining the slope direction, the slope data value, the slope type and the elevation of the sand dunes at each sample point according to the elevation data;
calculating the contribution rate of each sample point according to the coordinates and the elevation of each sample point;
And selecting the sample points with contribution rates, slope directions and slope types meeting preset conditions from all the sample points as three-dimensional seismic exploration physical points.
In one implementation of the embodiment of the present invention, the point location determining module 602 is configured to:
determining the downhill direction and the sand hill trend of the sand hill by utilizing the elevation data, and determining the included angle between the direction with the largest elevation change rate from each pixel to the adjacent pixels in each downhill direction of the sand hill and the designated direction as slope data;
Determining a reference angle range of each slope according to the sand dune trend, wherein the slope comprises a windward side and a leeward side;
Determining a windward slope and a leeward slope of the sand dune according to the reference angle range and the slope direction data;
and determining the slope direction of the sand dune at each sample point according to the coordinates of each sample point.
In one implementation of the embodiment of the present invention, the point location determining module 602 is configured to:
generating a Terrain Position Index (TPI) data body by using the elevation data;
generating a terrain slope SPC data body of the construction area according to the elevation data and the TPI data body;
according to the SPC data body and the coordinates of each sample point, determining the slope type of the sand dune at each sample point; the slope types of the sand dunes comprise sand valleys, slope feet, flat areas, medium slope areas, steep slope areas and sand ridges.
In one implementation of the embodiment of the present invention, the point location determining module 602 is configured to:
determining the height difference between each sample point and the theoretical point according to the elevation;
determining the plane distance between each sample point and the theoretical point according to the coordinates of each sample point;
And determining the ratio of the height difference to the plane distance as the contribution rate of the sample point.
In one implementation of the embodiment of the present invention, the point location determining module 602 is configured to:
selecting samples with contribution rate greater than 10% from all samples, slope direction being windward surface, and slope type being not sand valley, slope toe and steep slope area as primary selected samples;
and determining the sample point with the highest contribution rate in the initially selected sample points as a physical point of the three-dimensional seismic exploration.
In one implementation of an embodiment of the present invention, the apparatus 600 further includes:
the interpolation module 603 is configured to insert virtual shots between two adjacent three-dimensional seismic exploration physical points by using an interpolation method;
The point location determining module 602 is further configured to determine a virtual three-dimensional seismic exploration physical point location corresponding to each virtual shot point by using the virtual shot point as a theoretical point;
The curve determining module 604 is configured to perform curve fitting by using the three-dimensional seismic exploration physical point location and the virtual three-dimensional seismic exploration physical point location, and use the fitted curve as a construction route of the controllable seismic source.
In the embodiment of the invention, the slope direction, the slope data value, the slope type and the elevation of the sand dune at each sample point are obtained by processing the elevation data; after the data are further processed, digitized data, such as contribution rate, slope direction, slope position type and the like, which accord with the characteristics of the mobile desert area are obtained, the physical point positions of the three-dimensional seismic exploration of the area are selected by utilizing the digitized data, and the physical point positions are selected at positions which are most favorable for guaranteeing the quality of seismic data and reducing the construction difficulty, so that the requirements of reasonable, accurate and efficient optimization of the physical point positions in the three-dimensional seismic exploration of the mobile desert area are met.
The embodiment of the invention also provides a computer device, and fig. 7 is a schematic diagram of the computer device in the embodiment of the invention, where the computer device can implement all the steps in the method for automatically designing the three-dimensional seismic prospecting physical point location in the desert area, and the computer device specifically includes the following contents:
a processor (processor) 701, a memory (memory) 702, a communication interface (Communications Interface) 703, and a communication bus 704;
Wherein, the processor 701, the memory 702 and the communication interface 703 complete communication with each other through the communication bus 704; the communication interface 703 is used for implementing information transmission between related devices;
The processor 701 is configured to invoke a computer program in the memory 702, where the processor executes the computer program to implement the method for automatically designing a physical point location for three-dimensional seismic exploration in a desert area in the foregoing embodiment.
The embodiment of the invention also provides a computer readable storage medium which stores a computer program for executing the method for automatically designing the three-dimensional seismic exploration physical point location of the desert area.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the invention, and is not meant to limit the scope of the invention, but to limit the invention to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the invention are intended to be included within the scope of the invention.
Claims (12)
1. The automatic design method for the three-dimensional seismic exploration physical point positions of the desert area is characterized by comprising the following steps:
acquiring elevation data of a construction area and coordinates of preset theoretical points;
for each theoretical point, determining a three-dimensional seismic exploration physical point position corresponding to the theoretical point according to the following method:
Determining the coordinates of all sample points in the selected point range of the theoretical point according to the coordinates of the theoretical point; the selection range of the shot points is set to be a rectangular range with the position of the theoretical shot point as the center InLine m in the direction of 25m and the CrossLine direction of 120 m; the selection range of the detection point is set to be a rectangular range with the detection point theoretical point position as a center and the direction of InLine m and the direction of cross line of 25 m;
Determining the slope direction, the slope data value, the slope type and the elevation of the sand dunes at each sample point according to the elevation data;
calculating the contribution rate of each sample point according to the coordinates and the elevation of each sample point;
Selecting sample points with contribution rate, slope direction and slope position type meeting preset conditions from all the sample points as three-dimensional seismic exploration physical points;
After selecting the sample points with contribution rates, slope directions and slope types meeting preset conditions from all the sample points as three-dimensional seismic exploration physical points, the method further comprises the following steps: inserting virtual shots between every two adjacent three-dimensional seismic exploration physical points by using an interpolation method; taking the virtual shot points as theoretical points, and determining virtual three-dimensional seismic exploration physical points corresponding to each virtual shot point; performing curve fitting by utilizing the three-dimensional seismic exploration physical points and the virtual three-dimensional seismic exploration physical points, and taking the fitted curve as a construction route of the controllable seismic source;
The maximum offset distance of the virtual shot point in the InLine direction is allowed to be 2.5m, the maximum offset distance in the CrossLine direction is allowed to be 30m, and the selection point range of the virtual shot point is set to be a rectangular range of 5m in the direction InLine and 60m in the CrossLine direction with the virtual shot point position as the center.
2. The method of claim 1, wherein determining the slope of the dune at each sample point based on elevation data comprises:
determining the downhill direction and the sand hill trend of the sand hill by utilizing the elevation data, and determining the included angle between the direction with the largest elevation change rate from each pixel to the adjacent pixels in each downhill direction of the sand hill and the designated direction as slope data;
Determining a reference angle range of each slope according to the sand dune trend, wherein the slope comprises a windward side and a leeward side;
Determining a windward slope and a leeward slope of the sand dune according to the reference angle range and the slope direction data;
and determining the slope direction of the sand dune at each sample point according to the coordinates of each sample point.
3. The method of claim 1, wherein determining the slope type of the dune at each sample point based on the elevation data comprises:
generating a Terrain Position Index (TPI) data body by using the elevation data;
generating a terrain slope SPC data body of the construction area according to the elevation data and the TPI data body;
according to the SPC data body and the coordinates of each sample point, determining the slope type of the sand dune at each sample point; the slope types of the sand dunes comprise sand valleys, slope feet, flat areas, medium slope areas, steep slope areas and sand ridges.
4. The method of claim 1, wherein calculating the contribution rate of each sample based on the coordinates and elevation of each sample comprises:
determining the height difference between each sample point and the theoretical point according to the elevation;
determining the plane distance between each sample point and the theoretical point according to the coordinates of each sample point;
and determining the ratio of the height difference to the plane distance as the contribution rate of the sample point.
5. The method according to any one of claims 1 to 4, wherein selecting, as the three-dimensional seismic prospecting physical points, the points whose contribution rate, slope direction and slope position type meet the preset conditions from all the points, comprises:
selecting samples with contribution rate greater than 10% from all samples, slope direction being windward surface, and slope type being not sand valley, slope toe and steep slope area as primary selected samples;
and determining the sample point with the highest contribution rate in the initially selected sample points as a physical point of the three-dimensional seismic exploration.
6. An automated physical point location design device for three-dimensional seismic exploration in a desert area, which is characterized by comprising:
the acquisition module is used for acquiring elevation data of the construction area and coordinates of preset theoretical points;
For each theoretical point, the point position determining module is used for determining a three-dimensional seismic exploration physical point position corresponding to the theoretical point according to the following method:
Determining the coordinates of all sample points in the selected point range of the theoretical point according to the coordinates of the theoretical point; the selection range of the shot points is set to be a rectangular range with the position of the theoretical shot point as the center InLine m in the direction of 25m and the CrossLine direction of 120 m; the selection range of the detection point is set to be a rectangular range with the detection point theoretical point position as a center and the direction of InLine m and the direction of cross line of 25 m;
Determining the slope direction, the slope data value, the slope type and the elevation of the sand dunes at each sample point according to the elevation data;
calculating the contribution rate of each sample point according to the coordinates and the elevation of each sample point;
Selecting sample points with contribution rate, slope direction and slope position type meeting preset conditions from all the sample points as three-dimensional seismic exploration physical points;
the interpolation module is used for inserting virtual shot points between every two adjacent three-dimensional seismic exploration physical points by using an interpolation method;
The point position determining module is also used for determining virtual three-dimensional seismic exploration physical points corresponding to each virtual shot point by taking the virtual shot point as a theoretical point;
The curve determining module is used for performing curve fitting by utilizing the three-dimensional seismic exploration physical points and the virtual three-dimensional seismic exploration physical points, and taking the fitting curve as a construction route of the controllable seismic source;
The maximum offset distance of the virtual shot point in the InLine direction is allowed to be 2.5m, the maximum offset distance in the CrossLine direction is allowed to be 30m, and the selection point range of the virtual shot point is set to be a rectangular range of 5m in the direction InLine and 60m in the CrossLine direction with the virtual shot point position as the center.
7. The apparatus of claim 6, wherein the point location determination module is configured to:
determining the downhill direction and the sand hill trend of the sand hill by utilizing the elevation data, and determining the included angle between the direction with the largest elevation change rate from each pixel to the adjacent pixels in each downhill direction of the sand hill and the designated direction as slope data;
Determining a reference angle range of each slope according to the sand dune trend, wherein the slope comprises a windward side and a leeward side;
Determining a windward slope and a leeward slope of the sand dune according to the reference angle range and the slope direction data;
and determining the slope direction of the sand dune at each sample point according to the coordinates of each sample point.
8. The apparatus of claim 6, wherein the point location determination module is configured to:
generating a Terrain Position Index (TPI) data body by using the elevation data;
generating a terrain slope SPC data body of the construction area according to the elevation data and the TPI data body;
according to the SPC data body and the coordinates of each sample point, determining the slope type of the sand dune at each sample point; the slope types of the sand dunes comprise sand valleys, slope feet, flat areas, medium slope areas, steep slope areas and sand ridges.
9. The apparatus of claim 6, wherein the point location determination module is configured to:
determining the height difference between each sample point and the theoretical point according to the elevation;
determining the plane distance between each sample point and the theoretical point according to the coordinates of each sample point;
and determining the ratio of the height difference to the plane distance as the contribution rate of the sample point.
10. The apparatus according to any one of claims 6 to 9, wherein the point location determining module is configured to:
selecting samples with contribution rate greater than 10% from all samples, slope direction being windward surface, and slope type being not sand valley, slope toe and steep slope area as primary selected samples;
and determining the sample point with the highest contribution rate in the initially selected sample points as a physical point of the three-dimensional seismic exploration.
11. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the method of any of claims 1 to 5 when executing the computer program.
12. A computer readable storage medium, characterized in that the computer readable storage medium stores a computer program which, when executed by a processor, implements the method of any of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011202469.5A CN112462422B (en) | 2020-11-02 | 2020-11-02 | Automatic design method and device for three-dimensional seismic exploration physical point positions in desert area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011202469.5A CN112462422B (en) | 2020-11-02 | 2020-11-02 | Automatic design method and device for three-dimensional seismic exploration physical point positions in desert area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112462422A CN112462422A (en) | 2021-03-09 |
| CN112462422B true CN112462422B (en) | 2024-08-27 |
Family
ID=74834934
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011202469.5A Active CN112462422B (en) | 2020-11-02 | 2020-11-02 | Automatic design method and device for three-dimensional seismic exploration physical point positions in desert area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112462422B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112114354B (en) * | 2019-06-21 | 2024-03-05 | 中国石油天然气集团有限公司 | Mountain three-dimensional seismic exploration excitation point position optimization method and device |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107340534A (en) * | 2016-05-03 | 2017-11-10 | 中国石油化工股份有限公司 | A kind of observation system shot point Optimization Design based on dimensional topography gradient slope aspect |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005062074A1 (en) * | 2003-12-12 | 2005-07-07 | Exxonmobil Upstream Research Company | Method for seismic imaging in geologically complex formations |
| US8326537B2 (en) * | 2004-01-20 | 2012-12-04 | Westerngeco L.L.C. | Survey design using earth observation data |
| CN100549730C (en) * | 2006-07-12 | 2009-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | The ribbon earthquake collection method of a kind of intricately list area |
| US8615362B2 (en) * | 2008-10-10 | 2013-12-24 | Westerngeco L.L.C. | Near-surface geomorphological characterization based on remote sensing data |
| CN102043164A (en) * | 2009-10-26 | 2011-05-04 | 中国石油天然气股份有限公司 | Seismic survey line optimization design method utilizing remote sensing information |
| CN101900832B (en) * | 2010-06-25 | 2012-02-01 | 恒泰艾普石油天然气技术服务股份有限公司 | Ellipse expansion imaging method and device of seismic data processing under complicated ground surface condition |
| CN103091703B (en) * | 2013-01-17 | 2015-09-09 | 中国石油天然气股份有限公司 | Three-dimensional seismic excitation point optimization method for seismic exploration |
| US10481287B2 (en) * | 2016-08-05 | 2019-11-19 | Saudi Arabian Oil Company | Surface consistent statics solution and amplification correction |
| CN109581503B (en) * | 2017-09-28 | 2020-08-25 | 中国石油化工股份有限公司 | Method and system for solving problem of chromatography static correction boundary of desert exploration area |
| US10920585B2 (en) * | 2017-12-26 | 2021-02-16 | Saudi Arabian Oil Company | Determining sand-dune velocity variations |
| CN111353203B (en) * | 2018-12-20 | 2023-08-22 | 中国石油天然气集团有限公司 | Seismic acquisition arrangement layout optimization method, device and storage medium |
-
2020
- 2020-11-02 CN CN202011202469.5A patent/CN112462422B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107340534A (en) * | 2016-05-03 | 2017-11-10 | 中国石油化工股份有限公司 | A kind of observation system shot point Optimization Design based on dimensional topography gradient slope aspect |
Non-Patent Citations (3)
| Title |
|---|
| 张岩.流动沙漠区基于地表压实程度的物理点优选方法.2022年中国石油物探学术年会论文集（上册）.2022,第238-240页. * |
| 检波点位置偏移原则在准噶尔盆地沙漠区地震采集中的应用;王正军等;石油物探;20200131;第59卷(第01期);第32-37页 * |
| 流动沙漠区基于地表压实程度的物理点优选方法;张岩等;2022年中国石油物探学术年会论文集（上册）;20220930;第248-250页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112462422A (en) | 2021-03-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105783810B (en) | Engineering earthwork measuring method based on unmanned plane camera work | |
| CN110780307B (en) | Method for obtaining road cross section based on storage battery car-mounted laser point cloud mobile measurement system | |
| DE19516003A1 (en) | Terrain estimate based on an inclinometer | |
| CN112729130A (en) | Method for measuring height of tree canopy by satellite remote sensing | |
| CN105571571A (en) | Accumulated profile space structure information analysis method based on three-dimensional laser scanning | |
| CN110207676A (en) | The acquisition methods and device of a kind of field ditch pool parameter | |
| CN108427741B (en) | DEM relative error evaluation method based on large number of high-precision control points | |
| CN114295069A (en) | Side slope deformation monitoring method and system for unmanned aerial vehicle carrying three-dimensional laser scanner | |
| CN110411422A (en) | The planing method of builder's road based on BIM | |
| CN108520271B (en) | Submarine geomorphy type sorter design method based on factorial analysis | |
| CN112462422B (en) | Automatic design method and device for three-dimensional seismic exploration physical point positions in desert area | |
| Pardo-Pascual et al. | New methods and tools to analyze beach-dune system evolution using a Real-Time Kinematic Global Positioning System and Geographic Information Systems | |
| CN107860375A (en) | A kind of landslide disaster volume rapid extracting method based on three-dimensional laser scanning technique | |
| CN108460422A (en) | Submarine geomorphy kind identification method based on depth distribution feature | |
| CN116012613B (en) | Method and system for measuring and calculating earthwork variation of strip mine based on laser point cloud | |
| CN111080536A (en) | Self-adaptive filtering method for airborne laser radar point cloud | |
| Hapep et al. | Comparison of Different DEM Generation Methods based on Open Source Datasets. | |
| Zlatanović et al. | Comparison of an automated and manual method for calculating storm runoff response in ungauged catchments in Serbia | |
| Sakitai et al. | Method for extracting the source of falling rock from microtopography highlight map created by high-density aerial laser data | |
| Carlisle | Digital elevation model quality and uncertainty in DEM-based spatial modelling | |
| CN115018973A (en) | Low-altitude unmanned-machine point cloud modeling precision target-free evaluation method | |
| Gülci et al. | The evaluation of SfM technique in the determination of surface deformation on skidding roads following timber harvesting | |
| CN112308907A (en) | Air route planning method for close-range photogrammetry of slope body by using aircraft | |
| CN105737802B (en) | Accumulation section spatial structural form analysis method based on body-sensing camera work | |
| Gu et al. | Ground point extraction using self-adaptive-grid and point to surface comparison |
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 |