Local target's offset determination method and device, storage medium and electronic equipment

Abstract

The present application relates to a kind of local target's shot distance determination method, device, storage medium and electronic equipment, comprising: obtaining the grid parameter and speed of three-dimensional velocity grid model;The position parameter of local target is obtained;Based on the grid parameter and speed obtained, the average speed in each depth direction is calculated;Based on the average speed in each depth direction, the ground level offset list of the incidence angle of local target is obtained;From ground level offset list, the offset distance corresponding to position parameter is obtained, based on position parameter and offset distance, the ground exit position of local target is obtained;Based on grid parameter, the ground position point is determined, the reflection position of local target is obtained for this ground position point, based on ground position point, reflection position and ground exit position, the opening angle contour map and shot distance contour map of local target are constructed, based on opening angle contour map and shot distance contour map, the shot distance of local target is determined.Can improve the determination precision of shot distance.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geophysical exploration, and in particular to a method and device for determining shot-receiver distance of a local target, a storage medium and an electronic device. BACKGROUND

[0002] Seismic data acquisition is the first process in oil and gas seismic exploration engineering, and complete seismic data acquisition function is realized by using a seismic detector and a seismic exploration instrument. Determining appropriate shot-receiver distance is an extremely important process in seismic data acquisition. With the development of seismic exploration technology, the requirements for seismic exploration are also increasing. In recent years, pre-stack depth migration imaging technology has been widely used in the field of seismic imaging. Although seismic migration imaging calculation has developed rapidly, the accurate delineation of complex structures, especially salt dome boundaries, is still a difficult problem. Because the shape of the salt dome structure is no longer a simple fault fracture, therefore, how to determine the shot-receiver distance to improve the accuracy of the determined shot-receiver distance is particularly important for realizing full target imaging. SUMMARY

[0003] Therefore, the present application provides a method and device for determining shot-receiver distance of a local target, a storage medium and an electronic device.

[0004] Specifically, the present application is realized by the following technical solutions:

[0005] According to a first aspect of the present application, a method for determining shot-receiver distance of a local target is provided, which comprises:

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

[0007] obtaining position parameters of a pre-set local target in the three-dimensional velocity grid model;

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

[0009] based on the average velocities in each depth direction, obtaining a list of surface horizontal migration distances corresponding to an incident angle of the local target;

[0010] obtaining a migration distance corresponding to the position parameters from the list of surface horizontal migration distances according to the position parameters of the local target, and obtaining a surface exit position of the local target based on the position parameters of the local target and the migration distance;

[0011] Based on the grid parameters, a surface location point is determined, and the reflection position of the surface location point relative to the local target is obtained. Based on the surface location point, the reflection position, and the surface emission position, an angular contour map and a shot-receiver distance contour map of the local target are constructed. Based on the angular contour map and the shot-receiver distance contour map, the shot-receiver distance of the local target is determined.

[0012] The method for determining the shot-receiver distance of a local target in this technical solution involves acquiring grid parameters, the velocities corresponding to the grid parameters, and the position parameters of the local target. Based on the acquired grid parameters and the velocities corresponding to the grid parameters, the average velocity in each depth direction is calculated. Based on the average velocity in each depth direction, a list of horizontal offsets corresponding to each incident angle of the local target is obtained. According to the position parameters of the local target, the offset distance is obtained from the list of horizontal offsets, thereby obtaining the surface launch position of the local target based on the position parameters and offset distance of the local target. Then, based on the grid parameters, the surface location point is determined, and the reflection position of the surface location point relative to the local target is obtained. Based on the surface location point, the reflection position, and the surface launch position, an angular contour map and a shot-receiver distance contour map of the local target are constructed. Based on the angular contour map and the shot-receiver distance contour map, the shot-receiver distance of the local target is determined. Thus, a method for determining the gun-receiver distance based on the specific local target is proposed. By using the average velocity in each depth direction, the horizontal offset of the ground surface at each incident angle of the local target can be obtained. This allows for the accurate acquisition of the incident angle corresponding to different horizontal offsets of the ground surface. Therefore, the calculation of the gun-receiver distance based on the accurately determined incident angle can improve the accuracy of the determined gun-receiver distance, thereby enabling full-target imaging.

[0013] According to a second aspect of the present invention, a gun-receiver distance determination device for a local target is provided, the device comprising:

[0014] The model parameter acquisition module is used to acquire the mesh parameters of the pre-built 3D velocity mesh model and the velocities corresponding to the mesh parameters;

[0015] The target parameter acquisition module is used to acquire the position parameters of a pre-set local target in the three-dimensional velocity grid model.

[0016] The depth velocity acquisition module is used to calculate the average velocity in each depth direction in the three-dimensional velocity mesh model based on the acquired mesh parameters and the velocities corresponding to the mesh parameters.

[0017] The offset calculation module is used to obtain a list of horizontal offsets of the ground surface corresponding to the incident angle of the local target based on the average velocity in each depth direction.

[0018] The launch position determination module is used to obtain the offset distance corresponding to the position parameter from the list of horizontal offsets on the ground surface based on the position parameter of the local target, and to obtain the ground launch position of the local target based on the position parameter and offset distance of the local target.

[0019] The shot-receiver distance determination module is used to determine the surface location point based on the grid parameters, obtain the reflection position of the surface location point relative to the local target, construct the angular contour map and shot-receiver distance contour map of the local target based on the surface location point, the reflection position and the surface launch position, and determine the shot-receiver distance of the local target based on the angular contour map and the shot-receiver distance contour map.

[0020] According to a third aspect of the present invention, a storage medium is provided having a computer program stored thereon, wherein the program, when executed by a processor, implements the steps of the method for determining the gun-receiver distance of a local target in any possible implementation of the first aspect.

[0021] According to a fourth aspect of the present invention, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps of the method for determining the gun-receiver distance of a local target in any possible implementation of the first aspect. Attached Figure Description

[0022] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with the invention and, together with the description, serve to explain the principles of the invention.

[0023] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, those skilled in the art can obtain other drawings based on these drawings without creative effort.

[0024] Figure 1 This is a flowchart illustrating a method for determining the shot-receiver distance of a local target, provided in an embodiment of the present invention.

[0025] Figure 2 A schematic diagram of the horizontal offset of the ground surface at different incident angles in a method for determining the shot-receiver distance of a local target provided in an embodiment of the present invention;

[0026] Figure 3 This is a schematic diagram of the incident angles corresponding to different offset intervals in a method for determining the shot-receiver distance of a local target provided in an embodiment of the present invention.

[0027] Figure 4 This is a schematic diagram illustrating the determination of the reflection position in a method for determining the shot-receiver distance of a local target, provided in an embodiment of the present invention.

[0028] Figure 5 This is a schematic diagram of a local target gun-receiver distance determination device provided in an embodiment of the present invention;

[0029] Figure 6 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0030] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0031] For seismic data acquisition, determining the shot-receiver distance (SRT) to improve its accuracy and achieve full-target imaging is crucial for complex geological structures. This embodiment utilizes angular domain imaging, decomposing the complex geological area into different local targets. Based on the requirements of these local targets on the seismic acquisition system parameters, the optimal SRT parameters and the effective surface receiving range are derived. This allows for precise imaging of different sides, top surfaces, and other local targets within the complex geological area, achieving full-target imaging and providing methodological support for full-target seismic exploration. This approach can be used for the acquisition design of seismic acquisition systems in complex geological structures, optimizing the system and making it suitable for three-dimensional complex geological full-target seismic acquisition systems. It provides a method for optimizing the shot-receiver point locations in seismic acquisition systems and offers reasonable parameters for system design.

[0032] See Figure 1 This invention provides a method for determining the shot-receiver distance of a local target. This embodiment is applicable to determining the shot-receiver distance parameters of an existing three-dimensional mesh geological model (three-dimensional velocity mesh model) with velocity as its attribute. The shot-receiver distance is determined by the depth velocity of the three-dimensional velocity mesh model and the spatial position of the local target. The method may include the following steps:

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

[0034] In this embodiment, the three-dimensional velocity grid model is constructed based on the seismic data acquisition area to be acquired. As an optional embodiment, the grid parameters include, but are not limited to, the number of sampling points and the grid spacing between 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, the number of sampling points in the x-direction of the three-dimensional velocity mesh model is Nx, and the mesh spacing between the sampling points in the x-direction is dx; the number of sampling points in the y-direction is Ny, and the mesh spacing between the sampling points in the y-direction is dy; the number of sampling points in the z-direction is Nz, and the mesh spacing between the sampling points in the z-direction is dz, where the z-direction is the depth direction.

[0036] In this embodiment, taking a 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, and the grid spacing is dx = 30m; the number of sampling points in the y-direction is Ny = 571, and 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 targets in the three-dimensional velocity grid model;

[0038] In this embodiment, the local target is one of the sampling points. As an optional embodiment, it is a sampling point located on the side or top surface of a pre-laid three-dimensional velocity grid model below the ground. For the specific process of determining the local target, please refer to relevant technical documents, which will not be detailed here.

[0039] In this embodiment, as an optional implementation, the position parameters include, but are not limited to: position coordinates, azimuth along the tilt direction, and tilt angle along the tilt direction. The position coordinates of the local target are denoted as o(x,y,z) in meters, the azimuth along the tilt direction is denoted as Azimuth, and the tilt angle is denoted as Dip, both in degrees.

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

[0041] S103. Based on the acquired mesh parameters and the velocities corresponding to the mesh parameters, calculate the average velocity in each depth direction of the three-dimensional velocity mesh model;

[0042] In this embodiment, the average velocity at different depths of the three-dimensional velocity mesh model is obtained based on the velocity of the mesh parameters at each depth. The depth can be determined based on the number of faults contained in the three-dimensional velocity mesh model. As an optional embodiment, the average velocity in each depth direction of the three-dimensional velocity mesh model is calculated based on the obtained mesh parameters and the velocities corresponding to those parameters, including:

[0043] A11, for each target depth, obtain the velocity of each sampling point located at that target depth;

[0044] In this embodiment, multiple sampling points are set up at each target depth, and the number of target depths is equal to the total depth divided by the grid spacing in the z direction.

[0045] A12, calculate the velocity and value of each sampling point;

[0046] In this embodiment, the velocities of each sampling point at the target depth are summed and accumulated.

[0047] A13, obtain the number of sampling points at each sampling point located at the target depth;

[0048] A14, calculate the quotient of the velocity sum and the number of sampling points to obtain the average velocity at the target depth.

[0049] In this embodiment, as an optional implementation, the average velocity at the target depth is calculated 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) in meters per second, Sum represents summation, ix is ​​the index of the target depth in the x-direction with a starting index of 1 and an ending index of Nx, iy is the index of the target depth in the y-direction with a starting index of 1 and an ending index of Ny, iz is the target depth in the z-direction with a starting index of 1 and an ending index of Nz, / represents division, and * represents multiplication.

[0052] In this embodiment, the sampling points are arranged in a matrix, forming an Nx*Ny matrix.

[0053] In this embodiment, as an optional implementation, 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, 1057*571 sampling points are arranged, and 480 target depths are arranged in the depth direction.

[0054] By calculating the velocity at each sampling point at each target depth, the average velocity at different depths is 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. Based on the average velocity in each depth direction, obtain a list of horizontal offsets of the ground surface corresponding to the incident angle of the local target.

[0057] In this embodiment, for a local target, based on the average velocity in each depth direction, the horizontal offset of the ground surface corresponding to any incident angle (iangle) of the local target is calculated, and a list of horizontal offsets of the ground surface is generated based on the horizontal offsets of the ground surface corresponding to each incident angle.

[0058] In this embodiment, as an optional implementation, a list of horizontal surface offsets corresponding to the incident angle of the local target is obtained based on the average velocity in each depth direction, including:

[0059] B11, determine the starting index and maximum index of the incident angle, and, 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, determine the depth layer number of the local target;

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

[0061] In this embodiment, as an optional embodiment, the starting index of the incident angle iangle is 1, the maximum index is 89, the unit is degrees, 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. The grid spacing dz of the sampling points in the z direction is 30m. Then the depth layer of the local target is 4800 / 30=160.

[0063] B12, traverse in order of incident angle from the starting index to the maximum index, and for each incident angle, obtain the horizontal offset of the ground surface for that incident angle based on the average velocity and depth direction grid spacing corresponding to the depth layer from the ground surface to the local target.

[0064] In this embodiment, as an optional embodiment, the horizontal offset of the incident angle is obtained based on the average velocity and depth direction grid spacing corresponding to the depth layers from the ground surface to the local target, including:

[0065] C11, obtain the sine value of the 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 angle of incidence in degrees.

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

[0070] In this embodiment, the depth layer number is 160, indicating that there are 160 layers. From the ground surface to the local target, they are respectively denoted as the first layer, the second layer, ..., the 160th layer.

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

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

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

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

[0075] In this embodiment, as an optional implementation, the coefficients of the current layer are 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 coefficient of the current layer, V2 is the average velocity of the current layer, and v1 is the average velocity of the previous layer. For example, V_average[1] = 2098.8 and V_average[2] = 2098.8.

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

[0080] In this embodiment, when iz is not equal to 1, sintheata2 = V2 * sintheata1 / v1 is calculated, such that sintheata1 = sintheata2. The horizontal offset of the current layer is calculated using the following formula:

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

[0082] In this embodiment, for each incident angle, for example, iangle = 1, 2, ..., 89, steps C11 to C14 are executed respectively to obtain the horizontal offset of the ground surface for each incident angle. The horizontal offset of the ground surface of each layer is accumulated sequentially to obtain the horizontal offset of the ground surface of the last layer. The horizontal offset of the last layer is the horizontal offset of the ground surface of the incident angle.

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

[0084] In this embodiment, the incident angle index i represents the incident angle number i.

[0085] B13, Based on the horizontal offset of the ground surface for each incident angle, construct the list of horizontal offsets of the ground surface.

[0086] Figure 2 This diagram illustrates the horizontal surface offset at different incident angles in a method for determining the shot-receiver distance of a local target according to an embodiment of the present invention. In this embodiment, the horizontal direction is used as the horizontal surface offset, the vertical direction is used as the depth parameter value, and the area directly below where the horizontal surface offset is 0 is considered a local target, thus constructing a list of horizontal surface offsets.

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

[0088] The list of horizontal offsets of the ground surface is divided according to a preset offset interval to obtain multiple segment values. The incident angle corresponding to each segment value is obtained by interpolation to update the list of horizontal offsets of the ground surface.

[0089] In this embodiment, the horizontal offset of the ground surface corresponding to different incident angles is divided according to a preset offset interval, such as 1m, to obtain multiple segment values. The incident angle corresponding to each segment value is obtained by interpolation.

[0090] In this embodiment, as an optional embodiment, the preset interval is 1 meter. The horizontal offset of the ground surface with different incident angles is divided into 1-meter intervals. For example, the horizontal offset of the ground surface with the largest incident angle is divided into 1-meter intervals, and the incident angle corresponding to each 1-meter interval is calculated by interpolation. The list of horizontal offsets of the ground surface is updated to construct the list of incident angles: Angle(ioffset), in degrees, where ioffset is the offset interval.

[0091] In this embodiment, ioffset is represented by an accumulated index, where the starting index is 1 and the maximum index is Noffset. This effectively expands the number of incident angles and improves the accuracy of the incident angles.

[0092] In this embodiment, as an optional implementation, the maximum index of the offset interval is obtained using the following formula:

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

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

[0095] Figure 3 This diagram illustrates the incident angles corresponding to different offset intervals in a method for determining the gun-receiver distance of a local target, provided in an embodiment of the present invention. In this embodiment, the list of incident angles corresponding to different offset intervals, calculated using an 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. Here, Angle(j) represents the incident angle corresponding to the j-th offset interval. For example, with an offset interval of 1m, Angle(8650) represents an incident angle of 87.9957 degrees corresponding to a horizontal offset of 8650m from the ground surface.

[0097] S105. Based on the position parameters of the local target, obtain the offset distance corresponding to the position parameters from the list of horizontal offsets on the ground surface, and obtain the ground surface emission position of the local target based on the position parameters and offset distance of the local target.

[0098] In this embodiment, for the updated list of horizontal offsets on the ground surface, the offset distance corresponding to the tilt angle is obtained from the list of incident angles. Based on the tilt angle of the local target along the tilt direction, the horizontal offset (offset distance) L1 corresponding to the tilt angle is obtained from the list of incident angles. The ground surface position (x, y) of the local target is translated by a distance L1 along the azimuth direction to obtain the ground surface exit position T of the local target, that is, the position point where the local target is vertically exited on the ground. Here, the tilt angle of the local target along the tilt direction corresponds to the incident angle, and the magnitude of the tilt angle is the magnitude of the incident angle.

[0099] In this embodiment, as an optional embodiment, obtaining the surface emission position of the local target based on the position parameters and offset distance of the local target includes:

[0100] The local target's surface position is shifted by the offset distance along the azimuth direction to obtain the local target's surface launch position.

[0101] In this embodiment, as described above, based on the local target's tilt angle of 15 degrees, a search is performed from the incident angle list, and OFFset(15) = 971.28m is obtained. The offset distance L1 = 971.28 corresponding to this tilt angle is obtained. The local target's surface position (x, y) is translated by a distance of 971.28 along the azimuth direction, thus obtaining the local target's surface emission position T. The coordinates of point T are (6040.91m, 5685.5m, 0m).

[0102] S106. Based on the grid parameters, determine the surface location point, obtain the reflection position of the surface location point relative to the local target, construct the angular contour map and the shot-receiver distance contour map of the local target based on the surface location point, the reflection position and the surface emission position, and determine the shot-receiver distance of the local target based on the angular contour map and the shot-receiver distance contour map.

[0103] In this embodiment, based on the grid spacing of the sampling points in the x-direction and the grid spacing of the sampling points in the y-direction of the grid parameters, the surface location point S(ix,iy) of the three-dimensional velocity grid model is selected. The surface location point S(ix,iy) is the number of sampling points on the surface, i.e., Nx*Ny. For each surface location point S(ix,iy), the reflection position R(ixr,iyr) of the surface location point S(ix,iy) relative to the underground local target O is obtained, and then the local target opening angle of the local target is calculated.

[0104] In this embodiment, as an optional embodiment, obtaining the reflected position of the ground location point relative to the local target includes:

[0105] D11 connects the surface location point with the underground local target, and connects the surface location point with the surface launch point and extends it;

[0106] In this embodiment, the surface location point S is connected to the underground local target o, and the surface location point S is connected to the surface launch position T, and extended in the direction from the surface location point S to the surface launch position T.

[0107] In this embodiment, as an optional embodiment, ix = 101, iy = 101, and the corresponding location coordinates of the ground location point S(101,101) are (3000m, 3000m, 0m).

[0108] D12, in the triangle formed by the surface location point, the underground local target, and the surface launch point, obtain the first angle between the first straight line formed by the surface location point and the underground local target, and the second straight line formed by the underground local target and the surface launch point;

[0109] In this embodiment, the degree measure of angle SOT is calculated within the triangle SOT formed by the surface location point S, the underground local target o, and the surface emission point T. Given S(3000m, 3000m, 0m), o(5200m, 5200m, 4800m), and T(6040.91m, 5685.5m, 0m), the line segment lengths are: ST = 4056.64m, SO = 5720.14m, and TO = 4897.18m.

[0110] D13, on the extended line, determine the reflection position. The second angle between the third straight line formed by the reflection position and the underground local target, and the fourth straight line formed by the surface emission position and the underground local target, is equal to the first angle.

[0111] Figure 4 This is a schematic diagram illustrating the determination of the reflection position in a method for determining the shot-receiver distance of a local target according to an embodiment of the present invention. In this embodiment, a point R is taken on the extension line ST such that angle TOR equals angle SOT; this point R is the desired reflection position R. As described above, according to the law of cosines, the degree measure of angle SOT can be calculated as 44.08 degrees. Taking a point R on the extension line ST such that angle TOR equals angle SOT, point R is the desired reflection position.

[0112] In this embodiment, as an optional implementation, based on the surface location point, reflection location, and surface emission location, an angular contour map and a shot-receiver distance contour map of the local target are constructed, including:

[0113] E11, based on the second angle measure, obtain the local target subtended angle of the local target based on the surface location point;

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

[0115] E12, obtain the distance from the ground location point to the reflection location, and get the shot-receiver distance;

[0116] In this embodiment, the length of SR is calculated and recorded at the reflection position R(ixr, iyr) to form a shot-receiver distance file. In this embodiment, the length SR = 9235.83m is obtained, resulting in ixr = 330, iyr = 303, and the local target opening angle openAngle = 44.08 degrees. This is recorded at position R(330, 303) to form an opening angle file.

[0117] E13, construct the angle contour map based on the local target angles at various surface locations;

[0118] In this embodiment, the local target angle corresponding to the ground location point in all three-dimensional velocity grid models is calculated point by point to form an angle contour map.

[0119] E14, construct the shot-receiver distance contour map based on the shot-receiver distance of each surface location point.

[0120] In this embodiment, the SR is recorded at position R(330, 303) to form a shot-receiver distance file. The SR length corresponding to the ground position in all three-dimensional velocity grid models is calculated point by point to form a shot-receiver distance contour map.

[0121] In this embodiment, as an optional implementation, determining the shot-receiver distance of the local target based on the angular contour map and the shot-receiver distance contour map includes:

[0122] F11, based on the angle contour map, obtain the region where the angle of the local target is within the preset angle threshold, and obtain the optimal excitation and receiving region of the local target;

[0123] F12, in the contour map of the shot-receiver distance corresponding to the optimal excitation receiving area, obtain the maximum shot-receiver distance, and determine that the shot-receiver distance of the local target is greater than or equal to the maximum shot-receiver distance.

[0124] In this embodiment, based on the angular contour map and the shot-receiver distance contour map, the region where the angularity of the local target is within the preset angularity threshold is obtained, and the excitation and receiving region is obtained. The excitation and receiving region is mapped onto the shot-receiver distance contour map, and the maximum length between each surface position point and the corresponding reflection position in the mapped shot-receiver distance contour map is obtained, and the maximum shot-receiver distance is obtained. The shot-receiver distance of the local target is determined to be greater than or equal to the maximum shot-receiver distance.

[0125] In this embodiment, as an optional embodiment, the angular threshold is 60 degrees. In the angular contour map, the contour area with an angular threshold of 60 degrees is set, and the internal area within 60 degrees is the optimal excitation and reception area of ​​the local target.

[0126] In this embodiment, to ensure complete seismic data imaging of the local target, the determined shot-receiver distance must be greater than or equal to MaxOffset, where MaxOffset is the maximum length of each seismic SR within the 60-degree region. As an optional embodiment, the shot-receiver distance of 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.4m.

[0129] In this embodiment, based on seismic reflection imaging technology, the local geological target (local target) is analyzed. By using the principle of incident reflection, the reflection angle (local target angle) corresponding to different shot-receiver points (surface location points) is analyzed. The maximum value of the distance between the corresponding shot-receiver points on the surface is scanned, thereby obtaining the optimal shot-receiver distance.

[0130] Based on the same inventive concept, such as Figure 5 As shown, this embodiment of the invention also provides a device for determining the gun-receiver distance of a local target, the device comprising:

[0131] The model parameter acquisition module 501 is used to acquire the mesh parameters of the pre-built three-dimensional velocity mesh model and the velocity corresponding to the mesh parameters;

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

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

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

[0135] The depth velocity acquisition module 503 is used to calculate the average velocity in each depth direction in the three-dimensional velocity mesh model based on the acquired mesh parameters and the velocities corresponding to the mesh parameters.

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

[0137] A velocity acquisition unit (not shown in the figure) is used to acquire the velocity of each sampling point located at each target depth for each target depth.

[0138] The calculation unit is used to calculate the velocity and value of each sampling point.

[0139] The statistical unit is used to obtain the number of sampling points at each sampling point located at the target depth;

[0140] The mean acquisition unit is used to calculate the quotient of the velocity sum and the number of sampling points to obtain the average velocity at the target depth.

[0141] The offset calculation module 504 is used to obtain a list of horizontal offsets of the ground surface corresponding to the incident angle of the local target based on the average velocity in each depth direction.

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

[0143] The layer number determination unit is used to determine the starting index and the maximum index of the incident angle, and to determine 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] The traversal unit is used to traverse in order from the starting index to the maximum index according to the incident angle. For each incident angle, the horizontal offset of the ground surface is obtained based on the average velocity and depth direction grid spacing corresponding to the depth layer from the ground surface to the local target.

[0145] The offset list construction unit is used to construct the surface horizontal offset list based on the surface horizontal offset for each incident angle.

[0146] In this embodiment, as an optional implementation, the traversal unit is specifically used for:

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

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

[0149] For the layers after the first layer, the coefficient of the current layer is obtained based on the average velocity corresponding to the current layer, the average velocity corresponding to the previous layer, and the sine value of the incident angle.

[0150] Based on the horizontal offset of the previous layer, the grid spacing in the depth direction, and the coefficient of the current layer, the horizontal offset of the current layer is calculated until the current layer is the layer corresponding to the depth layer number, and the horizontal offset of the current layer is obtained for the incident angle.

[0151] The launch position determination module 505 is used to obtain the offset distance corresponding to the position parameter from the list of horizontal offsets on the ground surface according to the position parameter of the local target, and obtain the ground launch position of the local target based on the position parameter and offset distance of the local target.

[0152] In this embodiment, as an optional embodiment, obtaining the surface emission position of the local target based on the position parameters and offset distance of the local target includes:

[0153] The local target's surface position is shifted by the offset distance along the azimuth direction to obtain the local target's surface launch position.

[0154] The gun-receiver distance determination module 506 is used to determine the surface location point based on the grid parameters, obtain the reflection position of the surface location point relative to the local target, construct the angular contour map and the gun-receiver distance contour map of the local target based on the surface location point, the reflection position and the surface launch position, and determine the gun-receiver distance of the local target based on the angular contour map and the gun-receiver distance contour map.

[0155] In this embodiment, as an optional embodiment, the shot-receiver distance determination module 506 includes:

[0156] The excitation and receiving area determination unit is used to obtain the area where the local target's angle is within a preset angle threshold based on the angle contour map, and to obtain the optimal excitation and receiving area of ​​the local target.

[0157] The gun-receiver distance determination unit is used to obtain the maximum gun-receiver distance from the gun-receiver distance contour map corresponding to the optimal excitation and receiving area, and to determine that the gun-receiver distance of the local target is greater than or equal to the maximum gun-receiver distance.

[0158] In this embodiment, as an optional embodiment, the shot-receiver distance determination module 506 includes:

[0159] A connecting unit is used to connect a surface location point with an underground local target, and to connect a surface location point with a surface launch point and extend it.

[0160] The first angle calculation unit is used to obtain the first angle between the first straight line formed by the surface location point and the underground local target, and the second straight line formed by the underground local target and the surface emission point, in the triangle formed by the surface location point, the underground local target, and the surface emission point.

[0161] The second included angle calculation unit is used to determine the reflection position on the extended line. The second included angle between the third straight line formed by the reflection position and the underground local target and the fourth straight line formed by the surface emission position and the underground local target is equal to the first included angle.

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

[0163] Angle acquisition unit is used to acquire the local target angle based on the second angle value and the local target position point on the ground.

[0164] The shot-receiver distance acquisition unit is used to obtain the distance from the ground location point to the reflection location, thus obtaining the shot-receiver distance;

[0165] Angle contour map construction unit is used to construct the angle contour map based on the local target angles at various surface location points;

[0166] The shot-receiver distance contour map construction unit is used to construct the shot-receiver distance contour map based on the shot-receiver distances of various surface location points.

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

[0168] An incident angle segmentation module (not shown in the figure) is used to segment the list of horizontal offsets of the ground surface according to a preset offset interval to obtain multiple segmentation values. The incident angle corresponding to each segmentation value is obtained by interpolation to update the list of horizontal offsets of the ground surface.

[0169] In this embodiment, as an optional embodiment, the preset interval is 1 meter. The horizontal offset of the ground surface with the maximum incident angle is divided into 1-meter intervals, and the incident angle corresponding to each 1-meter interval is calculated by interpolation to construct an incident angle list.

[0170] Based on the same inventive concept, embodiments of the present invention also provide a storage medium storing a computer program thereon, which, when executed by a processor, implements the steps of the method for determining the gun-receiver distance of a local target in any of the above possible implementations.

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

[0172] Based on the same inventive concept, see [link to inventive concept] Figure 6This invention also provides an electronic device, including a memory 101 (e.g., non-volatile memory), a processor 102, and a computer program stored in the memory 101 and executable on the processor 102. When the processor 102 executes the program, it implements the steps of the local target gun-receiver distance determination method in any of the above possible implementations, which can be equivalent to the aforementioned local target gun-receiver distance determination device. Of course, the processor can also be used to process other data or perform calculations. This electronic device can be a PC, server, terminal, or other similar device.

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

[0174] It should be noted that the aforementioned local target gun-receiver distance determination device can be implemented by software. As a logical device, it is formed by the processor 102 of the electronic device in which it is located reading the computer program instructions stored in the non-volatile memory into the memory 103 for execution.

[0175] The embodiments of the subject matter and functional operation described in this specification can be implemented in the following ways: digital electronic circuits, tangibly embodied computer software or firmware, computer hardware including the structures disclosed in this specification and their structural equivalents, or combinations thereof. Embodiments of the subject matter described in this specification can be implemented as one or more computer programs, i.e., one or more modules of computer program instructions encoded on a tangible, non-transitory program carrier for execution by a data processing apparatus or for controlling the operation of a data processing apparatus. Alternatively or additionally, the program instructions may be encoded on artificially generated propagation signals, such as machine-generated electrical, optical, or electromagnetic signals, which are generated to encode information and transmit it to a suitable receiving device for execution by the data processing apparatus. The computer storage medium may be a machine-readable storage device, a machine-readable storage substrate, a random or serial access memory device, or combinations thereof.

[0176] The processing and logic flow described in this specification can be executed by one or more programmable computers that execute one or more computer programs to perform corresponding functions by operating on input data and generating output. The processing and logic flow can also be executed by special-purpose logic circuitry—such as FPGAs (Field-Programmable Gate Arrays) or ASICs (Application-Specific Integrated Circuits), and the device can also be implemented as special-purpose logic circuitry.

[0177] Suitable computers for executing computer programs include, for example, general-purpose and / or special-purpose microprocessors, or any other type of central processing unit. Typically, the central processing unit receives 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. Typically, a computer will also include one or more mass storage devices for storing data, such as disks, magneto-optical disks, or optical disks, or the computer will be operatively coupled to such mass storage devices to receive data from or transfer data to them, or both. However, a computer is not required to have such devices. Furthermore, a computer can 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 a few.

[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 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. Processors and memory may be supplemented by or incorporated into dedicated logic circuitry.

[0179] While this specification contains numerous specific implementation details, these should not be construed as limiting the scope of any invention or the scope of the claims, but rather are primarily used to describe features of specific embodiments of a particular invention. Certain features described in the various embodiments herein may also be implemented in combination in a single embodiment. Conversely, various features described in a single embodiment may also be implemented separately in various embodiments or in any suitable sub-combination. Furthermore, while features may function in certain combinations as described above and even initially claimed in this way, one or more features from a claimed combination may be removed from that combination in some cases, and a claimed combination may refer to a sub-combination or a variation thereof.

[0180] Similarly, although the operations are depicted in a specific order in the accompanying drawings, this should not be construed as requiring these operations to be performed in the specific order shown or sequentially, or requiring all illustrated operations to be performed to achieve the desired result. In some cases, multitasking and parallel processing may be advantageous. Furthermore, the separation of various system modules and components in the above embodiments should not be construed as requiring such separation 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 may be performed in a different order and still achieve the desired result. Furthermore, the processes depicted in the drawings are not necessarily shown in a specific order or sequence to achieve the desired result. In some implementations, multitasking and parallel processing may be advantageous.

[0182] It should be noted that, in this document, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes the element.

[0183] The above are merely specific embodiments of the present invention, enabling those skilled in the art to understand or implement the invention. Various modifications to these embodiments will be readily 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 invention. Therefore, the present invention is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features claimed herein.

Claims

1. A method for determining shot-to-target distances for a local target, characterized by, The method comprises the following steps: obtaining grid parameters of a pre-constructed three-dimensional velocity grid model and velocities corresponding to the grid parameters; obtaining position parameters of a local target pre-set in the three-dimensional velocity grid model; calculating average velocities in each depth direction of 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 offset values corresponding to an incident angle of the local target based on the average velocities in each depth direction; obtaining an offset distance corresponding to the position parameters of the local target from the list of surface horizontal offset values, and obtaining a surface exit position of the local target based on the position parameters of the local target and the offset distance; determining a surface position point based on the grid parameters, obtaining a reflection position of the local target relative to the surface position point, constructing an opening angle contour map and a shot-receiver distance contour map of the local target based on the surface position point, the reflection position and the surface exit position, and determining a shot-receiver distance of the local target based on the opening angle contour map and the shot-receiver distance contour map.

2. The local target's shot-to-hit determination method of claim 1, wherein, The method comprises the following steps: for each target depth, obtaining velocities of each sampling point located at the target depth; calculating a sum of velocities of the velocities of each sampling point; obtaining a number of sampling points of each sampling point located at the target depth; calculating a quotient value of the sum of velocities and the number of sampling points to obtain an average velocity of the target depth.

3. The local target shot gather determination method of claim 1, wherein, The method comprises the following steps: determining a starting index and a maximum index of the incident angle, and determining a depth layer number of the local target based on a depth direction parameter value in the position parameters of the local target and a depth direction grid interval in the grid parameters; traversing the incident angle from the starting index to the maximum index in sequence, and for each incident angle, obtaining a surface horizontal offset value of the incident angle based on average velocities and depth direction grid intervals corresponding to a depth layer number from the surface to the local target; constructing the list of surface horizontal offset values based on the surface horizontal offset value of each incident angle.

4. The local target's shot-to-hit determination method of claim 3, wherein, The method comprises the following steps: obtaining a sine value of the incident angle; calculating a surface horizontal offset value of a first layer in the depth layer number from the surface to the local target direction based on the sine value of the incident angle and the depth direction grid interval; for other layers after the first layer, obtaining a current layer coefficient according to an average velocity corresponding to the current layer, an average velocity corresponding to a previous layer of the current layer and the sine value of the incident angle; calculating a surface horizontal offset value of the current layer based on the surface horizontal offset value of the previous layer, the depth direction grid interval and the current layer coefficient until the current layer is the layer corresponding to the depth layer number, to obtain the surface horizontal offset value of the incident angle.

5. The local target's shot-distance determination method according to any one of claims 1 to 4, characterized in that, The method comprises the following steps: Linking the ground surface position point with the underground local target, and linking the ground surface position point with the ground surface exit position and extending; In the triangle formed by the ground surface position point, the underground local target and the ground surface exit position, the first included angle between the first straight line formed by the ground surface position point and the underground local target and the second straight line formed by the underground local target and the ground surface exit position is obtained; On the extension line, the reflection position is determined, and the second included angle between the third straight line formed by the reflection position and the underground local target and the fourth straight line formed by the ground surface exit position and the underground local target is equal to the first included angle.

6. The local target's shot-to-hit determination method of claim 5, wherein, The steps of constructing the local target opening angle contour map and the offset distance contour map of the local target based on the ground surface position point, the reflection position and the ground surface exit position include: Based on the second included angle, the local target opening angle of the local target based on the ground surface position point is obtained; The distance from the ground surface position point to the reflection position is obtained to obtain the offset distance; Based on the local target opening angle of each ground surface position point, the opening angle contour map is constructed; Based on the offset distance of each ground surface position point, the offset distance contour map is constructed.

7. The local target's shot-distance determination method according to any one of claims 1 to 4, characterized in that, The method further includes: The ground surface horizontal offset distance list is segmented according to a pre-set offset interval to obtain a plurality of segmentation values, and the incident angle corresponding to each segmentation value is obtained by using an interpolation method to update the ground surface horizontal offset distance list.

8. A local target shotpoint determination device characterized by, The local target offset distance determination device includes: A model parameter acquisition module is configured to acquire grid parameters of a pre-constructed three-dimensional velocity grid model and velocities corresponding to the grid parameters; A target parameter acquisition module is configured to acquire position parameters of a pre-set local target in the three-dimensional velocity grid model; A depth velocity acquisition module is configured to calculate average velocities in each depth direction of the three-dimensional velocity grid model based on the acquired grid parameters and velocities corresponding to the grid parameters; An offset distance calculation module is configured to acquire a ground surface horizontal offset distance list corresponding to an incident angle of the local target based on the average velocities in each depth direction; An exit position determination module is configured to obtain an offset distance corresponding to the position parameters from the ground surface horizontal offset distance list according to the position parameters of the local target, and acquire a ground surface exit position of the local target based on the position parameters of the local target and the offset distance; An offset distance determination module is configured to acquire a ground surface position point based on the grid parameters, acquire a reflection position of the local target for the ground surface position point, construct an opening angle contour map and an offset distance contour map of the local target based on the ground surface position point, the reflection position and the ground surface exit position, and determine the offset distance of the local target based on the opening angle contour map and the offset distance contour map.

9. A storage medium, characterized by A program or instruction is stored on a storage medium, and the program or instruction is executed by a processor to implement the steps of the local target offset distance determination method according to any one of claims 1 to 7.

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, The processor executes the program to implement the steps of the local target offset distance determination method according to any one of claims 1 to 7.

Images