Method and device for determining offset of local target, storage medium and electronic equipment

Abstract

The invention relates to a local target offset determination method and device, a storage medium and electronic equipment. The method comprises the following steps: acquiring grid parameters and speed of a three-dimensional speed grid model; acquiring a position parameter of the local target; calculating an average speed in each depth direction based on the acquired grid parameters and speeds; based on the average speed in each depth direction, acquiring a ground surface horizontal offset list of the incident angle of the local target; an offset corresponding to the position parameter is obtained from the earth surface horizontal offset list, and the earth surface emergent position of the local target is obtained based on the position parameter and the offset; the method comprises the steps of obtaining a grid parameter, determining a ground surface position point based on the grid parameter, obtaining a reflection position of the ground surface position point to a local target, constructing an opening angle isogram and a geophone offset isogram of the local target based on the ground surface position point, the reflection position and a ground surface emergent position, and determining the geophone offset of the local target based on the opening angle isogram and the geophone offset isogram. The accuracy of determining the offset can be improved.

Description

Technical Field

[0001] The present invention relates to the technical field of geophysical exploration, and particularly to a method, device, storage medium and electronic device for determining the shot-receiver offset of a local target. Background Art

[0002] Seismic data acquisition is the first process in the oil and gas seismic exploration project. Using seismic geophones and seismic exploration instruments, a complete seismic data acquisition function is realized. Among them, determining an appropriate shot-receiver offset is a very important process in seismic data acquisition. With the development of seismic exploration technology, the requirements for seismic exploration are also getting higher and higher. In recent years, the prestack depth migration imaging technology has been widely applied in the field of seismic imaging. Although the development of seismic migration imaging calculation is rapid, for complex structures, especially the accurate characterization of the salt dome boundary is still a difficult problem, because the shape of the salt dome structure is no longer a simple fault fracture. Therefore, specifically how to determine the shot-receiver offset to improve the accuracy of the determined shot-receiver offset is particularly important for achieving full-target imaging. Summary of the Invention

[0003] In view of this, the present invention provides a method, device, storage medium and electronic device for determining the shot-receiver offset of a local target.

[0004] Specifically, the present invention is realized through the following technical solutions:

[0005] According to the first aspect of the present invention, a method for determining the shot-receiver offset of a local target is provided. The method for determining the shot-receiver offset of a local target includes:

[0006] Obtain the grid parameters of the pre-constructed three-dimensional velocity grid model and the velocities corresponding to the grid parameters;

[0007] Obtain the position parameters of the local target preset in the three-dimensional velocity grid model;

[0008] Based on the obtained grid parameters and the velocities corresponding to the grid parameters, calculate the average velocity in each depth direction of the three-dimensional velocity grid model;

[0009] Based on the average velocity in each depth direction, obtain the list of surface horizontal offsets corresponding to the incident angle of the local target;

[0010] According to the position parameters of the local target, obtain the offset corresponding to the position parameters from the list of surface horizontal offsets, and based on the position parameters and the offset of the local target, obtain the surface emergence position of the local target;

[0011] Based on the grid parameters, determine the surface position points, obtain the reflection positions of the local target for these surface position points, and construct the angular aperture equal-value map and the offset equal-value map of the local target based on the surface position points, the reflection positions, and the surface emergence positions. Based on the angular aperture equal-value map and the offset equal-value map, determine the offset of the local target.

[0012] In the method for determining the offset of the local target in this technical solution, by obtaining the grid parameters, the velocities corresponding to the grid parameters, and the position parameters of the local target, based on the obtained grid parameters and the velocities corresponding to the grid parameters, calculate the average velocity in each depth direction; based on the average velocity in each depth direction, obtain the list of surface horizontal offsets corresponding to each incident angle of the local target; according to the position parameters of the local target, obtain the offset from the list of surface horizontal offsets, so as to obtain the surface emergence position of the local target based on the position parameters and the offset of the local target; then based on the grid parameters, determine the surface position points, obtain the reflection positions of the local target for these surface position points, construct the angular aperture equal-value map and the offset equal-value map of the local target based on the surface position points, the reflection positions, and the surface emergence positions, and based on the angular aperture equal-value map and the offset equal-value map, determine the offset of the local target. In this way, a method for specifically determining the offset based on the local target is proposed. By using the average velocity in each depth direction to obtain the surface horizontal offsets of each incident angle of the local target, the incident angles corresponding to different surface horizontal offsets can be accurately obtained, so that the offset calculation based on the accurately determined incident angles can improve the accuracy of the determined offset, and thus full-target imaging can be realized.

[0013] According to the second aspect of the present invention, there is provided an apparatus for determining the offset of a local target, and the apparatus for determining the offset of a local target includes:

[0014] A model parameter acquisition module, configured to acquire the grid parameters of a pre-constructed three-dimensional velocity grid model and the velocities corresponding to the grid parameters;

[0015] A target parameter acquisition module, configured to acquire the position parameters of a local target pre-set in the three-dimensional velocity grid model;

[0016] A depth velocity acquisition module, configured to calculate the average velocity in each depth direction in the three-dimensional velocity grid model based on the acquired grid parameters and the velocities corresponding to the grid parameters;

[0017] An offset calculation module, configured to obtain the list of surface horizontal offsets corresponding to the incident angles of the local target based on the average velocity in each depth direction;

[0018] An exit position determination module, configured to obtain an offset corresponding to the position parameter from a list of surface horizontal offsets according to the position parameter of the local target, and obtain the surface exit position of the local target based on the position parameter and the offset of the local target;

[0019] An offset determination module, configured to determine a surface position point based on the grid parameters, obtain a reflection position of the surface position point with respect to the local target, construct an aperture angle equal value map and an offset equal value map of the local target based on the surface position point, the reflection position and the surface exit position, and determine the offset of the local target based on the aperture angle equal value map and the offset equal value map.

[0020] According to a third aspect of the present invention, there is provided a storage medium, on which a computer program is stored, and when the program is executed by a processor, the steps of the method for determining the offset of a local target in any possible implementation manner of the first aspect are implemented.

[0021] According to a fourth aspect of the present invention, there is provided an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, and when the processor executes the program, the steps of the method for determining the offset of a local target in any possible implementation manner of the first aspect are implemented. Description of the Drawings

[0022] The drawings here are incorporated into the specification and form a part of the specification, showing embodiments consistent with the present invention, and are used together with the specification to explain the principles of the present invention.

[0023] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or related technologies. Obviously, for those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0024] Figure 1 It is a schematic flow chart of a method for determining the offset of a local target provided by an embodiment of the present invention;

[0025] Figure 2 It is a schematic diagram of surface horizontal offsets at different incident angles in a method for determining the offset of a local target provided by an embodiment of the present invention;

[0026] Figure 3 It is a schematic diagram of incident angles corresponding to different offset intervals in a method for determining the offset of a local target provided by an embodiment of the present invention;

[0027] Figure 4 It is a schematic diagram for determining the reflection position in a method for determining the offset of a local target provided by an embodiment of the present invention;

[0028] Figure 5 Schematic diagram of an offset determination device for local targets provided by an embodiment of the present invention;

[0029] Figure 6 Schematic diagram of the structure of an electronic device provided by an embodiment of the present invention. Detailed implementation manners

[0030] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the described embodiments are some but not all of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0031] For a work area with complex structures, specifically how to determine the offset to improve the accuracy of the determined offset and achieve full-target imaging plays an important role in seismic data acquisition. In this embodiment, based on angle-domain imaging, by decomposing the work area with complex structures into different local targets, according to the requirements of the local targets for the parameters of the seismic exploration acquisition and observation system, the optimized offset parameters and the effective surface receiving range are obtained, so as to perform accurate imaging on the local targets such as different sides and top surfaces of the work area with complex structures, achieve full-target imaging, and provide method and technical support for full-target seismic exploration. It can be used for the acquisition design of the seismic exploration acquisition and observation system for work areas with complex geological structure targets, achieve the purpose of optimizing the seismic exploration acquisition and observation system, be applicable to the seismic exploration acquisition and observation system for three-dimensional complex geology full targets, provide a method for optimizing the offset point positions of the seismic exploration acquisition and observation system, and provide reasonable parameters for the design of the seismic exploration acquisition and observation system.

[0032] See Figure 1 , an embodiment of the present invention provides a method for determining the offset of local targets. This embodiment is applicable to determining the offset parameters for a three-dimensional grid geological model (three-dimensional velocity grid model) with an existing attribute of velocity, and determining the offset based on the depth velocity of the three-dimensional velocity grid model and the spatial position of the local target. The method may include the following steps:

[0033] S101. Obtain the grid parameters of the pre-constructed three-dimensional velocity grid model and the velocities corresponding to the grid parameters;

[0034] In this embodiment, the three-dimensional velocity grid model is constructed based on the work area where seismic data is to be collected. As an optional embodiment, the grid parameters include but are not limited to: the number of sampling points, the grid spacing of the sampling points. Taking the number of sampling points as an example, the number of sampling points includes the number of sampling points in the x direction, the number of sampling points in the y direction, and the number of sampling points in the z direction.

[0035] In this embodiment, let the number of sampling points in the x - direction of the three - dimensional velocity grid model be Nx, and the grid spacing of the sampling points in the x - direction be dx; the number of sampling points in the y - direction be Ny, and the grid spacing of the sampling points in the y - direction be dy; the number of sampling points in the z - direction be Nz, and the grid spacing of the sampling points in the z - direction be dz, where the z - direction is the depth direction.

[0036] In this embodiment, taking the three - dimensional velocity grid model of a certain work area as an example, the number of sampling points in the x - direction of the constructed three - dimensional velocity grid model is set to Nx = 1057, the grid spacing is dx = 30m; the number of sampling points in the y - direction is Ny = 571, the grid spacing is dy = 30m; the number of sampling points in the z - direction is Nz = 480, and the grid spacing is dz = 30m.

[0037] S102. Obtain the position parameters of the pre - set local target in the three - dimensional velocity grid model;

[0038] In this embodiment, the local target is one of the sampling points. As an alternative embodiment, it is a sampling point on the side, top surface, etc. of the pre - deployed three - dimensional velocity grid model underground. For the specific process of determining the local target, reference can be made to relevant technical literature, which is not elaborated here.

[0039] In this embodiment, as an alternative embodiment, the position parameters include but are not limited to: position coordinates, azimuth along the inclined direction, and dip angle along the inclined direction. Among them, the position coordinates of the local target are denoted as o(x, y, z), with the unit of meter, the azimuth of the local target along the inclined direction is Azimuth, and the dip angle is Dip, both with the unit of degree.

[0040] In this embodiment, by setting the local target, the position of the local target in the three - dimensional velocity grid model is determined. For example, as an alternative embodiment, the determined position coordinates of the local target are (x, y, z)=(5200m, 5200m, 4800m), the azimuth Azimuth of the local target along the inclined direction is 60 degrees, and the dip angle Dip is 15 degrees.

[0041] S103. Calculate the average velocity in each depth direction of the three - dimensional velocity grid model based on the obtained grid parameters and the velocities corresponding to the grid parameters;

[0042] In this embodiment, based on the velocities of the grid parameters at each depth, the average velocities at different depths of the three - dimensional velocity grid model are obtained, where the depth can be determined according to the number of faults included in the three - dimensional velocity grid model. As an alternative embodiment, calculating the average velocity in each depth direction of the three - dimensional velocity grid model based on the obtained grid parameters and the velocities corresponding to the grid parameters includes:

[0043] A11. For each target depth, obtain the velocities of the sampling points located at that target depth.

[0044] In this embodiment, at each target depth, a plurality of sampling points are arranged. The number of target depths is equal to the quotient of the total depth divided by the grid spacing in the z direction.

[0045] A12. Calculate the sum value of the velocities of the sampling points.

[0046] In this embodiment, for the sampling points at the target depth, sum up the velocities of the sampling points.

[0047] A13. Obtain the number of sampling points located at that target depth.

[0048] A14. Calculate the quotient of the sum value and the number of sampling points to obtain the average velocity of the target depth.

[0049] In this embodiment, as an alternative embodiment, calculate the average velocity of the target depth using the following formula:

[0050] V_average[iz] = Sum(V[ix,iy,iz]) / (Nx * Ny)

[0051] In the formula, V[ix,iy,iz] represents the velocity of the sampling point (ix,iy,iz), with the unit of meters per second. Sum represents summation. ix is the index of the target depth in the x direction, with the starting index of 1 and the ending index of Nx. iy is the index of the target depth in the y direction, with the starting index of 1 and the ending index of Ny. iz is the target depth in the z direction, with the starting index of 1 and the ending index of Nz. / represents the division sign, and * represents the multiplication sign.

[0052] In this embodiment, the sampling points are arranged in a matrix to form a matrix of Nx * Ny.

[0053] In this embodiment, as an alternative embodiment, the starting index of ix is 1 and the ending index is 1057, the starting index of iy is 1 and the ending index is 571, and the starting index of iz is 1 and the ending index is 480. That is, for each target depth iz, there are 1057 * 571 sampling points arranged. In the depth direction, there are 480 target depths.

[0054] By calculating the velocities of the sampling points at each target depth, the average velocities at different depths are obtained as follows:

[0055] V_average[1] = 2098.8, V_average[2] = 2098.8, V_average[3] = 2098.8, …, V_average

[15] = 2099.3, …, V_average

[480] = 4123.0。

[0056] S104. Obtain a list of surface horizontal offsets corresponding to the incident angles of the local target based on the average velocities in each depth direction;

[0057] In this embodiment, for the local target, based on the average velocities in each depth direction, calculate the surface horizontal offsets corresponding to any incident angle (iangle) of the local target, and generate a list of surface horizontal offsets based on the surface horizontal offsets corresponding to each incident angle.

[0058] In this embodiment, as an alternative embodiment, obtaining a list of surface horizontal offsets corresponding to the incident angles of the local target based on the average velocities in each depth direction includes:

[0059] B11. Determine the starting index and the maximum index of the incident angle, and determine the number of depth layers of the local target based on the depth direction parameter value in the position parameters of the local target and the depth direction grid spacing in the grid parameters;

[0060] In this embodiment, each incident angle (iangle) corresponds to a surface horizontal offset, and the surface horizontal offsets of each incident angle form a list of surface horizontal offsets: OFFset(iangle).

[0061] In this embodiment, as an alternative embodiment, the starting index of the incident angle iangle is 1, the maximum index is 89, the unit is degree, and each index value corresponds to one degree.

[0062] In this embodiment, taking the above as an example, the depth direction parameter value in the position parameters of the local target is the z-direction coordinate value in the position coordinates (x, y, z) = (5200m, 5200m, 4800m), that is, 4800m, and the grid spacing dz of the sampling points in the z direction is 30m, then the number of depth layers of the local target is 4800 / 30 = 160.

[0063] B12. Traverse in the order of the incident angle from the starting index to the maximum index. For each incident angle, obtain the surface horizontal offset of the incident angle based on the average velocity and the depth direction grid spacing corresponding to the number of depth layers from the surface to the local target;

[0064] In this embodiment, as an alternative embodiment, based on the average velocity and the grid spacing in the depth direction respectively corresponding to the depth layers from the ground surface to the local target, obtaining the horizontal offset of the ground surface of this incident angle includes:

[0065] C11, obtaining the sine value of this incident angle;

[0066] In this embodiment, the sine value of the incident angle is calculated using the following formula:

[0067] sintheata1 = sin(iangle / 180.0 * PI);

[0068] In the formula, iangle is the degree of the incident angle.

[0069] C12, based on the sine value of the incident angle and the grid spacing in the depth direction, calculating the horizontal offset of the ground surface of the first layer among the depth layers from the ground surface to the local target direction;

[0070] In this embodiment, the depth layer number is 160, indicating 160 layers. From the ground surface to the local target direction, they are respectively denoted as the first layer, the second layer,..., the one hundred and sixtieth layer.

[0071] In this embodiment, as an alternative embodiment, the horizontal offset of the ground surface of the first layer is calculated using the following formula:

[0072] OFFSET = dz * sintheata1 / sqrt(1 - sintheata1 * sintheata1)

[0073] In the formula, sqrt represents taking the square root.

[0074] C13, for the other layers after the first layer, based on the average velocity corresponding to the current layer, the average velocity corresponding to the previous layer of the current layer, and the sine value of the incident angle, obtaining the current layer coefficient;

[0075] In this embodiment, as an alternative embodiment, the current layer coefficient is obtained using the following formula:

[0076] sintheata2 = V2 * sintheata1 / v1; where,

[0077] V2 = V_average[iz + 1]

[0078] In the formula, iz + 1 is the current layer, sintheata2 is the current layer coefficient, V2 is the average velocity of the current layer, v1 is the average velocity of the previous layer of the current layer. For example, V_average[1] = 2098.8, V_average[2] = 2098.8.

[0079] C14. Calculate the surface horizontal offset of the current layer based on the surface horizontal offset of the previous layer, the grid spacing in the depth direction, and the coefficient of the current layer, until the current layer is the layer corresponding to the depth layer number, to obtain the surface horizontal offset at this incident angle.

[0080] In this embodiment, when iz is not equal to 1, calculate sintheta2 = V2 * sintheta1 / v1, and make: sintheta1 = sintheta2. Use the following formula to calculate the surface horizontal offset of the current layer:

[0081] OFFSET = OFFSET + dz * sintheta1 / sqrt(1 - sintheta1 * sintheta1).

[0082] In this embodiment, for each incident angle, for example, iangle = 1, 2,..., 89, perform steps C11 to C14 respectively to obtain the surface horizontal offset at each incident angle. By successively accumulating the surface horizontal offsets of each layer, the surface horizontal offset of the last layer is obtained, and the surface horizontal offset of the last layer is the surface horizontal offset at this incident angle.

[0083] In this embodiment, through the above calculations, the surface horizontal offsets at each incident angle are respectively: OFFset(1) = 64.1m, OFFset(2) = 128.2m, OFFset(3) = 192.3m,..., OFFset(88) = 8288.6m, OFFset(89) = 8651.1m. Among them, OFFset(i) is the surface horizontal offset at the incident angle index i.

[0084] In this embodiment, the incident angle index i represents that the incident angle is i degrees.

[0085] B13. Construct the surface horizontal offset list based on the surface horizontal offsets at each incident angle.

[0086] Figure 2 This is a schematic diagram of the surface horizontal offsets at different incident angles in a method for determining the shot-receiver offset of a local target provided by an embodiment of the present invention. In this embodiment, with the horizontal direction as the surface horizontal offset and the vertical direction as the depth parameter value, directly below the vertical where the surface horizontal offset in the horizontal direction is 0 is the local target, and a surface horizontal offset list is constructed.

[0087] In this embodiment, the incident angles can also be refined to obtain a more refined surface horizontal offset list. Therefore, as an optional embodiment, the method further includes:

[0088] For the list of surface horizontal offset amounts, perform segmentation according to a preset offset interval to obtain a plurality of segmentation values, and use an interpolation method to obtain the incident angle corresponding to each segmentation value, so as to update the list of surface horizontal offset amounts.

[0089] In this embodiment, the surface horizontal offset amounts corresponding to different incident angles are segmented according to a preset offset interval, such as 1 m, to obtain a plurality of segmentation values, and an interpolation method is used to obtain the incident angle corresponding to each segmentation value.

[0090] In this embodiment, as an alternative embodiment, the preset interval is 1 meter. The surface horizontal offset amounts (offsets) corresponding to different incident angles are evenly divided into intervals of 1 m. For example, the surface horizontal offset amount of the maximum incident angle is segmented according to 1 m, and the incident angle corresponding to each 1 m interval is obtained by an interpolation method, and the list of surface horizontal offset amounts is updated to construct a list of incident angles: Angle(ioffset), with the unit of degree, where ioffset is the offset interval.

[0091] In this embodiment, ioffset is characterized by an accumulated index, where the starting index is 1 and the maximum index is Noffset. In this way, the number of incident angles can be effectively expanded and the accuracy of the incident angle can be improved.

[0092] In this embodiment, as an alternative embodiment, the following formula is used to obtain the maximum index of the offset interval:

[0093] Noffset = CELL(OFFset(89)), where CELL() represents rounding up.

[0094] In this embodiment, as an alternative embodiment, the starting index of ioffset is 1 and the maximum index is 8651.

[0095] Figure 3 It is a schematic diagram of the incident angles corresponding to different offset intervals in a method for determining the shot-receiver offset of a local target provided by an embodiment of the present invention. In this embodiment, the list of incident angles corresponding to different offset intervals calculated by the interpolation method is as follows:

[0096] Angle(1) = 0.0156, Angle(2) = 0.0312, Angle(3) = 0.0468, Angle(4) = 0.0624, Angle(5) = 0.0780,..., Angle(8650) = 87.9957, Angle(8651) = 87.9978. Among them, Angle(j) represents the incident angle corresponding to the jth offset interval. For example, when the offset interval is 1 m, Angle(8650) represents that the incident angle corresponding to the surface horizontal offset amount of 8650 m is 87.9957 degrees.

[0097] S105. According to the position parameter of the local target, obtain the offset corresponding to the position parameter from the list of surface horizontal offsets, and based on the position parameter and offset of the local target, obtain the surface emergence position of the local target;

[0098] In this embodiment, for the updated list of surface horizontal offsets, obtain the offset corresponding to the dip angle from the list of incident angles. According to the dip angle of the local target along the inclined direction, search in the list of incident angles to obtain the surface horizontal offset (offset) L1 corresponding to the dip angle, and translate the surface position (x, y) of the local target along the azimuth direction by a distance L1 to obtain the surface emergence position T of the local target, that is, the position point where the local target vertically emerges on the ground. Among them, the dip angle of the local target along the inclined direction corresponds to the incident angle, and the magnitude of the dip angle is the magnitude of the incident angle.

[0099] In this embodiment, as an alternative embodiment, obtaining the surface emergence position of the local target based on the position parameter and offset of the local target includes:

[0100] Translate the surface position of the local target along the azimuth direction by the offset to obtain the surface emergence position of the local target.

[0101] In this embodiment, as described above, according to the dip angle of the local target being 15 degrees, search in the list of incident angles, and there is OFFset(15)=971.28m, obtain the offset L1 corresponding to the dip angle as 971.28, and translate the surface position (x, y) of the local target along the azimuth direction by a distance of 971.28, that is, obtain the surface emergence position T of the local target, and the coordinates of point T are (6040.91m, 5685.5m, 0m).

[0102] S106. Based on the grid parameters, determine the surface position points, obtain the reflection position of the surface position points with respect to the local target, construct the opening angle contour map and offset contour map of the local target based on the surface position points, reflection positions, and surface emergence positions, and determine the offset of the local target based on the opening angle contour map and offset contour map.

[0103] In this embodiment, according to the grid spacing of the sampling points in the x direction and the grid spacing of the sampling points in the y direction in the grid parameters, select the surface position points S(ix, iy) of the three-dimensional velocity grid model. The surface position points S(ix, iy) are the number of sampling points on the surface, that is, Nx*Ny. For each surface position point S(ix, iy), obtain the reflection position R(ixr, iyr) of the surface position point S(ix, iy) with respect to the underground local target O, and then calculate the local opening angle openAngle of the local target.

[0104] In this embodiment, as an alternative embodiment, obtaining the reflection position of the surface location point with respect to the local target includes:

[0105] D11, connecting the surface location point and the subsurface local target, and connecting the surface location point and the surface emergence position and extending it;

[0106] In this embodiment, connect the surface location point S with the subsurface local target o, and connect the surface location point S with the surface emergence position T, and extend it in the direction from the surface location point S to the surface emergence position T.

[0107] In this embodiment, as an alternative embodiment, take ix = 101, iy = 101, and the position coordinates at the corresponding surface location point S(101, 101) are (3000m, 3000m, 0m).

[0108] D12, in the triangle formed by the surface location point, the subsurface local target, and the surface emergence position, obtain the degree of the first included angle between the first straight line formed by the surface location point and the subsurface local target and the second straight line formed by the subsurface local target and the surface emergence position;

[0109] In this embodiment, in the triangle SOT formed by the surface location point S, the subsurface local target o, and the surface emergence position T, find the degree of the angle SOT. Among them, S(3000m, 3000m, 0m), o(5200m, 5200m, 4800m), T(6040.91m, 5685.5m, 0m), and the line segment lengths can be obtained: ST = 4056.64m, SO = 5720.14m, TO = 4897.18m.

[0110] D13, on the extension line, determine the reflection position, and the degree of the second included angle between the third straight line formed by the reflection position and the subsurface local target and the fourth straight line formed by the surface emergence position and the subsurface local target is equal to the first included angle degree.

[0111] Figure 4 It is a schematic diagram for determining the reflection position in a method for determining the offset of a local target provided by an embodiment of the present invention. In this embodiment, take a point R on the extension line of ST such that the angle TOR is equal to the angle SOT, and this R point is the required reflection position R. As described above, according to the cosine theorem, the degree of the angle SOT can be obtained as 44.08 degrees. Take a point R on the extension line of ST such that the angle TOR is equal to the angle SOT, and at this time, the R point is the required reflection position.

[0112] In this embodiment, as an alternative embodiment, based on the surface location point, the reflection position, and the surface emergence position, constructing the angular aperture equal-value map and the offset equal-value map of the local target includes:

[0113] E11. Based on the second included angle degree, obtain the local target opening angle of the local target based on the surface position point.

[0114] In this embodiment, the local target opening angle openAngle is angle TOR, and record the TOR angle at the reflection position R(ixr, iyr) to form an opening angle file.

[0115] E12. Obtain the distance from the surface position point to the reflection position to get the shot-receiver offset.

[0116] In this embodiment, calculate the length of SR, record it at the reflection position R(ixr, iyr) to form a shot-receiver offset file. In this embodiment, the obtained length SR = 9235.83 m, get ixr = 330, iyr = 303, and the local target opening angle openAngle = 44.08 degrees, and record them at the position R(330, 303) to form an opening angle file.

[0117] E13. Based on the local target opening angles of each surface position point, construct the opening angle contour map.

[0118] In this embodiment, calculate the local target opening angles corresponding to the surface position points in all three-dimensional velocity grid models point by point to form an opening angle contour map.

[0119] E14. Based on the shot-receiver offsets of each surface position point, construct the shot-receiver offset contour map.

[0120] In this embodiment, record SR at the position R(330, 303) to form a shot-receiver offset file. Calculate the SR lengths corresponding to the surfaces in all three-dimensional velocity grid models point by point to form a shot-receiver offset contour map.

[0121] In this embodiment, as an alternative embodiment, based on the opening angle contour map and the shot-receiver offset contour map, determining the shot-receiver offset of the local target includes:

[0122] F11. Based on the opening angle contour map, obtain the region where the local target opening angle is within a preset opening angle threshold to get the optimal excitation and reception region of the local target.

[0123] F12. In the shot-receiver offset contour map corresponding to the optimal excitation and reception region, obtain the maximum shot-receiver offset, and determine that the shot-receiver offset of the local target is greater than or equal to the maximum shot-receiver offset.

[0124] In this embodiment, based on the opening angle contour map and the shot-receiver offset contour map, obtain the region where the local target opening angle is within a preset opening angle threshold to get the excitation and reception region, map the excitation and reception region to the shot-receiver offset contour map, obtain the maximum length between each surface position point and the corresponding reflection position in the mapped shot-receiver offset contour map to get the maximum shot-receiver offset, and determine the shot-receiver offset of the local target to be greater than or equal to the maximum shot-receiver offset.

[0125] In this embodiment, as an alternative embodiment, the included angle threshold is 60 degrees. In the included angle contour map, the contour line area with an included angle of 60 degrees is set, and the inner area within 60 degrees is the optimal excitation and reception area for the local target.

[0126] In this embodiment, in order to enable the local target to achieve complete seismic data imaging, it is required that the determined shot-receiver distance >= MaxOffset, where MaxOffset is the maximum length among the SR lengths within the 60-degree area. As an alternative embodiment, the shot-receiver distance for the local target is determined using the following formula:

[0127] MaxOffset = max{SR}, (openAngle < 60 degrees and openAngle > 0 degrees)

[0128] In this embodiment, MaxOffset = max{SR} = 9532.4 m.

[0129] In this embodiment, based on seismic reflection imaging technology, by analyzing the underground local geological target (local target), through the incident reflection principle, the reflection included angle (local target included angle) corresponding to different shot-receiver points (surface position points) is analyzed, and the maximum value of the distance corresponding to the surface shot-receiver points is scanned, thereby obtaining the optimal shot-receiver distance.

[0130] Based on the same inventive concept, as Figure 5 shown, the embodiment of the present invention further provides a device for determining the shot-receiver distance of a local target, the device comprising:

[0131] A model parameter acquisition module 501, configured to acquire grid parameters of a pre-constructed three-dimensional velocity grid model and the velocities corresponding to the grid parameters;

[0132] In this embodiment, as an alternative embodiment, the grid parameters include but are not limited to: the number of sampling points, the grid spacing of the sampling points.

[0133] A target parameter acquisition module 502, configured to acquire the position parameters of a local target pre-set in the three-dimensional velocity grid model;

[0134] In this embodiment, as an alternative embodiment, the position parameters include but are not limited to: position coordinates, azimuth angle along the inclined direction, and inclination angle along the inclined direction.

[0135] A depth velocity acquisition module 503, configured to calculate the average velocity in each depth direction of the three-dimensional velocity grid model based on the acquired grid parameters and the velocities corresponding to the grid parameters;

[0136] In this embodiment, as an alternative embodiment, the depth velocity acquisition module 503 includes:

[0137] A speed acquisition unit (not shown in the figure) for acquiring the speeds of each sampling point located at the target depth for each target depth;

[0138] A calculation unit for calculating the sum value of the speeds of each sampling point;

[0139] A statistics unit for obtaining the number of sampling points located at the target depth;

[0140] An average value acquisition unit for calculating the quotient of the sum value of the speeds and the number of sampling points to obtain the average speed of the target depth.

[0141] An offset calculation module 504 for obtaining a list of surface horizontal offsets corresponding to the incident angle of the local target based on the average speeds in each depth direction;

[0142] In this embodiment, as an alternative embodiment, the offset calculation module 504 includes:

[0143] A layer number determination unit for determining the starting index and the maximum index of the incident angle, and determining the depth layer number of the local target based on the depth direction parameter value in the position parameters of the local target and the depth direction grid spacing in the grid parameters;

[0144] A traversal unit for traversing in the order of the incident angle from the starting index to the maximum index. For each incident angle, based on the average speeds and the depth direction grid spacing corresponding to the depth layer numbers from the surface to the local target, obtain the surface horizontal offset of the incident angle;

[0145] An offset list construction unit for constructing the list of surface horizontal offsets based on the surface horizontal offsets of each incident angle.

[0146] In this embodiment, as an alternative embodiment, the traversal unit is specifically configured to:

[0147] Obtain the sine value of the incident angle;

[0148] Based on the sine value of the incident angle and the depth direction grid spacing, calculate the surface horizontal offset of the first layer in the depth layer numbers in the direction from the surface to the local target;

[0149] For the other layers after the first layer, obtain the current layer coefficient based on the average speed corresponding to the current layer, the average speed corresponding to the previous layer of the current layer, and the sine value of the incident angle;

[0150] Based on the surface horizontal offset of the previous layer, the depth direction grid spacing, and the current layer coefficient, calculate the surface horizontal offset of the current layer until the current layer is the layer corresponding to the depth layer number, to obtain the surface horizontal offset of the incident angle.

[0151] An exit position determination module 505, configured to obtain an offset corresponding to the position parameter from a list of surface horizontal offsets according to the position parameter of the local target, and obtain the surface exit position of the local target based on the position parameter and the offset of the local target;

[0152] In this embodiment, as an alternative embodiment, obtaining the surface exit position of the local target based on the position parameter and the offset of the local target includes:

[0153] Translating the surface position of the local target along the azimuth direction by the offset to obtain the surface exit position of the local target.

[0154] An offset determination module 506, configured to determine a surface position point based on the grid parameter, obtain the reflection position of the surface position point with respect to the local target, construct a dip angle equal-value map and an offset equal-value map of the local target based on the surface position point, the reflection position and the surface exit position, and determine the offset of the local target based on the dip angle equal-value map and the offset equal-value map.

[0155] In this embodiment, as an alternative embodiment, the offset determination module 506 includes:

[0156] An excitation and reception area determination unit, configured to obtain an area where the dip angle of the local target is within a preset dip angle threshold based on the dip angle equal-value map, and obtain the optimal excitation and reception area of the local target;

[0157] An offset determination unit, configured to obtain the maximum offset in the offset equal-value map corresponding to the optimal excitation and reception area, and determine that the offset of the local target is greater than or equal to the maximum offset.

[0158] In this embodiment, as an alternative embodiment, the offset determination module 506 includes:

[0159] A connection unit, configured to connect the surface position point to the underground local target, and connect and extend the surface position point to the surface exit position;

[0160] A first included angle calculation unit, configured to obtain the degree of a first included angle between a first straight line formed by the surface position point and the underground local target and a second straight line formed by the underground local target and the surface exit position in a triangle formed by the surface position point, the underground local target, and the surface exit position;

[0161] A second included angle calculation unit, configured to determine a reflection position on the extension line, and the degree of a second included angle between a third straight line formed by the reflection position and the underground local target and a fourth straight line formed by the surface exit position and the underground local target is equal to the degree of the first included angle.

[0162] In this embodiment, as another optional embodiment, the shot-receiver distance determination module 506 further includes:

[0163] An included angle acquisition unit, configured to acquire a local target included angle of the local target based on the surface position point based on the second included angle degree;

[0164] A shot-receiver distance acquisition unit, configured to acquire the distance from the surface position point to the reflection position to obtain the shot-receiver distance;

[0165] An included angle contour map construction unit, configured to construct the included angle contour map based on the local target included angles of each surface position point;

[0166] A shot-receiver distance contour map construction unit, configured to construct the shot-receiver distance contour map based on the shot-receiver distances of each surface position point.

[0167] In this embodiment, as an optional embodiment, the device further includes:

[0168] An incident angle segmentation module (not shown in the figure), configured to segment the surface horizontal offset amount list according to a preset offset interval to obtain a plurality of segmentation values, and use an interpolation method to obtain the incident angle corresponding to each segmentation value, so as to update the surface horizontal offset amount list.

[0169] In this embodiment, as an optional embodiment, the preset interval is 1 meter, the surface horizontal offset amount of the maximum incident angle is segmented according to 1 m, and the incident angle corresponding to each 1 m interval is obtained by an interpolation method to construct an incident angle list.

[0170] Based on the same inventive concept, an embodiment of the present invention further provides a storage medium, on which a computer program is stored, and when the program is executed by a processor, the steps of the method for determining the shot-receiver distance of the local target in any possible implementation manner described above are implemented.

[0171] Optionally, the storage medium may be a non-transitory computer-readable storage medium. For example, the non-transitory computer-readable storage medium may be a ROM, a random access memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, and an optical data storage device, etc.

[0172] Based on the same inventive concept, see Figure 6, an embodiment of the present invention further provides an electronic device, including a memory 101 (such as a non-volatile memory), a processor 102, and a computer program stored on the memory 101 and executable on the processor 102. When the processor 102 executes the program, it implements the steps of the method for determining the offset of the local target in any of the above possible implementation manners, which is equivalent to the device for determining the offset of the local target as described above. Of course, this processor can also be used to process other data or perform operations. The electronic device can be a device such as a PC, a server, or a terminal.

[0173] As Figure 6 shown, the electronic device generally may further include: a memory 103, a network interface 104, and an internal bus 105. In addition to these components, other hardware may also be included, which will not be elaborated here.

[0174] It should be noted that the above device for determining the offset of the local target can be implemented by software. As a logically meaningful device, it is formed by the processor 102 of the electronic device where it is located reading the computer program instructions stored in the non-volatile memory into the memory 103 and running them.

[0175] The embodiments of the subject matter and functional operations described in this specification can be implemented in the following: digital electronic circuits, tangible computer software or firmware, computer hardware including the structures disclosed in this specification and their structural equivalents, or a combination of one or more of them. The embodiments of the subject matter described in this specification can be implemented as one or more computer programs, that is, one or more modules in computer program instructions encoded on a tangible non-transitory program carrier to be executed by a data processing device or to control the operation of the data processing device. Alternatively or additionally, the program instructions can be encoded on an artificially generated propagated signal, such as a machine-generated electrical, optical, or electromagnetic signal, which is generated to encode information and transmit it to a suitable receiver device for execution by the data processing device. The computer storage medium can be a machine-readable storage device, a machine-readable storage substrate, a random or serial access memory device, or a combination of one or more of them.

[0176] The processes and logical flows described in this specification can be executed by one or more programmable computers executing one or more computer programs to perform corresponding functions by operating on input data and generating outputs. The processes and logical flows can also be executed by dedicated logic circuits - such as FPGAs (Field Programmable Gate Arrays) or ASICs (Application Specific Integrated Circuits), and the device can also be implemented as dedicated logic circuits.

[0177] Computers suitable for executing computer programs include, for example, general and / or special-purpose microprocessors, or any other type of central processing unit. Generally, the central processing unit will receive instructions and data from read-only memory and / or random access memory. The basic components of a computer include a central processing unit for implementing or executing instructions and one or more memory devices for storing instructions and data. Generally, a computer will also include one or more mass storage devices for storing data, such as magnetic disks, magneto-optical disks, or optical disks, etc., or the computer will be operably coupled to such mass storage devices to receive data therefrom or transfer data thereto, or both. However, a computer is not necessarily required to have such devices. In addition, a computer may be embedded in another device, such as a mobile phone, a personal digital assistant (PDA), a mobile audio or video player, a game console, a global positioning system (GPS) receiver, or a portable storage device such as a universal serial bus (USB) flash drive, to name just a few examples.

[0178] Computer-readable media suitable for storing computer program instructions and data include all forms of non-volatile memory, media, and memory devices, such as including semiconductor memory devices (e.g., EPROM, EEPROM, and flash memory devices), magnetic disks (e.g., internal hard disks or removable disks), magneto-optical disks, and CD-ROM and DVD-ROM disks. The processor and memory may be supplemented by, or incorporated in, special logic circuitry.

[0179] Although this specification contains many specific implementation details, these should not be construed as limiting the scope of any invention or the scope of what is claimed, but rather as mainly describing the features of specific embodiments of a particular invention. Certain features described in multiple embodiments in this specification may also be implemented in combination in a single embodiment. On the other hand, various features described in a single embodiment may also be implemented separately in multiple embodiments or in any suitable sub-combination. In addition, although features may act in certain combinations as described above and even be claimed as such initially, one or more features from a claimed combination may in some cases be removed from that combination, and the claimed combination may be directed to a sub-combination or a variation of a sub-combination.

[0180] Similarly, although operations are depicted in the drawings in a particular order, this should not be understood as requiring that the operations be performed in the particular order shown or sequentially, or that all illustrated operations be performed, to achieve the desired result. In some cases, multitasking and parallel processing may be advantageous. In addition, the separation of various system modules and components in the above embodiments should not be understood as required in all embodiments, and it should be understood that the described program components and systems can generally be integrated together in a single software product or packaged into multiple software products.

[0181] Thus, specific embodiments of the subject matter have been described. Other embodiments are within the scope of the appended claims. In some cases, the actions recited in the claims can be performed in a different order and still achieve the desired results. In addition, the processes depicted in the drawings do not necessarily require the particular order or sequential order shown to achieve the desired results. In some implementations, multitasking and parallel processing may be advantageous.

[0182] It should be noted that, in this article, relational terms such as "first" and "second" are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the existence of other identical elements in the process, method, article or device including the elements.

[0183] The above is only a specific embodiment of the present invention, so that those skilled in the art can understand or implement the present invention. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to the embodiments shown herein, but should conform to the widest scope consistent with the principles and novel features applied herein.

Claims

1. A method for determining the offset of a local target, characterized in that Including: Obtaining grid parameters of a pre - constructed three - dimensional velocity grid model and velocities corresponding to the grid parameters; Obtaining position parameters of a preset local target in the three - dimensional velocity grid model; Calculating the average velocity in each depth direction in the three - dimensional velocity grid model based on the obtained grid parameters and velocities corresponding to the grid parameters; Obtaining a list of surface horizontal offsets corresponding to the incident angle of the local target based on the average velocities in each depth direction; Obtaining the offset corresponding to the position parameter from the list of surface horizontal offsets according to the position parameter of the local target, and obtaining the surface emergence position of the local target based on the position parameter and offset of the local target; Determining surface position points based on the grid parameters, obtaining the reflection position of the surface position points with respect to the local target, constructing an opening - angle equal - value map and an offset - distance equal - value map of the local target based on the surface position points, reflection positions and surface emergence positions, and determining the offset - distance of the local target based on the opening - angle equal - value map and the offset - distance equal - value map.

2. The method for determining the offset of the shot point and the receiver for the local target according to claim 1, characterized in that, The calculating the average velocity in each depth direction in the three - dimensional velocity grid model based on the obtained grid parameters and velocities corresponding to the grid parameters includes: For each target depth, obtaining the velocities of each sampling point located at the target depth; Calculating the sum value of the velocities of each sampling point; Obtaining the number of sampling points located at the target depth; Calculating the quotient of the sum value of velocities and the number of sampling points to obtain the average velocity of the target depth.

3. The method for determining the offset of the local target according to claim 1, wherein The obtaining a list of surface horizontal offsets corresponding to the incident angle of the local target based on the average velocities in each depth direction includes: Determining the starting index and maximum index of the incident angle, and determining the number of depth layers of the local target based on the depth - direction parameter value in the position parameter of the local target and the depth - direction grid spacing in the grid parameters; Traversing in the order of the incident angle from the starting index to the maximum index. For each incident angle, obtaining the surface horizontal offset of the incident angle based on the average velocities and depth - direction grid spacing corresponding to the number of depth layers from the surface to the local target; Constructing the list of surface horizontal offsets based on the surface horizontal offsets of each incident angle.

4. The method for determining the offset of the local target according to claim 3, characterized in that The obtaining the surface horizontal offset of the incident angle based on the average velocities and depth - direction grid spacing corresponding to the number of depth layers from the surface to the local target includes: Obtaining the sine value of the incident angle; Calculating the surface horizontal offset of the first layer in the number of depth layers in the direction from the surface to the local target based on the sine value of the incident angle and the depth - direction grid spacing; For other layers after the first layer, obtaining the current - layer coefficient based on the average velocity corresponding to the current layer, the average velocity corresponding to the previous layer of the current layer and the sine value of the incident angle; Calculating the surface horizontal offset of the current layer based on the surface horizontal offset of the previous layer, the depth - direction grid spacing and the current - layer coefficient until the current layer is the layer corresponding to the number of depth layers, to obtain the surface horizontal offset of the incident angle.

5. The method for determining the offset of the local target according to any one of claims 1 to 4, characterized in that, The obtaining the reflection position of the surface position point with respect to the local target includes: Connect the surface position point to the underground local target, and connect the surface position point to the surface emergence position and extend it; In the triangle formed by the surface position point, the underground local target, and the surface emergence position, obtain the first angle degree between the first straight line formed by the surface position point and the underground local target and the second straight line formed by the underground local target and the surface emergence position; On the extension line, determine the reflection position, and the second angle degree between the third straight line formed by the reflection position and the underground local target and the fourth straight line formed by the surface emergence position and the underground local target is equal to the first angle degree.

6. The method for determining the offset of the shot-receiver for a local target according to claim 5, wherein Constructing the aperture angle equal value map and the shot-receiver offset equal value map of the local target based on the surface position point, the reflection position, and the surface emergence position includes: Based on the second angle degree, obtain the local target aperture angle of the local target based on the surface position point; Obtain the distance from the surface position point to the reflection position to get the shot-receiver offset; Based on the local target aperture angles of each surface position point, construct the aperture angle equal value map; Based on the shot-receiver offsets of each surface position point, construct the shot-receiver offset equal value map.

7. The method for determining the offset of the local target according to any one of claims 1 to 4, characterized in that The method further includes: Divide the surface horizontal offset list according to a preset offset interval to obtain multiple divided values, and use the interpolation method to obtain the incident angle corresponding to each divided value to update the surface horizontal offset list.

8. An offset determination device for a local target, characterized in that, The shot-receiver offset determination device for the local target includes: A model parameter acquisition module for acquiring the grid parameters of a pre-constructed three-dimensional velocity grid model and the velocities corresponding to the grid parameters; A target parameter acquisition module for acquiring the position parameters of a preset local target in the three-dimensional velocity grid model; A depth velocity acquisition module for calculating the average velocity in each depth direction in the three-dimensional velocity grid model based on the acquired grid parameters and the velocities corresponding to the grid parameters; An offset calculation module for obtaining a list of surface horizontal offsets corresponding to the incident angles of the local target based on the average velocities in each depth direction; An emergence position determination module for obtaining the offset corresponding to the position parameter from the list of surface horizontal offsets according to the position parameters of the local target, and obtaining the surface emergence position of the local target based on the position parameters and the offset of the local target; A shot-receiver offset determination module for determining a surface position point based on the grid parameters, obtaining the reflection position of the surface position point with respect to the local target, constructing the aperture angle equal value map and the shot-receiver offset equal value map of the local target based on the surface position point, the reflection position, and the surface emergence position, and determining the shot-receiver offset of the local target based on the aperture angle equal value map and the shot-receiver offset equal value map.

9. A storage medium, characterized in that, A program or instruction is stored on a storage medium, and when the program or instruction is run by a processor, the steps of the method for determining the shot-receiver offset of the local target according to any one of claims 1 to 7 are implemented.

10. An electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the program, the steps of the method for determining the shot-receiver offset of the local target according to any one of claims 1 to 7 are implemented.

Images