Method, device and equipment for establishing small-layer model on stratigraphic section
By calculating the central axis and depositional direction of stratigraphic block units, and performing sub-layer division and splicing under complex geological conditions, the problem of unsatisfactory sub-layer model morphology in traditional methods is solved, and refined sub-layer model construction is realized under conditions of fault development and severe stratigraphic folding is achieved.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- CHINA NAT PETROLEUM CORP
- Filing Date
- 2024-10-18
- Publication Date
- 2026-04-21
AI Technical Summary
In complex geological conditions with well-developed faults and severe stratigraphic folds, the traditional method of dividing strata along the vertical direction cannot accurately determine the top and bottom positions of the strata, resulting in suboptimal morphology of the strata models and limiting their application effectiveness.
By acquiring stratigraphic block units, calculating the central axis and depositional direction of the block, and using the projection of the inflection point onto the central axis as anchor points, sub-layers are divided along the depositional direction to establish sub-layer segments within the stratigraphic block unit. Sub-layer segments of the same stratigraphic unit are then spliced together to form a complete sub-layer.
Under complex geological conditions, local distortions were avoided, a good layer morphology was established, and the velocity model was established and iteratively updated during the drilling guidance process.
Smart Images
Figure CN121899892A_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of seismic signal processing and interpretation technology, and more specifically, to a method, apparatus, and equipment for establishing a small-layer model on a stratigraphic profile. Background Technology
[0002] In the field of geophysical exploration, the construction of layer models of geological bodies is crucial for accurately understanding subsurface geological structures. Related techniques primarily construct layer models by dividing the geological body into smaller layers along a vertical line. This method performs well in continuous, gently sloping strata, but when faced with complex geological conditions such as well-developed faults and severe folding, the difference between the true thickness of the strata and the vertical thickness becomes significant, resulting in less than ideal layer model morphology. In extreme cases, it may even be impossible to determine the top and bottom positions of the strata vertically, thus limiting its application effectiveness. Summary of the Invention
[0003] This application provides a method, apparatus, and equipment for establishing a sub-layer model on a stratigraphic profile, aiming to solve the problem that the sub-layer model obtained by dividing sub-layers in the vertical direction in related technologies has an unsatisfactory morphology in geological conditions with severe fault development and folding.
[0004] The first aspect of this application provides a method for establishing a sub-layer model on a stratigraphic profile, including: Obtain stratigraphic block units cut out by strata and faults; Based on the strike of the top and bottom layers of the stratigraphic block unit, the block center axis of the stratigraphic block unit is calculated, and the deposition direction of the stratigraphic block unit is simulated; the top and bottom layers include: the top layer and the bottom layer; Using the projection positions of each inflection point on the boundary of the stratigraphic block unit onto the central axis of the block as anchor points, subdivisions are performed along the depositional direction to obtain sublayer segments within the stratigraphic block unit. By splicing together the sub-layer segments of the stratigraphic block unit of the same stratum, a complete sub-layer in the stratum is established.
[0005] Optionally, the step of calculating the block center axis of the stratigraphic block unit based on the strike of the top and bottom layers of the stratigraphic block unit, and simulating the depositional direction of the stratigraphic block unit, includes: The boundaries of each of the stratigraphic block units are scanned to determine the top and bottom layers of the strata to which each of the stratigraphic block units belongs; Based on the position of the top stratum corresponding to each stratigraphic block unit in the overall top stratum of its respective stratum, the order of the stratigraphic block unit in the entire stratum is determined, and the stratigraphic block units in the same stratum are numbered to obtain the numbering order of each stratigraphic block unit. For each of the stratigraphic blocks, calculate the block's central axis and the depositional direction.
[0006] Optionally, calculating the central axis of the block for each of the formation blocks includes: The start and end positions of the top layer corresponding to the stratigraphic block unit are set as the first point and the second point, respectively; a first straight line is determined based on the first point and the second point, and the first midpoint between the first point and the second point is calculated; The start and end positions of the bottom stratum corresponding to the stratigraphic block unit are set as the third and fourth points, respectively; a second straight line is determined based on the third and fourth points, and the second midpoint between the third and fourth points is calculated; Calculate the third midpoint between the first midpoint and the second midpoint, and take the third midpoint as the center of the block; With the center of the block as the endpoint, draw a first perpendicular line to the first straight line to obtain the first foot of the perpendicular; With the center of the block as the endpoint, draw a second perpendicular line to the second straight line to obtain the second foot of the perpendicular; A first vector is defined with the center of the block as the starting point and the first perpendicular foot as the ending point; a second vector is defined with the center of the block as the starting point and the second perpendicular foot as the ending point. Rotate the first vector clockwise by 90° to obtain the third vector, and rotate the second vector counterclockwise by 90° to obtain the fourth vector; The third vector and the fourth vector are added together to obtain the fifth vector, and the direction of the fifth vector is taken as the direction of the block center axis of the stratigraphic block unit. Starting from the center of the block, and taking the direction of the central axis of the block as the direction, a central ray of the block is defined, and the straight line containing the central ray of the block is taken as the central axis of the block unit of the stratum.
[0007] Optionally, calculating the depositional direction for each of the stratigraphic block units includes: Calculate the projection of the first midpoint onto the central axis of the block to obtain the fifth point; The direction from the first midpoint to the fifth point is taken as the deposition direction of the stratigraphic block unit.
[0008] Optionally, the step of dividing the stratigraphic block unit into smaller layers along the depositional direction using the projection positions of each inflection point on the boundary of the stratigraphic block unit onto the central axis of the block as anchor points, to obtain smaller layer segments within the stratigraphic block unit, includes: For each of the stratigraphic block units, determine all inflection points on the boundary of the stratigraphic block unit, calculate the projection position of each inflection point on the block's central axis of the stratigraphic block unit, and sort the projection positions to obtain a sorted sequence. A boundary sampling ray sequence is defined with each projection position in the sequence as the passing point and the deposition direction as the direction; Calculate the set of intersection points between the sampling line containing each sampling ray in the boundary sampling ray sequence and the boundary of the stratigraphic block unit, and sort all the intersection points in the set of intersection points in the sampling ray direction to determine the target top intersection point and target bottom intersection point corresponding to the beginning and end of each sampling line; Based on the target top intersection point and the target bottom intersection point, the sub-layer division range of each sampling line is obtained; Divide the area of each sampling line into smaller layers to obtain the smaller layer segments within the stratigraphic block unit.
[0009] Optionally, determining the target top intersection point and target bottom intersection point corresponding to the beginning and end of each sampling line includes: If the top and bottom intersection points of the sampling line are found at the beginning and end, and the top and bottom intersection points are located at the top and bottom layers of the stratigraphic block unit, the top and bottom intersection points are respectively determined as the target top and bottom intersection points of the sampling line. If the top intersection point cannot be found on the sampling line or the top intersection point is not on the top stratum of the stratigraphic block unit, a first ray is defined with the opposite direction of the deposition direction as the direction and the corresponding projection position as the endpoint; the first intersection point of the first ray and the outer enclosing rectangle of the stratigraphic block unit is calculated, and the first intersection point is taken as the target top intersection point of the sampling line. If the sampling line cannot find the bottom intersection point or the bottom intersection point is not on the bottom stratum of the stratigraphic block unit, a second ray is defined with the direction of the deposition direction as the direction and the corresponding projection position as the endpoint; the second intersection point of the second ray and the outer enclosing rectangle of the stratigraphic block unit is calculated, and the second intersection point is taken as the target bottom intersection point of the sampling line.
[0010] Optionally, the step of dividing the stratigraphic block unit into smaller layers within the sub-layer division range of each sampling line to obtain sub-layer segments includes: The sedimentary model of the strata to which the stratigraphic block unit belongs is set, and the sedimentary model includes: parallel top model, parallel bottom model or equal scale model; Based on the depositional pattern, sub-layers are divided within the sub-layer division range of each sampling line to obtain sub-layer segments within the stratigraphic block unit.
[0011] Optionally, the step of splicing together the sublayer segments of the stratigraphic block unit within the same stratum to establish a complete sublayer in the stratum includes: According to the numbering order, the sub-layer segments of the stratigraphic block unit of the same stratum are spliced together to establish a complete sub-layer in the stratum.
[0012] The second aspect of this application provides a device for establishing a small-layer model on a river stratigraphic profile, comprising: The acquisition module is used to acquire stratigraphic block units cut out by strata and faults; The calculation module is used to simulate the deposition direction of the block based on the orientation of the top and bottom layers, and to calculate the central axis of the block. The partitioning module is used to divide the stratigraphic block into smaller layers along the depositional direction, using the projection of each inflection point on the block boundary onto the central axis of the block as anchor points, thus obtaining smaller layer segments within the stratigraphic block unit. A module is created to splice together small layers of the same stratigraphic block to establish a complete small layer within the stratigraphy.
[0013] A third aspect of this application provides an electronic device, including: a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method for establishing a small-layer model on a stratigraphic profile as described in the first aspect.
[0014] The method for establishing a small-layer model on a stratigraphic profile according to the embodiments of this application calculates the central axis of the stratigraphic block unit based on the strike of the top and bottom layers, and simulates the depositional direction of the stratigraphic block unit. Using the projection positions of each inflection point on the boundary of the stratigraphic block unit onto the central axis as anchor points, small-layer divisions are performed along the depositional direction to obtain small-layer segments within the stratigraphic block unit. These small-layer segments of the same stratigraphic block unit are then spliced together to establish a complete small layer within the stratigraphy. By starting from the stratigraphic block unit, dividing it into small layers along the simulated depositional direction, and splicing small layers within the same stratigraphic block, local distortions are avoided. This method can fully consider the morphology and strike of various parts of the overall stratigraphy, especially in cases where there are significant differences between the true thickness and vertical thickness of the strata, establishing a good small-layer morphology. It can be used in complex geological conditions with developed faults and severe stratigraphic folds, providing key technical support for the establishment and iterative updating of attribute models based on small-layer models, particularly velocity models during drilling guidance. Attached Figure Description
[0015] To more clearly illustrate the technical solutions of the embodiments of this application, the drawings used in the description of the embodiments of this application will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0016] Figure 1 This is a flowchart illustrating the steps of a method for establishing a sub-layer model on a stratigraphic profile according to an embodiment of this application. Figure 2 This is a schematic diagram illustrating the determination of the block center in a method for establishing a sub-layer model on a stratigraphic profile according to an embodiment of this application; Figure 3 This is a schematic diagram of the method for establishing a small-layer model on a stratigraphic profile according to an embodiment of this application, which determines the central axis of the block; Figure 4 This is a schematic diagram illustrating the method for determining the depositional direction in establishing a small-layer model on a stratigraphic profile according to an embodiment of this application; Figure 5 This is a schematic diagram of the projection of block boundary points onto the central axis of the block in a method for establishing a small-layer model on a stratigraphic profile according to an embodiment of this application. Figure 6 This is a schematic diagram of the method for determining the sampling rays in a method for establishing a sub-layer model on a stratigraphic profile according to an embodiment of this application; Figure 7 This is a schematic diagram of the effect of a method for establishing a sub-layer model on a stratigraphic profile according to an embodiment of this application; Figure 8 This is a schematic diagram of an electronic device according to this application. Detailed Implementation
[0017] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.
[0018] In the accompanying drawings, the size of constituent elements, the thickness of layers, or areas may sometimes be exaggerated for clarity. Therefore, any implementation of this disclosure is not necessarily limited to the dimensions shown in the drawings, and the shapes and sizes of the components in the drawings do not reflect true proportions. Furthermore, the drawings schematically illustrate ideal examples, and any implementation of this disclosure is not limited to the shapes or values shown in the drawings.
[0019] In the field of geophysical exploration, due to the complexity of geological conditions, geological structural models based on stratigraphic blocks derived from strata and faults often remain coarse. It is necessary to further subdivide stratigraphic blocks according to specific sedimentary conditions, establishing layer-by-layer models of geological bodies. Based on this, refined geological attribute models can be established using layers as units, taking into account the strike of the layer and its surrounding topology, and utilizing various geophysical and geological data, such as well data, as data sources. This facilitates the understanding of subsurface geological structures and the guidance of drilling operations. Therefore, the technology for constructing layer-by-layer models of geological bodies has significant research value in the fields of refined structural interpretation and integrated geoengineering.
[0020] In related technologies, the establishment of stratigraphic models is often based on individual strata. The entire geological profile is divided vertically into several stratigraphic units from top to bottom. Then, starting from the top of each traverse of the model or according to a given sampling step size, sampling points are determined, and the stratigraphic subdivision is performed along the vertical line from these sampling points. This method achieves good results when the strata are continuous and gently sloping. However, with the development of refined exploration and interpretation techniques, in geological conditions with well-developed faults and severe stratigraphic folds, the difference between the true thickness and vertical thickness of the strata is significant. The traditional method of subdividing strata vertically yields sub-strata models with less than ideal morphology. In extreme cases, it may even be impossible to determine the top and bottom positions of the strata vertically, rendering the method ineffective.
[0021] This application proposes a method for establishing layer models on stratigraphic profiles. It introduces a novel approach that starts from stratigraphic block units, simulates the depositional direction of the blocks, and then divides the blocks into layers along the simulated depositional direction. This method abandons the traditional approach of dividing layers vertically using strata as control units. Instead, it starts from stratigraphic blocks cut by strata and faults, and based on each stratigraphic block unit, simulates the depositional direction of the block according to the strike of its top and bottom layers. Layer division is then performed along the depositional direction using the projection of each inflection point on the block's central axis as anchor points. This method can establish layer models even in geological conditions with well-developed faults and severe stratigraphic folding, and effectively overcomes the limitations of conventional layer model establishment techniques when there is a significant difference between vertical thickness and true thickness. By distributing the task of layer division across individual stratigraphic block units, it avoids the local distortions that occur in traditional methods of establishing layer models on a single, unified basis.
[0022] Reference Figure 1 , Figure 1 This is a flowchart illustrating the steps of a method for establishing a sub-layer model on a stratigraphic profile according to an embodiment of this application. Figure 1 As shown, the method includes the following steps S11 to S14: Step S11: Obtain the stratigraphic block units cut by the strata and faults.
[0023] For complex geological formations characterized by well-developed faults and severe stratigraphic folds, data on subsurface geological structure and rock properties are obtained through seismic exploration, drilling, and geophysical exploration. Based on this data, stratigraphic interpretation and block unit division can be performed using geological modeling software, extracting stratigraphic blocks formed by strata and fault cutting. These stratigraphic blocks are three-dimensional rock masses with well-defined boundaries and geological properties.
[0024] Step S12: Calculate the block center axis of the stratigraphic block unit based on the strike of the top and bottom layers of the stratigraphic block unit, and simulate the deposition direction of the stratigraphic block unit; the top and bottom layers include: the top layer and the bottom layer.
[0025] The boundaries of each stratigraphic block unit are identified, and the strike data of the top and bottom layers of the stratigraphic block unit are extracted. The top and bottom layers include the top and bottom layers. Based on the strike of the top and bottom layers of the stratigraphic block unit, the block's central axis is calculated, located between the boundaries of the top and bottom layers. Based on the block's central axis, the depositional orientation of the stratigraphic block unit is simulated. Using the morphological information of the top and bottom layers of the stratigraphic block as a starting point, the depositional orientation of the stratigraphic block is simulated, providing guidance for subsequent sublayer division. This addresses the irrationality of sublayer division based solely on vertical thickness when there is a significant difference between vertical thickness and true stratigraphic thickness.
[0026] Further, the step of calculating the block center axis of the stratigraphic block unit based on the strike of the top and bottom layers of the stratigraphic block unit, and simulating the depositional direction of the stratigraphic block unit, includes: The boundaries of each of the stratigraphic block units are scanned to determine the top and bottom layers of the strata to which each of the stratigraphic block units belongs; Based on the position of the top stratum corresponding to each stratigraphic block unit in the overall top stratum of its respective stratum, the order of the stratigraphic block unit in the entire stratum is determined, and the stratigraphic block units in the same stratum are numbered to obtain the numbering order of each stratigraphic block unit. For each of the stratigraphic blocks, calculate the block's central axis and the depositional direction.
[0027] Using geological modeling software or image processing techniques, the boundaries of stratigraphic block units are precisely scanned to identify each unit. By comparing the lithology, thickness, and occurrence of each stratum, the portion of each stratigraphic block unit belonging to its overlying stratum can be further identified. Then, based on the position of the overlying stratum of a stratigraphic block unit within the overall top stratum, its order within the entire stratigraphic group is determined. Following this left-to-right order on the geological profile, the stratigraphic block units within the same stratum are numbered, resulting in their sequential numbering. For each stratigraphic block unit, its depositional orientation and central axis are calculated.
[0028] Further, for each of the formation block units, calculating the block's central axis includes: The start and end positions of the top layer corresponding to the stratigraphic block unit are set as the first point and the second point, respectively; a first straight line is determined based on the first point and the second point, and the first midpoint between the first point and the second point is calculated; The start and end positions of the bottom stratum corresponding to the stratigraphic block unit are set as the third and fourth points, respectively; a second straight line is determined based on the third and fourth points, and the second midpoint between the third and fourth points is calculated; Calculate the third midpoint between the first midpoint and the second midpoint, and take the third midpoint as the center of the block; With the center of the block as the endpoint, draw a first perpendicular line to the first straight line to obtain the first foot of the perpendicular; With the center of the block as the endpoint, draw a second perpendicular line to the second straight line to obtain the second foot of the perpendicular; A first vector is defined with the center of the block as the starting point and the first perpendicular foot as the ending point; a second vector is defined with the center of the block as the starting point and the second perpendicular foot as the ending point. Rotate the first vector clockwise by 90° to obtain the third vector, and rotate the second vector counterclockwise by 90° to obtain the fourth vector; The third vector and the fourth vector are added together to obtain the fifth vector, and the direction of the fifth vector is taken as the direction of the block center axis of the stratigraphic block unit. Starting from the center of the block, and taking the direction of the central axis of the block as the direction, a central ray of the block is defined, and the straight line containing the central ray of the block is taken as the central axis of the block unit of the stratum.
[0029] Reference Figure 2 The specific method for determining the center of the block is as follows: The starting and ending points of the top stratum to which this stratigraphic block unit belongs are designated as the first point (P_TopHorizon_Start) and the second point (P_TopHorizon_End), respectively. Based on the first and second points, a first straight line (L_TopHorizon) is determined, and the midpoint between the starting and ending points of the top stratum is calculated and designated as the first midpoint (P_TopHorizon_Center). The formula for calculating the first midpoint (P_TopHorizon_Center) is as follows: P_TopHorizon_Center = (P_TopHorizon_Start +P_TopHorizon_End) / 2.
[0030] The starting and ending points of the bottom stratum to which this stratigraphic block unit belongs are designated as the third point (P_BottomHorizon_Start) and the fourth point (P_BottomHorizon_End), respectively. Based on the third and fourth points, a second straight line (L_BottomHorizon) is determined, and the midpoint between the starting and ending points of the bottom stratum is calculated and designated as the second midpoint (P_BottomHorizon_Center). The formula for calculating the second midpoint (P_TopHorizon_Center) is as follows: P_BottomHorizon_Center=(P_BottomHorizon_Start +P_BottomHorizon_End) / 2.
[0031] Calculate the midpoint between the first and second midpoints to obtain the third midpoint. This third midpoint is the block center (P_Block_Center). The formula for calculating the block center (P_Block_Center) is: P_Block_Center = (P_TopHorizon_Center +P_BottomHorizon_Center) / 2.
[0032] Reference Figure 3 The specific method for determining the central axis of the block is as follows: With the block center (P_Block_Center) as the endpoint, draw the first perpendicular line to the first straight line (L_TopHorizon) to obtain the first foot of the perpendicular (P_TopHorizon_FP).
[0033] With the block center (P_Block_Center) as the endpoint, draw a second perpendicular line to the second straight line (L_BottomHorizon) to obtain the second foot of the perpendicular (P_BottomHorizon_FP).
[0034] Starting from the center of the block (P_Block_Center) and ending at the first perpendicular (P_TopHorizon_FP), define the first vector (V_TopHorizon). Starting from the center of the block (P_Block_Center) and ending at the second perpendicular (P_BottomHorizon_FP), define the second vector (V_BottomHorizon).
[0035] Rotate the first vector (V_TopHorizon) 90° clockwise to obtain the third vector (V_TopHorizon_Rotated), and rotate the second vector (V_BottomHorizon) 90° counterclockwise to obtain the fourth vector (V_BottomHorizon_Rotated).
[0036] Add the third vector (V_TopHorizon_Rotated) and the fourth vector (V_BottomHorizon_Rotated) to obtain the fifth vector. The direction of this fifth vector is taken as the direction of the block's central axis (V_BlockAxis_Direction) for this stratigraphic block unit. The formula for obtaining the direction of the block's central axis (V_BlockAxis_Direction) is: V_BlockAxis_Direction = V_TopHorizon_Rotated + V_BottomHorizon_Rotated.
[0037] Starting from the block center (P_Block_Center), and taking the direction of the block center axis (V_BlockAxis_Direction) as the direction, define the block center ray (R_BlockAxis_Direction), and take the straight line containing the block center ray (R_BlockAxis_Direction) as the block center axis (L_BlockAxis_Direction) of this stratigraphic block unit.
[0038] Reference Figure 4 The specific method for determining the deposition direction is as follows: Calculate the projection of the first midpoint (P_TopHorizon_Center) onto the block's central axis (L_BlockAxis_Direction) to obtain the fifth point (P_Pro_TopHorizon_Center). Define a vector using the fifth point (P_Pro_TopHorizon_Center) and the first midpoint (P_TopHorizon_Center) to obtain the depositional direction (Simulated_Deposit_Direction). The direction from the first midpoint to the fifth point is taken as the depositional direction (Simulated_Deposit_Direction) of this stratigraphic block unit. The formula for obtaining the depositional direction (Simulated_Deposit_Direction) is: Simulated_Deposit_Direction = P_Pro_TopHorizon_Center - P_TopHorizon_Center.
[0039] Step S13: Using the projection positions of each inflection point on the boundary of the stratigraphic block unit onto the central axis of the block as anchor points, subdivide the stratigraphic block unit into smaller layers along the depositional direction to obtain smaller layer segments within the stratigraphic block unit.
[0040] The inflection points on the boundary of each stratigraphic block are projected onto the calculated central axis of the block; these projection points serve as anchor points. Based on the projection locations and depositional directions, the stratigraphic block is subdivided into a series of sublayer segments along the depositional direction. This method facilitates the establishment of sublayer models in geological conditions with well-developed faults and severe stratigraphic folding, and effectively overcomes the limitations of conventional sublayer model establishment techniques when there is a significant difference between vertical thickness and true thickness.
[0041] Further, the step of dividing the stratigraphic block unit into smaller layers along the depositional direction using the projection positions of each inflection point on the boundary of the block onto the central axis of the block as anchor points, to obtain smaller layer segments within the stratigraphic block unit, includes: For each of the stratigraphic block units, determine all inflection points on the boundary of the stratigraphic block unit, calculate the projection position of each inflection point on the block's central axis of the stratigraphic block unit, and sort the projection positions to obtain a sorted sequence. A boundary sampling ray sequence is defined with each projection position in the sequence as the passing point and the deposition direction as the direction; Calculate the set of intersection points between the sampling line containing each sampling ray in the boundary sampling ray sequence and the boundary of the stratigraphic block unit, and sort all the intersection points in the set of intersection points in the sampling ray direction to determine the target top intersection point and target bottom intersection point corresponding to the beginning and end of each sampling line; Based on the target top intersection point and the target bottom intersection point, the sub-layer division range of each sampling line is obtained; Divide the area of each sampling line into smaller layers to obtain the smaller layer segments within the stratigraphic block unit.
[0042] Reference Figure 5 For each stratigraphic block unit, for all inflection points on the boundary, denoted as {P1, P2, P3, ..., Pn}, calculate their projection positions on the block's central axis (L_BlockAxis_Direction), and sort these positions based on their distances from the endpoints of the ray on the block's central ray (R_BlockAxis_Direction) to obtain the sorted sequence {O1, O2, O3, ..., On}.
[0043] Reference Figure 6 Based on the elements (projected positions) in the sequence {O1,O2,O3,,,On} as the passing points, and using the simulated deposition direction (Simulated_Deposit_Direction) as the direction, a sampling ray sequence {SampRay1,SampRay2,SampRay3,,,SampRayn} is defined, resulting in a boundary sampling ray sequence. The set of intersection points between the sampling line containing each sampling ray in the boundary sampling ray sequence and the boundary of the stratigraphic block unit is calculated. All intersection points in the set are sorted along the ray direction, and the first and last intersection points are found, thus determining the target top intersection point and target bottom intersection point (MaxDis_TopTerminal, MaxDis_BottomTerminal) corresponding to the first and last points of the sampling line. This yields the sublayer division range (MLT_Terminal, MLB_Terminal) of the sampling line. Furthermore, the sublayer division range on each sampling line is obtained, and the corresponding sublayer division range sequence is as follows: {(MLT_Terminal1, MLB_Terminal1), (MLT_Terminal2, MLB_Terminal2), (MLT_Terminal3, MLB_Terminal3), ... (MLT_Terminaln, MLB_Terminaln)}.
[0044] Further, determining the target top intersection point and target bottom intersection point corresponding to the beginning and end of each sampling line includes: If the top and bottom intersection points of the sampling line are found at the beginning and end, and the top and bottom intersection points are located at the top and bottom layers of the stratigraphic block unit, the top and bottom intersection points are respectively determined as the target top and bottom intersection points of the sampling line. If the top intersection point cannot be found on the sampling line or the top intersection point is not on the top stratum of the stratigraphic block unit, a first ray is defined with the opposite direction of the deposition direction as the direction and the corresponding projection position as the endpoint; the first intersection point of the first ray and the outer enclosing rectangle of the stratigraphic block unit is calculated, and the first intersection point is taken as the target top intersection point of the sampling line. If the sampling line cannot find the bottom intersection point or the bottom intersection point is not on the bottom stratum of the stratigraphic block unit, a second ray is defined with the direction of the deposition direction as the direction and the corresponding projection position as the endpoint; the second intersection point of the second ray and the outer enclosing rectangle of the stratigraphic block unit is calculated, and the second intersection point is taken as the target bottom intersection point of the sampling line.
[0045] Calculate the set of intersection points between the sampling line containing each sampling ray in the boundary sampling ray sequence and the boundary of the stratigraphic block unit. Sort all intersection points in the set along the ray direction, and find the first and last intersection points, which are designated as the top intersection point and the bottom intersection point. If the top and bottom intersection points of the sampling line are found, and these two points are located at the top and bottom layers of their respective stratigraphic block units, they are designated as the target top intersection point and target bottom intersection point (MaxDis_TopTerminal, MaxDis_BottomTerminal), thus obtaining the sublayer division range (MLT_Terminal, MLB_Terminal) of the sampling line.
[0046] If the sampling line cannot find a top intersection point or is not located at the top layer of its stratigraphic block unit, i.e., the target top intersection point MaxDis_TopTerminal cannot be determined, then the first ray is defined with the opposite direction of the depositional direction as the direction and the corresponding projection position (e.g., O1) of the sampling line as the endpoint. The first intersection point of this first ray with the rectangle surrounding the current stratigraphic block unit is calculated, and this first intersection point is taken as the target top intersection point of the sampling line. In this way, the MLT_Terminal of the sublayer division range of the sampling line can be obtained.
[0047] If the sampling line cannot find a bottom intersection point or is not located at the bottom layer of its stratigraphic block unit, i.e., the target bottom intersection point MaxDis_BottomTerminal cannot be determined, then a second ray is defined with the depositional direction as the direction and the projection position (e.g., O1) corresponding to the sampling line as the endpoint. The second intersection point of this second ray with the rectangle enclosing the current stratigraphic block unit is calculated, and this second intersection point is taken as the target bottom intersection point of the sampling line. In this way, the MLB_Terminal of the sublayer division range of the sampling line can be obtained.
[0048] Further, the step of dividing the stratigraphic block unit into smaller layers within the sub-layer division range of each sampling line to obtain sub-layer segments includes: The sedimentary model of the strata to which the stratigraphic block unit belongs is set, and the sedimentary model includes: parallel top model, parallel bottom model or equal scale model; Based on the depositional pattern, sub-layers are divided within the sub-layer division range of each sampling line to obtain sub-layer segments within the stratigraphic block unit.
[0049] Based on the sedimentary mode (parallel top, parallel bottom, equal scale) set by the user, sublayers are divided within the sublayer division range along each sampling line. For the parallel top-bottom mode, starting from the MicroLayer_TopTerminal (parallel top) or MicroLayer_BottomTerminal (parallel bottom), the position of the bottom (parallel top) or top (parallel bottom) of each sublayer is determined according to the set maximum sublayer thickness. If this position exceeds the MicroLayer_BottomTerminal (parallel top) or MicroLayer_TopTerminal (parallel bottom), then the MicroLayer_BottomTerminal (parallel top) or MicroLayer_TopTerminal (parallel bottom) is selected as the bottom (parallel top) or top (parallel bottom) position of the sublayer. If the stratum is in equal scale mode, then the maximum distance between the sublayer division ranges across all stratigraphic blocks within the same stratum is counted as the stratum thickness. Then, based on the set maximum sublayer thickness, the number of sublayers to be divided within that maximum thickness is determined. The sub-layers within all stratigraphic blocks of this formation are divided according to this number. This method ensures that the thickness of any sub-layer in this formation will not exceed the set maximum thickness.
[0050] Once all sub-layers within a stratigraphic block unit's sub-layer division sequence have been completed, the top positions of belonging to the same sub-layer within each sub-layer division are connected sequentially to obtain the top of a sub-layer; similarly, the bottom positions are connected sequentially to obtain the bottom of a sub-layer. Determining the top and bottom of a sub-layer allows for its definition, ultimately yielding the sub-layer fragments within that stratigraphic block unit.
[0051] Step S14: Splice together the sub-layer segments of the stratigraphic block unit in the same stratum to establish a complete sub-layer in the stratum.
[0052] By following the numbering order of stratigraphic block units within the same stratum, smaller layer fragments of the same stratigraphic block unit can be spliced together to establish a complete sub-layer within the stratum. By splicing together smaller layer fragments from different blocks, the established sub-layers within each stratigraphic block can be combined to form a complete sub-layer. This helps to address the problem of unsatisfactory sub-layer modeling in geological conditions where the overall strata are interrupted by faults or severely folded during refined geophysical interpretation.
[0053] Furthermore, the step of splicing together the sublayer segments of the stratigraphic block unit within the same stratum to establish a complete sublayer in the stratum includes: According to the numbering order, the sub-layer segments of the stratigraphic block unit of the same stratum are spliced together to establish a complete sub-layer in the stratum.
[0054] Based on the position of the top stratum to which a stratigraphic block unit belongs within the overall top stratum, the order of that stratigraphic block unit within the entire stratum is determined. Then, according to its left-to-right order on the geological profile, the stratigraphic block units within the same stratum are numbered, resulting in the numbering order of each stratigraphic block unit. Following the numbering order of the stratigraphic block units within the same stratum, the chronological order of the smaller layer fragments scattered within each block unit can be obtained, thus constructing the distribution of the entire layer, such as... Figure 7 As shown.
[0055] The method for establishing a small-layer model on a stratigraphic profile according to the embodiments of this application calculates the central axis of the stratigraphic block unit based on the strike of the top and bottom layers, and simulates the depositional direction of the stratigraphic block unit. Using the projection positions of each inflection point on the boundary of the stratigraphic block unit onto the central axis as anchor points, small-layer divisions are performed along the depositional direction to obtain small-layer segments within the stratigraphic block unit. These small-layer segments of the same stratigraphic block unit are then spliced together to establish a complete small layer within the stratigraphy. By starting from the stratigraphic block unit, dividing it into small layers along the simulated depositional direction, and splicing small layers within the same stratigraphic block, local distortions are avoided. This method can fully consider the morphology and strike of various parts of the overall stratigraphy, especially in cases where there are significant differences between the true thickness and vertical thickness of the strata, establishing a good small-layer morphology. It can be used in complex geological conditions with developed faults and severe stratigraphic folds, providing key technical support for the establishment and iterative updating of attribute models based on small-layer models, particularly velocity models during drilling guidance.
[0056] It should be noted that, for the sake of simplicity, the method embodiments are all described as a series of actions. However, those skilled in the art should understand that the embodiments of this disclosure are not limited to the described order of actions, because according to the embodiments of this disclosure, some steps can be performed in other orders or simultaneously. Secondly, those skilled in the art should also understand that the embodiments described in the specification are all preferred embodiments, and the actions involved are not necessarily required by the embodiments of this disclosure.
[0057] Based on the same inventive concept, embodiments of this application disclose an apparatus for establishing a small-layer model on a stratigraphic profile, comprising: The acquisition module is used to acquire stratigraphic block units cut out by strata and faults; The calculation module is used to simulate the deposition direction of the block based on the orientation of the top and bottom layers, and to calculate the central axis of the block. The partitioning module is used to divide the stratigraphic block into smaller layers along the depositional direction, using the projection of each inflection point on the block boundary onto the central axis of the block as anchor points, thus obtaining smaller layer segments within the stratigraphic block unit. A module is created to splice together small layers of the same stratigraphic block to establish a complete small layer within the stratigraphy.
[0058] Optionally, the calculation module is specifically used for: The boundaries of each of the stratigraphic block units are scanned to determine the top and bottom layers of the strata to which each of the stratigraphic block units belongs; Based on the position of the top stratum corresponding to each stratigraphic block unit in the overall top stratum of its respective stratum, the order of the stratigraphic block unit in the entire stratum is determined, and the stratigraphic block units in the same stratum are numbered to obtain the numbering order of each stratigraphic block unit. For each of the stratigraphic blocks, calculate the block's central axis and the depositional direction.
[0059] Optionally, the calculation module is specifically used for: The start and end positions of the top layer corresponding to the stratigraphic block unit are set as the first point and the second point, respectively; a first straight line is determined based on the first point and the second point, and the first midpoint between the first point and the second point is calculated; The start and end positions of the bottom stratum corresponding to the stratigraphic block unit are set as the third and fourth points, respectively; a second straight line is determined based on the third and fourth points, and the second midpoint between the third and fourth points is calculated; Calculate the third midpoint between the first midpoint and the second midpoint, and take the third midpoint as the center of the block; With the center of the block as the endpoint, draw a first perpendicular line to the first straight line to obtain the first foot of the perpendicular; With the center of the block as the endpoint, draw a second perpendicular line to the second straight line to obtain the second foot of the perpendicular; A first vector is defined with the center of the block as the starting point and the first perpendicular foot as the ending point; a second vector is defined with the center of the block as the starting point and the second perpendicular foot as the ending point. Rotate the first vector clockwise by 90° to obtain the third vector, and rotate the second vector counterclockwise by 90° to obtain the fourth vector; The third vector and the fourth vector are added together to obtain the fifth vector, and the direction of the fifth vector is taken as the direction of the block center axis of the stratigraphic block unit. Starting from the center of the block, and taking the direction of the central axis of the block as the direction, a central ray of the block is defined, and the straight line containing the central ray of the block is taken as the central axis of the block unit of the stratum.
[0060] Optionally, the calculation module is specifically used for: Calculate the projection of the first midpoint onto the central axis of the block to obtain the fifth point; The direction from the first midpoint to the fifth point is taken as the deposition direction of the stratigraphic block unit.
[0061] Optionally, the modules are specifically used for: For each of the stratigraphic block units, determine all inflection points on the boundary of the stratigraphic block unit, calculate the projection position of each inflection point on the block's central axis of the stratigraphic block unit, and sort the projection positions to obtain a sorted sequence. A boundary sampling ray sequence is defined with each projection position in the sequence as the passing point and the deposition direction as the direction; Calculate the set of intersection points between the sampling line containing each sampling ray in the boundary sampling ray sequence and the boundary of the stratigraphic block unit, and sort all the intersection points in the set of intersection points in the sampling ray direction to determine the target top intersection point and target bottom intersection point corresponding to the beginning and end of each sampling line; Based on the target top intersection point and the target bottom intersection point, the sub-layer division range of each sampling line is obtained; Divide the area of each sampling line into smaller layers to obtain the smaller layer segments within the stratigraphic block unit.
[0062] Optionally, the modules are specifically used for: If the top and bottom intersection points of the sampling line are found at the beginning and end, and the top and bottom intersection points are located at the top and bottom layers of the stratigraphic block unit, the top and bottom intersection points are respectively determined as the target top and bottom intersection points of the sampling line. If the top intersection point cannot be found on the sampling line or the top intersection point is not on the top stratum of the stratigraphic block unit, a first ray is defined with the opposite direction of the deposition direction as the direction and the corresponding projection position as the endpoint; the first intersection point of the first ray and the outer enclosing rectangle of the stratigraphic block unit is calculated, and the first intersection point is taken as the target top intersection point of the sampling line. If the sampling line cannot find the bottom intersection point or the bottom intersection point is not on the bottom stratum of the stratigraphic block unit, a second ray is defined with the direction of the deposition direction as the direction and the corresponding projection position as the endpoint; the second intersection point of the second ray and the outer enclosing rectangle of the stratigraphic block unit is calculated, and the second intersection point is taken as the target bottom intersection point of the sampling line.
[0063] Optionally, the modules are specifically used for: The sedimentary model of the strata to which the stratigraphic block unit belongs is set, and the sedimentary model includes: parallel top model, parallel bottom model or equal scale model; Based on the depositional pattern, sub-layers are divided within the sub-layer division range of each sampling line to obtain sub-layer segments within the stratigraphic block unit.
[0064] Optionally, the module is specifically used for: According to the numbering order, the sub-layer segments of the stratigraphic block unit of the same stratum are spliced together to establish a complete sub-layer in the stratum.
[0065] Based on the same inventive concept, this application discloses an electronic device, with reference to... Figure 8 , Figure 8 This is a schematic diagram of an electronic device illustrated in an embodiment of this disclosure. For example... Figure 8 As shown, the electronic device 100 includes a memory 110 and a processor 120. The memory 110 and the processor 120 are connected via a bus for communication. The memory 110 stores a computer program that can run on the processor 120 to implement the steps in the method for establishing a small-layer model on a stratigraphic profile disclosed in this embodiment.
[0066] As the device embodiment is basically similar to the method embodiment, the description is relatively simple, and relevant parts can be found in the description of the method embodiment.
[0067] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. The same or similar parts between the various embodiments can be referred to each other.
[0068] This application describes embodiments with reference to flowchart illustrations and / or block diagrams of methods, apparatuses, electronic devices, and computer program products according to embodiments of this application. It should be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing terminal device to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing terminal device, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.
[0069] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing terminal device to operate in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.
[0070] These computer program instructions can also be loaded onto a computer or other programmable data processing terminal equipment, causing a series of operational steps to be performed on the computer or other programmable terminal equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable terminal equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.
[0071] Although preferred embodiments of the present application have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the embodiments of the present application.
[0072] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only 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 terminal device 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 terminal device. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or terminal device that includes said element.
[0073] The above provides a detailed description of the method, apparatus, and equipment for establishing a small-layer model on a stratigraphic profile. Specific examples have been used to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of this application. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of this application. Therefore, the content of this specification should not be construed as a limitation of this application.
Claims
1. A method for establishing a sub-layer model on a stratigraphic profile, characterized in that, include: Obtain stratigraphic block units cut out by strata and faults; Based on the strike of the top and bottom layers of the stratigraphic block unit, the block center axis of the stratigraphic block unit is calculated, and the depositional direction of the stratigraphic block unit is simulated; the top and bottom layers include: the top layer and the bottom layer; Using the projection positions of each inflection point on the boundary of the stratigraphic block unit onto the central axis of the block as anchor points, subdivisions are performed along the depositional direction to obtain sublayer segments within the stratigraphic block unit. By splicing together the sub-layer segments of the stratigraphic block unit of the same stratum, a complete sub-layer in the stratum is established.
2. The method for establishing a sub-layer model on a stratigraphic profile according to claim 1, characterized in that, The step of calculating the block center axis of the stratigraphic block unit based on the strike of its top and bottom layers, and simulating the depositional direction of the stratigraphic block unit, includes: The boundaries of each of the stratigraphic block units are scanned to determine the top and bottom layers of the strata to which each of the stratigraphic block units belongs; Based on the position of the top stratum corresponding to each stratigraphic block unit in the overall top stratum of its respective stratum, the order of the stratigraphic block unit in the entire stratum is determined, and the stratigraphic block units in the same stratum are numbered to obtain the numbering order of each stratigraphic block unit. For each of the stratigraphic blocks, calculate the block's central axis and the depositional direction.
3. The method for establishing a sub-layer model on a stratigraphic profile according to claim 2, characterized in that, The calculation of the central axis of each of the formation blocks includes: The start and end positions of the top layer corresponding to the stratigraphic block unit are set as the first point and the second point, respectively; a first straight line is determined based on the first point and the second point, and the first midpoint between the first point and the second point is calculated; The start and end positions of the bottom stratum corresponding to the stratigraphic block unit are set as the third and fourth points, respectively; a second straight line is determined based on the third and fourth points, and the second midpoint between the third and fourth points is calculated; Calculate the third midpoint between the first midpoint and the second midpoint, and take the third midpoint as the center of the block; With the center of the block as the endpoint, draw a first perpendicular line to the first straight line to obtain the first foot of the perpendicular; With the center of the block as the endpoint, draw a second perpendicular line to the second straight line to obtain the second foot of the perpendicular; A first vector is defined with the center of the block as the starting point and the first perpendicular foot as the ending point; a second vector is defined with the center of the block as the starting point and the second perpendicular foot as the ending point. Rotate the first vector 90° clockwise to obtain the third vector, and rotate the second vector 90° counterclockwise to obtain the fourth vector; The third vector and the fourth vector are added together to obtain the fifth vector, and the direction of the fifth vector is taken as the direction of the block center axis of the stratigraphic block unit. Starting from the center of the block, and taking the direction of the central axis of the block as the direction, a central ray of the block is defined, and the straight line containing the central ray of the block is taken as the central axis of the block unit of the stratum.
4. The method for establishing a sub-layer model on a stratigraphic profile according to claim 3, characterized in that, The calculation of the depositional direction for each of the stratigraphic blocks includes: Calculate the projection of the first midpoint onto the central axis of the block to obtain the fifth point; The direction from the first midpoint to the fifth point is taken as the deposition direction of the stratigraphic block unit.
5. The method for establishing a sub-layer model on a stratigraphic profile according to claim 1, characterized in that, The process involves using the projection positions of each inflection point on the boundary of the stratigraphic block unit onto the central axis of the block as anchor points, and dividing the stratigraphic block unit into smaller layers along the depositional direction to obtain smaller layer segments within the stratigraphic block unit, including: For each of the stratigraphic block units, determine all inflection points on the boundary of the stratigraphic block unit, calculate the projection position of each inflection point on the block's central axis of the stratigraphic block unit, and sort the projection positions to obtain a sorted sequence. A boundary sampling ray sequence is defined with each projection position in the sequence as the passing point and the deposition direction as the direction; Calculate the set of intersection points between the sampling line containing each sampling ray in the boundary sampling ray sequence and the boundary of the stratigraphic block unit, and sort all the intersection points in the set of intersection points in the sampling ray direction to determine the target top intersection point and target bottom intersection point corresponding to the beginning and end of each sampling line; Based on the target top intersection point and the target bottom intersection point, the sub-layer division range of each sampling line is obtained; Divide the area of each sampling line into smaller layers to obtain the smaller layer segments within the stratigraphic block unit.
6. The method for establishing a sub-layer model on a stratigraphic profile according to claim 5, characterized in that, The determination of the target top intersection point and target bottom intersection point corresponding to the beginning and end of each sampling line includes: If the top and bottom intersection points of the sampling line are found at the beginning and end, and the top and bottom intersection points are located at the top and bottom layers of the stratigraphic block unit, the top and bottom intersection points are respectively determined as the target top and bottom intersection points of the sampling line. If the top intersection point cannot be found on the sampling line or the top intersection point is not on the top stratum of the stratigraphic block unit, a first ray is defined with the opposite direction of the deposition direction as the direction and the corresponding projection position as the endpoint; the first intersection point of the first ray and the outer enclosing rectangle of the stratigraphic block unit is calculated, and the first intersection point is taken as the target top intersection point of the sampling line. If the sampling line cannot find the bottom intersection point or the bottom intersection point is not on the bottom stratum of the stratigraphic block unit, a second ray is defined with the direction of the deposition direction as the direction and the corresponding projection position as the endpoint; the second intersection point of the second ray and the outer enclosing rectangle of the stratigraphic block unit is calculated, and the second intersection point is taken as the target bottom intersection point of the sampling line.
7. The method for establishing a sub-layer model on a stratigraphic profile according to claim 5, characterized in that, The process of dividing the stratigraphic block unit into smaller layers within the sub-layer division range of each sampling line to obtain sub-layer segments includes: The sedimentary model of the strata to which the stratigraphic block unit belongs is set, and the sedimentary model includes: parallel top model, parallel bottom model or equal scale model; Based on the depositional pattern, sub-layers are divided within the sub-layer division range of each sampling line to obtain sub-layer segments within the stratigraphic block unit.
8. The method for establishing a sub-layer model on a stratigraphic profile according to claim 2, characterized in that, The process of splicing together sub-layer segments of the same stratigraphic block unit to establish a complete sub-layer within the stratigraphy includes: According to the numbering order, the sub-layer segments of the stratigraphic block unit of the same stratum are spliced together to establish a complete sub-layer in the stratum.
9. A device for establishing a small-layer model on a stratigraphic profile, characterized in that, include: The acquisition module is used to acquire stratigraphic block units cut out by strata and faults; The calculation module is used to simulate the deposition direction of the block based on the orientation of the top and bottom layers, and to calculate the central axis of the block. The partitioning module is used to divide the stratigraphic block into smaller layers along the depositional direction, using the projection of each inflection point on the block boundary onto the block's central axis as anchor points, thus obtaining smaller layer segments within the stratigraphic block unit. A module is created to splice together small layers of the same stratigraphic block to establish a complete small layer within the stratigraphy.
10. An electronic device, characterized in that, include: A processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the program, implements the method for establishing a sub-layer model on a stratigraphic profile as described in any one of claims 1 to 8.