CN111273339A - 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 PDFInfo
- Publication number
- CN111273339A CN111273339A CN201811471763.9A CN201811471763A CN111273339A CN 111273339 A CN111273339 A CN 111273339A CN 201811471763 A CN201811471763 A CN 201811471763A CN 111273339 A CN111273339 A CN 111273339A
- 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.)
- Granted
Links
Images
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/02—Generating seismic energy
Abstract
The invention provides a shot point encryption method and system based on an obstacle target area, which 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 method has the beneficial effects of simple, quick and efficient shot point encryption.
Description
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 unexpected obstacle area after the completion of the reconnaissance. 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; (3) and (5) 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 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.
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 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 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 system based on an obstacle target area, which 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 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 description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flow chart of a shot 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 for importing preliminary layout in 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 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.
As used herein, the terms "first," "second," … …, etc. do not denote any order or order, nor are they used to limit the invention, but rather are used to distinguish one element from another element or operation described by the same technical terms.
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 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.
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 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 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 range of the preset observation system is as follows:
{(Geometryxmin，Geometryymin)，(Geometryxmax，Geometryymax)}
among them, GeometryxminIs the minimum value of the x coordinate of the observation system, GeometryxmaxIs the maximum value of the x coordinate of the observation system, GeometryyminBeing the minimum of the y-coordinate of the observation system, GeometryymaxIs 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，Cellhight}. Wherein, CellwidthCell, the width of the minimum grid of the observation system after griddinghightIs the height of the minimum grid of the gridded observation system.
S202: the obstacle region and the safe 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):
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, ObstaclenumNumber of obstacles, U is union operation, ObstacleiIs the ith obstacle, GridObstacle(x,y)∈ObstacleiThe grid (x, y) representing the barrier region after gridding belongs to the i-th barrier region.
For an obstacle with a 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 performing the gridding process on such an obstacle region, it is necessary to include a safe distance region, and the gridding process on the obstacle region and the safe distance region is as shown in formula (2):
wherein, ObsSafenumObsSafe as the number of obstacles having a safe distanceiFor the ith obstacle with a safe distance,
for an open operation, S is a structural operator,
is a barrier region, DsIs 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):
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＝{(TargetRectxmin，TargetRectymin)，
(TargetRectxmax，TargetRectymax)} (4)
wherein, TargetselectxminIs the minimum value of the x coordinate of the minimum rectangle surrounding the kernel, TargetRefxmaxTargetRef, the maximum x-coordinate value of the minimum rectangle bounding the kernelyminTargetRef, the minimum value of the y coordinate of the minimum rectangle bounding the kernelymaxThe 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 step S204 specifically 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-geophone point range (SR) responding 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 regionxmin，SRymin)，(SRxmax，SRymax) As shown in equation (5):
wherein，SRxminIs the minimum value of the x coordinate, SR, of the shot point rangeyminIs the minimum value of the y coordinate, SR, of the shot point rangexmaxMaximum value of x coordinate, SR, of shot point rangeymaxIs the maximum value of the y coordinate of the shot point range, { (SR)xmin,SRymin)，(SRxmax,SRymax) Is the shot check point range, CellwidthCell, being the width of the smallest grid within the smallest rectangular bounding kernelhightThe method comprises the steps of enclosing the height of a minimum grid in a kernel by a minimum rectangle, wherein PathWidth is the width of an arrangement slice of a target area, PathHeight is the height of the arrangement slice of the target area, LevelDept is the depth of a 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-checking pairs existing in the shot-checking point range to generate a shot-checking point range to be encrypted and a check point range corresponding to the shot-checking point range to be encrypted.
In specific implementation, finding shot-checking point range { (SR) for observation systemxmin，SRymin)，(SRxmax，SRymax) The shot-check pairs present in (1) } 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):
wherein, GridShot (x, y) represents the minimum grid unit corresponding to the shot point after gridding,
indicating that GridShot (x, y) is an added shot in the observation system,
denotes GridShot (x, y) as a checkpoint in the viewing system, GridObstacle(x,y)1 denotes that GridShot (x, y) is in the obstacle region, 1 denotes a shot to be added, and 0 denotes an unadditionable shot. x is the number ofAnd x-axis coordinate values representing GridShot (x, y), and y-axis coordinate values representing GridShot (x, y).
And gridding the detection points within the range of the shot detection points, as shown in formula (7):
wherein GridReceiver (x, y) represents the minimum grid unit corresponding to the checked point after gridding,
and if the minimum grid unit corresponding to the gridded detection point belongs to the detection point existing in the system, marking the grid receiver (x, y) corresponding to the detection point as 1, otherwise, marking 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 step S205 specifically executes 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 a reference point gridrreceiver (rx, ry) as a reference and gridtargetreion (x, y) as a center, calculating a shot GridShot (sx, sy) to be added:
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 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. 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 point GridShot (sx, sy) 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 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):
wherein, (sx, sy) is the position information of a shot to be added in the shot range to be encrypted,
for the bin attribute contribution, BinProperty, of the shot point GridShot (sx, sy) to be added corresponding to the target areanFor the nth bin attribute, BinPropertYFun, corresponding to the target region to be added shot GridShot (sx, sy)nIs the surface element attribute function of the target area, targetRegion is the target area, and ^ integralTargetRegionBinPropertyFunn(x, y) dxd 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 BinPropertynIf the number of covering the surface element is equal, the function BinPropertyFun is calculated by the surface elementnCalculating the bin attribute tributary value BinVal (sx, sy) of the shot point GridShot (sx, sy) to be added to the target area, and obtaining the bin attribute tributary value BinVal (sx, sy) as shown in the formula (10):
wherein BinVal (sx, sy) represents the surface element attribute contribution value of the shot point GridShot (sx, sy) to be added to the target area, and anIs 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 gray 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 diagram before shot point encryption, and fig. 10 is a bin calculation effect diagram after shot point encryption, and it can be known through comparison of the two drawings 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 barrier target area, the invention also provides a shot point encryption system based on the barrier target area, as described in the following embodiments. 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 repeated.
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 unit 101, surrounding core generation unit 102, discretization unit 103, contribution value generation unit 104, and encryption unit 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.
And the surrounding kernel generating unit 102 is configured to perform gridding on the target area generated by the limited offset 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 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.
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 obstacle 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 system based on an obstacle target area, which 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, there may be variations in the specific embodiments and the application scope, 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 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.
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-checking pairs existing in the shot-checking point range to generate the shot-checking point range to be encrypted and a check point range corresponding to the shot-checking 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:
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; 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:
wherein, SRxminIs the minimum value of the x coordinate, SR, of the shot point rangeyminIs the minimum value of the y coordinate, SR, of the shot point rangexmaxIs the maximum value of the x coordinate, SR, of the shot point rangeymaxFor the maximum y-coordinate value of the shot point range, { (SR)xmin，SRymin)，(SRxmax，SRymax) Is the shot check point range, CellwidthCell for the width of the minimum grid within the minimum rectangular bounding kernelhightThe height of the minimum grid in the minimum rectangle surrounding core is defined, PathWidth is the width of the arrangement piece of the target area, PathHeight is the height of the arrangement piece 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 x-axis coordinate value of the minimum grid in the minimum rectangle surrounding core, and y represents the y-axis coordinate value 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:
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, and (rx, ry) is a second grid index of any point in the point detection 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:
wherein, (Sx, sy) is the position information of a shot to be added in the range of the shot to be encrypted,
BinProperty corresponding to the bin attribute contribution degree of the target area for the shot point to be addednBinPropertyFun as the nth surface element attribute of the shot to be added corresponding to the target areanIs the surface element attribute function of the target area, targetRegion is the target area, and ^ integralTargetRegionBinPropertyFunn(x, y) dxd 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:
wherein BinVal (sx, sy) is a bin attribute contribution value of a shot point to be added corresponding to the target area, anIs 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 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 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 obstacle area in the shot point range to be encrypted according to the size of the attribute contribution value of each bin.
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 true CN111273339A (en) | 2020-06-12 |
| CN111273339B 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) |
Cited By (1)
| 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 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120020186A1 (en) * | 2009-04-16 | 2012-01-26 | Landmark Graphics Corporation | Seismic Imaging Systems and Methods Employing a Fast Target-Oriented Illumination Calculation |
| CN102540254A (en) * | 2010-12-29 | 2012-07-04 | 中国石油天然气集团公司 | Method for determining effective excitation region based on lighting analysis |
| CN103543465A (en) * | 2012-07-12 | 2014-01-29 | 中国石油天然气集团公司 | Method for determining effective excitation points on basis of illumination energy of target strata |
| US20140226438A1 (en) * | 2013-02-13 | 2014-08-14 | Wireless Seismic, Inc. | Assigned scheduled acquisition process in wireless exploration |
| CN104076387A (en) * | 2014-06-06 | 2014-10-01 | 中国石油集团东方地球物理勘探有限责任公司 | Shallow sea obstacle zone OBC seismic focus acquisition construction method |
| CN104502962A (en) * | 2014-12-16 | 2015-04-08 | 中国石油天然气集团公司 | Method and device for designing shot points |
| CN104502956A (en) * | 2014-12-16 | 2015-04-08 | 中国石油天然气集团公司 | Handling method for avoiding barrier |
| CN105487106A (en) * | 2014-09-18 | 2016-04-13 | 中国石油化工股份有限公司 | Method for supplementing shotpoint based on energy illumination on target stratum by Gaussian beams |
| CN105607115A (en) * | 2015-12-21 | 2016-05-25 | 中国石油天然气集团公司 | Shot point arrangement method and device |
-
2018
- 2018-12-04 CN CN201811471763.9A patent/CN111273339B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120020186A1 (en) * | 2009-04-16 | 2012-01-26 | Landmark Graphics Corporation | Seismic Imaging Systems and Methods Employing a Fast Target-Oriented Illumination Calculation |
| CN102540254A (en) * | 2010-12-29 | 2012-07-04 | 中国石油天然气集团公司 | Method for determining effective excitation region based on lighting analysis |
| CN103543465A (en) * | 2012-07-12 | 2014-01-29 | 中国石油天然气集团公司 | Method for determining effective excitation points on basis of illumination energy of target strata |
| US20140226438A1 (en) * | 2013-02-13 | 2014-08-14 | Wireless Seismic, Inc. | Assigned scheduled acquisition process in wireless exploration |
| CN104076387A (en) * | 2014-06-06 | 2014-10-01 | 中国石油集团东方地球物理勘探有限责任公司 | Shallow sea obstacle zone OBC seismic focus acquisition construction method |
| CN105487106A (en) * | 2014-09-18 | 2016-04-13 | 中国石油化工股份有限公司 | Method for supplementing shotpoint based on energy illumination on target stratum by Gaussian beams |
| CN104502962A (en) * | 2014-12-16 | 2015-04-08 | 中国石油天然气集团公司 | Method and device for designing shot points |
| CN104502956A (en) * | 2014-12-16 | 2015-04-08 | 中国石油天然气集团公司 | Handling method for avoiding barrier |
| CN105607115A (en) * | 2015-12-21 | 2016-05-25 | 中国石油天然气集团公司 | Shot point arrangement method and device |
Non-Patent Citations (6)
| Title |
|---|
| WEI TIE，等: "Some Technologies for Improving Seismic Acquisition Data Quality in Complicated Area", 《SPG/SEG BEIJING 2016 INTERNATIONAL GEOPHYSICAL CONFERENCE》 * |
| ZHAO HU，等: "An automatical infill shot method for uniform imaging of target layer", 《APPLIED GEOPHYSICS》 * |
| 倪宇东，等: "基于实际数据地震道局部相似属性的炮点设计新方法", 《CPSSEG北京2018国际地球物理会议暨展览电子论文集》 * |
| 刘远志，等: "城区特大型障碍物变观三维观测系统设计及应用-以德阳三维为例", 《复杂油气藏》 * |
| 杨光大，等: "乾安城区特殊观测系统设计及应用", 《石油地球物理勘探》 * |
| 董淑兰: "春光探区新层系地震采集观测系统设计研究", 《石油地质与工程》 * |
Cited By (1)
| 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111273339B (en) | 2022-12-02 |
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 | |
| Salvini et al. | The use of an unmanned aerial vehicle for fracture mapping within a marble quarry (Carrara, Italy): photogrammetry and discrete fracture network modelling | |
| CN104376595B (en) | A kind of three-dimensional road generation method cooperateed with based on airborne LiDAR and GIS | |
| Soliman et al. | Effects of vertical accuracy of digital elevation model (DEM) data on automatic lineaments extraction from shaded DEM | |
| Carrea et al. | Automatic rockfalls volume estimation based on terrestrial laser scanning data | |
| CN103941286A (en) | Rapid three-dimensional horizon interpretation method based on horizontal navigation | |
| CN111539904B (en) | Disaster vulnerability prediction method based on rainfall | |
| CN114547866B (en) | Prefabricated part intelligent detection method based on BIM-unmanned aerial vehicle-mechanical dog | |
| CN108647826B (en) | Construction method and device of soil heavy metal environmental risk prediction model | |
| Fanos et al. | Multi-scenario rockfall hazard assessment using LiDAR data and GIS | |
| Wang et al. | Fragmentation calculation method for blast muck piles in open-pit copper mines based on three-dimensional laser point cloud data | |
| Chen et al. | A novel image-based approach for interactive characterization of rock fracture spacing in a tunnel face | |
| Xu et al. | Comprehensive identification of statistical homogeneity of fractured rock masses for a candidate HLW repository site, China | |
| Zhao et al. | Combining ICESat-2 photons and Google Earth Satellite images for building height extraction | |
| CN111273339B (en) | Shot point encryption method and system based on barrier target area | |
| CN113537141B (en) | Method and system for rapidly detecting piping and landslide diseases of dykes and dams | |
| Wu et al. | Stability analysis of rock blocks around a cross-harbor tunnel using the improved morphological visualization method | |
| Xue et al. | Rock mass fracture maps prediction based on spatiotemporal image sequence modeling | |
| Miky et al. | A combined contour lines iteration algorithm and Delaunay triangulation for terrain modeling enhancement | |
| Zhan et al. | A stepwise approach for 3D fracture intersection analysis and application to a hydropower station in Southwest China | |
| Sharifi et al. | Efficiency evaluating of automatic lineament extraction by means of remote sensing (Case study: Venarch, Iran) | |
| CN116226640A (en) | Bad geological feature extraction method for tunnel geophysical prospecting inversion data | |
| Zhang et al. | A novel approach to investigating 3d fracture Connectivity in Ultrahigh steep rock slopes | |
| Xiao et al. | Visualisation of the characteristics of mining subsidence with the consideration of topography: a case study in China | |
| 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 |