Method and apparatus for constructing borehole comparison profile model constrained by formation distribution pattern

Abstract

The present invention discloses a method and device for constructing a borehole correlation profile model constrained by a formation distribution pattern. The method includes: (1) reading a borehole set D of measured boreholes on a profile line and a borehole formation thickness set T; (2) reading any adjacent combination of measured boreholes from the borehole set D, analyzing the formation distribution pattern, and forming a formation distribution pattern flag list; (3) interpolating to generate virtual boreholes between the adjacent combinations of measured boreholes, and forming a borehole set PD with the adjacent combinations of measured boreholes; (4) based on the formation pinch-out rule, initially inferring the thicknesses of each formation of each borehole in the borehole set PD according to the formation distribution pattern flag list M, and storing them in a formation thickness matrix PT; (5) optimizing the formation thickness and constructing all formation surfaces between the adjacent combinations of measured boreholes; (6) generating a correlation profile according to the constructed formation surfaces and the borehole set PD. The modeling of the present invention has higher quality, higher efficiency, and higher automation degree.

Description

Technical Field

[0001] The present invention relates to geographic information technology, and particularly to a method and device for constructing a borehole comparison profile model constrained by a formation distribution pattern. Background Art

[0002] Constructing a formation correlation profile based on boreholes is the main means for stratigraphers to correlate formations. The main construction method is to convert the attribute data of boreholes and formations into geometric elements for drawing, thereby realizing the visualization of profile data. The discontinuous distribution or uneven thickness of formations in space often increases the complexity of constructing a formation profile. Therefore, the scientific understanding of the formation distribution pattern and the correct application of inference rules are the basis and prerequisite for reasonably constructing a formation correlation profile in complex situations.

[0003] Existing methods mostly adopt simple formation pinch-out processing rules, resulting in unsatisfactory profile modeling quality in complex situations such as continuous pinch-out of multiple formations and discontinuous missing of multiple formations, which conflicts with the formation development law. In addition, although the generated profile can be further edited and optimized manually by experts, the efficiency is low and the profile modeling quality is not stable enough. Summary of the Invention

[0004] Object of the Invention: Aiming at the problems existing in the prior art, the present invention provides a method and device for constructing a borehole comparison profile model constrained by a formation distribution pattern with higher efficiency and better quality.

[0005] Technical Solution: The method for constructing a borehole comparison profile model constrained by a formation distribution pattern according to the present invention includes:

[0006] (1) Reading a borehole set D of measured boreholes on a profile line and a borehole formation thickness set T;

[0007] (2) Reading any adjacent combination of measured boreholes from the borehole set D, establishing a missing flag list SL of the missing extracted formations between the adjacent combinations of measured boreholes based on the borehole formation thickness set T, constructing a formation missing interval set DS according to SL, and then parsing the formation distribution pattern based on DS to form a formation distribution pattern flag list M;

[0008] (3) Interpolating to generate virtual boreholes between the adjacent combinations of measured boreholes, inferring the bottom elevation and ground elevation of each virtual borehole, and storing the adjacent combinations of measured boreholes and the virtual boreholes therebetween into a borehole set PD;

[0009] (4) Based on the formation pinch-out law, initially inferring the thickness of each formation of each borehole in the borehole set PD according to the formation distribution pattern flag list M, and storing it into a formation thickness matrix PT;

[0010] (5) Optimize the formation thickness according to the bottom elevation of each virtual borehole and the ground elevation, and construct all formation surfaces between the adjacent measured borehole combinations;

[0011] (6) Generate a comparison profile between the adjacent measured borehole combinations according to the constructed formation surfaces and the borehole set PD;

[0012] (7) Loop through steps (2)-(6) until all adjacent measured borehole combinations in the borehole set D are traversed to obtain all comparison profiles.

[0013] Further, step (1) includes;

[0014] (1-1) Read the borehole information of the measured boreholes on the profile line to obtain the borehole set D = {d i | i = 1, 2,..., dn}; where, d i represents the i-th borehole, and dn represents the number of boreholes;

[0015] (1-2) Read the borehole formation information of the measured boreholes on the profile line to obtain the formation thickness set T = {t i，j | i = 1, 2,..., dn; j = 1, 2,..., sni}; where, t i，j represents the thickness of the j-th formation of the borehole d i , and sni represents the total number of formations drilled by the borehole d i .

[0016] Further, step (2) includes;

[0017] (2-1) Read any group of adjacent measured boreholes d i and d i+1 from the borehole set D;

[0018] (2-2) Traverse the formation thickness set T, and match the missing flag list SL = {sl i | j = 1, 2,..., sni} of the missing extracted formations between the borehole d i+1 according to the following formula; j | j = 1, 2,..., sni};

[0019]

[0020] Where, sl j represents the formation missing flag of the formation j, t i，j , t i+1，j represent the thicknesses of the formation j in the boreholes d i , d i+1 , and sni represents the total number of formations drilled by the borehole d i ;

[0021] (2 - 3) Based on the list of formation missing flags SL i , construct borehole d i and the set of formation missing intervals DS between borehole d i+1 and borehole d u as DS = {ds u | u = 1, 2,..., dsn}; where ds u = (start, end, count) is a triple representing the u-th formation missing interval, dsn is the number of missing intervals, start is the starting formation number of interval ds u , end is the ending formation number of interval ds u , and count is the number of formations in interval ds

[0022] (2 - 4) For each formation missing interval ds in the set of formation missing intervals DS u , traverse the missing situation of all formations within interval ds u , and calculate the formation distribution pattern flag m u within interval ds according to the following formula, and store it in the list of formation distribution pattern flags M = {m u | u = 1, 2,..., dsn}; u | u = 1, 2,..., dsn};

[0023]

[0024] Furthermore, step (2 - 3) includes;

[0025] (2 - 3 - 1) Sequentially traverse the list of missing flags SL until sl j is not equal to 0, assign the value of j to start, and assign the value of count as 1;

[0026] (2 - 3 - 2) Continue to traverse the list of missing flags SL. When j is less than or equal to sni and |sl j | is equal to |sl start |, increment the value of count by 1; otherwise, assign the value of j - 1 to end, create a new interval element ds u = (start, end, count), and add it to the set of formation missing intervals DS, then return to step (2 - 3 - 2); where || is the absolute value symbol;

[0027] (2 - 3 - 3) Loop through steps (2 - 3 - 1) to (2 - 3 - 2) until the list of missing flags SL is traversed completely to obtain the set of formation missing intervals DS.

[0028] Furthermore, step (3) includes;

[0029] (3 - 1) Traverse the set DS to obtain all ds in DSu (count) is maximized, denoted as the maximum number of missing strata in the interval max, which is used as the number of virtual boreholes to be inserted. Among them, ds u (count) is the number of strata in the interval ds u ;

[0030] (3-2) Create an empty borehole set PD = {pd k |k = 1, 2,..., max + 2}, which stores virtual boreholes and the adjacent measured boreholes d i and d i+1 ; Among them, the value of pd1 is assigned to the measured borehole d i , and the value of pd max+2 is assigned to the measured borehole d i+1 , and pd2 to pd max+1 are the 1st to max virtual boreholes in sequence;

[0031] (3-3) Traverse the borehole set PD in sequence, and calculate the coordinates and bottom elevation pd k of the virtual borehole to be inserted through linear interpolation according to the following formula k (bottom);

[0032]

[0033] Among them, pd k (x) and pd k (y) are the abscissa and ordinate of the virtual borehole pd k , d i (x) and d i (y) are the abscissa and ordinate of the borehole d i , d i+1 (x) and d i+1 (y) are the abscissa and ordinate of the borehole d i+1 , d i (bottom), d i+1 (bottom) and pd k (bottom) represent the bottom elevation of the boreholes d i , d i+1 and pd k respectively;

[0034] (3-4) Extract the ground elevation pd k of the virtual borehole pd k from the DEM data according to the coordinates of the virtual borehole pd k (h), k = 2,..., max + 1, and complete the assignment of each borehole in the borehole set PD.

[0035] Furthermore, step (4) includes;

[0036] (4-1) Construct the formation thickness matrix PT of the borehole set PD = {pt k，j | k = 1, 2,..., max + 2; j = 1, 2,..., sni}; where pt k，j represents the thickness of the j-th formation of the k-th borehole in the borehole set PD, and the initial value is assigned as 0, and sni represents the total number of formations encountered by borehole d i ;

[0037] (4-2) Read any formation interval ds from the set DS u , when m u equals 2, the formations within the interval ds u are continuous, then calculate the formation thickness according to the following formula and execute step (4-7); otherwise, execute step (4-3);

[0038]

[0039] In the formula, t i，j , t i+1，j represent the thickness of formation j in boreholes d i , d i+1 ;

[0040] (4-3) Initialize the number of formations with pinched-out in borehole d i count1 to 0, initialize the number of formations with pinched-out in borehole d i+1 count2 to 0, initialize formation j to ds u (start), ds u (start) is the starting formation number of ds u ;

[0041] (4-4) According to the following formula, update the number of boreholes with pinched-out formations count1 and count2, and infer the virtual borehole subscript PN where the formation is pinched out;

[0042]

[0043] In the formula, sl j represents the formation missing flag of formation j in the missing flag list SL;

[0044] (4-5) According to the following formula, infer the thickness pt of the formation in each virtual borehole k，j ;

[0045]

[0046] (4-6) Judge whether j is equal to ds u (end), ds u (end) is the ending formation number of ds uThe termination formation number. If so, end the iteration and execute step (4-7). Otherwise, set j = j + 1 and execute step (4-4);

[0047] (4-7) Loop and execute steps (4-2)-(4-6) until all formation intervals in the set DS are traversed, and complete the calculation of the formation thickness matrix PT.

[0048] Furthermore, step (5) includes;

[0049] (5-1) Optimize the formation thickness matrix PT based on the inferred bottom elevation and ground elevation of the virtual borehole;

[0050] (5-2) Calculate the top elevation of each formation of each borehole according to the following formula and store it in the top elevation matrix UP = {up k，j |k = 1, 2,..., max + 2; j = 1, 2,..., sni};

[0051]

[0052] In the formula, up k，j represents the top elevation of formation j of borehole pd in the borehole set PD k , pd k (h) represents the ground elevation of borehole pd k , pt k，j represents the thickness of the jth formation of borehole pd in the formation thickness matrix PT k , sni represents the total number of formations drilled by borehole d i , and max + 2 represents the number of boreholes in the borehole set PD;

[0053] (5-3) Calculate the bottom elevation of each formation of each borehole according to the following formula and store it in the bottom elevation matrix DOWN = {down k，j |k = 1, 2,..., max + 2; j = 1, 2,..., sni};

[0054] down k，j = up k，j - pt k，j , k = 1, 2,..., max + 2, j = 1, 2,..., sni

[0055] In the formula, down k，j represents the bottom elevation of formation j of borehole pd in the borehole set PD; k

[0056] (5-4) Construct borehole d according to the missing flag list SL, the optimized formation thickness matrix PT, the top elevation matrix UP, and the bottom elevation matrix DOWN(5-4) Construct borehole d according to the missing flag list SL, the optimized formation thickness matrix PT, the top elevation matrix UP, and the bottom elevation matrix DOWN i and di+1 All the strata between them.

[0057] Further, step (5-1) includes;

[0058] (5-1-1) Calculate the depth pd k of each borehole pd in the borehole set PD based on the formation thickness matrix PT according to the following formula k (depth);

[0059] pd k (depth) = ∑ j pt k，j , k = 1, 2,..., max + 2

[0060] (5-1-2) Based on the ground elevation and bottom elevation of the borehole, calculate the optimized bottom elevation pd k (depth)' of the borehole pd according to the following formula, and infer the formation thickness pt' k , and update the formation thickness matrix PT; k，j

[0061]

[0062] In the formula, pd k (h) represents the ground elevation of the borehole pd k , and pd k (bottom) represents the bottom elevation of the borehole pd k .

[0063] Further, step (5-4) includes;

[0064] (5-4-1) Set j = 1;

[0065] (5-4-2) Construct an empty point set P;

[0066] (5-4-3) Traverse each element in the j-th column of the formation thickness matrix PT. When pt k，j is greater than 0, add the coordinates (k, up k，j ) and (k, down k，j ) to the point set P;

[0067] (5-4-4) When the j-th element sl j of the missing flag list SL = -1, sequentially traverse the j-th column of the thickness matrix PT, read the row number fk of the first zero value, and add the coordinate point (k, up fk，j ) to the point set P; when sl j = 1, sequentially traverse the j-th column of the thickness matrix PT, read the row subscript lk of the last zero value, and add the coordinate point (pd k ​(x), pd k (y), up lk，j ) are added to the point set P;

[0068] (5 - 4 - 5) If the point set P is not empty, a polygon is constructed based on the point set P as the ground layer sp i，j ;

[0069] (5 - 4 - 6) Determine whether j is equal to the number of columns of the formation thickness matrix PT. If so, end the iteration; otherwise, set j = j + 1 and return to execute step (5 - 4 - 2).

[0070] The device for constructing a borehole correlation profile model constrained by the formation distribution pattern of the present invention includes a processor and a computer program stored on a memory and executable on the processor. When the processor executes the program, the above - mentioned method is implemented.

[0071] Beneficial effects: Compared with the prior art, the significant advantages of the present invention are as follows: The present invention can better apply to relatively complex formation distribution patterns such as continuous pinch - out and continuous absence of formations, which not only improves the modeling quality of borehole correlation profiles but also has a high degree of automation and higher efficiency. Description of the Drawings

[0072] Figure 1 are the DEM data and borehole data used in this embodiment;

[0073] Figure 2 is a screenshot of the borehole table used in this embodiment;

[0074] Figure 3 is a screenshot of the formation table used in this embodiment;

[0075] Figure 4 is a flowchart of the method for constructing a borehole correlation profile model constrained by the formation distribution pattern provided by the present invention;

[0076] Figure 5 is a schematic diagram of the formation thickness inference of the present invention;

[0077] Figure 6 is d 33 and d 34 comparative profile diagrams before thickness optimization;

[0078] Figure 7 is d 33 and d 34 comparative profile diagrams after thickness optimization. Detailed Embodiments

[0079] The technical solution of the present invention will be further described in detail below. In this embodiment, the Quaternary boreholes and DEM data in a certain experimental area are selected ( Figure 1 ,Figure 2 and Figure 3 ) As experimental data, the projection coordinate system of the DEM and borehole vector data used in the experiment is the WGS84 coordinate system. This will be further illustrated below with reference to the accompanying drawings and by describing a specific embodiment.

[0080] As Figure 4 shown, this embodiment provides a method for modeling a borehole comparison profile constrained by a formation distribution pattern, which specifically includes the following steps.

[0081] (1) Read the borehole set D of the measured boreholes on the profile line and the borehole formation thickness set T.

[0082] This step specifically includes:

[0083] (1-1) Read the borehole information of the measured boreholes on the profile line to obtain the borehole set D = {d i | i = 1, 2,..., dn}; where d i represents the i-th borehole, and dn represents the number of boreholes; in this embodiment, dn = 69;

[0084] (1-2) Read the borehole formation information of the measured boreholes on the profile line to obtain the formation thickness set T = {t i，j | i = 1, 2,..., dh; j = 1, 2,..., sni}; where t i，j represents the thickness of the j-th formation of borehole d i , and sni represents the total number of formations drilled by borehole d i ; in this embodiment, sni = 35..

[0085] (2) Read any adjacent measured borehole combination from the borehole set D, establish a missing flag list SL of the missing extracted formations between the adjacent measured borehole combinations based on the borehole formation thickness set T, construct a formation missing interval set DS according to SL, and then analyze the formation distribution pattern based on DS to form a formation distribution pattern flag list M.

[0086] This step specifically includes:

[0087] (2-1) Read any group of adjacent measured boreholes d i and d i+1 from the borehole set D; in this embodiment, taking i = 33 as an example;

[0088] (2-2) Traverse the borehole formation thickness set T, and match the missing flag list SL = {sl i | j = 1, 2,..., sni} of the missing extracted formations between borehole d i+1 and borehole d j according to the following formula;

[0089]

[0090] Among them, sl j represents the formation missing flag of formation j, and t i，j , t i+1，j represent the thickness of formation j in borehole d i , d i+1 . sni represents the total number of formations encountered in borehole d i ;

[0091] (2-3) Based on the formation missing flag list SL i , construct the formation missing interval set DS = {ds i} between borehole d i+1 and borehole d u | u = 1, 2,..., dsn}; where ds u = (start, end, count) is a triple, representing the u-th formation missing interval, dsn is the number of missing intervals, start is the starting formation number of interval ds u , end is the ending formation number of interval ds u , and count is the number of formations in interval ds u ; in this embodiment, when i = 33, ds1 = (1, 1, 0), dsn = 9; the specific steps are as follows:

[0092] (2-3-1) Traverse the missing flag list SL sequentially until sl j is not equal to 0, assign the value of j to start, and assign the value of count to 1;

[0093] (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 tart |, increment the value of count by 1; otherwise, assign the value of j-1 to end, create a new interval element ds u = (start, end, count), and add it to the formation missing interval set DS, then return to step (2-3-2); where || is the absolute value symbol;

[0094] (2-3-3) Loop through steps (2-3-1) to (2-3-2) until the missing flag list SL is traversed, and obtain the formation missing interval set DS;

[0095] (2-4) For each formation missing interval ds u in the formation missing interval set DS, traverse the missing situation of all formations in interval ds u , and calculate the interval ds according to the following formulau The formation distribution pattern flag m within u is stored in the formation distribution pattern flag list M = {m u | u = 1, 2,..., dsn};

[0096]

[0097] In this embodiment, when i = 33, M = {0, 1, 0, 1, 0, 1, 0, 2, 0}.

[0098] (3) Interpolate to generate virtual boreholes between the adjacent measured borehole combinations, and infer the bottom elevation and ground elevation of each virtual borehole, and store the adjacent measured borehole combinations and the virtual boreholes between them in the borehole set PD.

[0099] This step specifically includes:

[0100] (3 - 1) Traverse the set DS to obtain the maximum value of all ds u (count) in DS, denoted as the maximum number of missing formation layers max, which is used as the number of virtual boreholes to be inserted. Among them, ds u (count) is the number of formation layers in the interval ds u ; In this embodiment, when i = 33, max = 11;

[0101] (3 - 2) Create an empty borehole set PD = {pd k | k = 1, 2,..., max + 2}, which stores the virtual boreholes and the read adjacent measured boreholes d i and d i+1 ; Among them, the value of pd1 is assigned to the measured borehole d i , the value of pd max+2 is assigned to the measured borehole d i+1 , and pd2 to pd max+1 are the 1st to max virtual boreholes in sequence;

[0102] (3 - 3) Sequentially traverse the borehole set PD, and calculate the coordinates and bottom elevation pd k (bottom) of the virtual borehole pd to be inserted through linear interpolation according to the following formula; k (bottom);

[0103]

[0104] Among them, pd k (x) and pd k (y) are the abscissa and ordinate of the virtual borehole pd k , and d i (x) and d i (y) are the abscissa and ordinate of the borehole d iThe abscissa and ordinate of d i+1 d(x) and d i+1 d(y) are the abscissa and ordinate of borehole d i+1 The abscissa and ordinate of d i d(bottom), d i+1 d(bottom) and pd k d(bottom) respectively represent the bottom elevation of boreholes d i , d i+1 and pd k ;

[0105] (3 - 4) Extract the ground elevation pd k of the virtual borehole pd from the DEM data according to the coordinates of the virtual borehole pd k pd(h), k = 2,..., max + 1, and complete the assignment of each borehole in the borehole set PD, where the ground elevation pd k pd(h), the calculation formula is as follows: k

[0106]

[0107] where x0 is the abscissa of the DEM origin, y0 is the ordinate of the DEM origin, C is the DEM pixel size, is the floor function symbol. In this embodiment, C = 100 meters.

[0108] (4) Based on the law of stratigraphic pinch-out, preliminarily infer the thickness of each stratum of each borehole in the borehole set PD according to the list M of stratigraphic distribution pattern flags, and store it in the stratum thickness matrix PT.

[0109] This step specifically includes;

[0110] (4 - 1) Construct the stratum thickness matrix PT of the borehole set PD = {pt k，j |k = 1, 2,..., max + 2; j = 1, 2,..., sni}; where pt k，j represents the thickness of the jth stratum of the kth borehole in the borehole set PD, and the initial value is assigned to 0, and sni represents the total number of strata drilled by borehole d i ;

[0111] (4 - 2) Read any stratum interval ds u from the set DS. When m u is equal to 2, the strata in the interval ds u are continuous, then calculate the stratum thickness according to the following formula and execute step (4 - 7); otherwise, execute step (4 - 3);

[0112]

[0113] ​Where t i，j and t i+1，j represent the thickness of formation j in boreholes d i and d i+1 ;

[0114] (4 - 3) Initialize the number of formations in the extinguished borehole count1 to 0 for borehole d i , initialize the number of formations in the extinguished borehole count2 to 0 for borehole d i+1 , and initialize formation j to ds u (start), where ds u (start) is the starting formation number of ds u ;

[0115] (4 - 4) Update the number of extinguished boreholes count1 and count2 according to the following formula, and infer the virtual borehole subscript PN where the formation is extinguished;

[0116]

[0117] Where sl j represents the formation missing flag of formation j in the missing flag list SL;

[0118] (4 - 5) Infer the thickness pt k，j of the formation in each virtual borehole according to the following formula, as shown in Figure 5 ;

[0119]

[0120] (4 - 6) Judge whether j is equal to ds u (end), where ds u (end) is the ending formation number of ds u . If so, end the iteration and execute step (4 - 7), otherwise set j = j + 1 and execute step (4 - 4);

[0121] (4 - 7) Loop and execute steps (4 - 2) - (4 - 6) until all formation intervals in the set DS are traversed, and complete the calculation of the formation thickness matrix PT, as shown in Figure 6 ;

[0122] (5) Optimize the formation thickness based on the bottom elevation and ground elevation of each virtual borehole, and construct all formation surfaces between the adjacent measured borehole combinations.

[0123] This step specifically includes;

[0124] (5 - 1) Optimize the formation thickness matrix PT based on the inferred bottom elevation and ground elevation of the virtual borehole; the specific steps are:

[0125] (5-1-1) Calculate the depth pd of each borehole pd in the borehole set PD based on the formation thickness matrix PT according to the following formula k of the borehole pd k (depth);

[0126] pd k (depth) = ∑ j pt k，j , k = 1, 2,..., max + 2

[0127] (5-1-2) Based on the ground elevation and bottom elevation of the borehole, calculate the optimized bottom elevation pd k (depth)' of the borehole pd, and infer the formation thickness pt' k , and update the formation thickness matrix PT; k，j

[0128]

[0129] In the formula, pd k (h) represents the ground elevation of the borehole pd k , pd k (bottom) represents the bottom elevation of the borehole pd k ;

[0130] (5-2) Calculate the top elevation of each formation of each borehole according to the following formula and store it in the top elevation matrix UP = {up k，j |k = 1, 2,..., max + 2; j = 1, 2,..., sni};

[0131]

[0132] In the formula, up k，j represents the top elevation of the j-th formation of the borehole pd in the borehole set PD k , pd k (h) represents the ground elevation of the borehole pd k , pt k，j represents the thickness of the j-th formation of the borehole pd in the formation thickness matrix PT k , sni represents the total number of formations encountered by the borehole d i , max + 2 represents the number of boreholes in the borehole set PD;

[0133] (5-3) Calculate the bottom elevation of each formation of each borehole according to the following formula and store it in the bottom elevation matrix DOWN = {down k，j |k = 1, 2,..., max + 2; j = 1, 2,..., sni};

[0134] down k，j= up k，j -pt k，j , where \(k = 1, 2, \cdots, max + 2\), \(j = 1, 2, \cdots, sni\)

[0135] In the formula, down k，j represents the bottom elevation of the formation \(j\) of the borehole \(pd\) in the borehole set \(PD\); k

[0136] (5 - 4) Construct all the formation surfaces between boreholes \(d\) i and \(d\) i+1 according to the missing flag list \(SL\), the optimized formation thickness matrix \(PT\), the top elevation matrix \(UP\), and the bottom elevation matrix \(DOWN\). The specific implementation steps are as follows:

[0137] (5 - 4 - 1) Set \(j = 1\);

[0138] (5 - 4 - 2) Construct an empty point set \(P\);

[0139] (5 - 4 - 3) Traverse each element in the \(j\)-th column of the formation thickness matrix \(PT\). When \(pt\) k，j is greater than 0, add the coordinates \((k, up\) k，j ) and \((k, down\) k，j ) to the point set \(P\);

[0140] (5 - 4 - 4) When the \(j\)-th element \(sl\) of the missing flag list \(SL\) j = -1, sequentially traverse the \(j\)-th column of the thickness matrix \(PT\), read the row number \(fk\) of the first zero value, and add the coordinate point \((k, up\) fk，j ) to the point set \(P\); when \(sl\) j = 1, sequentially traverse the \(j\)-th column of the thickness matrix \(PT\), read the row subscript \(lk\) of the last zero value, and add the coordinate point \((pd\) k (x), \(pd\) k (y), \(up\) lk，j ) to the point set \(P\);

[0141] (5 - 4 - 5) If the point set \(P\) is not empty, construct a polygon based on the point set \(P\) as the formation surface \(sp\) i，j ;

[0142] (5 - 4 - 6) Determine whether \(j\) is equal to the number of columns of the formation thickness matrix \(PT\). If so, end the iteration; otherwise, set \(j = j + 1\) and return to execute step (5 - 4 - 2).

[0143] (6) Generate the comparison profile between the adjacent measured borehole combinations according to the constructed formation surfaces and the borehole set \(PD\). In this embodiment, when \(i = 33\), the comparison profile between \(d\) 33 and \(d\) 34 is as follows​Figure 7 as shown

[0144] (7) Repeat steps (2)-(6) until all adjacent measured borehole combinations in the borehole set D are traversed, and all comparison profiles are obtained.

[0145] The device for constructing a borehole comparison profile model constrained by the formation distribution pattern described in the present invention includes a processor and a computer program stored in a memory and executable on the processor. When the processor executes the program, the above method is implemented.

[0146] In this embodiment, only the winCharts plug-in under the.NET 4.5 environment is used to visually display the comparison profile and borehole data. This method can also be applied to visualization tools such as echarts.

[0147] The above-disclosed is only a preferred embodiment of the present invention, and it cannot be used to limit the scope of the rights of the present invention. Therefore, equivalent changes made according to the claims of the present invention still fall within the scope covered by the present invention.

Claims

1. A method for constructing a borehole comparison profile model constrained by a stratigraphic distribution pattern, characterized in that The method specifically includes: (1) Reading a borehole set D of measured boreholes on a profile line and a borehole formation thickness set T; (2) Reading any adjacent combination of measured boreholes from the borehole set D, establishing a missing flag list SL of the missing extracted formation between adjacent combinations of measured boreholes based on the borehole formation thickness set T, constructing a formation missing interval set DS according to SL, and then parsing the formation distribution pattern based on DS to form a formation distribution pattern flag list M; (3) Interpolating and generating virtual boreholes between the adjacent combinations of measured boreholes, inferring the bottom elevation and ground elevation of each virtual borehole, and storing the adjacent combinations of measured boreholes and the virtual boreholes between them into a borehole set PD; (4) Based on the formation pinch-out rule, initially inferring the formation thickness of each borehole in the borehole set PD according to the formation distribution pattern flag list M and storing it into a formation thickness matrix PT; Step (4) includes: (4-1) Construct the formation thickness matrix PT of the borehole set PD = {pt k，j | k = 1, 2,..., max + 2; j = 1, 2,..., sni}; where pt k，j represents the thickness of the j-th formation of the k-th borehole in the borehole set PD, and the initial value is assigned as 0 for all, and sni represents the total number of formations encountered by borehole d i drilled through; (4-2) Read any formation interval ds from the set DS u When m u equals 2, and the formations within the interval ds u are continuous, calculate the formation thickness according to the following formula and execute step (4-7); otherwise, execute step (4-3); where t i，j and t i+1，j represent the thickness of formation j in boreholes d i and d i+1 ; (4-3) Drill hole d i Initialize the number of strata in the drilled hole that has thinned out, count1, to 0, and drill hole d i+1 Initialize the number of strata in the drilled hole that has thinned out, count2, to 0, and initialize stratum j to ds u (start), ds u (start) is the starting stratum number of ds u ; (4-4) According to the following formula, updating the number of boreholes with pinched-out formations count1 and count2, and inferring the subscript PN of the virtual borehole where the formation pinches out; where sl j represents the formation missing flag of formation j in the missing flag list SL; (4-5) Infer the thickness pt of the formation in each virtual borehole according to the following formula k，j ; (4-6) Determine whether j is equal to ds u (end), ds u (end) is ds u The termination formation number of ds. If so, end the iteration and execute step (4-7). Otherwise, set j = j + 1 and execute step (4-4); (4-7) Repeatedly executing steps (4-2)-(4-6) until all formation intervals in the set DS are traversed to complete the calculation of the formation thickness matrix PT; (5) Optimizing the formation thickness according to the bottom elevation and ground elevation of each virtual borehole, and constructing all formation surfaces between the adjacent combinations of measured boreholes; (6) Generating a comparison profile between the adjacent combinations of measured boreholes according to the constructed formation surfaces and the borehole set PD; (7) Repeatedly executing steps (2)-(6) until all adjacent combinations of measured boreholes in the borehole set D are traversed to obtain all comparison profiles.

2. The method for constructing a borehole comparison profile model constrained by a stratigraphic distribution pattern according to claim 1, characterized in that: Step (1) includes: (1-1) Read the borehole information of the measured boreholes on the profile line to obtain the borehole set D = {d i | i = 1, 2,..., dn}; where d i represents the i-th borehole, and dn represents the number of boreholes; (1-2) Read the borehole formation information of the actual boreholes on the profile line to obtain the formation thickness set T = {t i，j | i = 1, 2,..., dn; j = 1, 2,..., sni}; where, t i，j represents the thickness of the j-th formation of borehole d i , and sni represents the total number of formations encountered by borehole d i .

3. The method for constructing a borehole comparison profile model constrained by a stratigraphic distribution pattern according to claim 1, characterized in that: Step (2) includes: (2-1) Read any set of adjacent measured boreholes d from the borehole set D i and d i+1 ; (2-2) Traverse the set T of borehole formation thicknesses, and match borehole d according to the following formula i and borehole d i+1 to extract the missing flag list SL of the formation between them = {sl j | j = 1, 2,..., sni}; Among them, sl j represents the formation missing flag of formation j, t i，j and t i+1，j represent the thickness of formation j in boreholes d i and d i+1 , and sni represents the total number of formations encountered in borehole d i ; (2-3) Based on the list of formation missing markers SL i , construct borehole d i and borehole d i+1 to form a set of formation missing intervals DS = {ds u | u = 1, 2,..., dsn}; where ds u = (start, end, count) is a triple representing the u-th formation missing interval, dsn is the number of missing intervals, start is the starting formation number of interval ds u , end is the ending formation number of interval ds u , and count is the number of formations in interval ds u ; (2-4) For each stratum missing interval ds in the stratum missing interval set DS u , traverse the interval ds u The missing conditions of all strata within the interval ds are calculated according to the following formula u The stratigraphic distribution pattern within m u , stored in the stratum distribution mode flag list M = {m u |u=1,2,...,dsn}; 4. The method for constructing a borehole comparison profile model constrained by a stratigraphic distribution pattern according to claim 3, characterized in that: Step (2-3) includes: (2-3-1)Traverse the missing flag list SL sequentially until sl j is not equal to 0, assign the value of j to start, and assign 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 start |, increment the count value by 1; otherwise, assign the value of j - 1 to end, create a new interval element ds u = (start, end, count), and add it to the formation missing interval set DS, then return to step (2-3-2); where, || is the absolute value symbol; (2-3-3) Repeatedly executing steps (2-3-1) to (2-3-2) until the missing flag list SL is traversed to obtain the formation missing interval set DS.

5. The method for constructing a borehole comparison profile model constrained by the formation distribution pattern according to claim 1, wherein: Step (3) includes: (3-1) Traverse the set DS to obtain all ds in DS u (count), and record the maximum value as the maximum number of missing strata in the interval, denoted as max, which is the number of virtual boreholes to be inserted. Among them, ds u (countt) is the number of strata in the interval ds u ; (3-2) Create a new empty set of boreholes PD = {pd k | k = 1, 2,..., max + 2}, store the virtual boreholes and the adjacent measured boreholes d i and d i+1 ; among them, the value of pd1 is assigned to the measured borehole d i , pd max+2 value is assigned to the measured borehole d i+1 , pd2 to pd max+1 are the 1st to max virtual boreholes in sequence; (3-3) Traverse the drilling set PD in sequence, and calculate the coordinates and bottom elevation pd of the virtual drilling hole pd to be inserted through linear interpolation according to the following formula k ; k (bottom); Among them, pd k (x) and pd k (y) are the abscissa and ordinate of the virtual borehole pd k . d i (x) and d i (y) are the abscissa and ordinate of the borehole d i . d i+1 (x) and d i+1 (y) are the abscissa and ordinate of the borehole d i+1 . d i (bottom), d i+1 (bottom) and pd k (bottom) respectively represent the bottom elevations of the boreholes d i , d i+1 and pd k ; (3-4) According to the virtual borehole pd k Extract the coordinates of the virtual borehole pd from the DEM data k Ground elevation pd of k (h), k = 2,..., max + 1, complete the assignment of each borehole in the borehole set PD.

6. The method for constructing a borehole comparison profile model constrained by the formation distribution pattern according to claim 1, wherein: Step (5) includes: (5-1) Optimizing the formation thickness matrix PT based on the inferred bottom elevation and ground elevation of the virtual boreholes; (5-2) Calculate the elevation of the top surface of each formation for each borehole according to the following formula and store it in the top surface elevation matrix UP = {up k，j | k = 1, 2,..., max + 2; j = 1, 2,..., sni}; where, up k，j represents the elevation of the top surface of the formation j of the borehole pd in the borehole set PD k , pd k (h) represents the ground elevation of the borehole pd k , pt k，j represents the thickness of the j-th formation of the borehole pd in the formation thickness matrix PT k , sni represents the total number of formations encountered by the borehole d i , max + 2 represents the number of boreholes in the borehole set PD; (5-3) Calculate the elevation of the bottom surface of each stratum for each borehole according to the following formula and store it in 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 where, down k，j represents the bottom elevation of formation j of borehole pd in borehole set PD k ; (5-4) Construct all the formation surfaces between boreholes d i and d i+1 according to the missing flag list SL, the optimized formation thickness matrix PT, the top elevation matrix UP, and the bottom elevation matrix DOWN.

7. The method for constructing a borehole comparison profile model constrained by the formation distribution pattern according to claim 6, wherein: Step (5-1) includes: (5-1-1) The depth pd of each borehole pd in the borehole set PD is calculated according to the following formula based on the formation thickness matrix PT k ; k (depth); pd k (depth) = ∑ j pt k，j , k = 1, 2,..., max + 2 (5-1-2) Calculate the borehole pd based on the ground elevation and bottom elevation of the borehole according to the following formula k Optimized bottom elevation pd k (depth)′, and infer the thickness pt′ of each stratum k，j , and update the stratum thickness matrix PT; where pd k (h) represents the ground elevation of borehole pd k and pd k (bottom) represents the bottom elevation of borehole pd k ​ 8. The method for constructing a borehole comparison profile model constrained by the formation distribution pattern according to claim 6, wherein: Step (5-4) includes: (5-4-1) Setting j = 1; (5-4-2) Constructing an empty point set P; (5-4-3) Traverse each element in the j-th column of the formation thickness matrix PT. When pt k，j is greater than 0, add the coordinates (k, up k，j ) and (k, down k，j ) to the point set P; (5-4-4) When the j-th element sl of the missing flag list SL j = -1, sequentially traverse the j-th column of the thickness matrix PT, read the row number fk of the first zero value, and add the coordinate point (k, up fk，j ) to the point set P; when sl j = 1, sequentially traverse the j-th column of the thickness matrix PT, read the row subscript lk of the last zero value, and add the coordinate point (pd k (x), pd k (y), up lk,j ) to the point set P; (5-4-5) If the point set P is not empty, construct a polygon based on the point set P as the ground surface layer sp i，j ; (5-4-6) Judging whether j is equal to the number of columns of the formation thickness matrix PT. If so, ending the iteration; otherwise, setting j = j + 1 and returning to execute step (5-4-2).

9. A device for constructing a borehole comparison profile model constrained by a formation distribution pattern, comprising a processor and a computer program stored on a memory and executable on the processor, characterized in that: When the processor executes the program, it implements the method described in any one of claims 1-8.

Images