CN114647927A - Method and device for constructing drilling contrast section model constrained by stratum distribution mode - Google Patents

Abstract

The invention discloses a method and a device for constructing a drilling contrast profile model constrained by a stratum distribution mode, wherein the method comprises the following steps: (1) reading a drilling set D and a drilling stratum thickness set T of actually measured drilling holes on the section line; (2) reading any adjacent actual measurement drilling hole combination from the drilling hole set D, analyzing the stratum distribution mode, and forming a stratum distribution mode mark list; (3) interpolating to generate virtual drill holes between the adjacent measured drill hole combinations, and forming a drill hole set PD with the adjacent measured drill hole combinations; (4) preliminarily deducing the thickness of each stratum of each borehole in the borehole set PD according to a stratum pinch-out rule according to a stratum distribution mode mark list M, and storing the thickness of each stratum into a stratum thickness matrix PT; (5) optimizing the thickness of the stratum and constructing all stratum surfaces between the adjacent actually measured drilling hole combinations; (6) a contrast profile is generated from the constructed formation face and the borehole set PD. The invention has higher modeling quality, higher efficiency and higher automation degree.

Description

Method and device for constructing drilling contrast profile model constrained by stratum distribution mode

Technical Field

The invention relates to a geographic information technology, in particular to a method and a device for constructing a drilling contrast profile model constrained by a stratum distribution mode.

Background

The construction of stratigraphic correlation sections based on boreholes is the primary means by which a stratigrapher contrasts the strata. The method mainly comprises the steps of converting attribute data of the drill hole and the stratum into geometric elements for drawing, and further realizing visualization of profile data. The formation is spatially discontinuous or has a non-uniform thickness that tends to increase the complexity of constructing a profile of the formation. Therefore, the correct application of the scientific cognition and inference rules of the stratum distribution mode is the basis and the premise for reasonably constructing the stratum contrast section under the complex condition.

The existing method mostly adopts a simple stratum pinch-out treatment rule, so that the profile modeling quality under complex conditions of continuous pinch-out of multiple stratums, intermittent deletion of multiple stratums and the like is not ideal enough, and conflicts with the stratum development rule. In addition, the generated profile can be further edited and optimized manually by experts, but the efficiency is low, and the quality of the profile modeling is not stable enough.

Disclosure of Invention

The purpose of the invention is as follows: aiming at the problems in the prior art, the invention provides a method and a device for constructing a drilling contrast section model constrained by a stratum distribution mode, which have higher efficiency and better quality.

The technical scheme is as follows: the method for constructing the drilling contrast profile model constrained by the stratum distribution mode comprises the following steps:

(1) reading a drilling set D and a drilling stratum thickness set T of actually measured drilling holes on the section line;

(2) reading any adjacent actual measurement drilling combination from the drilling set D, establishing a missing mark list SL of a missing extraction stratum between the adjacent actual measurement drilling combinations based on the drilling stratum thickness set T, establishing a stratum missing interval set DS according to the SL, and analyzing a stratum distribution mode based on the DS to form a stratum distribution mode mark list M;

(3) interpolating to generate virtual drill holes between the adjacent actually measured drill hole combinations, deducing the bottom surface elevation and the ground elevation of each virtual drill hole, and storing the virtual drill holes between the adjacent actually measured drill hole combinations into a drill hole set PD;

(4) preliminarily deducing the thickness of each stratum of each borehole in the borehole set PD according to a stratum pinch-out rule according to a stratum distribution mode mark list M, and storing the thickness of each stratum into a stratum thickness matrix PT;

(5) optimizing the stratum thickness according to the bottom surface elevation and the ground elevation of each virtual drilling hole, and constructing all stratum surfaces between the adjacent actually-measured drilling hole combinations;

(6) generating a comparison profile between the adjacent actually measured drilling hole combinations according to the constructed ground surface and the drilling hole set PD;

(7) and (4) circularly executing the steps (2) to (6) until all the adjacent measured borehole combinations in the borehole set D are traversed to obtain all the comparison section maps.

Further, the step (1) comprises;

(1-1) reading the drilling information of the actually measured drilling holes on the section line to obtain a drilling hole set D ═ D _i1, · dn } for 1, 2, · dn; wherein d is_iRepresenting the ith borehole, dn representing the number of boreholes;

(1-2) reading the borehole stratum information of the actually measured borehole on the section line to obtain a stratum thickness set T ═ T _i，j1, · dn, | i ═ 1, 2; j ═ 1, 2., sni }; wherein, t_i，jIndicating a borehole d_iSni denotes the borehole d_iThe total number of strata encountered.

Further, the step (2) comprises;

(2-1) reading any group of adjacent measured boreholes D from the borehole set D_iAnd d_i+1；

(2-2) traversing the borehole stratum thickness set T, matching the borehole d according to the following formula_iAnd drilling a hole d_i+1The list SL ═ { SL ═ of missing flags of the inter-missing extraction formation_j|j＝1，2，...，sni}；

Wherein, sl is_jFormation lost flag, t, representing formation j_i，j、t_i+1，jIndicating a borehole d_i、d_i+1Thickness of middle formation j, sni denotes borehole d_iThe total number of strata encountered by drilling;

(2-3) stratum-missing flag list SL_iBuilding a borehole d_iAnd drilling a hole d_i+1Inter-strata missing interval set DS ═ DS_u1, 2., dsn }; wherein ds_uThe triplet represents the u-th formation missing interval, dsn represents the number of missing intervals, and start represents the interval ds_uIs the starting formation number, end is the interval ds_uThe end stratum number of (1), count being the interval ds_uThe number of strata in (1);

(2-4) for each stratigraphic missing interval DS in the stratigraphic missing interval set DS_uOver the interval ds_uAll strata missing in the interval and the interval ds is calculated according to the following formula_uInner stratigraphic distribution pattern mark m_uAnd storing the stratum distribution mode mark list M ═ M_u|u＝1，2，...，dsn}；

Further, the step (2-3) comprises;

(2-3-1) sequentially traversing the missing flag list SL until SL_jThe value of j is given as start and the count value is given as 1, which are not equal to 0;

(2-3-2) continue to traverse the missing flag list SL when j is less than or equal to sni and | SL_j| is equal to | sl_startWhen |, the count value is increased by 1; otherwise, assigning the end value as j-1, and newly building an interval element ds_u(start, end, count), and add to the missing interval set DS, and then return to step (2-3-2); wherein, | | is an absolute value symbol;

and (2-3-3) circularly executing the steps (2-3-1) to (2-3-2) until the missing mark list SL is traversed to obtain a stratum missing interval set DS.

Further, the step (3) comprises;

(3-1) traversing the set DS, and acquiring all the DS in the DS_u(count) as the number of virtual boreholes to be inserted, the maximum value, denoted as the maximum interval missing formation number max, wherein ds_u(count) is the interval ds_uThe number of strata in (1);

(3-2) creating a set of empty boreholes PD ═ { PD ═ 2_k1, 2, max +2}, storing the virtual borehole and the read adjacent measured borehole d_iAnd d_i+1(ii) a Wherein pd is₁Assigning a value to the measured borehole d_i，pd_max+2Assigning a value to the measured borehole d_i+1，pd₂To pd _max+11 st to max virtual drill holes are formed in sequence;

(3-3) sequentially traversing the borehole set PD, and calculating a virtual borehole PD to be inserted by linear interpolation according to the following formula_kCoordinates and bottom elevation pd of_k(bottom)；

Wherein pd is_k(x) And pd_k(y) is a virtual borehole pd_kAbscissa and ordinate of (a), d_i(x) And d_i(y) is a bore hole d_iAbscissa and ordinate of (d)_i+1(x) And d_i+1(y) is a bore hole d_i+1Abscissa and ordinate of (a), d_i(bottom)、d_i+1(bottom) and pd_k(bottom) respectively represents the drilled holes d_i、d_i+1And pd_kThe hole bottom elevation;

(3-4) according to the virtual borehole pd_kCoordinate extraction of virtual borehole pd from DEM data_kGround elevation pd_k(h) Max +1, completes each borehole assignment in the set of boreholes PD.

Further, the step (4) comprises;

(4-1) constructing a formation thickness matrix PT ═ PT for the borehole set PD _k，j1, · 2, · max + 2; j ═ 1, 2., sni }; wherein pt is_k，jRepresents the thickness of the jth formation of the kth borehole in the borehole set PD and assigns an initial value of 0And sni denotes a borehole d_iThe total number of strata encountered by drilling;

(4-2) reading any stratum interval DS from the set DS_uWhen m is_uWhen 2, interval ds_uIf the inner stratum is continuous, calculating the thickness of the stratum according to the following formula, and executing the step (4-7); otherwise, executing the step (4-3);

in the formula, t_i，j、t_i+1，jIndicating a borehole d_i、d_i+1The thickness of the middle formation j;

(4-3) drilling a hole d_iNumber of drilled formations quenched count₁Setting the initial value to 0, drilling a hole d_i+1Number of drilled formations quenched count₂An initial value of 0 is assigned to the formation j and an initial value of ds is assigned to the formation j_u(start)，ds_u(start) is ds_uThe starting formation number of (a);

(4-4) updating the number of holes drilled in the ground which have been extinguished, count, according to the following formula₁And count₂And deducing a virtual borehole subscript PN of the formation pinch-out;

in the formula, sl_jA formation missing flag representing a formation j in the missing flag list SL;

(4-5) estimating the thickness pt of the formation in each of the pseudo drilled holes according to the following equation_k，j；

(4-6) judging whether j is equal to ds_u(end)，ds_u(end) is ds_uIf yes, ending iteration, executing step (4-7), otherwise, setting j to j +1, and executing step (4-4);

and (4-7) circularly executing the steps (4-2) - (4-6) until all stratum intervals in the set DS are traversed, and finishing the calculation of the stratum thickness matrix PT.

Further, the step (5) comprises;

(5-1) optimizing a stratum thickness matrix PT based on the inferred bottom surface elevation and ground elevation of the virtual borehole;

(5-2) calculating the top surface elevation of each stratum of each borehole according to the following formula, and storing a top surface elevation matrix UP ═ { UP ═ UP }_k，j|k＝1，2，...，max+2；j＝1，2，...，sni}；

In the formula, up_k，jRepresenting a borehole PD in a borehole set PD_kTop surface elevation, pd, of formation j_k(h) Representing a borehole pd_kGround elevation of pt_k，jRepresenting a borehole pd in a formation thickness matrix PT_kSni denotes the borehole d_iThe total number of the stratums encountered by drilling, and max +2 represents the number of the drilling holes in the drilling hole set PD;

(5-3) calculating the elevation of the bottom surface of each stratum of each borehole according to the following formula, and storing the bottom surface elevation matrix DOWN ═ DOWN_k，j|k＝1，2，...，max+2；j＝1，2，...，sni}；

down_k，j＝up_k，j-pt_k，j，k＝1，2，...，max+2，j＝1，2，...，sni

In the formula, down_k，jRepresenting a borehole PD in a borehole set PD_kThe elevation of the bottom surface of the stratum j;

(5-4) constructing a borehole d according to the missing mark list SL, the optimized formation thickness matrix PT, the top elevation matrix UP and the bottom elevation matrix DOWN_iAnd d_i+1All the terrain levels in between.

Further, the step (5-1) comprises;

(5-1-1) calculating each borehole PD in the borehole set PD according to the following formula based on the formation thickness matrix PT_kIs pd of_k(depth)；

pd_k(depth)＝∑_jpt_k，j，k＝1，2，...，max+2

(5-1-2) calculating the borehole pd according to the following equation based on the surface elevation and the bottom surface elevation of the borehole_kOptimized bottom surface elevation pd_k(depth) ' and deducing the thickness pt ' of each stratum '_k，jUpdating a stratum thickness matrix PT;

in the formula, pd_k(h) Representing a borehole pd_kGround elevation, pd_k(bottom) denotes the borehole pd_kThe bottom surface elevation.

Further, the step (5-4) comprises;

(5-4-1) setting j ═ 1;

(5-4-2) constructing an empty point set P;

(5-4-3) traversing each element of jth column of the formation thickness matrix PT when PT_k，jWhen the coordinate is greater than 0, the coordinate (k, up) is set_k，j) And (k, down)_k，j) Adding to the point set P;

(5-4-4) when the jth element SL of the flag list SL is missing_jWhen the value is-1, sequentially traversing jth column of the thickness matrix PT, reading row number fk of the first zero value, and combining the coordinate points (k, up)_fk，j) Adding to the point set P; when sl is_jWhen 1, sequentially traversing jth column of the thickness matrix PT, reading the last row index lk of zero value, and dividing the coordinate point (pd)_k(x)，pd_k(y)，up_lk，j) Adding to the point set P;

(5-4-5) if the point set P is not empty, constructing a polygon based on the point set P as a ground-level sp_i，j；

(5-4-6) judging whether j is equal to the number of columns of the stratum thickness matrix PT, if so, ending the iteration, otherwise, changing j to j +1, and returning to execute the step (5-4-2).

The device for constructing the stratum distribution mode constrained borehole contrast profile model comprises a processor and a computer program which is stored on a memory and can run on the processor, wherein the processor realizes the method when executing the program.

Has the advantages that: compared with the prior art, the invention has the following remarkable advantages: the method can be better suitable for more complex stratum distribution modes such as continuous pinch-out and continuous deletion of the stratum, not only improves the modeling quality of the drilling contrast section, but also has higher automation degree and higher efficiency.

Drawings

FIG. 1 is DEM data and borehole data as used in the present embodiment;

FIG. 2 is a grab of the present embodiment using a drill table;

FIG. 3 is a snapshot of the present embodiment using a formation table;

FIG. 4 is a flow chart of a method of modeling a borehole contrast profile constrained by formation distribution patterns provided by the present invention;

FIG. 5 is a schematic representation of formation thickness extrapolation of the present invention;

FIG. 6 is d₃₃And d₃₄A comparison cross-sectional view before thickness optimization;

FIG. 7 is d₃₃And d₃₄Cross-sectional comparison of the thickness optimized.

Detailed Description

As will be described in further detail below, in this embodiment, a fourth system of drilling and DEM data (fig. 1, 2, and 3) in a certain experimental area is selected as experimental data, and a projection coordinate system of the DEM and drilling vector data adopted in the experiment is a WGS84 coordinate system. The following further description is provided by describing a specific embodiment in conjunction with the accompanying drawings.

As shown in fig. 4, the present embodiment provides a method for modeling a contrast profile of a borehole with a formation distribution pattern constraint, which specifically includes the following steps.

(1) And reading a borehole set D and a borehole stratum thickness set T of the actual measurement borehole on the section line.

The method specifically comprises the following steps:

(1-1) reading the drilling information of the actually measured drilling holes on the section line to obtain a drilling hole set D ═ D _i1, · dn } for 1, 2, · dn; wherein d is_iRepresenting the ith borehole, dn representing the number of boreholes; in the present embodiment, dn — 69;

(1-2) reading the borehole stratum information of the actually measured borehole on the section line to obtain a stratum thickness set T ═ T_i，jDh, | i ═ 1, 2, · dh; j ═ 1, 2, ·, sni }; wherein, t_i，jIndicating a borehole d_iSni denotes the borehole d_iThe total number of strata encountered is sni ═ 35 in this embodiment. .

(2) Reading any adjacent actual measurement drilling combination from the drilling set D, establishing a missing mark list SL of a missing extraction stratum between the adjacent actual measurement drilling combinations based on the drilling stratum thickness set T, establishing a stratum missing interval set DS according to the SL, and analyzing a stratum distribution mode based on the DS to form a stratum distribution mode mark list M.

The method specifically comprises the following steps:

(2-1) reading any group of adjacent measured boreholes D from the borehole set D_iAnd d_i+1(ii) a In the present embodiment, i — 33 is taken as an example;

(2-2) traversing the borehole stratum thickness set T, matching the borehole d according to the following formula_iAnd drilling a hole d_i+1The list SL ═ { SL ═ of missing flags of the inter-missing extraction formation_j|j＝1，2，...，sni}；

Wherein, sl is_jFormation lost flag, t, representing formation j_i，j、t_i+1，jIndicating a borehole d_i、d_i+1Thickness of middle formation j, sni denotes borehole d_iThe total number of strata encountered by drilling;

(2-3) stratum-missing flag list SL_iBuilding a borehole d_iAnd drilling a hole d_i+1Inter-strata missing interval set DS ═ DS _u1, 2.., dsn }; wherein ds_uThe triplet represents the u-th formation missing interval, dsn represents the number of missing intervals, and start represents the interval ds_uInitial formation ofNumber end is the interval ds_uThe end stratum number of (1), count being the interval ds_uThe number of strata in (1); in the present embodiment, ds when i is 33₁(1, 1, 0), dsn-9; the method comprises the following specific steps:

(2-3-1) sequentially traversing the missing flag list SL until SL_jThe value of j is given as start and the count value is given as 1, which are not equal to 0;

(2-3-2) continue to traverse the missing flag list SL when j is less than or equal to sni and | SL_j| is equal to | sls_tartWhen |, the count value is increased by 1; otherwise, assigning the end value as j-1, and newly building an interval element ds_u(start, end, count), and add to the missing interval set DS, and then return to step (2-3-2); wherein, | | is an absolute value symbol;

(2-3-3) circularly executing the steps (2-3-1) to (2-3-2) until the missing mark list SL is traversed to obtain a stratum missing interval set DS;

(2-4) for each stratigraphic missing interval DS in the stratigraphic missing interval set DS_uOver the interval ds_uAll formations in the interval are missing and the interval ds is calculated as follows_uInner stratigraphic distribution mode mark m_uAnd storing the stratum distribution pattern mark list M ═ M_u|u＝1，2，...，dsn}；

In the present embodiment, when i is 33, M is {0, 1, 0, 1, 0, 1, 0, 2, 0 }.

(3) And generating virtual drill holes between the adjacent actually measured drill hole combinations through interpolation, deducing the bottom surface elevation and the ground surface elevation of each virtual drill hole, and storing the virtual drill holes between the adjacent actually measured drill hole combinations into a drill hole set PD.

The method comprises the following steps:

(3-1) traversing the set DS, and acquiring all the DS in the DS_u(count) as the number of virtual boreholes to be inserted, the maximum value, denoted as the maximum interval missing formation number max, wherein ds_u(count) is the interval ds_uThe number of strata in (1); in the present embodiment, when i is 33, max is 11;

(3-2) creating a set of empty boreholes PD ═ { PD ═ 2_k1, 2, max +2}, storing the virtual borehole and the read adjacent measured borehole d_iAnd d_i+1(ii) a Wherein pd is₁Assigning a value to the measured borehole d_i，pd_max+2Assigning the value to the measured borehole d_i+1，pd₂To pd _max+11 st to max virtual drill holes are formed in sequence;

(3-3) sequentially traversing the borehole set PD, and calculating a virtual borehole PD to be inserted according to the following formula through linear interpolation_kCoordinates and bottom elevation pd of_k(bottom)；

Wherein pd is_k(x) And pd_k(y) is a virtual borehole pd_kAbscissa and ordinate of (d)_i(x) And d_i(y) is a bore hole d_iAbscissa and ordinate of (a), d_i+1(x) And d_i+1(y) is a bore hole d_i+1Abscissa and ordinate of (d)_i(bottom)、d_i+1(bottom) and pd_k(bottom) respectively represents the drilled holes d_i、d_i+1And pd_kThe hole bottom elevation;

(3-4) according to the virtual borehole pd_kCoordinate extraction of virtual borehole pd from DEM data_kGround elevation pd_k(h) K 2, max +1, completing each borehole assignment in the borehole set PD, wherein the surface elevation PD_k(h) The calculation formula is as follows:

wherein x is₀For DEM origin abscissa, y₀Is the vertical coordinate of the DEM origin, C is the DEM pixel size,

to round the symbol down. In this embodiment, C is 100 meters.

(4) And preliminarily deducing the thickness of each stratum of each borehole in the borehole set PD according to the stratum pinch-out rule according to the stratum distribution mode mark list M, and storing the thickness of each stratum into a stratum thickness matrix PT.

The method specifically comprises the following steps;

(4-1) constructing a formation thickness matrix PT ═ { PT) for the borehole set PD _k，j1, 2,. said, max + 2; j ═ 1, 2, ·, sni }; wherein pt is_k，jThickness of jth formation of kth borehole in borehole set PD, and initial values of 0, sni for borehole d_iThe total number of strata encountered by drilling;

(4-2) reading any stratum interval DS from the set DS_uWhen m is_uWhen 2, interval ds_uIf the inner stratum is continuous, calculating the thickness of the stratum according to the following formula, and executing the step (4-7); otherwise, executing the step (4-3);

in the formula, t_i，j、t_i+1，jIndicating a borehole d_i、d_i+1The thickness of the middle formation j;

(4-3) drilling a hole d_iNumber of drilled formations quenched count₁Setting the initial value to 0, drilling a hole d_i+1Number of drilled formations quenched count₂An initial value of 0 is assigned to the formation j and an initial value of ds is assigned to the formation j_u(start)，ds_u(start) is ds_uThe starting formation number of (a);

(4-4) updating the number of holes drilled in the ground which have been extinguished, count, according to the following formula₁And count₂And deducing a virtual borehole subscript PN of the formation pinch-out;

in the formula, sl_jA stratigraphic missing flag representing a stratigraphic j in the missing flag list SL;

(4-5) estimating the thickness pt of the formation in each of the pseudo drilled holes according to the following equation_k，jAs shown in fig. 5;

(4-6) judging whether j is equal to ds_u(end)，ds_u(end) is ds_uIf yes, ending iteration, executing step (4-7), otherwise, setting j to j +1, and executing step (4-4);

and (4-7) circularly executing the steps (4-2) - (4-6) until all stratum intervals in the set DS are traversed, and completing the calculation of the stratum thickness matrix PT, as shown in FIG. 6.

(5) And optimizing the thickness of the stratum according to the bottom surface elevation and the ground elevation of each virtual drilling hole, and constructing all the stratum surfaces between the adjacent actually-measured drilling hole combinations.

The method specifically comprises the following steps;

(5-1) optimizing a formation thickness matrix PT based on the inferred bottom surface elevation and ground elevation of the virtual borehole; the method comprises the following specific steps:

(5-1-1) calculating each borehole PD in the borehole set PD according to the following formula based on the formation thickness matrix PT_kIs pd of_k(depth)；

pd_k(depth)＝∑_jpt_k，j，k＝1，2，...，max+2

(5-1-2) calculating the borehole pd according to the following formula based on the surface level and the bottom level of the borehole_kBottom surface elevation pd after optimization_k(depth) ' and deducing the thickness pt ' of each stratum '_k，jUpdating a stratum thickness matrix PT;

in the formula, pd_k(h) Representing a borehole pd_kGround elevation, pd_k(bottom) denotes the borehole pd_kThe bottom surface elevation of the plate;

(5-2) calculating the top surface elevation of each stratum of each borehole according to the following formula, and storing a top surface elevation matrix UP ═ { UP ═ UP }_k，j|k＝1，2，...，max+2；j＝1，2，...，sni}；

In the formula, up_k，jRepresenting a borehole PD in a borehole set PD_kTop surface elevation, pd, of formation j_k(h) Representing a borehole pd_kGround elevation of pt_k，jRepresenting a borehole pd in a formation thickness matrix PT_kSni denotes the borehole d_iThe total number of the stratums encountered by drilling, and max +2 represents the number of the drilling holes in the drilling hole set PD;

(5-3) calculating the elevation of the bottom surface of each stratum of each borehole according to the following formula, and storing the bottom surface elevation matrix DOWN ═ DOWN_k，j|k＝1，2，...，max+2；j＝1，2，...，sni}；

down_k，j＝up_k，j-pt_k，j，k＝1，2，...，max+2，j＝1，2，...，sni

In the formula, down_k，jRepresenting a borehole PD in a borehole set PD_kThe elevation of the bottom surface of the stratum j;

(5-4) constructing a borehole d according to the missing mark list SL, the optimized formation thickness matrix PT, the top elevation matrix UP and the bottom elevation matrix DOWN_iAnd d_i+1The method comprises the following specific steps of:

(5-4-1) setting j ═ 1;

(5-4-2) constructing an empty point set P;

(5-4-3) traversing each element of jth column of the formation thickness matrix PT when PT_k，jWhen the coordinate is greater than 0, the coordinate (k, up) is set_k，j) And (k, down)_k，j) Adding to the point set P;

(5-4-4) when the jth element SL of the flag list SL is missing_jWhen the coordinate point (k, up) is determined as-1, the jth column of the thickness matrix PT is sequentially traversed, the row number fk of the first zero value is read, and the coordinate point (k, up) is determined_fk，j) Adding to the point set P;when sl is_jWhen 1, sequentially traversing jth column of the thickness matrix PT, reading the last row index lk of zero value, and dividing the coordinate point (pd)_k(x)，pd_k(y)，up_lk，j) Adding to the point set P;

(5-4-5) if the point set P is not empty, constructing a polygon based on the point set P as a ground-level sp_i，j；

(5-4-6) judging whether j is equal to the number of columns of the stratum thickness matrix PT, if so, ending the iteration, otherwise, changing j to j +1, and returning to execute the step (5-4-2).

(6) And generating a comparison section between the adjacent actually measured drilling hole combinations according to the constructed stratum surface and the drilling hole set PD. In this embodiment, when i is 33, d₃₃And d₃₄The cross section of the comparison between them is shown in FIG. 7.

(7) And (4) circularly executing the steps (2) to (6) until all the adjacent measured borehole combinations in the borehole set D are traversed to obtain all the comparison section maps.

The device for constructing the stratum distribution mode constrained borehole contrast profile model comprises a processor and a computer program which is stored on a memory and can run on the processor, wherein the processor realizes the method when executing the program.

In the embodiment, the method is only based on visualization display of the contrasted section and the drilling data by the winCharts plug-in the NET 4.5 environment, and the method can also be applied to visualization tools such as echarts and the like.

While the invention has been described in connection with what is presently considered to be the most practical and preferred embodiment, it is to be understood that the invention is not to be limited to the disclosed embodiment, but on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the appended claims.

Claims (10)

1. A method for constructing a drilling contrast profile model constrained by a stratum distribution mode is characterized by specifically comprising the following steps of:

(1) reading a drilling set D and a drilling stratum thickness set T of actually measured drilling holes on the section line;

(2) reading any adjacent actual measurement drilling combination from the drilling set D, establishing a missing mark list SL of a missing extraction stratum between the adjacent actual measurement drilling combinations based on the drilling stratum thickness set T, establishing a stratum missing interval set DS according to the SL, and analyzing a stratum distribution mode based on the DS to form a stratum distribution mode mark list M;

(3) interpolating to generate virtual drill holes between the adjacent actually measured drill hole combinations, deducing the bottom surface elevation and the ground elevation of each virtual drill hole, and storing the virtual drill holes between the adjacent actually measured drill hole combinations into a drill hole set PD;

(4) preliminarily deducing the thickness of each stratum of each borehole in the borehole set PD according to a stratum pinch-out rule according to a stratum distribution mode mark list M, and storing the thickness of each stratum into a stratum thickness matrix PT;

(5) optimizing the stratum thickness according to the bottom surface elevation and the ground elevation of each virtual drilling hole, and constructing all stratum surfaces between the adjacent actually-measured drilling hole combinations;

(6) generating a comparison profile between the adjacent actually measured drilling hole combinations according to the constructed ground surface and the drilling hole set PD;

(7) and (5) circularly executing the steps (2) to (6) until all the adjacent measured borehole combinations in the borehole set D are traversed to obtain all the comparison section maps.

2. The method of constructing a formation distribution pattern-constrained borehole contrast profile model according to claim 1, wherein: the step (1) comprises the following steps;

(1-1) reading the drilling information of the actually measured drilling holes on the section line to obtain a drilling hole set D ═ D_i1, · dn } for 1, 2, · dn; wherein d is_iRepresenting the ith borehole, dn representing the number of boreholes;

(1-2) reading the borehole stratum information of the actually measured borehole on the section line to obtain a stratum thickness set T ═ T_i，j1, · dn, | i ═ 1, 2; j ═ 1, 2, ·, sni }; wherein, t_i，jIndicating a borehole d_iSni denotes the borehole d_iThe total number of strata encountered.

3. The method of constructing a formation distribution pattern-constrained borehole contrast profile model according to claim 1, wherein: the step (2) comprises the following steps;

(2-1) reading any group of adjacent measured boreholes D from the borehole set D_iAnd d_i+1；

(2-2) traversing the borehole stratum thickness set T, matching the borehole d according to the following formula_iAnd drilling a hole d_i+1The list SL ═ { SL ═ of missing flags of the inter-missing extraction formation_j|j＝1，2，...，sni}；

Wherein, sl is_jFormation lost flag, t, representing formation j_i，j、t_i+1，jIndicating a borehole d_i、d_i+1Thickness of middle formation j, sni denotes borehole d_iThe total number of strata encountered by drilling;

(2-3) stratum-missing flag list SL_iBuilding a borehole d_iAnd drilling a hole d_i+1Inter-strata missing interval set DS ═ DS_u1, 2.., dsn }; wherein ds_uThe triplet represents the u-th formation missing interval, dsn represents the number of missing intervals, and start represents the interval ds_uIs the starting formation number, end is the interval ds_uThe end stratum number of (1), count being the interval ds_uThe number of strata in (1);

(2-4) for each stratigraphic missing interval DS in the stratigraphic missing interval set DS_uOver the interval ds_uAll formations in the interval are missing and the interval ds is calculated as follows_uInner stratigraphic distribution pattern mark m_uAnd storing the stratum distribution mode mark list M ═ M_u|u＝1，2，...，dsn}；

4. The method of constructing a formation distribution pattern-constrained borehole contrast profile model according to claim 3, characterized by: the step (2-3) comprises the following steps;

(2-3-1) sequentially traversing the missing flag list SL until SL_jIf not, assigning the value of j to start and the value of count to 1;

(2-3-2) continue to traverse the missing flag list SL when j is less than or equal to sni and | SL_j| is equal to | sl_startWhen |, the count value is increased by 1; otherwise, assigning the end value as j-1, and newly building an interval element ds_u(start, end, count), and add to the missing interval set DS, and then return to step (2-3-2); wherein, | | is an absolute value symbol;

and (2-3-3) circularly executing the steps (2-3-1) to (2-3-2) until the missing mark list SL is traversed to obtain a stratum missing interval set DS.

5. The method of constructing a formation distribution pattern-constrained borehole contrast profile model according to claim 1, wherein: the step (3) comprises the following steps;

(3-1) traversing the set DS, and acquiring all the DS in the DS_u(count) as the number of virtual boreholes to be inserted, the maximum value, denoted as the maximum interval missing formation number max, wherein ds_u(count) is the interval ds_uThe number of strata in (1);

(3-2) creating a set of empty boreholes PD ═ PD { PD ═ PD }_k1, 2, max +2}, storing the virtual borehole and the read adjacent measured borehole d_iAnd d_i+1(ii) a Wherein pd is₁Assigning a value to the measured borehole d_i，pd_max+2Assigning a value to the measured borehole d_i+1，pd₂To pd_max+11 st to max virtual drill holes are formed in sequence;

(3-3) sequentially traversing the borehole set PD, and calculating a virtual borehole PD to be inserted according to the following formula through linear interpolation_kCoordinates and bottom elevation pd of_k(bottom)；

Wherein pd is_k(x)And pd_k(y) is a virtual borehole pd_kAbscissa and ordinate of (d)_i(x) And d_i(y) is a bore hole d_iAbscissa and ordinate of (d)_i+1(x) And d_i+1(y) is a bore hole d_i+1Abscissa and ordinate of (d)_i(bottom)、d_i+1(bottom) and pd_k(bottom) respectively represents the drilled holes d_i、d_i+1And pd_kThe hole bottom elevation;

(3-4) according to the virtual borehole pd_kCoordinate extraction of virtual borehole pd from DEM data_kGround elevation pd_k(h) Max +1, completes each borehole assignment in the set of boreholes PD.

6. The method of constructing a formation distribution pattern-constrained borehole contrast profile model according to claim 1, wherein: the step (4) comprises the following steps;

(4-1) constructing a formation thickness matrix PT ═ { PT) for the borehole set PD_k，j1, · 2, · max + 2; j ═ 1, 2, ·, sni }; wherein pt is_k，jThickness of jth formation of kth borehole in borehole set PD, and initial values of 0, sni for borehole d_iThe total number of strata encountered by drilling;

(4-2) reading any stratum interval DS from the set DS_uWhen m is_uWhen 2, interval ds_uIf the inner stratum is continuous, calculating the thickness of the stratum according to the following formula, and executing the step (4-7); otherwise, executing the step (4-3);

in the formula, t_i，j、t_i+1，jIndicating a borehole d_i、d_i+1The thickness of the middle formation j;

(4-3) drilling a hole d_iNumber of drilled formations quenched count₁Setting the initial value to 0, drilling a hole d_i+1Number of drilled formations quenched count₂An initial value of 0 is assigned to the formation j and an initial value of ds is assigned to the formation j_u(start)，ds_u(start) is ds_uThe starting formation number of (a);

(4-4) updating the number of holes drilled in the ground which have been extinguished, count, according to the following formula₁And count₂And deducing a virtual borehole subscript PN of the formation pinch-out;

in the formula, sl_jA formation missing flag representing a formation j in the missing flag list SL;

(4-5) estimating the thickness pt of the formation in each of the pseudo drilled holes according to the following equation_k，j；

(4-6) judging whether j is equal to ds_u(end)，ds_u(end) is ds_uIf yes, ending iteration, executing step (4-7), otherwise, setting j to j +1, and executing step (4-4);

and (4-7) circularly executing the steps (4-2) - (4-6) until all stratum intervals in the set DS are traversed, and finishing the calculation of the stratum thickness matrix PT.

7. The method of constructing a formation distribution pattern constrained borehole contrast profile model according to claim 1, wherein: the step (5) comprises the following steps;

(5-1) optimizing a formation thickness matrix PT based on the inferred bottom surface elevation and ground elevation of the virtual borehole;

(5-2) calculating the top surface elevation of each formation for each borehole according to the following equation, and storing the top surface elevation matrix UP ═_k，j|k＝1，2，...，max+2；j＝1，2，...，sni}；

In the formula, up_k，jTo representDrilling PD in drilling set PD_kTop surface elevation, pd, of formation j_k(h) Representing a borehole pd_kGround elevation of (pt)_k，jRepresenting a borehole pd in a formation thickness matrix PT_kSni denotes the borehole d_iThe total number of the stratums encountered by drilling, and max +2 represents the number of the drilling holes in the drilling hole set PD;

(5-3) calculating the elevation of the bottom surface of each stratum of each borehole according to the following formula, and storing the bottom surface elevation matrix DOWN ═ DOWN_k，j|k＝1，2，...，max+2；j＝1，2，...，sni}；

down_k，j＝up_k，j-pt_k，j，k＝1，2，...，max+2，j＝1，2，...，sni

In the formula, down_k，jRepresenting a borehole PD in a borehole set PD_kThe elevation of the bottom surface of the stratum j;

(5-4) constructing a borehole d according to the missing mark list SL, the optimized formation thickness matrix PT, the top elevation matrix UP and the bottom elevation matrix DOWN_iAnd d_i+1All the terrain levels in between.

8. The method of constructing a stratigraphic distribution pattern-constrained borehole contrasted profile model according to claim 7, characterized by: the step (5-1) comprises;

(5-1-1) calculating each borehole PD in the borehole set PD according to the following formula based on the formation thickness matrix PT_kIs pd of_k(depth)；

pd_k(depth)＝∑_jpt_k，j，k＝1，2，...，max+2

(5-1-2) calculating the borehole pd according to the following formula based on the surface level and the bottom level of the borehole_kBottom surface elevation pd after optimization_k(depth) ' and deducing the thickness pt ' of each stratum '_k，jUpdating a stratum thickness matrix PT;

in the formula, pd_k(h) Indicating a boreholepd_kGround elevation, pd_k(bottom) denotes the borehole pd_kThe bottom surface elevation of (a).

9. The method of constructing a formation distribution pattern-constrained borehole contrast profile model according to claim 7, wherein: the step (5-4) comprises the following steps;

(5-4-1) setting j ═ 1;

(5-4-2) constructing an empty point set P;

(5-4-3) traversing each element of jth column of the formation thickness matrix PT when PT_k，jWhen the coordinate is greater than 0, the coordinate (k, up) is set_k，j) And (k, down)_k，j) Adding to the point set P;

(5-4-4) when the jth element SL of the flag list SL is missing_jWhen the coordinate point (k, up) is determined as-1, the jth column of the thickness matrix PT is sequentially traversed, the row number fk of the first zero value is read, and the coordinate point (k, up) is determined_fk，j) Adding to the point set P; when sl is_jWhen 1, sequentially traversing jth column of the thickness matrix PT, reading the last row index lk of zero value, and dividing the coordinate point (pd)_k(x)，pd_k(y)，up_lk，j) Adding to the point set P;

(5-4-5) if the point set P is not empty, constructing a polygon based on the point set P as a ground-level sp_i，j；

(5-4-6) judging whether j is equal to the number of columns of the stratum thickness matrix PT, if so, ending the iteration, otherwise, changing j to j +1, and returning to execute the step (5-4-2).

10. A formation distribution pattern constrained borehole contrast profile model construction apparatus comprising a processor and a computer program stored on a memory and executable on the processor, wherein: the processor, when executing the program, implements the method of any of claims 1-9.

Images