CN111273339B - Shot point encryption method and system based on barrier target area - Google Patents

Shot point encryption method and system based on barrier target area Download PDF

Info

Publication number
CN111273339B
CN111273339B CN201811471763.9A CN201811471763A CN111273339B CN 111273339 B CN111273339 B CN 111273339B CN 201811471763 A CN201811471763 A CN 201811471763A CN 111273339 B CN111273339 B CN 111273339B
Authority
CN
China
Prior art keywords
shot
target area
point
encrypted
range
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
Application number
CN201811471763.9A
Other languages
Chinese (zh)
Other versions
CN111273339A (en
Inventor
潘英杰
李伟波
田磊
白志宏
赵君
许银坡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201811471763.9A priority Critical patent/CN111273339B/en
Publication of CN111273339A publication Critical patent/CN111273339A/en
Application granted granted Critical
Publication of CN111273339B publication Critical patent/CN111273339B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/02Generating seismic energy

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)
  • Radar Systems Or Details Thereof (AREA)
  • Storage Device Security (AREA)

Abstract

The invention provides a shot point encryption method and a shot point encryption system based on an obstacle target area, wherein the shot point encryption method comprises the following steps: importing the acquired boundary data of the obstacle area into a preset observation system and gridding the observation system; gridding a target area generated by using a limited offset distance and determining a minimum rectangular surrounding core according to the target area; generating a shot point range to be encrypted and a check point range in the observation system according to the minimum rectangular surrounding core, and respectively meshing the shot point range to be encrypted and the check point range; generating bin attribute contribution values of the shot points to be added in the shot point range to be encrypted to the target area according to the minimum grids of the target area and the check points in the check point range; and selecting the encrypted shot points corresponding to the barrier areas in the shot point range to be encrypted according to the attribute contribution values of the surface elements. The application has the beneficial effects of simply, quickly and efficiently carrying out shot point encryption.

Description

Shot point encryption method and system based on barrier target area
Technical Field
The invention relates to the technical field of seismic acquisition, in particular to a shot point encryption method and system for an obstacle target area.
Background
Today, seismic exploration areas are mostly areas with complex surface conditions, in which various obstacles exist, and the forms of the obstacles include the following: surface and ground objects which can not be normally arranged at a shot point, such as regions with frequent human activities in town regions, dense building regions, houses, factories, bridges and the like, infrastructure places, roads, railways, dams, oil pipelines and the like, and water areas such as rivers, aquaculture ponds and the like; important facilities with certain safety distance need to be arranged, such as cultural relics and other protection areas, railway tunnels and other, military pipe areas and the like; some artificially set areas that cannot be blasted; an obstacle area that has occurred unexpectedly after the survey is completed. The obstacles in the exploration area bring great difficulty to field seismic data acquisition, so that cannons and inspection points cannot be arranged according to the designed observation system rule, shallow layers in the seismic data are lost, and the overall understanding of the structural blocks and the resource evaluation of the pits are influenced.
Several methods are commonly used to treat the obstacle area: (1) By adopting large-distance recovery shot points or direct blank shots, the construction design is unreasonable, the covering times are uneven, the data of the shallow and middle layers of the collected data are lost, and the geological exploration task is influenced; (2) And carrying out encrypted shot point processing on the obstacle area by using the empirical value. If bilateral encryption shot points are carried out around the barrier, the encryption quantity is carried out according to 1.2-1.5 times of the quantity of the shot points in the original barrier, and the like, the seismic data acquisition quality is improved to a certain extent by carrying out encryption shot point processing on the barrier area, but the encryption shot points are often designed by experience, the result varies from person to person, the error is large, and the size and the imaging effect of a seismic section gap cannot be really predicted; and (3) simulating a calculation method. When an observation system is designed, firstly, encrypted shot point processing is carried out on an obstacle area, then surface element coverage data, illumination and the like of a target area are simulated and calculated, finally, the designed observation system parameters are optimized through analyzing and calculating results, the surface element coverage data, the illumination and the like of the target area are recalculated, repeated correction is carried out until the result is satisfied, repeated trial and error are also needed according to experience by using the method, the pertinence is poor, the result cannot be directly quantitatively displayed, the shot point encryption process is complicated, repeated iteration is needed, the calculation amount is large, the consumed time is long, and the timeliness is poor.
Therefore, how to provide a simple, fast and efficient shot encryption technology is a technical problem to be solved urgently at present.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention provides a shot point encryption method and system based on a barrier target area, which aim to encrypt shot points to the target area in a targeted and quantitative manner by quickly calculating the bin attribute contribution of the shot points to be added to the target area, reduce the blindness of shot point encryption, and improve the scientificity and rationality of shot point encryption, thereby improving the coverage times and coverage uniformity of the target area, reducing the loss of shallow data, improving the quality of imaging data, having good application prospects in places with dense large barrier areas such as urban areas and the like, and having the beneficial effects of simply, quickly and efficiently carrying out shot point encryption.
In order to achieve the above object, the present invention provides a shot encryption method based on an obstacle target area, including:
importing the acquired boundary data of the obstacle area into a preset observation system, and gridding the observation system;
gridding a target area generated by using a limited offset distance, and determining a minimum rectangular surrounding core according to the target area;
generating a shot point range to be encrypted and a point detection range in the observation system according to the minimum rectangular surrounding core, and respectively meshing the shot point range to be encrypted and the point detection range; each minimum grid in the shot point range to be encrypted corresponds to a shot point to be added; each minimum grid of the detection point range corresponds to one detection point; each inspection point at least corresponds to one shot point to be added; the shot point range to be encrypted corresponds to the inspection point range one by one;
generating a surface element attribute contribution value of each shot to be added in the shot range to be encrypted to the target area according to each minimum grid of the target area and each detection point in the detection point range;
and selecting the encrypted shot points corresponding to the barrier area in the shot point range to be encrypted according to the size of the attribute contribution value of each surface element.
The invention also provides a shot point encryption system based on the barrier target area, which comprises the following components:
the guiding-in unit is used for guiding the acquired boundary data of the obstacle area into a preset observation system and meshing the observation system;
the surrounding kernel generating unit is used for gridding a target area generated by using a limited offset distance and determining a minimum rectangular surrounding kernel according to the target area;
the discrete unit is used for generating a shot point range to be encrypted and a point detection range in the observation system according to the minimum rectangular surrounding core, and respectively meshing the shot point range to be encrypted and the point detection range; each minimum grid in the shot point range to be encrypted corresponds to a shot point to be added; each minimum grid of the detection point range corresponds to one detection point; each inspection point at least corresponds to one shot point to be added; the shot point range to be encrypted corresponds to the inspection point range one by one;
a contribution value generating unit, configured to generate, according to each minimum grid of the target area and each check point in the check point range, a bin attribute contribution value of each shot point to be added in the shot point range to be encrypted to the target area;
and the encryption unit is used for selecting the encrypted shot points corresponding to the obstacle area in the shot point range to be encrypted according to the size of the attribute contribution value of each bin.
The invention provides a shot point encryption method and a shot point encryption system based on an obstacle target area, which comprise the following steps: importing the acquired boundary data of the obstacle area into a preset observation system, and gridding the observation system; gridding a target area generated by using a limited offset distance, and determining a minimum rectangular surrounding core according to the target area; generating a shot point range to be encrypted and a point detection range in the observation system according to the minimum rectangular surrounding core, and respectively meshing the shot point range to be encrypted and the point detection range; each minimum grid in the shot point range to be encrypted corresponds to a shot point to be added; each minimum grid of the detection point range corresponds to one detection point; each inspection point at least corresponds to one shot point to be added; the shot point range to be encrypted corresponds to the inspection point range one by one; generating bin attribute contribution values of the shot points to be added in the shot point range to be encrypted to the target area according to the minimum grids of the target area and the check points in the check point range; and selecting the encrypted shot points corresponding to the obstacle area in the shot point range to be encrypted according to the size of the attribute contribution value of each surface element. This application is through the quantitative encrypted shot point to the target area, has the scientificity and the rationality that improve the shot point and encrypt, improves the regional coverage number of times of target and covers even degree, reduces shallow data disappearance and simple, quick, high-efficient beneficial effect who carries out the shot point and encrypt.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the prior art descriptions will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flow chart of a shot point encryption method based on an obstacle target area of the present application;
FIG. 2 is a flowchart of a shot encryption method based on an obstacle target area according to an embodiment of the present application;
FIG. 3 is a diagram of an observation system importing preliminary layout according to an embodiment of the present application;
FIG. 4 is a diagram of an observation system for importing obstacle data according to an embodiment of the present application;
FIG. 5 is a flowchart illustrating a method for implementing step S204 according to an embodiment of the present application;
fig. 6 is a flowchart illustrating a specific implementation method of step S205 in an embodiment of the present application;
FIG. 7 is a diagram illustrating contribution values of shots to be added according to an embodiment of the present application;
FIG. 8 is a diagram of an observation system to which shot contributions are to be added according to an embodiment of the present application;
FIG. 9 is a plot of the effect of bin calculations on shot encryption in an embodiment of the present application;
FIG. 10 is a chart illustrating the effect of bin calculations after shot encryption in an embodiment of the present application;
fig. 11 is a schematic structural diagram of a shot point encryption system based on an obstacle target area according to the present application.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The terms "first," "second," "8230," "8230," and the like as used herein do not particularly denote any order or sequence, nor are they intended to limit the invention, but rather are used to distinguish one element from another or from another element described in the same technical term.
As used herein, the terms "comprising," "including," "having," "containing," and the like are open-ended terms that mean including, but not limited to.
As used herein, "and/or" includes any and all combinations of the described items.
Aiming at the defects in the prior art, the invention provides a shot point encryption method based on an obstacle target area, the flow chart of which is shown in figure 1, and the method comprises the following steps:
s101: and importing the acquired boundary data of the obstacle area into a preset observation system, and gridding the observation system.
S102: and gridding the target area generated by using the limited offset distance, and determining a minimum rectangular surrounding core according to the target area.
S103: and generating a shot point range to be encrypted and a point detection range in the observation system according to the minimum rectangular surrounding core, and respectively meshing the shot point range to be encrypted and the point detection range. Each minimum grid in the shot point range to be encrypted corresponds to a shot point to be added; each minimum grid of the detection point range corresponds to one detection point; each inspection point at least corresponds to one shot point to be added; the shot point range to be encrypted corresponds to the inspection point range one by one.
S104: and generating a bin attribute contribution value of each shot to be added in the shot range to be encrypted to the target area according to each minimum grid of the target area and each detection point in the detection point range.
S105: and selecting the encrypted shot points corresponding to the barrier areas in the shot point range to be encrypted according to the attribute contribution values of the surface elements.
As can be seen from the process shown in fig. 1, the method guides the acquired boundary data of the obstacle area into a preset observation system, and performs gridding on the observation system; gridding a target area generated by using the limited offset distance, and determining a minimum rectangular surrounding core according to the target area; generating a shot point range to be encrypted and a point detection range in the observation system according to the minimum rectangular surrounding core, and respectively meshing the shot point range to be encrypted and the point detection range; each minimum grid in the shot point range to be encrypted corresponds to a shot point to be added; each minimum grid of the detection point range corresponds to one detection point; each inspection point at least corresponds to one shot point to be added; the shot point range to be encrypted corresponds to the inspection point range one by one; generating a surface element attribute contribution value of each shot point to be added in the shot point range to be encrypted to the target area according to each minimum grid of the target area and each detection point in the detection point range; and selecting the encrypted shot points corresponding to the barrier areas in the shot point range to be encrypted according to the attribute contribution values of the surface elements. This application is through the quantitative encrypted shot point to the target area, has the scientificity and the rationality that improve the shot point and encrypt, improves the regional coverage number of times of target and covers even degree, reduces shallow data disappearance and simple, quick, high-efficient beneficial effect who carries out the shot point and encrypt.
In order to make those skilled in the art better understand the present invention, a more detailed embodiment is listed below, as shown in fig. 2, an embodiment of the present invention provides a shot encryption method based on an obstacle target area, the method includes the following steps:
s201: and importing the acquired boundary data of the obstacle area into a preset observation system, and gridding the observation system.
In specific implementation, as shown in fig. 3, an observation system is first designed according to actual production requirements of a exploration area. As shown in fig. 4, the acquired boundary data of the obstacle area is imported into a preset observation system, and obstacle avoidance design for the obstacle is implemented according to the boundary data of the obstacle area and the safety distance of the obstacle area. The boundary data of the obstacle area is obtained by using a high-precision digital satellite photo and seismic measurement, and the boundary data of the obstacle area is generated according to the description of the obstacle area. The boundary data of the obstacle region includes: coordinates of the obstacle region, safety distance, and the like. The safe distance of the barrier area is that the barrier area is expanded outwards according to a certain distance, and shot point excitation cannot be carried out within the safe distance of the barrier.
Assume that the scope of the preset observation system is as follows:
{(Geometry xmin ,Geometry ymin ),(Geometry xmax ,Geometry ymax )}
among them, geometry xmin Is the minimum value of the x coordinate of the observation system, geometry xmax Is the maximum value of the x coordinate of the observation system, geometry ymin Being the minimum of the y-coordinate of the observation system, geometry ymax Is the maximum value of the y coordinate of the viewing system.
Gridding the observation system into discrete grid data, each minimum grid having a size of { Cell } width ,Cell hight }. Wherein, cell width Cell, the width of the minimum grid of the observation system after meshing hight Is the height of the minimum grid of the gridded observation system.
S202: the obstacle region and the safety distance of the obstacle region are gridded.
In specific implementation, the obstacle area and the safe distance of the obstacle area are gridded, and the gridded obstacle area is used for determining which positions can not be added with shot points. And marking the obstacle on the minimum grid unit after gridding.
The obstacle area gridding formula is shown as formula (1):
Figure GDA0003739544720000061
wherein, gridObstacle (x,y) A grid (x, y) of the gridded Obstacle region, x being the x-coordinate of the gridded Obstacle region, y being the y-coordinate of the gridded Obstacle region, obstacle num Number of obstacles, U is union operation, obstacle i Is the ith obstacle, gridObstacle (x,y) ∈Obstacle i The grid (x, y) indicating the gridded obstacle region belongs to the i-th obstacle region.
For the obstacle with safe distance, when the shot point is encrypted, the shot point to be added needs to be added beyond the safe distance of the obstacle area. When the grid processing is performed on such an obstacle region, it is necessary to include a safe distance region, and the grid processing is performed on the obstacle region and the safe distance region as shown in formula (2):
Figure GDA0003739544720000062
wherein, the ObsSafe num ObsSafe as the number of obstacles having a safe distance i For the ith obstacle with a safe distance,
Figure GDA0003739544720000063
for an open operation, S is a structural operator,
Figure GDA0003739544720000064
is a barrier region, D s Is an operator region.
S203: and gridding the target area generated by using the limited offset distance, and determining a minimum rectangular surrounding core according to the target area.
Wherein the target area comprises at least a part of the obstacle area. The area information of the target area includes: the width of the alignment sheet, the height of the alignment sheet, the depth of the target layer, the maximum reflection angle, etc., but the invention is not limited thereto.
In specific implementation, the target area generated by using the limited offset distance is gridded, as shown in formula (3):
Figure GDA0003739544720000071
wherein gridtargeregion (x, y) is a minimum grid coordinate (i.e. a first index) of the gridded target region, and TargetRegion is a target region range. A minimum grid of the target area is marked by 0 or 1, a 0 on the minimum grid indicates that the grid coordinate does not belong to the target area, and a 1 on the minimum grid indicates that the grid coordinate belongs to the target area.
Determining a minimum rectangular surrounding kernel according to the gridded target area, as shown in formula (4):
TargetRect={(TargetRect xmin ,TargetRect ymin ),(TargetRect xmax ,TargetRect ymax )} (4)
wherein, targetlect xmin Is the minimum value of the x coordinate of the minimum rectangle surrounding the kernel, targetRef xmax TargetRef, the maximum x-coordinate value of the minimum rectangle bounding the kernel ymin Minimum value of y coordinate, targetRef, for the minimum rectangle bounding the kernel ymax The minimum rectangle encloses the maximum value of the Y coordinate of the kernel.
S204: and generating a shot point range to be encrypted and a point detection range in the observation system according to the minimum rectangular surrounding core, and respectively meshing the shot point range to be encrypted and the point detection range.
Each minimum grid in the shot point range to be encrypted corresponds to a shot point to be added; each minimum grid of the detection point range corresponds to one detection point; each inspection point at least corresponds to one shot point to be added; the shot point range to be encrypted corresponds to the inspection point range one by one.
As shown in fig. 5, the specific execution of step S204 includes the following steps:
s301: judging whether each minimum grid in the minimum rectangular surrounding core belongs to a target area; if so, generating a shot point range which responds to each minimum grid in the minimum rectangular surrounding core according to the width of the arrangement piece, the height of the arrangement piece, the depth of the target layer and the maximum reflection angle of the target area.
Specifically, it is determined whether each minimum mesh in the minimum rectangular bounding kernel belongs to the target region, that is, whether the flag of each minimum mesh gridtargeregion (x, y) in the minimum rectangular bounding kernel is 1, as shown in formula (3).
If the minimum grid GridTargetRegion (x, y) in the minimum rectangular bounding kernel belongs to the target region, generating a shot point range { (SR) giving response to each minimum grid in the minimum rectangular bounding kernel according to the width of the permutation piece, the height of the permutation piece, the depth of the target layer and the maximum reflection angle of the target region xmin ,SR ymi n),(SR xmax ,SR ymax ) As shown in equation (5):
Figure GDA0003739544720000081
wherein SRxmin is the minimum value of x coordinate of shot point range, SR ymin Is the minimum value of the y coordinate, SR, of the shot point range xmax Maximum value of x coordinate, SR, of shot point range ymax Maximum value of y coordinate of shot-to-shot point range { (SR) xmin ,SR ymin ),(SR xmax ,SR ymax ) Is the shot check point range, cell width Cell, being the width of the smallest grid within the smallest rectangular bounding kernel hight The minimum rectangle surrounds the height of the minimum grid in the kernel, pathWidth is the width of the arrangement slice of the target area, pathheight is the height of the arrangement slice of the target area, levelDept is the depth of the target layer of the target area, and theta is the maximum reflection angle of the target area. x represents the x-axis coordinate value of GridTargetRegion (x, y) and y represents the y-axis coordinate value of GridTargetRegion (x, y).
S302: and gridding shot-test pairs existing in the shot-test point range to generate a shot point range to be encrypted and a test point range corresponding to the shot point range to be encrypted.
In specific implementation, searching shot and check point range { (SR) for observation system xmin ,SR ymin ),(SR xmax ,SR ymax ) The shot pairs present in the } are gridded.
And gridding the shot points within the shot detection point range to generate shot points to be added, wherein the gridding formula is shown as a formula (6):
Figure GDA0003739544720000082
wherein, gridShot (x, y) represents the minimum grid unit corresponding to the shot point after gridding,
Figure GDA0003739544720000083
indicating that GridShot (x, y) is the added shot in the observation system,
Figure GDA0003739544720000084
denotes GridShot (x, y) as a checkpoint in the viewing system, gridObstacle (x,y) =1 indicates that GridShot (x, y) is in an obstacle region, 1 represents a shot to be added, and 0 represents a shot that cannot be added. x represents the x-axis coordinate value of GridShot (x, y), and y represents the y-axis coordinate value of GridShot (x, y).
And gridding the detection points within the range of the shot detection points, as shown in formula (7):
Figure GDA0003739544720000091
wherein GridReceiver (x, y) represents the minimum grid unit corresponding to the checked point after gridding,
Figure GDA0003739544720000093
indicating that the minimum grid unit corresponding to the gridded detection point belongs to the detection point existing in the system of the observation systemAnd marking the grid receiver (x, y) corresponding to the detection point as 1, otherwise, marking the grid receiver (x, y) as 0.x represents the x-axis coordinate value of GridReceiver (x, y), and y represents the y-axis coordinate value of GridReceiver (x, y).
And generating a shot point range to be encrypted according to each shot point GridShot (x, y) to be added, and generating a detection point range according to each gridded detection point GridReceiver (x, y).
S205: and generating the bin attribute contribution value of each shot point to be added in the shot point range to be encrypted to the target area according to each minimum grid of the target area and each check point in the check point range.
As shown in fig. 6, the specific execution of step S205 includes the following steps:
s401: and respectively calculating the shot points to be added corresponding to each inspection point in the inspection point range by taking each minimum grid of the target area as a central point. Wherein, the shot information of the shot to be added comprises: an x-coordinate and a y-coordinate; each minimum grid of the target area includes: a first grid index; the point detection comprises the following steps: a second grid index.
Specifically, with the detected point GridReceiver (rx, ry) as a reference and GridTargetRegion (x, y) as a center, calculating the GridShot (sx, sy) to be added:
Figure GDA0003739544720000092
wherein sx is an x coordinate of the shot point GridShot (sx, sy) to be added, sy is a y coordinate of the shot point GridShot (sx, sy) to be added, position information of the shot point GridShot (sx, sy) to be added is (sx, sy), (x, y) is a first grid index of any minimum grid gridtargetregetion (x, y) of the gridded target region, and (rx, ry) is a second grid index of any inspection point GridReceiver (rx, ry) in the inspection point range.
S402: and judging whether the shot to be added corresponding to each shot point information belongs to the shot point range to be encrypted and is in a state to be encrypted or not according to each shot point information. And if so, integrating the bin attribute function of the target area to generate the contribution of each shot point to be added to the bin attribute of the target area.
Specifically, whether the shot GridShot (sx, sy) to be added corresponding to each shot point information belongs to the range of the shot points to be encrypted and is in a state to be encrypted is sequentially judged, that is, the identifier of the minimum grid unit corresponding to the GridShot (sx, sy) is 1.
If the shot point GridShot (sx, sy) to be added corresponding to the shot point information (sx, sy) belongs to the shot point range to be encrypted and the identifier of the minimum grid unit corresponding to the GridShot (sx, sy) is 1, integrating the bin attribute function of the target area to generate the bin attribute contribution of the shot point GridShot (sx, sy) to the target area, wherein a specific calculation formula is shown as a formula (9):
Figure GDA0003739544720000101
wherein, (sx, sy) is the position information of a shot point to be added in the range of the shot point to be encrypted,
Figure GDA0003739544720000102
for the bin attribute contribution degree, binProperty, of the shot GridShot (sx, sy) to be added corresponding to the target area n For the nth bin attribute, binPropertYFun, corresponding to the target region to be added shot GridShot (sx, sy) n Is the bin attribute function of the target region, targetRegion is the target region, - [ integral ] integral [ integral ] factor TargetRegion BinPropertyFun n (x, y) dxdy is the integral operation of the bin attribute function for the target area.
S403: and generating a bin attribute contribution value of each shot point to be added to the target area according to the bin attribute contribution degree of each shot point to be added to the target area.
Specifically, assume that each bin attribute BinProperty n If the covering times of the bin are equal, the bin calculates the function BinPropertyFun n Calculating the bin attribute tributary value BinVal (sx, sy) of the shot point GridShot (sx, sy) to be added to the target area, as shown in a formula (10):
Figure GDA0003739544720000103
wherein BinVal (sx, sy) represents the bin attribute contribution value of the shot GridShot (sx, sy) to be added to the target area, a n Is a predetermined constant coefficient, and c is a constant.
As shown in fig. 7, the triangular region in fig. 7 is the target region, the left side of the triangle is marked as the shot GridShot (sx, sy) to be added, and the right side is the inspection point. The number marked on each shot to be added is the surface element attribute membership value of each shot to be added GridShot (sx, sy) to the target area, and the larger the numerical value is, the larger the surface element attribute membership value is. Wherein the grey shots represent shots that have been added that cannot be shot-encrypted. Fig. 8 is a view showing the contribution value of the shot to be added in the present embodiment.
S206: and selecting the encrypted shot points corresponding to the barrier areas in the shot point range to be encrypted according to the attribute contribution values of the surface elements.
Specifically, the bin attribute contribution values of the shot points to be added to the target area are sorted from large to small, a plurality of shot points to be added with the largest bin attribute contribution values of the target area are selected as encrypted shot points corresponding to the barrier area according to design requirements, shot point encryption is carried out, and shot point encryption of the target area with all shallow and middle layer data missing is completed.
Fig. 9 is a bin calculation effect graph before shot point encryption, and fig. 10 is a bin calculation effect graph after shot point encryption, and it can be known through comparison of the two figures that the coverage times and the coverage uniformity of the target area are improved, shallow data loss is reduced, and the quality of imaging data is improved.
Based on the same application concept as the shot point encryption method based on the obstacle target area, the invention also provides a shot point encryption system based on the obstacle target area, which is described in the following embodiment. Because the principle of solving the problem of the shot point encryption system based on the obstacle target area is similar to that of the shot point encryption method based on the obstacle target area, the implementation of the shot point encryption system based on the obstacle target area can refer to the implementation of the shot point encryption method based on the obstacle target area, and repeated parts are not described again.
Fig. 11 is a schematic structural diagram of a shot point encryption system based on an obstacle target area provided in the present application, and as shown in fig. 11, the shot point encryption system based on the obstacle target area includes: import section 101, surrounding core generation section 102, discrete section 103, contribution value generation section 104, and encryption section 105.
An importing unit 101 is configured to import the acquired boundary data of the obstacle area into a preset observation system, and perform meshing on the observation system.
A surrounding kernel generating unit 102, configured to perform gridding on the target area generated by using the shift-limited distance, and determine a minimum rectangular surrounding kernel according to the target area.
And the discrete unit 103 is configured to generate a shot point range to be encrypted and a point detection range in the observation system according to the minimum rectangular surrounding kernel, and perform meshing on the shot point range to be encrypted and the point detection range respectively. Each minimum grid in the shot point range to be encrypted corresponds to one shot point to be added; each minimum grid of the detection point range corresponds to one detection point; each inspection point at least corresponds to one shot point to be added; the shot point range to be encrypted corresponds to the inspection point range one by one.
And a contribution value generating unit 104, configured to generate, according to each minimum grid of the target area and each point within the point detection range, a bin attribute contribution value of each to-be-added shot within the shot range to be encrypted to the target area.
And the encryption unit 105 is used for selecting the encrypted shot points corresponding to the barrier area within the shot point range to be encrypted according to the size of the attribute contribution value of each surface element.
The invention provides a shot point encryption method and a shot point encryption system based on an obstacle target area, wherein the shot point encryption method comprises the following steps: importing the acquired boundary data of the obstacle area into a preset observation system, and gridding the observation system; gridding a target area generated by using the limited offset distance, and determining a minimum rectangular surrounding core according to the target area; generating a shot point range to be encrypted and a point detection range in the observation system according to the minimum rectangular surrounding core, and respectively meshing the shot point range to be encrypted and the point detection range; each minimum grid in the shot point range to be encrypted corresponds to a shot point to be added; each minimum grid of the detection point range corresponds to one detection point; each inspection point at least corresponds to one shot point to be added; the shot point range to be encrypted corresponds to the inspection point range one by one; generating bin attribute contribution values of the shot points to be added in the shot point range to be encrypted to the target area according to the minimum grids of the target area and the check points in the check point range; and selecting the encrypted shot points corresponding to the barrier areas in the shot point range to be encrypted according to the attribute contribution values of the surface elements. This application has the quantitative shot point of encrypting to the target area, has the scientificity and the rationality that improve shot point encryption, improves target area's cover number of times and cover even degree, reduces shallow layer data disappearance and simple, quick, high-efficient beneficial effect who carries out shot point encryption.
As will be appreciated by one skilled in the art, 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 principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed, and in summary, the content of the present specification should not be construed as a limitation to the present invention.

Claims (10)

1. A shot point encryption method based on an obstacle target area is characterized by comprising the following steps:
importing the acquired boundary data of the obstacle area into a preset observation system, and gridding the observation system;
gridding a target area generated by using a limited offset distance, and determining a minimum rectangular surrounding core according to the target area;
generating a shot point range to be encrypted and a point detection range in the observation system according to the minimum rectangular surrounding core, and respectively meshing the shot point range to be encrypted and the point detection range; each minimum grid in the shot point range to be encrypted corresponds to a shot point to be added; each minimum grid of the detection point range corresponds to one detection point; each inspection point at least corresponds to one shot point to be added; the shot point range to be encrypted corresponds to the inspection point range one by one;
generating a surface element attribute contribution value of each shot to be added in the shot range to be encrypted to the target area according to each minimum grid of the target area and each detection point in the detection point range;
and selecting the encrypted shot points corresponding to the obstacle area in the shot point range to be encrypted according to the size of the attribute contribution value of each surface element.
2. The shot encryption method based on an obstacle target area according to claim 1, further comprising: and gridding the obstacle area and the safe distance of the obstacle area.
3. The shot point encryption method based on an obstacle target area according to claim 2, wherein the target area comprises at least a portion of the obstacle area; the region information of the target region includes: the arrangement sheet width, the arrangement sheet height, the target layer depth and the maximum reflection angle.
4. The shot point encryption method based on the obstacle target area according to claim 3, wherein the generating a shot point range to be encrypted and a point detection range in the observation system according to the minimum rectangular surrounding kernel, and respectively gridding the shot point range to be encrypted and the point detection range comprises:
judging whether each minimum grid in the minimum rectangular surrounding core belongs to the target area; if so, generating a shot point range which responds to each minimum grid in the minimum rectangular surrounding core according to the width of the arrangement piece, the height of the arrangement piece, the depth of the target layer and the maximum reflection angle of the target area;
and gridding shot-test pairs existing in the shot-test point range to generate the shot-test point range to be encrypted and a test point range corresponding to the shot-test point range to be encrypted.
5. The method according to claim 1, wherein the generating a bin attribute contribution value of each shot to be added in the shot range to be encrypted to the target area according to each minimum grid of the target area and each check point in the check point range comprises:
calculating the shot points to be added corresponding to each detection point in the detection point range by taking each minimum grid of the target area as a central point; the shot point information of the shot points to be added comprises the following steps: an x-coordinate and a y-coordinate; each minimum grid of the target region includes: a first grid index; the point detection comprises the following steps: a second grid index;
judging whether the shot points to be added corresponding to the shot point information belong to the shot point range to be encrypted and are in a state to be encrypted according to the shot point information; if so, integrating the bin attribute function of the target area to generate the bin attribute contribution of each shot to be added to the target area;
and generating a bin attribute contribution value of each shot point to be added to the target area according to the bin attribute contribution degree of each shot point to be added to the target area.
6. The method of claim 4, wherein the calculation formula for generating a shot point range responsive to each minimum grid in the minimum rectangular bounding kernel according to the width of the permutation piece, the height of the permutation piece, the depth of the target layer and the maximum reflection angle of the target region is as follows:
Figure FDA0003739544710000021
wherein, SR xmin Is the minimum value of the x coordinate, SR, of the shot point range ymin Is the minimum value of the y coordinate, SR, of the shot-check point range xmax Is the maximum value of the x coordinate, SR, of the shot-inspection point range ymax Is the maximum value of the y coordinate of the shot point range { (SR) xmin ,SR ymin ),(SR xmax ,SR ymax ) Is the shot check point range, cell width Cell for the width of the minimum grid within the minimum rectangular bounding kernel hight The height of the minimum grid in the minimum rectangle surrounding core is defined, pathWidth is the width of the arrangement slice of the target area, pathheight is the height of the arrangement slice of the target area, levelDept is the target layer depth of the target area, theta is the maximum reflection angle of the target area, x represents the coordinate value of the x axis of the minimum grid in the minimum rectangle surrounding core, and y represents the coordinate value of the y axis of the minimum grid in the minimum rectangle surrounding core.
7. The shot encryption method based on the obstacle target area according to claim 5, wherein the calculation formula for calculating the shot to be added corresponding to each inspection point in the inspection point range by taking each minimum grid of the target area as a central point is as follows:
Figure FDA0003739544710000031
wherein sx is an x coordinate of the shot to be added, sy is a y coordinate of the shot to be added, (sx, sy) is position information of the shot to be added, (x, y) is a first grid index of any minimum grid of the target area after gridding, gridReceiver (rx, ry) is any detected point in the detected point range, (rx, ry) is a second grid index of any detected point in the detected point range, and if GridReceiver (rx, ry) =1 indicates that a detected point exists in the detected point range.
8. The shot encryption method based on the obstacle target area according to claim 5, wherein the calculation formula for generating the contribution degree of each shot to be added to the bin attribute of the target area according to the integration of the bin attribute function of the target area is as follows:
Figure FDA0003739544710000032
wherein, (sx, sy) is the position information of a shot to be added in the range of the shot to be encrypted,
Figure FDA0003739544710000033
BinProperty corresponding to the bin attribute contribution degree of the target area for the shot point to be added n BinPropertyFun as the nth surface element attribute of the shot to be added corresponding to the target area n Is the surface element attribute function of the target area, targetRegion is the target area, and ^ integral TargetRegion BinPropertyFun n (x, y) dxdy is an integral operation of the bin attribute function of the target area.
9. The shot encryption method based on the obstacle target area according to claim 8, wherein the calculation formula for generating the bin attribute contribution value of each shot to be added to the target area according to the bin attribute contribution degree of each shot to be added to the target area is as follows:
Figure FDA0003739544710000034
wherein BinVal (sx, sy) is a bin attribute contribution value of a shot point to be added corresponding to the target area, a n Is a predetermined constant coefficient, and c is a constant.
10. A shot encryption system based on an obstacle target area, comprising:
the guiding-in unit is used for guiding the acquired boundary data of the obstacle area into a preset observation system and meshing the observation system;
the surrounding kernel generating unit is used for gridding a target area generated by using the limited offset distance and determining a minimum rectangular surrounding kernel according to the target area;
the dispersion unit is used for generating a shot point range to be encrypted and a point detection range in the observation system according to the minimum rectangular surrounding core, and respectively gridding the shot point range to be encrypted and the point detection range; each minimum grid in the shot point range to be encrypted corresponds to a shot point to be added; each minimum grid of the detection point range corresponds to one detection point; each inspection point at least corresponds to one shot point to be added; the shot point range to be encrypted corresponds to the inspection point range one by one;
a contribution value generating unit, configured to generate, according to each minimum grid of the target area and each inspection point in the inspection point range, a bin attribute contribution value of each to-be-added shot point in the to-be-encrypted shot point range to the target area;
and the encryption unit is used for selecting the encrypted shot points corresponding to the barrier area in the shot point range to be encrypted according to the size of each surface element attribute contribution value.
CN201811471763.9A 2018-12-04 2018-12-04 Shot point encryption method and system based on barrier target area Active CN111273339B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811471763.9A CN111273339B (en) 2018-12-04 2018-12-04 Shot point encryption method and system based on barrier target area

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811471763.9A CN111273339B (en) 2018-12-04 2018-12-04 Shot point encryption method and system based on barrier target area

Publications (2)

Publication Number Publication Date
CN111273339A CN111273339A (en) 2020-06-12
CN111273339B true CN111273339B (en) 2022-12-02

Family

ID=70996568

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811471763.9A Active CN111273339B (en) 2018-12-04 2018-12-04 Shot point encryption method and system based on barrier target area

Country Status (1)

Country Link
CN (1) CN111273339B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112415573A (en) * 2020-10-30 2021-02-26 中国石油天然气集团有限公司 Shot and geophone point arrangement method and device based on barrier

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010120301A1 (en) * 2009-04-16 2010-10-21 Landmark Graphics Corporation Seismic imaging systems and methods employing a fast target-oriented illumination calculation
CN102540254B (en) * 2010-12-29 2014-11-26 中国石油天然气集团公司 Method for determining effective excitation region based on lighting analysis
CN103543465B (en) * 2012-07-12 2016-06-08 中国石油天然气集团公司 The method of effective shot point is determined based on zone of interest illumination energy
US20140226438A1 (en) * 2013-02-13 2014-08-14 Wireless Seismic, Inc. Assigned scheduled acquisition process in wireless exploration
CN104076387B (en) * 2014-06-06 2017-07-07 中国石油集团东方地球物理勘探有限责任公司 A kind of shallow sea barrier OBC earthquake-capturing focus construction methods
CN105487106B (en) * 2014-09-18 2018-03-09 中国石油化工股份有限公司 A kind of benefit big gun method based on the illumination of Gaussian ray bundle target zone energy
CN104502956A (en) * 2014-12-16 2015-04-08 中国石油天然气集团公司 Handling method for avoiding barrier
CN104502962B (en) * 2014-12-16 2017-03-08 中国石油天然气集团公司 A kind of method and device of design shot point
CN105607115B (en) * 2015-12-21 2018-02-02 中国石油天然气集团公司 The method and apparatus that shot point is laid

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
春光探区新层系地震采集观测系统设计研究;董淑兰;《石油地质与工程》;20180525(第03期);第56-58页 *

Also Published As

Publication number Publication date
CN111273339A (en) 2020-06-12

Similar Documents

Publication Publication Date Title
Pan et al. A deterministic-stochastic identification and modelling method of discrete fracture networks using laser scanning: Development and case study
Gargoum et al. Assessing stopping and passing sight distance on highways using mobile LiDAR data
Soliman et al. Effects of vertical accuracy of digital elevation model (DEM) data on automatic lineaments extraction from shaded DEM
Wang et al. Refined micro-scale geological disaster susceptibility evaluation based on UAV tilt photography data and weighted certainty factor method in Mountainous Area
CN103941286A (en) Rapid three-dimensional horizon interpretation method based on horizontal navigation
CN111539904B (en) Disaster vulnerability prediction method based on rainfall
Suh et al. Subsidence hazard assessment at the Samcheok coalfield, South Korea: a case study using GIS
CN114547866B (en) Prefabricated part intelligent detection method based on BIM-unmanned aerial vehicle-mechanical dog
CN113537141B (en) Method and system for rapidly detecting piping and landslide diseases of dykes and dams
CN106291671A (en) A kind of automatic troubleshooting method of stereo observing system based on satellite image data
Zhao et al. Combining ICESat-2 photons and Google Earth Satellite images for building height extraction
Fanos et al. Multi-scenario rockfall hazard assessment using LiDAR data and GIS
Chen et al. A novel image-based approach for interactive characterization of rock fracture spacing in a tunnel face
Fanos et al. Assessment of multi-scenario rockfall hazard based on mechanical parameters using high-resolution airborne laser scanning data and GIS in a tropical area
Wang et al. Fragmentation calculation method for blast muck piles in open-pit copper mines based on three-dimensional laser point cloud data
Xu et al. Comprehensive identification of statistical homogeneity of fractured rock masses for a candidate HLW repository site, China
CN111273339B (en) Shot point encryption method and system based on barrier target area
Agrawal et al. Topographic classification of North Eastern Region of India using geospatial technique and following seismic code provisions
Tarquini et al. Uncertainties in lava flow hazard maps derived from numerical simulations: The case study of Mount Etna
Zhan et al. A stepwise approach for 3D fracture intersection analysis and application to a hydropower station in Southwest China
Ji et al. An automated method to build 3D multi-scale geological models for engineering sedimentary layers with stratum lenses
Wu et al. Stability analysis of rock blocks around a cross-harbor tunnel using the improved morphological visualization method
Richardson Modeling the construction and evolution of distributed volcanic fields on Earth and Mars
Zhang et al. A novel approach to investigating 3D fracture connectivity in ultrahigh steep rock slopes
Rákay et al. Surveying of inaccessible rock faces and volume calculation of the irregular solids using robotic total station.

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