Methods, apparatus, devices, storage media, and products for determining formation dip

Abstract

The application discloses a method, device, equipment, storage medium and product for determining a stratum dip angle, and belongs to the technical field of geophysical seismic exploration. The method comprises the following steps: determining a horizontal well design track to be drilled; determining a seismic arbitrary line based on the horizontal well design track, and dividing the seismic arbitrary line into a plurality of straight line segments; determining intersection information of the plurality of straight line segments and a seismic grid, wherein the intersection information of each straight line segment and the seismic grid comprises intersection point coordinates and an intersection type; determining seismic trace data in a preset spatial window with the intersection point coordinates as the center based on the intersection type; determining seismic trace data of the intersection point coordinates based on the seismic trace data in the preset spatial window; drawing a seismic profile based on the seismic trace data of a plurality of intersection point coordinates included in the seismic arbitrary line; and determining a stratum dip angle based on the seismic profile. The application can improve the accuracy of the determined stratum dip angle.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geophysical seismic exploration, and in particular to a method and device for determining a stratum dip angle, equipment, a storage medium and a product. BACKGROUND

[0002] With the continuous progress of three-dimensional seismic exploration technology and the continuous deepening of the application of three-dimensional seismic data in the field of oilfield exploration and development, the importance of seismic data is increasing. In the implementation of horizontal well drilling in oilfields, predicting the stratum dip angle in front of the drill bit based on seismic data has become a means to guide drilling.

[0003] In the related art, a seismic arbitrary line is determined, seismic trace data observed along the seismic arbitrary line is plotted to obtain a seismic profile of the seismic arbitrary line, and a stratum dip angle is determined based on the dip angle of the seismic profile. However, when plotting the seismic profile of the seismic arbitrary line, a nearest borrowing drawing method is often used, which will project seismic trace data not on the seismic arbitrary line onto the seismic profile, resulting in a "sawtooth" phenomenon on the seismic profile, that is, the seismic profile is not accurate, and thus the stratum dip angle determined based on the dip angle of the seismic profile will not be accurate, that is, the accuracy of the stratum dip angle determined in the related art is poor. SUMMARY

[0004] The embodiments of the present application provide a method, device, equipment, storage medium and product for determining a stratum dip angle, which can improve the accuracy of the determined stratum dip angle. The technical solution is as follows:

[0005] In one aspect, a method for determining a stratum dip angle is provided, and the method comprises:

[0006] determining a horizontal well design trajectory to be drilled;

[0007] based on the horizontal well design trajectory, determining a seismic arbitrary line, and dividing the seismic arbitrary line into a plurality of straight line segments;

[0008] determining intersection information of the plurality of straight line segments and a seismic grid, and the intersection information of any straight line segment and the seismic grid includes intersection point coordinates and an intersection type;

[0009] based on the intersection type, determining seismic trace data in a preset spatial window centered on the intersection point coordinates;

[0010] based on the seismic trace data in the preset spatial window, determining seismic trace data of the intersection point coordinates;

[0011] based on the seismic trace data of a plurality of intersection point coordinates included in the seismic arbitrary line, plotting a seismic profile;

[0012] based on the seismic profile, determining a stratum dip angle.

[0013] In some embodiments, the determining the intersection information of the plurality of straight line segments and the seismic grid comprises:

[0014] For any straight line segment, determining a coordinate of an intersection point of the straight line segment and the seismic grid;

[0015] determining a first number and a second number, the first number being a number of intersection points of the straight line segment and a main line included in the seismic grid, and the second number being a number of intersection points of the straight line segment and a tie line included in the seismic grid;

[0016] based on the first number and the second number, determining an intersection type of the straight line segment and the seismic grid.

[0017] In some embodiments, the determining the intersection type of the straight line segment and the seismic grid based on the first number and the second number comprises:

[0018] in a case where the first number is greater than the second number, determining that the intersection type of the straight line segment and the seismic grid is a first oblique intersection;

[0019] in a case where the first number is less than the second number, determining that the intersection type of the straight line segment and the seismic grid is a second oblique intersection.

[0020] In some embodiments, the determining seismic trace data in a preset spatial window centered on the coordinate of the intersection point based on the intersection type comprises:

[0021] in a case where the intersection type is the first oblique intersection, searching for seismic trace data in a preset spatial window to both sides in a direction of a main line of the seismic grid centered on the coordinate of the intersection point;

[0022] in a case where the intersection type is the second oblique intersection, searching for seismic trace data in a preset spatial window to both sides in a direction of a tie line of the seismic grid centered on the coordinate of the intersection point.

[0023] In some embodiments, the preset spatial window includes a plurality of seismic traces;

[0024] the determining seismic trace data of the coordinate of the intersection point based on the seismic trace data in the preset spatial window comprises:

[0025] encrypting seismic trace data of the plurality of seismic traces in a F-K domain to obtain seismic trace data of the coordinate of the intersection point.

[0026] In some embodiments, the preset spatial window includes a plurality of seismic traces;

[0027] determine the seismic trace data of the intersection point coordinate based on the seismic trace data in the preset spatial window.

[0028] determine the weight of the plurality of seismic traces;

[0029] perform weighted summation on the seismic trace data of the plurality of seismic traces based on the weight of the plurality of seismic traces to obtain the seismic trace data of the intersection point coordinate.

[0030] In another aspect, an apparatus for determining a formation dip is provided, and the apparatus comprises:

[0031] a first determining module configured to determine a horizontal well design trajectory to be drilled;

[0032] a second determining module configured to determine a seismic arbitrary line based on the horizontal well design trajectory, and divide the seismic arbitrary line into a plurality of straight line segments;

[0033] a third determining module configured to determine intersection information of the plurality of straight line segments and a seismic grid, the intersection information of any straight line segment and the seismic grid comprising an intersection point coordinate and an intersection type;

[0034] a fourth determining module configured to determine seismic trace data in a preset spatial window centered on the intersection point coordinate based on the intersection type;

[0035] a fifth determining module configured to determine the seismic trace data of the intersection point coordinate based on the seismic trace data in the preset spatial window;

[0036] a drawing module configured to draw a seismic profile based on the seismic trace data of the plurality of intersection point coordinates included in the seismic arbitrary line;

[0037] a sixth determining module configured to determine a formation dip based on the seismic profile.

[0038] In some embodiments, the third determining module is configured to, for any straight line segment, determine intersection point coordinates of the straight line segment and the seismic grid; determine a first number and a second number, the first number being a number of intersection points of the straight line segment and a main survey line included in the seismic grid, and the second number being a number of intersection points of the straight line segment and a tie line included in the seismic grid; and determine the intersection type of the straight line segment and the seismic grid based on the first number and the second number.

[0039] In some embodiments, the third determining module is configured to, in a case where the first number is greater than the second number, determine the intersection type of the straight line segment and the seismic grid to be a first oblique intersection; and in a case where the first number is less than the second number, determine the intersection type of the straight line segment and the seismic grid to be a second oblique intersection.

[0040] In some embodiments, the fourth determining module is configured to, in the case of the first oblique intersection, search for seismic trace data within a preset spatial window to both sides in the main line direction of the seismic grid with the intersection point coordinate as the center; and in the case of the second oblique intersection, search for seismic trace data within a preset spatial window to both sides in the connecting line direction of the seismic grid with the intersection point coordinate as the center.

[0041] In some embodiments, the preset spatial window includes a plurality of seismic traces.

[0042] The fifth determining module is configured to encrypt seismic trace data of the plurality of seismic traces in the F-K domain to obtain the seismic trace data of the intersection point coordinate.

[0043] In some embodiments, the preset spatial window includes a plurality of seismic traces.

[0044] The fifth determining module is configured to determine weights of the plurality of seismic traces, and perform weighted summation on seismic trace data of the plurality of seismic traces based on the weights of the plurality of seismic traces to obtain the seismic trace data of the intersection point coordinate.

[0045] In another aspect, a computer device is provided, which includes one or more processors and one or more memories, and the one or more memories store at least one piece of program code, which is loaded and executed by the one or more processors to implement the method for determining a stratigraphic dip angle according to any of the possible implementation manners.

[0046] In another aspect, a computer readable storage medium is provided, which stores at least one piece of program code, which is loaded and executed by a processor to implement the method for determining a stratigraphic dip angle according to any of the possible implementation manners.

[0047] In another aspect, a computer program product is provided, which stores at least one piece of program code, which is loaded and executed by a processor to implement the method for determining a stratigraphic dip angle according to any of the possible implementation manners.

[0048] The technical scheme provided by the embodiments of the present application has at least the following beneficial effects:

[0049] In an embodiment of the present application, by determining the intersection coordinates and intersection type of the seismic arbitrary line and the seismic grid, the seismic trace data within a preset spatial window centered on the intersection coordinates can be determined based on the intersection coordinates and the intersection type. Based on these seismic trace data, the seismic trace data of the intersection coordinates can be accurately calculated. It can be seen that the present application accurately calculates the seismic trace data of the intersection coordinates, which can eliminate the "sawtooth" phenomenon caused by the seismic profile drawn using the nearest borrowing method. Therefore, the present application can improve the accuracy of the seismic profile determined by the seismic trace data based on the intersection coordinates, and thereby improve the accuracy of the stratum dip determined based on the seismic profile. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0051] Figure 1 This is a flow chart of a method for determining formation dip provided in an embodiment of the present application;

[0052] Figure 2 This is a flow chart of another method for determining formation dip provided by an embodiment of the present application;

[0053] Figure 3 This is a schematic diagram of the intersection of a designed trajectory of a horizontal well to be drilled and a seismic grid and the target layer structure provided by an embodiment of the present application;

[0054] Figure 4 is a schematic diagram of a seismic profile before optimization provided in an embodiment of the present application;

[0055] Figure 5 This is a schematic diagram of the intersection of a straight line segment and a seismic grid and the location points close to the intersection, provided in an embodiment of the present application;

[0056] Figure 6 This is a schematic diagram of the location points of the actual removal of seismic trace data after optimization provided by an embodiment of the present application;

[0057] Figure 7 is a comparative schematic diagram of a seismic profile provided in an embodiment of the present application;

[0058] Figure 8 is a schematic diagram of determining a seismic profile provided in an embodiment of the present application;

[0059] Figure 9 is a schematic diagram of determining a formation dip provided in an embodiment of the present application;

[0060] Figure 10 is a structural schematic diagram of a device for determining a formation dip angle provided by an embodiment of the present application;

[0061] Figure 11 is a structural schematic diagram of a server provided by an embodiment of the present application;

[0062] Figure 12 is a structural schematic diagram of a terminal provided by an embodiment of the present application. DETAILED DESCRIPTION

[0063] In order to make the objects, technical solutions and advantages of the present application clearer, the embodiments of the present application will be further described in detail below with reference to the drawings.

[0064] The terms "first", "second", "third", and "fourth" and the like in the specification and claims of the present application and the drawings are used to distinguish different objects, and are not used to describe a particular order. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product or device including a series of steps or units is not limited to the listed steps or units, but can optionally include other steps or units not listed or can optionally include other steps or units inherent to the process, method, product or device.

[0065] It should be noted that the information (including but not limited to user device information, user personal information, etc.), data (including but not limited to data for analysis, stored data, displayed data, etc.) and signals involved in the present application are all authorized by the user or fully authorized by all parties, and the collection, use and processing of related data need to comply with relevant laws, regulations and standards of relevant countries and regions. For example, the seismic data involved in the present application is obtained under full authorization.

[0066] The execution subject of the method for determining a formation dip angle provided by the embodiments of the present application is a computer device; the computer device can be a terminal or a server. In some embodiments, the terminal can be a smart watch, a smart phone, a tablet computer, a notebook computer, a desktop computer, etc., but is not limited thereto. The server is an independent physical server, can also be a server cluster or a distributed system composed of multiple physical servers, can also be a cloud server providing cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, CDN (Content Delivery Network), and basic cloud computing services such as big data and artificial intelligence platforms.

[0067] Figure 1This is a flow chart of a method for determining formation dip provided by an embodiment of the present application. The execution subject of this method may be a computer device; see Figure 1 , the method comprising:

[0068] Step 101: Determine the design trajectory of the horizontal well to be drilled.

[0069] In some embodiments, the computer device determines exploration data for the target horizontal well, the exploration data including a planned trajectory for the horizontal well to be drilled. In this step, the computer device obtains the planned trajectory for the horizontal well to be drilled from the exploration data. In this embodiment of the present application, formation dip is predicted to provide a basis for drilling.

[0070] Step 102: Based on the designed horizontal well trajectory, determine a seismic arbitrary line, and divide the seismic arbitrary line into multiple straight line segments.

[0071] In some embodiments, the computer device extracts a seismic arbitrary line along the designed horizontal well trajectory; for example, the computer device determines a well connection line based on the designed horizontal well trajectory, and determines the well connection line as the seismic arbitrary line.

[0072] In some embodiments, the computer device can divide the seismic arbitrary line into multiple uniform straight line segments, that is, the length of each straight line segment is equal; or, the computer device can divide the seismic arbitrary line into multiple uneven straight line segments, that is, the length of each straight line segment is not equal.

[0073] Step 103: Determine intersection information between a plurality of straight line segments and the seismic grid. The intersection information between any straight line segment and the seismic grid includes intersection coordinates and intersection type.

[0074] The seismic grid includes a main survey line and a tie line, and the main survey line is perpendicular to the tie line; for example, if the main survey line is a horizontal line, the tie line is a vertical line. The intersection information of any straight line segment and the seismic grid includes the intersection coordinates of the intersection of any straight line segment and the main survey line and / or the intersection coordinates of any straight line segment and the tie line; and the intersection of any straight line segment and the seismic grid can be one or more. The intersection type can be a first oblique intersection or a second oblique intersection. The first oblique intersection is used to indicate that the number of intersections between the straight line segment and the main survey line is greater than the number of intersections between the straight line segment and the tie line, and the second oblique intersection is used to indicate that the number of intersections between the straight line segment and the tie line is greater than the number of intersections between the straight line segment and the main survey line.

[0075] Step 104: Based on the intersection type, determine the seismic trace data within a preset spatial window centered on the intersection coordinates.

[0076] The computer device determines a target direction matching the intersection type; and searches for seismic trace data in a preset space centered on the intersection point coordinate and in the target direction. Different intersection types correspond to different search directions; for example, the first oblique intersection corresponds to the main line direction; and the second oblique intersection corresponds to the connecting line direction. When the intersection type is the first oblique intersection, the computer device searches for seismic trace data in a preset space window in the main line direction of the seismic grid and on both sides of the intersection point coordinate; and when the intersection type is the second oblique intersection, the computer device searches for seismic trace data in a preset space window in the connecting line direction of the seismic grid and on both sides of the intersection point coordinate. The seismic trace data in the preset space window includes seismic trace data of at least one seismic trace. The size of the preset space window can be set and changed as needed, and in the embodiments of the present application, the size of the preset space window is not specifically limited.

[0077] It should be noted that when there are multiple intersection point coordinates, for each intersection point coordinate, the seismic trace data in a preset space window centered on the intersection point coordinate is determined based on the intersection type.

[0078] Step 105: Based on the seismic trace data in the preset space window, the seismic trace data of the intersection point coordinate is determined.

[0079] The seismic trace data in the preset space window includes seismic trace data of at least one seismic trace, and in this step, the seismic trace data of at least one seismic trace is operated to obtain one seismic trace data, which is the seismic trace data of the intersection point coordinate.

[0080] Step 106: Based on the seismic trace data of multiple intersection point coordinates included in the seismic arbitrary line, a seismic profile is drawn.

[0081] The seismic profile can be a depth profile; and the seismic profile can be obtained in two ways, one way is to obtain it from a time profile by time-depth conversion, and the other way is to obtain it by depth migration.

[0082] Step 107: Based on the seismic profile, the stratigraphic dip angle is determined.

[0083] The computer device determines the height difference h of the target layer and the actual planar distance l of the adjacent two seismic traces based on the seismic profile, and determines the stratigraphic dip angle a based on the height difference h and the actual planar distance l; for example, the computer device determines the stratigraphic dip angle a based on h and l by the following formula one;

[0084]

[0085] The computer device determines the height difference h of the target layer based on the seismic profile, which can include: the computer device determines the height difference of the seismic similar waveform based on the seismic profile, and the height difference is the height difference h of the target layer; the height difference of the seismic similar waveform can be the height difference of two wave peaks, the height difference of two wave troughs, or the height difference between two zero-crossing points.

[0086] In the embodiments of the present application, by determining the intersection coordinates and intersection types of the seismic arbitrary line and the seismic grid, the seismic trace data in the preset spatial window centered on the intersection coordinates can be determined according to the intersection coordinates and the intersection types, and the seismic trace data of the intersection coordinates can be accurately calculated according to the seismic trace data. As can be seen, the present application accurately calculates the seismic trace data of the intersection coordinates, which can eliminate the "jagged" phenomenon caused by the nearest borrowing method of the seismic profile. Therefore, the present application can improve the accuracy of the seismic profile determined based on the seismic trace data of the intersection coordinates, and further improve the accuracy of the stratigraphic dip determined based on the seismic profile.

[0087] Figure 2 is a flowchart of a method for determining a stratigraphic dip provided by the embodiments of the present application, referring to Figure 2 , the method includes:

[0088] Step 201: The computer device determines a horizontal well design trajectory to be drilled.

[0089] In some embodiments, this step is the same as step 101, which will not be described here.

[0090] For example, the dashed line in Figure 3 , Figure 3 is the horizontal well design trajectory, Figure 3 the grid in is the seismic grid. Figure 3 the solid line in is the target layer, which is to determine the stratigraphic dip of the target layer in the embodiments of the present application. As can be seen in Figure 3 , the target layer structure tends to be parallel to the contact line direction of the seismic grid, and the horizontal well to be drilled is projected on the plane and is oblique to the direction of the seismic grid. If the borrowing method in the related art is used to extract the seismic trace data, the seismic trace on the seismic arbitrary line will be projected onto the seismic profile, and the stratigraphic structure contour direction is oblique to the seismic arbitrary line direction, which causes the "jagged" phenomenon of the seismic profile, as shown in Figure 4 , and Figure 4 cannot provide the function of predicting the stratigraphic dip in the horizontal well drilling.

[0091] Step 202: The computer device determines a seismic arbitrary line based on the horizontal well design trajectory, and divides the seismic arbitrary line into a plurality of straight line segments.

[0092] In some embodiments, this step is the same as step 102, which will not be repeated here.

[0093] For example, the computer device divides the seismic arbitrary line into multiple straight line segments, takes out one straight line segment, such as Figure 5 , and draws it as shown in Figure 5 . The black dots in Figure 5 are the positions of the actual extracted seismic trace data before optimization. The black dots in Figure 6 are the positions of the seismic trace data we want to extract.

[0094] Step 203: For any straight line segment, the computer device determines the intersection coordinates of the straight line segment and the seismic grid.

[0095] The seismic grid is a depth domain seismic grid, and the seismic grid includes a main line and a tie line, and the main line is perpendicular to the tie line; for example, the main line is a horizontal line, and the tie line is a vertical line. Any straight line segment can have one or more intersections with the seismic grid; accordingly, the computer device will determine one or more intersection coordinates in this step.

[0096] Step 204: The computer device determines a first number and a second number, the first number being the number of intersections of the straight line segment with the main line included in the seismic grid, and the second number being the number of intersections of the straight line segment with the tie line included in the seismic grid.

[0097] Step 205: The computer device determines the intersection type of the straight line segment and the seismic grid based on the first number and the second number.

[0098] In the case where the first number is greater than the second number, the computer device determines that the intersection type of the straight line segment and the seismic grid is the first oblique intersection; in the case where the first number is less than the second number, the computer device determines that the intersection type of the straight line segment and the seismic grid is the second oblique intersection.

[0099] In some embodiments, in the case where the first number is equal to the second number, the computer device determines that the straight line segment is along the diagonal direction; in the case where the first number is 0 and the second number is not 0, the computer device determines that the straight line segment is along the main line direction; in the case where the second number is 0 and the first number is not 0, the computer device determines that the straight line segment is along the tie line direction. In addition, in the case where the computer device determines that the straight line segment is along the diagonal direction, or in the case where the computer determines that the straight line segment is along the main line direction, or in the case where the computer device determines that the straight line segment is along the tie line direction, the computer device does not perform special processing, but in the case where the computer device determines that the intersection type of the straight line segment and the seismic grid is the first oblique intersection or the second oblique intersection, the computer device executes step 206.

[0100] Step 206: The computer device determines the seismic trace data in a preset spatial window centered on the intersection coordinates based on the intersection type.

[0101] In the case of the first oblique intersection, the computer device searches for seismic trace data within a preset spatial window to the left and right of the intersection point coordinate in the main line direction of the seismic grid. In the case of the second oblique intersection, the computer device searches for seismic trace data within a preset spatial window to the left and right of the intersection point coordinate in the connecting line direction of the seismic grid.

[0102] For example, continuing to refer to Figure 4 , the computer device searches for seismic trace data within a preset spatial window to the left and right of the intersection point coordinate in the connecting line direction of the seismic grid, Figure 4 The block in

[0103] Step 207: The computer device determines the seismic trace data of the intersection point coordinate based on the seismic trace data within the preset spatial window.

[0104] The preset spatial window includes multiple seismic traces. In some embodiments, the seismic trace data of the intersection point is obtained by encrypting the searched seismic trace data within the preset spatial window in the F-K domain, and the coordinates of the data are set as the coordinates of the intersection point. Correspondingly, step 207 can be: the computer device encrypts the seismic trace data of the multiple seismic traces in the F-K domain to obtain the seismic trace data of the intersection point coordinate.

[0105] The computer device obtains the seismic trace data of the interpolation seismic trace by interpolating the seismic trace data of the adjacent two seismic traces between the adjacent two seismic traces, and then obtains the seismic trace data of the further interpolation seismic trace by interpolating the seismic trace data of the adjacent two interpolation seismic traces between the adjacent two interpolation seismic traces, and then obtains the seismic trace data of the still further interpolation seismic trace by interpolating the seismic trace data of the adjacent two still further interpolation seismic traces between the adjacent two still further interpolation seismic traces, and so on, until the distance between the final interpolation position and the intersection position is less than a preset distance position, at which time the seismic trace data of the interpolation seismic trace is the seismic trace data of the intersection point coordinate.

[0106] For example, the preset window includes 4 seismic traces, namely seismic trace 1, seismic trace 2, seismic trace 3 and seismic trace 4, and the intersection coordinate is located between seismic trace 2 and seismic trace 3; in this step, the computer device interpolates seismic trace data between seismic trace 1 and seismic trace 2 based on the seismic trace data of seismic trace 1 and the seismic trace data of seismic trace 2; for example, the seismic trace data of the interpolated seismic trace 0.5; then interpolates seismic trace data between seismic trace 3 and seismic trace 4 based on the seismic trace data of seismic trace 3 and the seismic trace data of seismic trace 4; for example, the seismic trace data of the interpolated seismic trace 3.5; then interpolates seismic trace data between seismic trace 0.5 and seismic trace 3.5 based on the seismic trace data of seismic trace 0.5 and the seismic trace data of seismic trace 3.5; for example, the seismic trace data of the interpolated seismic trace 2.5, which is the seismic trace data of the intersection coordinate.

[0107] The application provides a method for optimizing the use of seismic prediction of formation dip in horizontal well drilling.

[0108] In some embodiments, the seismic trace data of the intersection coordinate is obtained by weighted summation of the seismic trace data of the plurality of seismic traces; correspondingly, step 207 can be implemented by the following steps (1) to (2), comprising:

[0109] (1) The computer device determines the weights of the plurality of seismic traces.

[0110] The sum of the weights of the plurality of seismic traces is 1, and the weights of the plurality of seismic traces are not completely the same. In some embodiments, the computer device determines the perpendicular distances of the plurality of seismic traces from the intersection coordinate, and for each seismic trace, determines a weight matching the perpendicular distance of the seismic trace from the intersection coordinate, and the perpendicular distance and the weight are negatively correlated, that is, the farther the perpendicular distance, the lower the weight.

[0111] For example, the preset window includes 4 seismic traces, namely seismic trace 1, seismic trace 2, seismic trace 3 and seismic trace 4; the perpendicular distances of seismic trace 1, seismic trace 2, seismic trace 3 and seismic trace 4 from the intersection coordinate are 20 meters, 15 meters, 10 meters and 5 meters respectively, and the weights of seismic trace 1, seismic trace 2, seismic trace 3 and seismic trace 4 are 0.1, 0.2, 0.3 and 0.4 respectively.

[0112] (2) The computer device performs weighted summation on the seismic trace data of the plurality of seismic traces based on the weights of the plurality of seismic traces to obtain the seismic trace data of the intersection coordinate.

[0113] For each seismic trace, the computer device determines the product of the weight of the seismic trace and the seismic trace data of the seismic trace, obtains the weighted seismic trace data of the seismic trace, sums the weighted seismic trace data of multiple seismic traces, and obtains the seismic trace data of the intersection point coordinates.

[0114] In the embodiments of the present application, the computer device determines the seismic trace data of the intersection point coordinates by weighted summation of the seismic trace data of multiple seismic traces, can determine the seismic trace data of the intersection point coordinates in combination with the influence degree of each seismic trace on the seismic trace data of the intersection point coordinates, and thus improves the accuracy of the determined seismic trace data of the intersection point coordinates.

[0115] Step 208: The computer device draws a seismic profile based on the seismic trace data of the intersection point coordinates included in the seismic arbitrary line.

[0116] In some embodiments, the present step is the same as step 106, which will not be described here.

[0117] For example, referring to Figure 7 , the computer device determines the seismic profile based on the seismic trace data extracted by the borrowing method in the related art. Figure 7 The left side figure is the seismic profile determined based on the way of extracting seismic trace data by the borrowing method in the related art, from which it can be seen that the seismic profile has obvious discontinuity, that is, the seismic profile has a "sawtooth" phenomenon. Figure 7 The right side figure is based on the method of the present application. Since the computer device encrypts the seismic trace data of multiple seismic traces in the F-K domain to obtain the seismic trace data of the intersection point coordinates, and uses an accurate interpolation algorithm, the accuracy is improved, so the phenomenon of stratigraphic discontinuity no longer appears, that is, the "sawtooth" phenomenon of the seismic profile no longer appears. Therefore, the seismic profile can be calculated for the stratigraphic dip to provide drilling reference for horizontal well implementation.

[0118] It should be noted that the computer device divides the seismic arbitrary line into multiple straight line segments, and performs the above steps 203-208 for each straight line segment, so that the above steps 203-208 are repeatedly performed for the multiple straight line segments along the horizontal well design trajectory. Figure 3 The multiple straight line segments on the horizontal well design trajectory are processed in this way to obtain the seismic profile of the seismic arbitrary line along the horizontal well design trajectory, which can eliminate the "sawtooth" phenomenon in the seismic profile. Figure 4 , so that the stratigraphic occurrence is more reasonable and tends to be the real underground condition, and the obtained seismic profile is shown in Figure 8 .

[0119] Step 209: The computer device determines the stratigraphic dip based on the seismic profile.

[0120] In some embodiments, the present step is the same as step 107, which will not be described here.

[0121] For example, referring to Figure 9 , the computer device determines the stratigraphic dip based on the seismic profileFigure 8 In the seismic profile, along the target layer ( Figure 9 The white line in FIG1 ), between each two adjacent seismic traces, the height difference between the two peaks and the actual plane distance l between the adjacent seismic traces are used to determine the dip angle of the formation through the trigonometric formula (Formula 1 in step 107).

[0122] In an embodiment of the present application, by determining the intersection coordinates and intersection type of the seismic arbitrary line and the seismic grid, the seismic trace data within a preset spatial window centered on the intersection coordinates can be determined based on the intersection coordinates and the intersection type. Based on these seismic trace data, the seismic trace data of the intersection coordinates can be accurately calculated. It can be seen that the present application accurately calculates the seismic trace data of the intersection coordinates, which can eliminate the "sawtooth" phenomenon caused by the seismic profile drawn using the nearest borrowing method. Therefore, the present application can improve the accuracy of the seismic profile determined by the seismic trace data based on the intersection coordinates, and thereby improve the accuracy of the stratum dip determined based on the seismic profile.

[0123] Figure 10 This is a schematic diagram of the structure of a device for determining the dip angle of a formation provided in an embodiment of the present application, see Figure 10 , the device comprises:

[0124] The first determination module 1001 is used to determine the design trajectory of the horizontal well to be drilled;

[0125] The second determining module 1002 is configured to determine a seismic arbitrary line based on the horizontal well design trajectory, and divide the seismic arbitrary line into a plurality of straight line segments;

[0126] The third determining module 1003 is used to determine the intersection information of multiple straight line segments and the seismic grid, where the intersection information of any straight line segment and the seismic grid includes the intersection coordinates and the intersection type;

[0127] A fourth determining module 1004 is configured to determine, based on the intersection type, seismic trace data within a preset spatial window centered on the intersection coordinates;

[0128] A fifth determining module 1005 is configured to determine the seismic trace data of the intersection coordinates based on the seismic trace data within the preset spatial window;

[0129] A drawing module 1006 is used to draw a seismic profile based on the seismic trace data of the coordinates of the multiple intersection points included in the seismic arbitrary line;

[0130] The sixth determination module 1007 is used to determine the formation dip based on the seismic profile.

[0131] In some embodiments, the third determination module 1003 is configured to determine, for any straight line segment, intersection point coordinates of the straight line segment and the seismic grid; determine a first quantity and a second quantity, the first quantity being a number of intersection points of the straight line segment and a main survey line included in the seismic grid, and the second quantity being a number of intersection points of the straight line segment and a tie line included in the seismic grid; and determine, based on the first quantity and the second quantity, an intersection type of the straight line segment and the seismic grid.

[0132] In some embodiments, the third determination module 1003 is configured to determine, in a case where the first quantity is greater than the second quantity, that the intersection type of the straight line segment and the seismic grid is a first oblique intersection; and determine, in a case where the first quantity is less than the second quantity, that the intersection type of the straight line segment and the seismic grid is a second oblique intersection.

[0133] In some embodiments, the fourth determination module 1004 is configured to search, in a case where the intersection type is the first oblique intersection, for seismic trace data within a preset spatial window to both sides in a direction of the main survey line of the seismic grid, with the intersection point coordinates as a center; and search, in a case where the intersection type is the second oblique intersection, for seismic trace data within the preset spatial window to both sides in a direction of the tie line of the seismic grid, with the intersection point coordinates as the center.

[0134] In some embodiments, the preset spatial window includes a plurality of seismic traces.

[0135] The fifth determination module 1005 is configured to encrypt seismic trace data of the plurality of seismic traces in an F-K domain to obtain the intersection point coordinates of the seismic trace data.

[0136] In some embodiments, the preset spatial window includes a plurality of seismic traces.

[0137] The fifth determination module 1005 is configured to determine weights of the plurality of seismic traces; and perform weighted summation on seismic trace data of the plurality of seismic traces based on the weights of the plurality of seismic traces to obtain the intersection point coordinates of the seismic trace data.

[0138] In the embodiments of the present application, by determining intersection point coordinates and an intersection type of a seismic arbitrary line and a seismic grid, the seismic trace data within a preset spatial window with the intersection point coordinates as a center can be determined according to the intersection point coordinates and the intersection type, and the intersection point coordinates of the seismic trace data can be accurately calculated according to the seismic trace data. It can be seen that the present application is to accurately calculate the intersection point coordinates of the seismic trace data, which can eliminate the "sawtooth" phenomenon caused by the drawing mode of borrowing nearby traces of the seismic profile. Therefore, the present application can improve the accuracy of the seismic profile determined based on the intersection point coordinates of the seismic trace data, and further improve the accuracy of the stratum dip angle determined based on the seismic profile.

[0139] It should be noted that the device for determining the dip angle of the stratum provided in the above embodiment is only used for determining the dip angle of the stratum by way of example, and in actual application, the above functions can be completed by different functional modules according to needs, that is, the internal structure of the computer device is divided into different functional modules to complete all or part of the above functions. In addition, the device for determining the dip angle of the stratum and the method for determining the dip angle of the stratum provided in the above embodiment belong to the same concept, and the specific implementation process is described in the method embodiment, which will not be repeated here.

[0140] In some embodiments, in the case of the computer device being a server, referring to Figure 11 , Figure 11 is a structural schematic diagram of a server provided by an embodiment of the present application. The server 1100 can have great differences due to different configurations or performances, and can include one or more processors (central processing units, CPUs) 1101 and one or more memories 1102, wherein the memory 1102 stores at least one program code, the at least one program code is loaded and executed by the processor 1101 to realize the above method for determining the dip angle of the stratum. Of course, the server 1100 can also have a wired or wireless network interface, a keyboard, and an input and output interface and other components for realizing the functions of the device, so as to perform input and output. The server 1100 can also include other components for realizing the functions of the device, which will not be described here.

[0141] In some embodiments, in the case of the computer device being a terminal, Figure 12 A structural block diagram of a terminal 1200 provided by an exemplary embodiment of the present application is shown. The terminal 1200 can be a smart phone, a tablet computer, a VR device, or a vehicle-mounted terminal, etc. Generally, the terminal 1200 includes a processor 1201 and a memory 1202.

[0142] The processor 1201 can include one or more processing cores, such as a 4-core processor, an 8-core processor, and the like. The processor 1201 can be implemented in at least one of a hardware form of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), a PLA (Programmable Logic Array). The processor 1201 can also include a main processor and a co-processor, the main processor being a processor for processing data in an awake state, also referred to as a CPU (Central Processing Unit), and the co-processor being a low-power processor for processing data in a standby state. In some embodiments, the processor 1201 can be integrated with a GPU (Graphics Processing Unit) for rendering and drawing content required to be displayed by the display screen. In some embodiments, the processor 1201 can further include an AI (Artificial Intelligence) processor for processing computing operations related to machine learning.

[0143] The memory 1202 can include one or more computer-readable storage media that can be non-transitory. The memory 1202 can also include a high-speed random access memory, and a nonvolatile memory such as one or more disk storage devices, flash storage devices. In some embodiments, the non-transitory computer-readable storage medium in the memory 1202 is used to store at least one program code for being executed by the processor 1201 to implement the method of determining the dip angle of the stratum provided by the method embodiments in the present application.

[0144] In some embodiments, the terminal 1200 can also optionally include a peripheral device interface 1203 and at least one peripheral device. The processor 1201, the memory 1202, and the peripheral device interface 1203 can be connected through a bus or a signal line. Each peripheral device can be connected to the peripheral device interface 1203 through a bus, a signal line, or a circuit board. Specifically, the peripheral device includes at least one of a radio frequency circuit 1204, a display screen 1205, a camera assembly 1206, an audio circuit 1207, and a power supply 1208.

[0145] The peripheral interface 1203 can be used to connect at least one I / O (Input / Output) related peripheral device to the processor 1201 and the memory 1202. In some embodiments, the processor 1201, the memory 1202 and the peripheral interface 1203 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 1201, the memory 1202 and the peripheral interface 1203 can be implemented on a separate chip or circuit board, and the present embodiments are not limited in this regard.

[0146] The radio frequency circuit 1204 is configured to receive and send RF (Radio Frequency) signals, also known as electromagnetic signals. The radio frequency circuit 1204 communicates with communication networks and other communication devices through electromagnetic signals. The radio frequency circuit 1204 converts electric signals into electromagnetic signals for transmission, or converts electromagnetic signals received into electric signals. Optionally, the radio frequency circuit 1204 includes an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a subscriber identity module card, and the like. The radio frequency circuit 1204 can communicate with other terminals through at least one wireless communication protocol. The wireless communication protocol includes, but is not limited to, a metropolitan area network, various generations of mobile communication networks (2G, 3G, 4G and 5G), a wireless local area network and / or a WiFi (Wireless Fidelity) network. In some embodiments, the radio frequency circuit 1204 can also include NFC (Near Field Communication) related circuitry, and the present application is not limited in this regard.

[0147] The display screen 1205 is configured to display a UI (User Interface). The UI can include graphics, text, icons, video, and any combination thereof. When the display screen 1205 is a touch display screen, the display screen 1205 is further configured to capture touch signals on or above the surface of the display screen 1205. The touch signals can be input to the processor 1201 as control signals for processing. In this case, the display screen 1205 can also be configured to provide virtual buttons and / or virtual keyboard, also known as soft buttons and / or soft keyboard. In some embodiments, the display screen 1205 can be one, arranged on the front panel of the terminal 1200; in other embodiments, the display screen 1205 can be at least two, arranged on different surfaces of the terminal 1200 or in a folding design; in still other embodiments, the display screen 1205 can be a flexible display screen, arranged on a curved surface or a folding surface of the terminal 1200. Even, the display screen 1205 can also be arranged in an irregular shape, i.e., a special-shaped screen. The display screen 1205 can be made of LCD (Liquid Crystal Display), OLED (Organic Light-Emitting Diode), etc.

[0148] The camera assembly 1206 is configured to capture images or videos. Optionally, the camera assembly 1206 includes a front-facing camera and a rear-facing camera. Typically, the front-facing camera is arranged on the front panel of the terminal, and the rear-facing camera is arranged on the back of the terminal. In some embodiments, the rear-facing camera is at least two, which is any one of a main camera, a depth-of-field camera, a wide-angle camera, and a telephoto camera, to realize the background blur function by fusing the main camera and the depth-of-field camera, the panoramic shooting and VR (Virtual Reality) shooting function by fusing the main camera and the wide-angle camera, or other fusion shooting functions. In some embodiments, the camera assembly 1206 can further include a flash. The flash can be a single-color-temperature flash or a dual-color-temperature flash. The dual-color-temperature flash refers to the combination of a warm light flash and a cold light flash, which can be used for light compensation under different color temperatures.

[0149] The audio circuit 1207 can include a microphone and a speaker. The microphone is used to collect sound waves of the user and the environment, and convert the sound waves into an electrical signal input to the processor 1201 for processing, or input to the radio frequency circuit 1204 to realize voice communication. For the purpose of stereo sound collection or noise reduction, the microphone can be multiple, respectively arranged at different parts of the terminal 1200. The microphone can also be an array microphone or an omnidirectional collection type microphone. The speaker is used to convert the electrical signal from the processor 1201 or the radio frequency circuit 1204 into sound waves. The speaker can be a traditional diaphragm speaker, or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, not only can the electrical signal be converted into a sound wave that humans can hear, but also can be converted into a sound wave that humans cannot hear for ranging purposes. In some embodiments, the audio circuit 1207 can also include a headphone jack.

[0150] The power supply 1208 is used to supply power to each component in the terminal 1200. The power supply 1208 can be alternating current, direct current, disposable battery or rechargeable battery. When the power supply 1208 includes a rechargeable battery, the rechargeable battery can support wired charging or wireless charging. The rechargeable battery can also be used to support fast charging technology.

[0151] In some embodiments, the terminal 1200 also includes one or more sensors 1209. The one or more sensors 1209 include, but are not limited to, an acceleration sensor 1210, a gyroscope sensor 1211, a pressure sensor 1212, an optical sensor 1213, and a proximity sensor 1214.

[0152] The acceleration sensor 1210 can detect the acceleration magnitude in three coordinate axes of the coordinate system established by the terminal 1200. For example, the acceleration sensor 1210 can be used to detect the components of the gravitational acceleration in three coordinate axes. The processor 1201 can control the display screen 1205 to display the user interface in a landscape view or a portrait view according to the gravitational acceleration signal collected by the acceleration sensor 1210. The acceleration sensor 1210 can also be used for game or user motion data collection.

[0153] The gyroscope sensor 1211 can detect the body direction and rotation angle of the terminal 1200, and the gyroscope sensor 1211 can collect 3D actions of the user on the terminal 1200 in cooperation with the acceleration sensor 1210. The processor 1201 can realize the following functions according to the data collected by the gyroscope sensor 1211: motion sensing (such as changing the UI according to the user's tilt operation), image stabilization when shooting, game control, and inertial navigation.

[0154] The pressure sensor 1212 can be disposed at the side frame of the terminal 1200 and / or the lower layer of the display screen 1205. When the pressure sensor 1212 is disposed at the side frame of the terminal 1200, the holding signal of the user to the terminal 1200 can be detected, and the left-hand or right-hand recognition or the shortcut operation can be performed by the processor 1201 according to the holding signal collected by the pressure sensor 1212. When the pressure sensor 1212 is disposed at the lower layer of the display screen 1205, the operability control on the UI interface can be controlled by the processor 1201 according to the pressure operation of the user to the display screen 1205. The operability control includes at least one of the button control, the scroll bar control, the icon control, and the menu control.

[0155] The optical sensor 1213 is used to collect the ambient light intensity. In one embodiment, the processor 1201 can control the display brightness of the display screen 1205 according to the ambient light intensity collected by the optical sensor 1213. Specifically, when the ambient light intensity is high, the display brightness of the display screen 1205 is increased; when the ambient light intensity is low, the display brightness of the display screen 1205 is decreased. In another embodiment, the processor 1201 can also dynamically adjust the shooting parameter of the camera assembly 1206 according to the ambient light intensity collected by the optical sensor 1213.

[0156] The proximity sensor 1214, also called the distance sensor, is usually disposed at the front panel of the terminal 1200. The proximity sensor 1214 is used to collect the distance between the user and the front of the terminal 1200. In one embodiment, when the proximity sensor 1214 detects that the distance between the user and the front of the terminal 1200 gradually decreases, the display screen 1205 is switched from the bright screen state to the off-screen state by the processor 1201; when the proximity sensor 1214 detects that the distance between the user and the front of the terminal 1200 gradually increases, the display screen 1205 is switched from the off-screen state to the bright screen state by the processor 1201.

[0157] Those skilled in the art can understand that the structures shown in the above embodiments are not a limitation on the terminal 1200, and the terminal 1200 can include more or fewer components than those shown in the figures, or combine certain components, or adopt a different component arrangement. Figure 12

[0158] ​In the example embodiments, the application further provides a computer readable storage medium, which stores at least one program code. The at least one program code is loaded and executed by a processor to implement the method for determining the dip angle of the stratum according to any of the above implementation manners. Optionally, the storage medium can be a non-transitory computer readable storage medium, for example, the non-transitory computer readable storage medium can be a ROM (Read-Only Memory), a RAM (Random Access Memory), a CD-ROM (Compact Disc Read-Only Memory), a magnetic tape, a floppy disk, an optical data storage device, etc.

[0159] In the example embodiments, the application further provides a computer program product, which stores at least one program code. The at least one program code is loaded and executed by a processor to implement the method for determining the dip angle of the stratum according to the application.

[0160] In some embodiments, the computer program related to the application can be deployed to be executed on one computer device, or on multiple computer devices located in one place, or on multiple computer devices distributed in multiple places and interconnected through a communication network, which can constitute a blockchain system.

[0161] Those of ordinary skill in the art can understand that all or part of the steps of the above-mentioned embodiments can be completed by hardware, or by a program instructing relevant hardware, and the program can be stored in a computer readable storage medium. The storage medium mentioned above can be a ROM, a magnetic disk, or an optical disk, etc.

[0162] The above description is only optional embodiments of the application and does not limit the application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the application shall be included in the protection scope of the application.

Claims

1. A method for determining formation dip, characterized in that: The method comprises: Determine the design trajectory of the horizontal well to be drilled; Determine a seismic arbitrary line based on the horizontal well design trajectory, and divide the seismic arbitrary line into a plurality of straight line segments; For any straight line segment, determining the coordinates of the intersection of the straight line segment and the seismic grid; Determining a first number and a second number, wherein the first number is the number of intersections between the straight line segment and the main survey line included in the seismic grid, and the second number is the number of intersections between the straight line segment and the tie line included in the seismic grid; When the first number is greater than the second number, determining that the intersection type of the straight line segment and the seismic grid is a first oblique intersection; When the first number is less than the second number, determining that the intersection type of the straight line segment and the seismic grid is a second oblique intersection; When the intersection type is the first oblique intersection, searching for seismic trace data within a preset spatial window in both sides in the main survey line direction of the seismic grid with the intersection coordinates as the center; When the intersection type is the second oblique intersection, taking the intersection coordinates as the center, searching for seismic trace data within a preset spatial window in both sides in the direction of the tie line of the seismic grid; Determining the seismic trace data of the intersection coordinates based on the seismic trace data within the preset spatial window; Drawing a seismic profile based on seismic trace data of a plurality of intersection coordinates included in the seismic arbitrary line; Based on the seismic profile, the dip of the formation is determined.

2. The method according to claim 1, characterized in that The preset spatial window includes a plurality of seismic traces; The determining of the seismic trace data of the intersection coordinates based on the seismic trace data within the preset spatial window includes: The seismic trace data of the plurality of seismic traces are encrypted in the FK domain to obtain the seismic trace data of the intersection coordinates.

3. The method according to claim 1, characterized in that The preset spatial window includes a plurality of seismic traces; The determining of the seismic trace data of the intersection coordinates based on the seismic trace data within the preset spatial window includes: determining weights of the plurality of seismic traces; Based on the weights of the multiple seismic traces, weighted summation is performed on the seismic trace data of the multiple seismic traces to obtain the seismic trace data of the intersection coordinates.

4. A device for determining the dip of a formation, characterized in that: The device comprises: The first determination module is used to determine the design trajectory of the horizontal well to be drilled; A second determining module is configured to determine a seismic arbitrary line based on the horizontal well design trajectory, and divide the seismic arbitrary line into a plurality of straight line segments; a third determining module configured to determine, for any straight line segment, the coordinates of the intersection of the straight line segment and the seismic grid; determine a first number and a second number, wherein the first number is the number of intersections between the straight line segment and the main survey line included in the seismic grid, and the second number is the number of intersections between the straight line segment and the tie lines included in the seismic grid; if the first number is greater than the second number, determine the intersection type of the straight line segment and the seismic grid to be a first oblique intersection; if the first number is less than the second number, determine the intersection type of the straight line segment and the seismic grid to be a second oblique intersection; A fourth determining module is configured to, when the intersection type is a first oblique intersection, search for seismic trace data within a preset spatial window in both directions in the direction of the main survey line of the seismic grid with the intersection coordinates as the center; and when the intersection type is a second oblique intersection, search for seismic trace data within a preset spatial window in both directions in the direction of the tie line of the seismic grid with the intersection coordinates as the center; a fifth determining module, configured to determine the seismic trace data of the intersection coordinates based on the seismic trace data within the preset spatial window; A drawing module, configured to draw a seismic profile based on seismic trace data of a plurality of intersection coordinates included in the seismic arbitrary line; The sixth determination module is used to determine the formation dip based on the seismic profile.

5. A computer device, characterized in that: The computer device includes one or more processors and one or more memories, wherein at least one program code is stored in the one or more memories, and the at least one program code is loaded and executed by the one or more processors to implement the method for determining the formation inclination as described in any one of claims 1 to 3.

6. A computer-readable storage medium, characterized in that The storage medium stores at least one program code, which is loaded and executed by a processor to implement the method for determining formation inclination according to any one of claims 1 to 3.

7. A computer program product, characterized in that The computer program product stores at least one program code, and the at least one program code is loaded and executed by a processor to implement the method for determining formation inclination according to any one of claims 1 to 3.

Images