A method and device for determining fracture characterization information and a computer device

Abstract

The present specification relates to the technical field of fracture parameter estimation, and particularly relates to a fracture characterization information determination method and device and computer equipment. The fracture characterization information determination method comprises: determining fracture dip angle average value, fracture dip angle standard deviation and fracture dip angle variation range data based on a formation microresistivity scan image corresponding to a fracture; updating a preset equivalent elastic stiffness matrix of organic matter-rich shale without a fracture to obtain a target equivalent elastic stiffness matrix by using the fracture dip angle average value, the fracture dip angle standard deviation and the fracture dip angle variation range data; and processing the target equivalent elastic stiffness matrix based on an orthogonal anisotropic parameter formula to obtain target orthogonal anisotropic parameters, so as to determine anisotropic characterization information of the fracture. By using the embodiment of the present specification, the quantification of fracture information in an orthogonal anisotropic medium is realized, and the accuracy of the quantification result is improved compared with the quantification of the fracture in an isotropic medium.

Description

TECHNICAL FIELD

[0001] The present specification relates to the technical field of fracture parameter estimation, and particularly relates to a fracture characterization information determination method and device and computer equipment. BACKGROUND

[0002] At present, a model of equivalent elastic properties of fractured rock is usually used to determine fracture characterization information, and the model of equivalent elastic properties of fractured rock is constructed under the assumption that the background medium of fractured rock is isotropic. However, in many fractured formations, the background medium is formed by thin layers, that is, the fractures develop in an inherent anisotropic medium. At this time, the model of equivalent elastic properties of rock considering only the isotropic medium has poor applicability. Moreover, when the model is used to analyze the rock physical properties and seismic response mechanism of rock, the accuracy of the results is also low.

[0003] How to extend the model of equivalent elastic properties of fractured rock to an anisotropic background medium to improve the accuracy of fracture characterization information is a problem to be solved in the prior art. SUMMARY

[0004] To solve the problems in the prior art, the embodiments of the present specification provide a fracture characterization information determination method and device and computer equipment, which realize the construction of a model of equivalent elastic properties of fractured organic-rich shale, determine a target equivalent elastic stiffness matrix based on the model of equivalent elastic properties of fractured rock, and further determine fracture characterization information quantified in an anisotropic medium based on the target equivalent elastic stiffness matrix and an anisotropic parameter formula of orthogonal medium, so as to realize the extension of the model of equivalent elastic properties of fractured rock to an anisotropic background medium and improve the accuracy of quantified fracture characterization information.

[0005] To solve the above technical problems, the specific technical solutions of the present specification are as follows:

[0006] On one hand, the embodiments of the present specification provide a fracture characterization information determination method, which comprises,

[0007] determining fracture dip angle average value, fracture dip angle standard deviation and fracture dip angle variation range data based on a formation microresistivity scan image corresponding to the fractures;

[0008] updating a preset equivalent elastic stiffness matrix of fractured organic-rich shale without fractures to contain fractures by using the fracture dip angle average value, the fracture dip angle standard deviation and the fracture dip angle variation range data, to obtain a target equivalent elastic stiffness matrix; and

[0009] processing the target equivalent elastic stiffness matrix based on an anisotropic parameter formula of orthogonal medium to obtain target anisotropic parameters of orthogonal medium, so as to determine anisotropic characterization information of the fractures.

[0010] Further, the data of the average of the fracture dip angle, the standard deviation of the fracture dip angle and the range of the fracture dip angle are used to update the preset equivalent elastic stiffness matrix of the shale without fractures to obtain a target equivalent elastic stiffness matrix, and the target equivalent elastic stiffness matrix further includes,

[0011] Based on the data of the average of the fracture dip angle, the standard deviation of the fracture dip angle and the range of the fracture dip angle, determine the flexibility change data; and

[0012] Based on the preset equivalent elastic stiffness matrix of the shale without fractures and the flexibility change data, determine the target equivalent elastic stiffness matrix.

[0013] Further, the flexibility change data further includes,

[0014]

[0015] Wherein, the θ represents the fracture dip angle, the σ represents the standard deviation of the fracture dip angle, and the θ0 represents the average of the fracture dip angle. max The θ represents the maximum fracture dip angle in the fracture dip angle, the θ represents the minimum fracture dip angle in the fracture dip angle, the θ represents the fracture dip angle, the θ0 represents the average of the fracture dip angle, the σ represents the standard deviation of the fracture dip angle, and the N represents the number of the fracture dip angle. min The θ represents the fracture dip angle, the σ represents the standard deviation of the fracture dip angle, and the θ0 represents the average of the fracture dip angle. θ The ΔS represents the preset flexibility change data, and the G(θ) represents the fracture dip angle distribution function.

[0016] Further, the fracture dip angle distribution function further includes,

[0017]

[0018] Wherein, the θ represents the fracture dip angle, the σ represents the standard deviation of the fracture dip angle, and the θ0 represents the average of the fracture dip angle.

[0019] Further, the determination of the target equivalent elastic stiffness matrix based on the preset equivalent elastic stiffness matrix of the shale without fractures and the flexibility change data further includes,

[0020] Invert the preset equivalent elastic stiffness matrix of the shale without fractures to obtain a preset flexibility matrix;

[0021] Calculate the sum of the flexibility matrix and the flexibility change data to obtain a target flexibility matrix; and

[0022] Invert the target flexibility matrix to obtain the target equivalent elastic stiffness matrix.

[0023] Further, the target equivalent elastic stiffness matrix further includes,

[0024]

[0025] wherein, the I represents a unit matrix, the C b represents the preset equivalent elastic stiffness matrix of the shale without the fracture, the p represents the density, and the ΔS' is the flexibility change data.

[0026] Further, the orthogonal medium anisotropy parameter formula includes a first plane orthogonal medium anisotropy parameter formula, a second plane orthogonal medium anisotropy parameter formula and a third plane orthogonal medium anisotropy parameter formula, and the processing of the target equivalent elastic stiffness matrix based on the orthogonal medium anisotropy parameter formula to obtain the target orthogonal medium anisotropy parameter further includes,

[0027] processing the target equivalent elastic stiffness matrix based on the first plane orthogonal medium anisotropy parameter formula to obtain a first plane target orthogonal medium anisotropy parameter;

[0028] processing the target equivalent elastic stiffness matrix based on the second target orthogonal medium anisotropy parameter formula to obtain a second plane target orthogonal medium anisotropy parameter; and

[0029] processing the target equivalent elastic stiffness matrix based on the third target orthogonal medium anisotropy parameter formula to obtain a third plane target orthogonal medium anisotropy parameter.

[0030] In another aspect, the embodiments of the present specification also provide a fracture representation information determination method, including,

[0031] a first determination unit configured to determine fracture dip angle average value, fracture dip angle standard deviation and fracture dip angle variation range data based on a formation microresistivity scan image corresponding to the fracture;

[0032] an updating unit configured to update a preset equivalent elastic stiffness matrix of the shale without the fracture by using the fracture dip angle average value, the fracture dip angle standard deviation and the fracture dip angle variation range data to obtain a target equivalent elastic stiffness matrix; and

[0033] a processing unit configured to process the target equivalent elastic stiffness matrix based on an orthogonal medium anisotropy parameter formula to obtain a target orthogonal medium anisotropy parameter, so as to determine anisotropy representation information of the fracture.

[0034] In another aspect, the embodiments of the present specification also provide a computer device, including a memory, a processor and a computer program stored in the memory and executable on the processor, and the processor executes the computer program to realize the above method.

[0035] In another aspect, the embodiments of the present specification also provide a computer readable storage medium, having stored thereon computer instructions, which, when executed by a processor, implement the method described above.

[0036] With the embodiments of the present specification, the equivalent elastic stiffness matrix of the preset organic-rich shale without cracks is updated to obtain a target equivalent elastic stiffness matrix according to the average value of the crack dip angle, the standard deviation of the crack dip angle and the crack dip angle variation range data determined from the micro-resistivity scan image of the formation corresponding to the cracks. Then, the target equivalent elastic stiffness matrix is processed based on the orthogonal anisotropic parameter formula to obtain target orthogonal anisotropic parameters, so as to determine the anisotropy characterization information of the cracks. Thus, the model of the rock equivalent elastic properties of the organic-rich shale with cracks is constructed. Then, the target equivalent elastic stiffness matrix and the orthogonal anisotropic parameter formula are used to determine the crack characterization information quantified in the anisotropic medium, so that the anisotropic background medium is expanded for the model of the rock equivalent elastic properties of the cracks, and the accuracy of the quantified crack characterization information is improved. BRIEF DESCRIPTION OF DRAWINGS

[0037] In order to more clearly illustrate the technical solutions in the embodiments of the present specification or the prior art, the drawings needed to be used in the embodiments or the prior art description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present specification, and other drawings can be obtained by those skilled in the art without creative labor.

[0038] Figure 1 An implementation system schematic diagram of a crack characterization information determination method according to an embodiment of the present specification is shown;

[0039] Figure 2 A flowchart of a crack characterization information determination method according to an embodiment of the present specification is shown;

[0040] Figure 3 A flowchart of a target equivalent elastic stiffness matrix determination method according to an embodiment of the present specification is shown;

[0041] Figure 4 A flowchart of a target equivalent elastic stiffness matrix determination method according to another embodiment of the present specification is shown;

[0042] Figure 5A A schematic diagram of a crack distribution according to an embodiment of the present specification is shown;

[0043] Figure 5BA principle diagram of a method for constructing a preset uncracked organic-rich shale model is shown in the embodiment of the present specification.

[0044] Figure 5C A schematic diagram of first elastic stiffness data is shown in the embodiment of the present specification.

[0045] Figure 5D A schematic diagram of second elastic stiffness data is shown in the embodiment of the present specification.

[0046] Figure 5E A schematic diagram of third elastic stiffness data is shown in the embodiment of the present specification.

[0047] Figure 6 A structural schematic diagram of a fracture characterization information determination device is shown in the embodiment of the present specification.

[0048] Figure 7 A structural schematic diagram of a computer device is shown in the embodiment of the present specification.

[0049]

Explanation of reference numerals

[0050] 101, user terminal;

[0051] 102, server;

[0052] 501, first sub-target equivalent elastic stiffness data;

[0053] 502, second sub-target equivalent elastic stiffness data;

[0054] 503, third sub-target equivalent elastic stiffness data;

[0055] 504, first sub-target orthogonal medium anisotropy parameter;

[0056] 505, second sub-target orthogonal medium anisotropy parameter;

[0057] 506, third sub-target orthogonal medium anisotropy parameter;

[0058] 511, fourth sub-target equivalent elastic stiffness data;

[0059] 512, fifth sub-target equivalent elastic stiffness data;

[0060] 513, sixth sub-target equivalent elastic stiffness data;

[0061] 514, fourth sub-target orthogonal medium anisotropy parameter;

[0062] 515, fifth sub-target orthogonal medium anisotropy parameter;

[0063] 516, sixth sub-target orthogonal medium anisotropy parameter;

[0064] 521, seventh sub-target equivalent elastic stiffness data;

[0065] 522, eighth sub-target equivalent angle elastic stiffness data;

[0066] 523, ninth sub-target equivalent angle elastic stiffness data;

[0067] 524, seventh sub-target orthogonal medium anisotropy parameter;

[0068] 525, eighth sub-target orthogonal medium anisotropy parameter;

[0069] 526, ninth sub-target orthogonal medium anisotropy parameter;

[0070] 610, first determination unit;

[0071] 620, update unit;

[0072] 630, processing unit;

[0073] 702, computer device;

[0074] 704, processing device;

[0075] 706, storage resource;

[0076] 708, driving mechanism;

[0077] 710, input / output module;

[0078] 712, input device;

[0079] 714, output device;

[0080] 716, presentation device;

[0081] 718, graphical user interface;

[0082] 720, network interface;

[0083] 722, communication link;

[0084] 724, communication bus. DETAILED DESCRIPTION

[0085] The technical solutions in the embodiments of the present specification will be described clearly and completely in the following combined with the drawings in the embodiments of the present specification. Obviously, the described embodiments are only part of the embodiments of the present specification, not all the embodiments. Based on the embodiments in the present specification, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present specification.

[0086] It should be noted that the terms "first", "second", and the like in the description and in the claims of the present specification and the above drawings are used to distinguish similar objects and are not necessarily used to describe a particular sequential or chronological order. It should be understood that the use of such terms as "first", "second", and the like, if used herein, are used to distinguish one object from another, rather than to denote a particular temporal or chronological order. It is to be understood that the data thus designated can be interchanged, where appropriate, so that the embodiments of the present specification described herein can be carried out in other than the order shown or described herein. Furthermore, the terms "comprise" and "have" and any variations thereof, are intended to cover non-exclusive inclusion, for example, processes, methods, devices, products, or equipment that include a list of steps or units as processes, methods, devices, products, or equipment are not necessarily limited to those steps or units clearly listed, but can include other steps or units not clearly listed or inherent to such processes, methods, products, or equipment.

[0087] It should be noted that the steps shown in the flowchart of the drawings can be executed in a computer system such as a set of computer executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described herein can be executed in a different order than that shown herein.

[0088] Figure 1 An implementation system schematic diagram of a fracture characterization information determination method according to an embodiment of the present specification is shown, which can include a user terminal 101 and a server 102, the user terminal 101 and the server 102 communicate with each other through a network, the network can include a local area network (LAN), a wide area network (WAN), the Internet or a combination thereof, and is connected to a website, a user device (such as a computing device) and a backend system. After receiving the formation micro-resistivity scan image, the server 102 determines the fracture dip angle average value, the fracture dip angle standard deviation and the fracture dip angle variation range data based on the formation micro-resistivity scan image; uses the fracture dip angle average value, the fracture dip angle standard deviation and the fracture dip angle variation range data to update the preset equivalent elastic stiffness matrix of the shale containing no fractures, to obtain a target equivalent elastic stiffness matrix; processes the target equivalent elastic stiffness matrix based on an orthogonal anisotropic parameter formula, to obtain a target orthogonal anisotropic parameter, and sends the target orthogonal anisotropic parameter to the user terminal 101. In addition, after determining the target orthogonal anisotropic parameter, the server 102 can also determine the anisotropy characterization information of the fracture based on the target orthogonal anisotropic parameter, and send the anisotropy characterization information of the fracture to the user terminal 101.

[0089] It should be noted that the server can also be network connected with the acquisition device to receive the formation micro-resistivity scan image collected and sent by the acquisition device.

[0090] Optionally, the server 102 can be a node of a cloud computing system (not shown in the figure), or each server 102 can be a separate cloud computing system, including multiple computers interconnected by a network and working as a distributed processing system.

[0091] In an optional embodiment, the user terminal 101 can include electronic devices such as but not limited to smart phones, acquisition devices, desktop computers, tablets, notebooks, smart speakers, digital assistants, augmented reality (AR) / virtual reality (VR) devices, smart wearable devices, etc. Optionally, the operating system running on the electronic device can include but is not limited to Android system, IOS system, Linux, Windows, etc.

[0092] In addition, it needs to be explained that, Figure 1 The illustrated is only one application environment provided by the present specification, and in actual application, a plurality of user terminals 101 can also be included, which is not limited by the present specification.

[0093] Figure 2 The illustrated is a flowchart of a fracture characterization information determination method according to an embodiment of the present specification. The determination process of the fracture characterization information is described in the figure, but more or fewer operation steps can be included based on conventional or non-inventive labor. The order of steps listed in the embodiment is only one of the many execution orders of the steps, and does not represent the only execution order. In actual system or device product execution, the method order shown in the embodiment or the figure can be executed in sequence or in parallel. Specifically, as Figure 2 As shown, the method can include:

[0094] S210, determining the average value of the fracture dip angle, the standard deviation of the fracture dip angle, and the range data of the fracture dip angle based on the formation microresistivity scan image corresponding to the fracture;

[0095] S220, updating the preset equivalent elastic stiffness matrix of the organic-rich shale without fracture by using the average value of the fracture dip angle, the standard deviation of the fracture dip angle, and the range data of the fracture dip angle, to obtain a target equivalent elastic stiffness matrix;

[0096] S230, processing the target equivalent elastic stiffness matrix based on the orthogonal anisotropy parameter formula to obtain a target orthogonal anisotropy parameter, so as to determine the anisotropy characterization information of the fracture.

[0097] According to the embodiments of the present specification, the average value of the fracture dip, the standard deviation of the fracture dip and the range of the fracture dip data are determined according to the micro-resistivity scanning image of the formation corresponding to the fracture, the equivalent elastic stiffness matrix of the preset shale without fracture is updated to contain the fracture, and the target equivalent elastic stiffness matrix is obtained. Then, the target equivalent elastic stiffness matrix is processed based on the orthogonal anisotropic parameter formula of the medium to obtain the target orthogonal anisotropic parameter of the medium, so as to determine the anisotropy characterization information of the fracture. Thus, the model of the rock equivalent elastic property of the shale containing the fracture is constructed; then, based on the target equivalent elastic stiffness matrix obtained from the model of the rock equivalent elastic property of the shale containing the fracture and the orthogonal anisotropic parameter formula of the medium, the fracture characterization information quantified in the anisotropic medium is determined, thereby the anisotropic background medium expansion is carried out for the model of the rock equivalent elastic property of the fracture, and the accuracy of the quantified fracture characterization information is improved.

[0098] According to one embodiment of the present specification, the average value of the fracture dip, the standard deviation of the fracture dip and the range of the fracture dip data can be determined based on the received micro-resistivity scanning image (FMI) of the formation corresponding to the fracture, the fracture information is determined by analyzing the received micro-resistivity scanning image (FMI) of the formation corresponding to the fracture, and the average value of the fracture dip, the standard deviation of the fracture dip and the range of the fracture dip data are determined by Gaussian function fitting for the fracture information. The fracture information may, for example, be any information representing the fracture, such as the number of fractures and the fracture dip.

[0099] The preset equivalent elastic stiffness matrix of the shale without fracture is obtained. The preset equivalent elastic stiffness matrix of the shale without fracture can be a preset equivalent elastic stiffness matrix of the shale without fracture, or can be obtained by processing the shale rock based on the existing method for determining the equivalent elastic stiffness matrix of the shale without fracture.

[0100] After the preset equivalent elastic stiffness matrix of the shale without fracture is determined, the compliance matrix is obtained by performing compliance processing on the preset equivalent elastic stiffness matrix of the shale without fracture. Then, the average value of the fracture dip, the standard deviation of the fracture dip and the range of the fracture dip data are used to correct the compliance matrix to obtain the target equivalent elastic stiffness matrix.

[0101] Specifically, the average value of the fracture dip, the standard deviation of the fracture dip and the range of the fracture dip data are used to correct the compliance matrix to obtain the target equivalent elastic stiffness matrix, for example, the average value of the fracture dip, the standard deviation of the fracture dip and the range of the fracture dip data are used to determine the correction number of each data in the compliance matrix, and then the sum of each data in the compliance matrix and the corresponding correction number is determined to obtain the target equivalent elastic stiffness matrix.

[0102] The target equivalent elastic stiffness matrix may be, for example, a 6x6 matrix, as shown in the following equation (1).

[0103]

[0104] wherein C ij are constants.

[0105] The orthogonal medium anisotropy parameter formula includes a first plane orthogonal medium anisotropy parameter formula, a second plane orthogonal medium anisotropy parameter formula, and a third plane orthogonal medium anisotropy parameter formula. The target equivalent elastic stiffness matrix is processed based on the orthogonal medium anisotropy parameter formula to obtain a target orthogonal medium anisotropy parameter. For example, the first plane target orthogonal medium anisotropy parameter is obtained by processing the target equivalent elastic stiffness matrix based on the first plane orthogonal medium anisotropy parameter formula; the second plane target orthogonal medium anisotropy parameter is obtained by processing the target equivalent elastic stiffness matrix based on the second target orthogonal medium anisotropy parameter formula; and the third plane target orthogonal medium anisotropy parameter is obtained by processing the target equivalent elastic stiffness matrix based on the third target orthogonal medium anisotropy parameter formula.

[0106] More specifically, the first plane is, for example, an [x2, x3] symmetry plane, and the first plane orthogonal medium anisotropy formula is as shown in the following equation (2).

[0107]

[0108] wherein ε, δ, and γ are target orthogonal medium anisotropy parameters, respectively.

[0109] The second plane is, for example, an [x1, x3] symmetry plane, and the second plane orthogonal medium anisotropy formula is as shown in the following equation (3).

[0110]

[0111] wherein ε, δ, and γ are target orthogonal medium anisotropy parameters, respectively.

[0112] The third plane is, for example, an [x1, x2] symmetry plane, and the third plane orthogonal medium anisotropy formula is as shown in the following equation (4).

[0113]

[0114] wherein ε, δ, and γ are target orthogonal medium anisotropy parameters, respectively.

[0115] The data in formula (1) is processed by formulas (2), (3) and (4) to obtain the anisotropic parameters of the nine target orthogonal media in the three planes.

[0116] After determining the anisotropy parameters of the target orthogonal medium, the anisotropy characterization information of the crack can be determined based on these parameters. This anisotropy characterization information can include, for example, information on the degree of anisotropy of the orthogonal medium and anisotropic velocity information.

[0117] Figure 3 The diagram shows a flowchart of a method for determining a target equivalent elastic stiffness matrix according to an embodiment of this specification. This figure describes a process for determining a target equivalent elastic stiffness matrix, but based on conventional or non-creative work, it may include more or fewer operational steps. Specifically, as shown... Figure 3 As shown, the method may include:

[0118] S321, Based on the average crack inclination angle, standard deviation of crack inclination angle, and range of crack inclination angle variation data, determine the compliance change data;

[0119] S322, based on the preset equivalent elastic stiffness matrix and flexibility change data of organic-rich shale without cracks, determine the target equivalent elastic stiffness matrix.

[0120] According to another embodiment of this specification, when determining the change in compliance data based on the average crack inclination angle, the standard deviation of crack inclination angle, and the range of crack inclination angle variation, the following formula (5) can be used, for example.

[0121]

[0122] Where, θ max The maximum crack dip angle, θ, represents the crack dip angle. min The minimum crack inclination angle is represented by θ, the crack inclination angle is represented by θ0, the average crack inclination angle is represented by σ, and the standard deviation of the crack inclination angle is represented by N. θ The transformation matrix is ​​represented by ΔS, which represents the preset compliance change data, and G(θ) represents the crack dip angle distribution function.

[0123] Specifically, the transformation matrix is ​​shown in the following formula (6).

[0124]

[0125] Where θ represents the crack dip angle.

[0126] The crack dip angle distribution function can be, for example, as shown in the following formula (7).

[0127]

[0128] wherein, θ represents the fracture dip angle, σ represents the fracture dip angle standard deviation, and θ0 represents the fracture dip angle average value.

[0129] It should be noted that, due to the introduction of the Gaussian function in the present specification to describe the distribution characteristics of the fracture dip angle, the information represented by the determined target equivalent elastic stiffness matrix is more abundant, and the anisotropy representation information of the fracture is more abundant.

[0130] Based on the preset equivalent elastic stiffness matrix of the shale without fractures and the flexibility change data, the target equivalent elastic stiffness matrix can be determined, for example, the flexibility of the preset equivalent elastic stiffness matrix of the shale without fractures is processed to obtain a preset flexibility matrix, and the flexibility matrix is corrected by using the flexibility change data to obtain the target equivalent elastic stiffness matrix as shown in the above formula (1).

[0131] Figure 4 The flowchart of another embodiment of a method for determining a target equivalent elastic stiffness matrix is shown. In this figure, another determination process of the target equivalent elastic stiffness matrix is described, but more or less operation steps can be included based on conventional or non-creative labor. Specifically, as shown in the above formula (1), the method can include: Figure 4

[0132] S4221, inverse of the preset equivalent elastic stiffness matrix of the shale without fractures is calculated to obtain a preset flexibility matrix;

[0133] S4222, the sum of the flexibility matrix and the flexibility change data is calculated to obtain a target flexibility matrix;

[0134] S4223, inverse of the target flexibility matrix is calculated to obtain a target equivalent elastic stiffness matrix.

[0135] According to another embodiment of the present specification, the inverse of the preset equivalent elastic stiffness matrix of the shale without fractures is calculated to obtain a preset flexibility matrix, which can be specifically shown in the following formula (8).

[0136]

[0137] wherein, C b is the preset equivalent elastic stiffness matrix of the shale without fractures, S b is the preset flexibility matrix, and C ib is a constant. After the inverse operation, the obtained preset flexibility matrix is shown in the following formula (9).

[0138]

[0139] wherein, S b ​S is a preset compliance matrix, and ΔS' is compliance change data. bi is a constant.

[0140] The preset compliance matrix and the compliance change data obtained via Figure 3 The sum of the compliance change data obtained via the steps shown in the figure obtains a target compliance matrix, specifically, as shown in the following formula (10).

[0141] S = S b + ΔS' formula (10)

[0142] Wherein, S is a target compliance matrix, S b is a preset compliance matrix, and ΔS' is compliance change data.

[0143] After determining the target compliance matrix, the target equivalent elastic stiffness matrix is obtained by inverting the target compliance matrix in the manner shown in the following formula (11).

[0144]

[0145] Wherein, C is a target equivalent elastic stiffness matrix, I represents a unit matrix, C b represents a preset uncracked organic-rich shale equivalent elastic stiffness matrix, ρ represents density, and ΔS' is compliance change data.

[0146] Figure 5A Fig. 1 shows a schematic diagram of a fracture distribution according to an embodiment of the present specification; Figure 5B Fig. 2 shows a schematic diagram of a method for constructing a preset uncracked organic-rich shale model according to an embodiment of the present specification; Figure 5C Fig. 3 shows a schematic diagram of first elastic stiffness data according to an embodiment of the present specification; Figure 5D Fig. 4 shows a schematic diagram of second elastic stiffness data according to an embodiment of the present specification; Figure 5E Fig. 5 shows a schematic diagram of third elastic stiffness data according to an embodiment of the present specification.

[0147] Taking an A shale work area as an example, well logging data and formation microresistivity scanning images corresponding to fractures are collected to determine fracture information, and then statistical analysis is performed on the fracture information, and the results are shown in Fig. 6. Figure 5A Fig. 6 shows a fracture dip angle histogram when the average value of the fracture dip angle is π / 2, and specifically Figure 5A Fig. 7 shows two kinds of fracture distributions when the standard deviation of the fracture dip angle is 1 and 4, respectively. Figure 5A Based on the statistical analysis of the fracture information, the average value of the fracture dip angle, the standard deviation of the fracture dip angle, and the range of the fracture dip angle data are obtained.

[0148] The preset uncracked organic-rich shale equivalent elastic stiffness matrix may be obtained, for example, as shown in the following formula (12). Figure 5BThe shown construction preset does not contain a crack rich organic shale model determined, specifically, first, the equivalent average (Backus average) is used to construct a clay mixture containing different clay minerals (clay fast); combined with the self-consistent approximation (SCA) theory and the differential equivalent medium theory (DEM), the saturated bound water pores are added to the obtained clay block to form a completely horizontally arranged clay-water mixture; the local clay block is rotated and superimposed to simulate the stratification of shale; the upper and lower limits (Hashin-Shtrikman-Walpole) are used to average the organic matter and siliceous, calcareous and other minerals (quartz, feldspar, pyrite, calcite and dolomite, etc.) to construct a homogeneous mixture, and the clay-water mixture is added using the anisotropic SCA model; finally, the anisotropic DEM model is used to add free gas / water pores to obtain shale, and based on the shale made, the equivalent elastic stiffness matrix of the shale without cracks is obtained.

[0149] Based on the equivalent elastic stiffness matrix of the shale without cracks Figure 2 The method shown is used to modify and update the equivalent elastic stiffness matrix of the shale without cracks to obtain a target equivalent elastic stiffness matrix, and then the target equivalent elastic stiffness matrix is processed using formulas (2), (3) and (4) to obtain target orthorhombic anisotropy parameters, as shown in details in Figure 5C 、 Figure 5D and Figure 5E

[0150] Specifically, when the standard deviation of the fracture inclination angle is 1 and the range of the fracture inclination angle is π / 2, the obtained target equivalent elastic stiffness matrix and target orthorhombic anisotropy parameters are as shown in Figure 5C As shown in Figure 5C , the first sub-target equivalent elastic stiffness data 501 shows the changes of C 11 , C 22 and C 33 in formula (1) with the standard deviation of the fracture inclination angle. The second sub-target equivalent elastic stiffness data 502 shows the changes of C 44 , C 55 and C 66 in formula (1) with the standard deviation of the fracture inclination angle. The third sub-target equivalent elastic stiffness data 503 shows the changes of C 13 , C 23 and C 12 in formula (1) with the standard deviation of the fracture inclination angle. The first sub-target orthorhombic anisotropy parameter 504 shows the changes of ε (1) , ε (2) and ε (3) determined by formulas (2), (3) and (4) with the standard deviation of the fracture inclination angle. The second sub-target orthorhombic anisotropy parameter 505 shows the changes of γ​(1) γ (2) and γ (3) The variation with the standard deviation of the crack dip angle. The anisotropy parameter 506 of the orthogonal medium of the third sub-target is shown as δ determined by formulas (2), (3) and (4). (1) δ (2) and δ (3) The change with the standard deviation of the crack inclination angle.

[0151] Similarly, the average crack dip angle is π / 2. When the crack dip angle varies within the range of π / 2, the obtained target equivalent elastic stiffness matrix and target orthogonal medium anisotropy parameters are as follows: Figure 5D As shown. Figure 5D As shown, the equivalent elastic stiffness data 511 of the fourth sub-target illustrates C in formula (1). 11 C 22 and C 33 The variation with the standard deviation of the crack inclination angle. The fifth sub-target's isoangular elastic stiffness data 512 shows C in formula (1). 44 C 55 and C 66 The variation with the standard deviation of the crack inclination angle. The sixth sub-target's isoangular elastic stiffness data 513 shows C in formula (1). 13 C 23 and C 12 The standard deviation of the crack dip angle varies. The anisotropy parameter 514 of the orthogonal medium of the fourth sub-target is shown as ε determined by formulas (2), (3) and (4). (1) ε (2) and ε (3) The variation with the standard deviation of the crack dip angle. The anisotropy parameter 515 of the orthogonal medium of the fifth sub-target is shown as γ determined by formulas (2), (3) and (4). (1) γ (2) and γ (3) The standard deviation of the crack inclination angle varies. The anisotropy parameter 516 of the orthogonal medium of the sixth sub-target is shown as δ determined by formulas (2), (3) and (4). (1) δ (2) and δ (3) The change with the standard deviation of the crack inclination angle.

[0152] Similarly, when the average crack inclination angle is π / 2 and the standard deviation of the crack inclination angle is 1, the obtained target equivalent elastic stiffness matrix and target orthogonal medium anisotropy parameters are as follows: Figure 5E As shown. Figure 5E As shown, the equivalent elastic stiffness data 521 of the seventh sub-target illustrates C in formula (1). 11 C 22 and C 33with the standard deviation of the fracture dip. The eighth sub-target equiangular elastic stiffness data 522 shows C 44 , C 5s , and C 66 with the standard deviation of the fracture dip. The ninth sub-target equiangular elastic stiffness data 523 shows C 13 , C 23 , and C 12 with the standard deviation of the fracture dip. The seventh sub-target orthorhombic anisotropy parameter 524 shows ε (1) , ε (2) , and ε (3) with the standard deviation of the fracture dip. The eighth sub-target orthorhombic anisotropy parameter 525 shows γ (1) , γ (2) , and γ (3) with the standard deviation of the fracture dip. The ninth sub-target orthorhombic anisotropy parameter 526 shows ε (1) , δ (2) , and δ (3) with the standard deviation of the fracture dip.

[0153] Figure 6 Fig. 1 shows a structural schematic diagram of a fracture characterization information determination apparatus according to an embodiment of the present specification. As shown in Fig. 1, the apparatus comprises, Figure 6

[0154] A first determination unit 610 is configured to determine fracture dip mean value, fracture dip standard deviation, and fracture dip variation range data based on a formation microresistivity scan image corresponding to the fracture.

[0155] An updating unit 620 is configured to update a preset equivalent elastic stiffness matrix of shale without fracture by using the fracture dip mean value, the fracture dip standard deviation, and the fracture dip variation range data, to obtain a target equivalent elastic stiffness matrix.

[0156] A processing unit 630 is configured to process the target equivalent elastic stiffness matrix based on an orthorhombic anisotropy parameter formula, to obtain a target orthorhombic anisotropy parameter, so as to determine anisotropy characterization information of the fracture.

[0157] Since the principle of solving the problem of the apparatus is similar to that of the method, the implementation of the apparatus can be referred to the implementation of the method, and the repeated parts will not be described herein.

[0158] As shown in Fig. 1, the apparatus comprises, Figure 7 ​Fig. 7 shows a schematic diagram of a computer device 702 according to an embodiment of the present specification. The apparatuses in the present specification can be the computer device in the present embodiment, which executes the method of the present specification. The computer device 702 can include one or more processing devices 704, such as one or more central processing units (CPUs), each of which can implement one or more hardware threads. The computer device 702 can also include any storage resources 706 for storing any kind of information, such as code, settings, data, and the like. Without limitation, for example, the storage resources 706 can include any one or combination of the following: any type of RAM, any type of ROM, a flash memory device, a hard disk, an optical disk, and the like. More generally, any storage resource can store information using any technology. Further, any storage resource can provide volatile or non-volatile retention of information. Further, any storage resource can represent a fixed or removable component of the computer device 702. In one case, the computer device 702 can perform any operation of the associated instructions when the processing device 704 executes the associated instructions stored in any storage resource or combination of storage resources. The computer device 702 also includes one or more drive mechanisms 708, such as a hard disk drive mechanism, an optical disk drive mechanism, and the like, for interacting with any storage resources.

[0159] The computer device 702 can also include an input / output module 710 (I / O) for receiving various inputs (via input devices 712) and for providing various outputs (via output devices 714). One particular output mechanism can include a presentation device 716 and an associated graphical user interface (GUI) 718. In other embodiments, the input / output module 710 (I / O), the input devices 712, and the output devices 714 can also not be included, just as a computer device in a network. The computer device 702 can also include one or more network interfaces 720 for exchanging data with other devices via one or more communication links 722. One or more communication buses 724 couple the above-described components together.

[0160] The communication links 722 can be implemented in any manner, for example, through a local area network, a wide area network (e.g., the Internet), a point-to-point connection, and the like, or any combination thereof. The communication links 722 can include any combination of hardwired links, wireless links, routers, gateway functionality, name servers, and the like governed by any protocol or combination of protocols.

[0161] The present specification also provides a computer-readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the above method.

[0162] The embodiment of the present specification further provides a computer program product, the computer program product comprising a computer program, the computer program being executed by a processor to implement the method described above.

[0163] Those skilled in the art should understand that the embodiment of the present specification can be provided as a method, a system, or a computer program product. Therefore, the present specification can adopt a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present specification can adopt the form of a computer program product implemented on one or more computer usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) containing computer usable program codes.

[0164] The present specification is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to the embodiments of the present specification. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of flows and / or blocks in the flowcharts 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, a special-purpose computer, an embedded processor, or other programmable data processing devices to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate a means for implementing the functions specified in the flowcharts and / or block diagrams. Figure One one or more flows and / or blocks Figure One an apparatus that carries out the functions specified in one or more flows and / or blocks.

[0165] These computer program instructions can also be stored in a computer readable storage medium capable of guiding a computer or other programmable data processing devices to work in a specific way, so that the instructions stored in the computer readable storage medium produce a product including instruction apparatus, which implements the functions specified in the flowcharts and / or block diagrams. Figure One one or more flows and / or blocks Figure One an apparatus that carries out the functions specified in one or more flows and / or blocks.

[0166] These computer program instructions can also be loaded into a computer or other programmable data processing device, so that a series of operation steps are executed on the computer or other programmable device to produce a computer implemented process, so that the instructions executed on the computer or other programmable device provide a means for implementing the functions specified in the flowcharts and / or block diagrams. Figure One one or more flows and / or blocks Figure One Figure One an apparatus that carries out the functions specified in one or more flows and / or blocks.

[0167] The above specific embodiments further illustrate the purpose, technical solutions and beneficial effects of the present specification, and it should be understood that the above are only specific embodiments of the present specification and are not used to limit the protection scope of the present specification, and any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present specification shall be included in the protection scope of the present specification.

Claims

1. A method for determining crack characterization information, characterized in that, include: Based on the formation microresistivity scanning images corresponding to the fractures, the average fracture dip angle, standard deviation of fracture dip angle, and range of fracture dip angle variation data were determined. Using the average value of the crack dip angle, the standard deviation of the crack dip angle, and the range of crack dip angle variation data, the preset equivalent elastic stiffness matrix of organic-rich shale without cracks is updated with cracks to obtain the target equivalent elastic stiffness matrix. as well as The target equivalent elastic stiffness matrix is ​​processed based on the orthogonal medium anisotropic parameter formula to obtain the target orthogonal medium anisotropic parameters, so as to determine the anisotropic characterization information of the crack. The step of updating the preset equivalent elastic stiffness matrix of organic-rich shale without fractures by using the average value of the fracture dip angle, the standard deviation of the fracture dip angle, and the range of fracture dip angle variation data, to obtain the target equivalent elastic stiffness matrix, includes: Based on the average crack inclination angle, the standard deviation of the crack inclination angle, and the range of crack inclination angle variation data, the compliance change data are determined; and Based on the preset equivalent elastic stiffness matrix of organic-rich shale without cracks and the flexibility change data, the target equivalent elastic stiffness matrix is ​​determined. The compliance change data includes: Wherein, the θ max The maximum crack inclination angle, θ, represents the crack inclination angle. min The minimum crack inclination angle is characterized by θ, the crack inclination angle is characterized by θ0, the average crack inclination angle is characterized by σ, and the standard deviation of the crack inclination angle is characterized by N. θ The transformation matrix is ​​represented by ΔS, which represents the preset compliance change data, and G(θ) represents the crack dip angle distribution function. The crack dip angle distribution function includes: Wherein, θ represents the crack inclination angle, σ represents the standard deviation of the crack inclination angle, and θ0 represents the average value of the crack inclination angle; The determination of the target equivalent elastic stiffness matrix based on the preset equivalent elastic stiffness matrix of organic-rich shale without cracks and the flexibility change data includes: Inverting the preset equivalent elastic stiffness matrix of the organic-rich shale without cracks yields the preset flexibility matrix. The target compliance matrix is ​​obtained by summing the compliance matrix and the compliance change data; and Inverting the target flexibility matrix yields the target equivalent elastic stiffness matrix; The target equivalent elastic stiffness matrix includes: Wherein, I represents the identity matrix, and C b The preset equivalent elastic stiffness matrix of the organic-rich shale without cracks is used to characterize the density, and ΔS' is the data on the change in flexibility.

2. The method according to claim 1, characterized in that, The orthogonal medium anisotropy parameter formula includes a first-plane orthogonal medium anisotropy parameter formula, a second-plane orthogonal medium anisotropy parameter formula, and a third-plane orthogonal medium anisotropy parameter formula. The target equivalent elastic stiffness matrix is ​​processed based on the orthogonal medium anisotropy parameter formula to obtain the target orthogonal medium anisotropy parameters, including: The equivalent elastic stiffness matrix of the target is processed based on the anisotropic parameter formula of the first planar orthogonal medium to obtain the anisotropic parameters of the first planar target orthogonal medium. Based on the formula for the anisotropic parameters of the second target orthogonal medium, the equivalent elastic stiffness matrix of the target is processed to obtain the anisotropic parameters of the second planar target orthogonal medium; and The equivalent elastic stiffness matrix of the target is processed based on the formula for the anisotropic parameters of the third target orthogonal medium to obtain the anisotropic parameters of the third plane target orthogonal medium.

3. A device for determining crack characterization information, characterized in that, include: The first determining unit is used to determine the average fracture dip angle, the standard deviation of the fracture dip angle, and the range of fracture dip angle variation based on the formation microresistivity scanning image corresponding to the fracture. The update unit is used to update the preset equivalent elastic stiffness matrix of organic shale without cracks by using the average value of the crack dip angle, the standard deviation of the crack dip angle, and the range of crack dip angle variation data, so as to obtain the target equivalent elastic stiffness matrix. as well as The processing unit is used to process the target equivalent elastic stiffness matrix based on the orthogonal medium anisotropic parameter formula to obtain the target orthogonal medium anisotropic parameters, so as to determine the anisotropic characterization information of the crack. The updating unit is further configured to determine the compliance change data based on the average crack inclination angle, the standard deviation of the crack inclination angle, and the range of crack inclination angle variation. as well as Based on the preset equivalent elastic stiffness matrix of organic-rich shale without cracks and the flexibility change data, the target equivalent elastic stiffness matrix is ​​determined. The compliance change data includes: Wherein, the θ max The maximum crack inclination angle, θ, represents the crack inclination angle. min The minimum crack inclination angle is characterized by θ, the crack inclination angle is characterized by θ0, the average crack inclination angle is characterized by σ, and the standard deviation of the crack inclination angle is characterized by N. θ The transformation matrix is ​​represented by ΔS, which represents the preset compliance change data, and G(θ) represents the crack dip angle distribution function. The crack dip angle distribution function includes: Wherein, θ represents the crack inclination angle, σ represents the standard deviation of the crack inclination angle, and θ0 represents the average value of the crack inclination angle; The determination of the target equivalent elastic stiffness matrix based on the preset equivalent elastic stiffness matrix of organic-rich shale without cracks and the flexibility change data includes: Inverting the preset equivalent elastic stiffness matrix of the organic-rich shale without cracks yields the preset flexibility matrix. The target compliance matrix is ​​obtained by summing the compliance matrix and the compliance change data; and Inverting the target flexibility matrix yields the target equivalent elastic stiffness matrix; The target equivalent elastic stiffness matrix includes: Wherein, I represents the identity matrix, and C b The preset equivalent elastic stiffness matrix of the organic-rich shale without cracks is used to characterize the density, and ΔS' is the data on the change in flexibility.

4. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method of any one of claims 1-2.

5. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, performs the method of any one of claims 1-2.

Images