Low-angle seam recognition method and device
By combining the resistivity intrusion correction difference ratio method and the three-porosity ratio method, the low-angle seams are comprehensively identified, which solves the problem of identifying low-angle seams in the existing technology, and accurately identifying and predicting low-angle seams, which improves the oil field development effect.
Patent Information
- Application Number
- CN202011072077.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2020-10-09
- Publication Date
- 2025-07-29
- Estimated Expiration
- 2040-10-09
AI Technical Summary
The problem of identifying and predicting low-angle seams in the prior art, especially in carbonate reservoirs, is mainly focused on the identification of high-angle tectonic seams, while relatively little research on low-angle seams has led to poor water injection and development results in oil field development.
Combined with the resistivity intrusion correction difference ratio method and the three-porosity ratio method, two methods are used to identify low-angle seams by obtaining resistivity data, porosity data and imaging logging data, and the resistivity intrusion correction difference ratio method is used to identify the crack inclination angle, and the three-porosity ratio method is used to identify the crack development intensity.
It has achieved good identification of low-angle seams, can accurately predict the development intensity of cracks, provides guidance for oil field development plans, and improves the oil field development effect.
Smart Images

Figure CN114325869B_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to geological exploration technology, and more specifically, to a method and device for identifying low-angle fractures. Background Art
[0002] Fractures are widely developed in carbonate reservoirs. Among the 35 large carbonate reservoirs discovered in the world so far, 12 have developed natural fractures, accounting for 34%. Fractures have an important impact on the oil-water seepage direction and the seepage velocity of underground fluids. Therefore, fracture prediction is one of the important topics in oilfield development research.
[0003] In the prior art, the research on fracture problems mainly focuses on the prediction of high-angle structural fractures (the angle of structural fractures is greater than 30°). The identification and prediction of high-angle structural fractures by predecessors can be summarized into 5 methods: (1) predicting fractures by logging and seismic information methods; (2) predicting fractures by numerical simulation methods; (3) predicting fractures by dynamic data methods; (4) predicting fractures by other qualitative methods (such as the distance from the fault and the structural position, etc.); (5) predicting fractures by non-linear theory methods. Among the above methods, the logging and seismic information method is the method with the best application effect in the mine field at present. The above methods mainly focus on high-angle structural fractures, and the research on low-angle fractures is relatively less. At present, the identification and prediction of low-angle fractures are still difficult problems. Summary of the Invention
[0004] In order to identify low-angle fractures and at least solve one defect in the prior art, the present invention provides a method for identifying low-angle fractures, including:
[0005] Obtaining resistivity data, porosity data, and imaging logging data of a target area;
[0006] Generating a resistivity fracture identification result according to the resistivity data;
[0007] Generating a porosity fracture identification result according to the porosity data;
[0008] Identifying low-angle fractures according to the resistivity fracture identification result, the porosity fracture identification result, and the imaging logging data.
[0009] In an embodiment of the present invention, the resistivity data includes: shallow lateral resistivity, deep lateral resistivity, and the true resistivity value of the formation corrected for invasion;
[0010] The porosity data includes: neutron porosity, density porosity, and acoustic porosity.
[0011] In an embodiment of the present invention, the generating a resistivity fracture identification result according to the resistivity data includes:
[0012] Generate a resistivity fracture identification result using the resistivity invasion correction difference ratio method based on the resistivity data described above.
[0013] In an embodiment of the present invention, the generating a porosity fracture identification result based on the porosity data described above includes:
[0014] Generate a porosity fracture identification result using the three-porosity identification method based on the porosity data described above.
[0015] In an embodiment of the present invention, the generating a resistivity fracture identification result based on the resistivity data described above includes:
[0016] Determine the deep and shallow dual laterolog difference ratio value according to the shallow laterolog resistivity value, the formation true resistivity value after invasion correction in the resistivity data described above, and the following formula;
[0017]
[0018] Generate a resistivity fracture identification result according to the determined deep and shallow dual laterolog difference ratio value;
[0019] Wherein, RTC is the deep and shallow dual laterolog difference ratio value, R lls is the shallow laterolog resistivity value, R t is the formation true resistivity value after invasion correction;
[0020] Wherein, R t = 2.589R lld - 1.589R lls , R lld is the deep laterolog resistivity value.
[0021] In an embodiment of the present invention, the generating a porosity fracture identification result based on the porosity data described above includes:
[0022] Determine the total porosity according to the neutron porosity, density porosity in the porosity data described above, and the following formula;
[0023]
[0024] Determine the secondary porosity according to the determined total porosity, acoustic porosity, and the following formula;
[0025]
[0026] Generate a porosity fracture identification result according to the determined secondary porosity;
[0027] Wherein, Φ N is the neutron porosity, Φ D is the density porosity, Φ S is the acoustic porosity, Φ T is the total porosity;
[0028] wherein, R p is the secondary porosity.
[0029] In an embodiment of the present invention, the identification of low-angle fractures based on the resistivity fracture identification result, the porosity fracture identification result, and the imaging logging data includes:
[0030] Comparing the resistivity fracture identification result and the porosity fracture identification result with the imaging logging data respectively to determine whether the identification rate of low-angle fractures by the resistivity fracture identification result is higher than a preset threshold;
[0031] If it is determined that the identification rate of low-angle fractures by the resistivity fracture identification result is higher than the preset threshold, then the porosity fracture identification result corresponding to the well section where the ratio of deep and shallow dual laterolog differences is less than zero is used as the low-angle fracture identification result;
[0032] If it is determined that the identification rate of low-angle fractures by the resistivity fracture identification result is not higher than the preset threshold, then the porosity fracture identification result of the well section where the ratio of deep and shallow dual laterolog differences is greater than or equal to zero is set to zero, and a part of the porosity fracture identification result after being set to zero is used as the low-angle fracture identification result.
[0033] In an embodiment of the present invention, the preset threshold is 70%. Meanwhile, the present invention also provides a low-angle fracture identification device, including:
[0034] A data acquisition module, configured to acquire resistivity data, porosity data, and imaging logging data of a target area;
[0035] A resistivity fracture identification module, configured to generate a resistivity fracture identification result according to the resistivity data;
[0036] A porosity fracture identification module, configured to generate a porosity fracture identification result according to the porosity data;
[0037] A low-angle fracture identification module, configured to identify low-angle fractures according to the resistivity fracture identification result, the porosity fracture identification result, and the imaging logging data.
[0038] In an embodiment of the present invention, the resistivity fracture identification module generating a resistivity fracture identification result according to the resistivity data includes:
[0039] Generating a resistivity fracture identification result according to the resistivity data by using the resistivity invasion correction difference ratio method.
[0040] In an embodiment of the present invention, the porosity fracture identification module generating a porosity fracture identification result according to the porosity data includes:
[0041] Generate a porosity fracture identification result using the three-porosity identification method based on the porosity data described above.
[0042] In an embodiment of the present invention, the resistivity fracture identification module includes:
[0043] A differential ratio determination unit for generating a determination of the deep and shallow dual lateral differential ratio according to the shallow lateral resistivity value, the invaded-corrected formation true resistivity value in the resistivity data, and the following formula;
[0044]
[0045] A resistivity fracture identification unit for generating a resistivity fracture identification result according to the determined deep and shallow dual lateral differential ratio;
[0046] Wherein, RTC is the deep and shallow dual lateral differential ratio, R lls is the shallow lateral resistivity value, and R t is the invaded-corrected formation true resistivity value;
[0047] Wherein, R t = 2.589R lld - 1.589R lls and R lld is the deep lateral resistivity value.
[0048] In an embodiment of the present invention, the porosity fracture identification module includes:
[0049] A total porosity determination unit for determining the total porosity according to the neutron porosity and density porosity in the porosity data and the following formula;
[0050]
[0051] Secondary porosity, for determining the secondary porosity according to the determined total porosity, acoustic porosity, and the following formula;
[0052]
[0053] A pore fracture identification unit for generating a porosity fracture identification result according to the determined secondary porosity;
[0054] Wherein, Φ N is the neutron porosity, Φ D is the density porosity, Φ S is the acoustic porosity, and Φ T is the total porosity;
[0055] Wherein, R p is the secondary porosity.
[0056] In an embodiment of the present invention, the low-angle fracture identification module includes:
[0057] A comparison unit for comparing the resistivity fracture identification result and the porosity fracture identification result with the imaging logging data respectively to determine whether the identification rate of low-angle fractures in the resistivity fracture identification result is higher than a preset threshold;
[0058] An identification unit for determining that the identification rate of low-angle fractures in the resistivity fracture identification result is higher than the preset threshold, and taking the porosity fracture identification result corresponding to the well section where the deep and shallow dual laterolog difference ratio is less than zero as the low-angle fracture identification result;
[0059] Determining that the identification rate of low-angle fractures in the resistivity fracture identification result is not higher than the preset threshold, zeroing the porosity fracture identification result of the well section where the deep and shallow dual laterolog difference ratio is greater than or equal to zero, and taking part of the zeroed porosity fracture identification result as the low-angle fracture identification result.
[0060] Meanwhile, the present invention also provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the computer program, the above method is implemented.
[0061] Meanwhile, the present invention also provides a computer-readable storage medium storing a computer program for executing the above method.
[0062] In the low-angle fracture identification scheme provided by the present invention, through processing, the fractures identified by the three-porosity ratio method are only low-angle fractures. Through this method, the goal of comprehensively applying the resistivity invasion correction difference ratio method and the three-porosity ratio method to identify low-angle fractures is achieved, giving full play to the advantage of the resistivity invasion correction difference ratio method in being able to identify the fracture dip angle. At the same time, the advantages of the three-porosity ratio method in having a good identification effect on low-angle fractures and being able to predict the fracture development intensity are exerted, realizing a better identification of low-angle fractures, providing guidance for subsequent oilfield development plan deployment and oilfield water injection development.
[0063] To make the above and other objects, features, and advantages of the present invention more obvious and understandable, the following specific preferred embodiments are given in conjunction with the accompanying drawings and are described in detail as follows. BRIEF DESCRIPTION OF THE DRAWINGS
[0064] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for the description of the embodiments or the prior art. Obviously, the following drawings are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.
[0065] Figure 1Flow chart of the low-angle fracture identification method provided by the present invention;
[0066] Figure 2 Block diagram of the low-angle fracture identification device provided by the present invention;
[0067] Figure 3 Schematic diagram in the embodiment of the present invention;
[0068] Figure 4 Schematic diagram in the embodiment of the present invention;
[0069] Figure 5 Schematic diagram in the embodiment of the present invention;
[0070] Figure 6 Schematic diagram of the electronic device provided in the embodiment of the present invention. Detailed implementation manners
[0071] Next, the technical solutions in the embodiments of the present invention will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.
[0072] Taking the logging and seismic information method as an example, due to the limitation of the resolution of seismic data, it is difficult to effectively distinguish the formation from the low-angle fractures. Logging is currently a method that can meet the accuracy requirements for low-angle fracture prediction. Based on conventional logging data, methods for fracture identification such as the resistivity invasion correction difference ratio method based on deep and shallow dual laterologs and the three-porosity ratio method based on density logging, neutron logging, and acoustic logging have been mainly proposed. The resistivity invasion correction difference ratio method proposed based on the conventional logging deep and shallow dual laterolog curves is the only conventional logging method in the prior art that can identify the fracture dip angle. However, due to the influence of the fluid properties in the fractures and the filling degree of the fractures, the reliability of low-angle fracture prediction is currently poor, and the prediction effect for high-angle fractures and oblique fractures (>30°) is better. In addition, the resistivity invasion correction difference ratio method can only identify whether fractures are developed and is difficult to predict the development intensity of fractures. The three-porosity ratio method can better identify fractures with different dip angles, especially for low-angle fracture prediction.
[0073] The present invention combines the resistivity invasion correction difference ratio method and the three-porosity ratio method to provide a low-angle fracture identification method, as Figure 1 shown, including:
[0074] Step S101, obtaining resistivity data, porosity data, and imaging logging data of the target area;
[0075] Step S102: Generate a resistivity fracture identification result based on the resistivity data;
[0076] Step S103: Generate a porosity fracture identification result based on the porosity data;
[0077] Step S104: Conduct low-angle fracture identification based on the resistivity fracture identification result, porosity fracture identification result, and the imaging logging data.
[0078] In the prior art, there is a phenomenon of a large number of low-angle fractures developing in some oilfields. Therefore, it is necessary to accurately predict the distribution of such fractures to provide guidance for the deployment of oilfield development plans. In the present invention, resistivity invasion correction difference ratio method and triple porosity ratio method are used to identify fractures respectively. Then, the advantage of the resistivity invasion correction difference ratio method in identifying fracture dip angles is utilized. Since in this example, the resistivity invasion correction difference ratio method has a poor effect on identifying low-angle fractures (the identification rate is lower than 70%), the high-angle fractures and oblique fractures identified by the resistivity invasion correction difference ratio method are used as the basis for judgment to screen out the high-angle fractures and oblique fractures identified by the triple porosity ratio method. Then, operations are performed on these identified high-angle fractures and oblique fractures (such as using the existing Petrel software), to eliminate the high-angle fractures and oblique fractures identified on the logging curves of the triple porosity ratio method, so that the fractures identified by the triple porosity ratio method are only low-angle fractures. Through this method, the goal of comprehensively applying the resistivity invasion correction difference ratio method and the triple porosity ratio method to identify low-angle fractures is achieved, giving full play to the advantage of the resistivity invasion correction difference ratio method in being able to identify fracture dip angles, and at the same time giving full play to the advantage of the triple porosity ratio method in having a good effect on identifying low-angle fractures and being able to predict the fracture development intensity, achieving a better identification of low-angle fractures. It provides guidance for subsequent oilfield development plan adjustment and oilfield water injection development.
[0079] In the embodiment of the present invention, the resistivity data includes: shallow lateral resistivity, deep lateral resistivity;
[0080] In this embodiment, based on the obtained resistivity data, the resistivity invasion correction difference ratio method is used to identify fractures with different dip angles in the study area. The steps for the resistivity invasion correction difference ratio method to achieve fracture prediction are as follows:
[0081]
[0082] In the formula, RTC is the difference ratio of the deep and shallow dual laterals;
[0083] R lls is the shallow lateral resistivity value;
[0084] R t is the true formation resistivity value after invasion correction, and its calculation formula is R t = 2.589R lld - 1.589Rlls ;
[0085] Wherein, R lld is the deep lateral resistivity value.
[0086] When RTC > 0, the dual laterolog curve shows a positive anomaly, indicating a high-angle fracture; when RTC = 0, it is an oblique fracture; when RTC < 0, the dual laterolog curve shows a negative anomaly, indicating a low-angle fracture.
[0087] The porosity data described above includes: neutron porosity, density porosity, and acoustic porosity;
[0088] In this embodiment, based on the obtained porosity data, fractures with different dip angles in the study area are identified using the three-porosity ratio method. The steps for identifying fractures in the study area using the three-porosity ratio method are as follows:
[0089]
[0090]
[0091] In the formula, Φ N is the neutron porosity;
[0092] Φ D is the density porosity;
[0093] Φ S is the acoustic porosity;
[0094] Φ T is the total porosity.
[0095] Among them, neutron logging and density logging represent the magnitude of the total porosity, and acoustic logging mainly reflects the porosity of the rock matrix. Therefore, constructing R p mainly reflects the magnitude of the secondary porosity.
[0096] When R p > 0, it indicates the development of secondary reservoir spaces, such as fractures, dissolution pores, and microfractures. The larger R p , the more developed the fractures are.
[0097] In the embodiment of the present invention, the identification of low-angle fractures based on the resistivity fracture identification result, porosity fracture identification result, and the imaging logging data includes:
[0098] Comparing the resistivity fracture identification result and the porosity fracture identification result with the imaging logging data respectively to determine whether the identification rate of the resistivity fracture identification result for low-angle fractures is higher than a preset threshold;
[0099] If it is determined that the recognition rate of the resistivity fracture recognition result for low-angle fractures is higher than the preset threshold, then the porosity fracture recognition result corresponding to the well section with a deep-shallow dual laterolog difference ratio less than zero is used as the low-angle fracture recognition result;
[0100] If it is determined that the recognition rate of the resistivity fracture recognition result for low-angle fractures is not higher than the preset threshold, then the porosity fracture recognition result of the well section with a deep-shallow dual laterolog difference ratio greater than or equal to zero is reset to zero, and a part of the reset porosity fracture recognition result is used as the low-angle fracture recognition result.
[0101] Specifically, if it is determined that the recognition rate of the resistivity fracture recognition result for low-angle fractures is higher than the preset threshold, that is, if the recognition rate of the deep-shallow dual laterolog difference ratio for low-angle fractures is high (determined by a preset threshold in this embodiment, the threshold is 70%), then the recognition is carried out using the deep-shallow dual laterolog difference ratio less than zero, and the porosity fracture recognition result corresponding to the well section with a deep-shallow dual laterolog difference ratio less than zero is used as the low-angle fracture recognition result;
[0102] If the recognition rate of the deep-shallow dual laterolog difference ratio for low-angle fractures is low (less than 70%), then it is determined to use the deep-shallow dual laterolog difference ratio greater than or equal to zero for fracture dip determination, and the triple porosity method curve values of the well section with a deep-shallow dual laterolog difference ratio greater than or equal to zero are set to zero.
[0103] Specifically, in the embodiment of the present invention, by comparing the fracture results identified by two methods, namely the resistivity invasion correction difference ratio method and the triple porosity ratio method, with the imaging logging, in the embodiment of the present invention, the fracture recognition results of the conventional logging are statistically analyzed respectively with a fracture dip of 30° as the boundary, that is, the recognition results of high-angle fractures, oblique fractures and low-angle fractures in the study area by the resistivity invasion correction difference ratio method and the triple porosity ratio method.
[0104] If the resistivity invasion correction difference ratio method has a good recognition effect on low-angle fractures, then the result identified by the triple porosity ratio method is determined using RTC < 0. If RTC is less than zero, it is considered that the triple porosity ratio method identifies low-angle fractures and the result of the triple porosity ratio method remains unchanged; if the condition is not met, it means that the fractures identified by the triple porosity ratio method are high-angle fractures, and the result identified by the triple porosity ratio method is set to 0. Thus, the goal of identifying and predicting the development intensity of low-angle fractures based on two technologies is achieved.
[0105] If the resistivity invasion correction difference ratio method has poor recognition effect on low-angle fractures, then use RTC≥0 to judge the results identified by the triple porosity ratio method. If RTC is greater than or equal to zero, it is considered that the triple porosity ratio method identifies high-angle fractures or oblique fractures, and the results of the triple porosity ratio method are classified as zero; if the condition that RTC is greater than or equal to zero is not satisfied, it means that the fractures identified by the triple porosity ratio method are low-angle fractures or there are no fractures, then the results identified by the triple porosity ratio method remain unchanged. Thus, the goal of identifying and predicting the development intensity of low-angle fractures based on the two technologies is achieved.
[0106] The resistivity invasion correction difference ratio method has extremely poor recognition effect on low-angle fractures in some areas, which is mainly related to the influence of the fluid properties in the fractures and the filling degree of the fractures. In the embodiments of the present invention, for the extremely poor recognition effect of the resistivity invasion correction difference ratio method on low-angle fractures in some areas, another idea can be adopted, that is, use RTC≥0 to judge the results identified by the triple porosity ratio method:
[0107] If RTC is greater than or equal to zero, it is considered that the triple porosity ratio method identifies high-angle fractures, and the results identified by the triple porosity ratio method are classified as 0. If the condition is not satisfied, it means that the fractures identified by the triple porosity ratio method are low-angle fractures, and the results of the triple porosity ratio method remain unchanged, and the goal of identifying and predicting the development intensity of low-angle fractures based on the two technologies is achieved.
[0108] At present, there is relatively little research on the prediction of low-angle fracture problems. Due to the unclear understanding of the fracture distribution, the injected water breaks through quickly after the oilfield is developed by water injection, and the oil wells are extremely prone to water flooding and shut-in or early conversion to water injection, resulting in extremely poor oilfield development effects. The present invention comprehensively uses the resistivity invasion correction difference ratio method and the triple porosity ratio method, gives play to the advantage of the resistivity invasion correction difference ratio method in judging the fracture dip angle, and at the same time utilizes the advantages of the triple porosity ratio method in high recognition rate of low-angle fractures and ability to predict fracture strength. Thus, not only can low-angle fractures be identified with a high recognition rate, but also the development intensity of low-angle fractures can be predicted, laying a foundation for the prediction of low-angle fractures between wells. Through the effective identification and prediction of low-angle fractures, it provides guidance for the deployment of oilfield development plans and the optimization of well patterns.
[0109] Meanwhile, the present invention also provides a low-angle fracture identification device, as Figure 2 shown, including:
[0110] A data acquisition module 201 for acquiring resistivity data, porosity data and imaging logging data of a target area;
[0111] A resistivity fracture identification module 202 for generating a resistivity fracture identification result according to the resistivity data;
[0112] The porosity fracture identification module 203 is configured to generate a porosity fracture identification result based on the porosity data;
[0113] The low-angle fracture identification module 204 is configured to identify low-angle fractures based on the resistivity fracture identification result, the porosity fracture identification result, and the imaging logging data.
[0114] In an embodiment of the present invention, the low-angle fracture identification module includes:
[0115] A comparison unit configured to compare the resistivity fracture identification result and the porosity fracture identification result with the imaging logging data respectively to determine whether the identification rate of the resistivity fracture identification result for low-angle fractures is higher than a preset threshold;
[0116] That is, to determine the quality of the identification effect of the resistivity fracture identification result for low-angle fractures, so as to judge whether to use the deep-shallow dual laterolog difference ratio less than zero or the deep-shallow dual laterolog difference ratio greater than or equal to zero to judge the fracture dip angle;
[0117] An identification unit configured to, if it is determined that the identification rate of the resistivity fracture identification result for low-angle fractures is higher than the preset threshold, use the porosity fracture identification result corresponding to the well section where the deep-shallow dual laterolog difference ratio is less than zero as the low-angle fracture identification result;
[0118] If it is determined that the identification rate of the resistivity fracture identification result for low-angle fractures is not higher than the preset threshold, zero the porosity fracture identification result corresponding to the well section where the deep-shallow dual laterolog difference ratio is greater than or equal to zero, and use a part of the zeroed porosity fracture identification result as the low-angle fracture identification result.
[0119] That is, if the resistivity invasion correction difference ratio method has a good identification effect on low-angle fractures, then use the deep-shallow dual laterolog difference ratio less than zero, and the result of the triple porosity ratio method remains unchanged. In the well sections where the deep-shallow dual laterolog difference ratio is not less than zero, zero the result of the triple porosity ratio method, and use the processed porosity fracture identification result as the low-angle fracture identification;
[0120] If the resistivity invasion correction difference ratio method has a poor identification effect on low-angle fractures, then use RTC≥0 to judge the result identified by the triple porosity ratio method. If RTC is greater than or equal to zero, it is considered that the triple porosity ratio method identifies high-angle fractures or oblique fractures, and zero the result of the triple porosity ratio method; if the condition that RTC is greater than or equal to zero is not met, it means that the fractures identified by the triple porosity ratio method are low-angle fractures or there are no fractures, then the result identified by the triple porosity ratio method remains unchanged.
[0121] For those skilled in the art, according to the description of the foregoing embodiments, the specific implementation manners of the low-angle fracture identification device of the present invention can be clearly obtained, and will not be elaborated herein.
[0122] The present invention borrows the identification results of the resistivity invasion correction difference ratio method and the triple porosity ratio method for fractures. In the example of the present invention, based on the identification result of the high-angle fractures by the resistivity invasion correction difference ratio method, the high-angle fractures and the oblique fractures identified by the triple porosity ratio method are eliminated, achieving the goal of identifying low-angle fractures in the study area by the triple porosity ratio method.
[0123] The present invention successfully identifies high-angle fractures and oblique fractures by using the resistivity invasion correction difference ratio method, thereby eliminating the high-angle fractures and oblique fractures among all the fractures identified by the triple porosity ratio method, achieving the goal of identifying low-angle fractures by the triple porosity ratio method and reflecting the development intensity of low-angle fractures, giving play to the advantage of the high identification rate of the triple porosity ratio method for low-angle fractures, and thus realizing the identification and prediction of low-angle fractures with a high identification rate.
[0124] Figure 3 It is a schematic flow chart of this embodiment. Figure 4 It is a template diagram for comprehensively identifying low-angle fractures by the resistivity invasion method and the triple porosity ratio method in an embodiment of the present invention.
[0125] Figure 5 It is a schematic diagram of the result of predicting low-angle fractures in a certain area by using the technical solution of the present invention.
[0126] This embodiment also provides an electronic device, which can be a desktop computer, a tablet computer, a mobile terminal, etc., and this embodiment is not limited thereto. In this embodiment, the electronic device can refer to the embodiments of the foregoing method and device, the content of which is incorporated herein, and the repeated parts will not be described again.
[0127] Figure 6 It is a schematic block diagram of the system composition of the electronic device 600 according to an embodiment of the present invention. As Figure 6 shown, the electronic device 600 may include a central processing unit 100 and a memory 140; the memory 140 is coupled to the central processing unit 100. It should be noted that this figure is exemplary; other types of structures can also be used to supplement or replace this structure to achieve telecommunication functions or other functions.
[0128] In one embodiment, the low-angle fracture identification function can be integrated into the central processing unit 100. Among them, the central processing unit 100 can be configured to perform the following controls:
[0129] Obtain resistivity data, porosity data, and imaging logging data of the target area;
[0130] Generate a resistivity fracture identification result according to the resistivity data;
[0131] Generate a porosity fracture identification result according to the porosity data;
[0132] Identify low-angle fractures based on the resistivity fracture identification results, porosity fracture identification results, and the imaging logging data.
[0133] In an embodiment of the present invention, the resistivity data includes: shallow lateral resistivity, deep lateral resistivity, and the true formation resistivity value after invasion correction;
[0134] The porosity data includes: neutron porosity, density porosity, and acoustic porosity.
[0135] In an embodiment of the present invention, generating the resistivity fracture identification result according to the resistivity data includes:
[0136] Generate the resistivity fracture identification result using the resistivity invasion correction difference ratio method based on the resistivity data.
[0137] In an embodiment of the present invention, generating the porosity fracture identification result according to the porosity data includes:
[0138] Generate the porosity fracture identification result using the three-porosity identification method based on the porosity data.
[0139] In an embodiment of the present invention, generating the resistivity fracture identification result according to the resistivity data includes:
[0140] Determine the deep and shallow dual lateral difference ratio according to the shallow lateral resistivity value and the true formation resistivity value after invasion correction in the resistivity data and the following formula;
[0141]
[0142] Generate the resistivity fracture identification result according to the determined deep and shallow dual lateral difference ratio;
[0143] where RTC is the deep and shallow dual lateral difference ratio, R lls is the shallow lateral resistivity value, and R t is the true formation resistivity value after invasion correction;
[0144] where R t = 2.589R lld - 1.589R lls and R lld is the deep lateral resistivity value.
[0145] In an embodiment of the present invention, generating the porosity fracture identification result according to the porosity data includes:
[0146] Determine the total porosity according to the neutron porosity, density porosity in the porosity data and the following formula;
[0147]
[0148] Determine the secondary porosity based on the determined total porosity, acoustic porosity and the following formula;
[0149]
[0150] Generate a porosity fracture identification result based on the determined secondary porosity;
[0151] Wherein, Φ N is the neutron porosity, Φ D is the density porosity, Φ S is the acoustic porosity, Φ T is the total porosity;
[0152] Wherein, R p is the secondary porosity.
[0153] In an embodiment of the present invention, the identification of low-angle fractures based on the resistivity fracture identification result, the porosity fracture identification result and the imaging logging data includes:
[0154] Compare the resistivity fracture identification result and the porosity fracture identification result with the imaging logging data respectively to determine whether the identification rate of the resistivity fracture identification result for low-angle fractures is higher than a preset threshold;
[0155] If it is determined that the identification rate of the resistivity fracture identification result for low-angle fractures is higher than the preset threshold, then use the porosity fracture identification result corresponding to the well section where the ratio of deep and shallow dual laterolog differences is less than zero as the low-angle fracture identification result;
[0156] If it is determined that the identification rate of the resistivity fracture identification result for low-angle fractures is not higher than the preset threshold, then set to zero the porosity fracture identification result corresponding to the well section where the ratio of deep and shallow dual laterolog differences is greater than or equal to zero, and use a part of the porosity fracture identification result after zeroing as the low-angle fracture identification result.
[0157] In an embodiment of the present invention, the preset threshold is 70%.
[0158] As Figure 6 shown, the electronic device 600 may further include: a communication module 110, an input unit 120, an audio processing unit 130, a display 160, and a power supply 170. It should be noted that the electronic device 600 does not necessarily have to include all the components shown in Figure 6 ; in addition, the electronic device 600 may further include components not shown in Figure 6 , and reference may be made to the prior art.
[0159] As Figure 6As shown, the central processing unit 100, sometimes also referred to as a controller or operation control, may include a microprocessor or other processor device and / or logic device. The central processing unit 100 receives inputs and controls the operations of the various components of the electronic device 600.
[0160] Among them, the memory 140 can be, for example, one or more of a buffer, a flash memory, a hard drive, a removable medium, a volatile memory, a non-volatile memory, or other suitable devices. It can store the above information related to failures, and can also store programs for executing relevant information. And the central processing unit 100 can execute the programs stored in the memory 140 to achieve information storage or processing, etc.
[0161] The input unit 120 provides inputs to the central processing unit 100. The input unit 120 is, for example, a key or a touch input device. The power supply 170 is used to supply power to the electronic device 600. The display 160 is used to display display objects such as images and texts. The display can be, for example, an LCD display, but is not limited thereto.
[0162] The memory 140 can be a solid-state memory. For example, it can be a read-only memory (ROM), a random access memory (RAM), a SIM card, etc. It can also be a memory that stores information even when powered off, can be selectively erased, and has more data. Examples of such a memory are sometimes referred to as EPROMs, etc. The memory 140 can also be some other type of device. The memory 140 includes a buffer memory 141 (sometimes referred to as a buffer). The memory 140 can include an application / function storage unit 142, which is used to store application programs and function programs or the processes for operating the electronic device 600 through the central processing unit 100.
[0163] The memory 140 can also include a data storage unit 143, which is used to store data, such as contacts, digital data, pictures, sounds, and / or any other data used by the electronic device. The driver storage unit 144 of the memory 140 can include various drivers for the communication functions of the electronic device and / or for executing other functions of the electronic device (such as a messaging application, an address book application, etc.).
[0164] The communication module 110 is a transmitter / receiver 110 that transmits and receives signals via the antenna 111. The communication module (transmitter / receiver) 110 is coupled to the central processing unit 100 to provide input signals and receive output signals, which can be the same as in the case of a conventional mobile communication terminal.
[0165] Based on different communication technologies, in the same electronic device, multiple communication modules 110 can be provided, such as a cellular network module, a Bluetooth module, and / or a wireless local area network module, etc. The communication module (transmitter / receiver) 110 is also coupled to a speaker 131 and a microphone 132 via an audio processor 130 to provide an audio output via the speaker 131 and receive an audio input from the microphone 132, so as to implement normal telecommunication functions. The audio processor 130 may include any suitable buffers, decoders, amplifiers, etc. Additionally, the audio processor 130 is also coupled to a central processor 100, so that recording can be performed on the local machine through the microphone 132, and the sound stored on the local machine can be played through the speaker 131.
[0166] An embodiment of the present invention also provides a computer-readable program, wherein when the program is executed in an electronic device, the program causes the computer to execute the low-angle seam recognition method as described in the above embodiment in the electronic device.
[0167] An embodiment of the present invention also provides a storage medium storing a computer-readable program, wherein the computer-readable program causes the computer to execute the low-angle seam recognition described in the above embodiment in an electronic device.
[0168] The preferred embodiments of the present invention have been described above with reference to the accompanying drawings. Many features and advantages of these embodiments are clear from this detailed description, and thus the appended claims are intended to cover all such features and advantages that fall within the true spirit and scope of these embodiments. In addition, since many modifications and changes are readily envisioned by those skilled in the art, the embodiments of the present invention are not to be limited to the exact structures and operations illustrated and described, but may cover all suitable modifications and equivalents falling within their scope.
[0169] Those skilled in the art should understand that the embodiments of the present invention may be provided as a method, a system, or a computer program product. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.
[0170] The present invention is described with reference to the flowcharts and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and combinations 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 the processors of general purpose computers, special purpose computers, embedded processors, or other programmable data processing devices to produce a machine, such that the instructions executed by the processors of the computer or other programmable data processing devices produce means for implementing the functions specified in the flow Figure 1 one or more flows and / or blocks Figure 1 or means for implementing the functions specified in a block or multiple blocks.
[0171] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, such that the instructions stored in the computer-readable memory produce a manufactured article including instruction means for implementing the functions specified in the flow Figure 1 one or more flows and / or blocks Figure 1 or a block or multiple blocks.
[0172] These computer program instructions can also be loaded onto a computer or other programmable data processing device, such that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, and thus the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in the flow Figure 1 one or more flows and / or blocks Figure 1 or a block or multiple blocks.
[0173] Specific embodiments are used in the present invention to elaborate on the principles and implementation manners of the present invention. The descriptions of the above embodiments are only for helping to understand the method and its core idea of the present invention; at the same time, for those of ordinary skill in the art, according to the idea of the present invention, there will be changes in the specific implementation manners and application scopes. In summary, the content of this specification should not be construed as a limitation to the present invention.
Claims
1. A method for identifying low-angle seams, characterized in that, The method described above includes: Obtaining resistivity data, porosity data, and imaging logging data of the target area; Generating a resistivity fracture identification result based on the resistivity data; Generating a porosity fracture identification result based on the porosity data; Performing low-angle fracture identification based on the resistivity fracture identification result, porosity fracture identification result, and the imaging logging data; The performing low-angle fracture identification based on the resistivity fracture identification result, porosity fracture identification result, and the imaging logging data includes: Comparing the resistivity fracture identification result and the porosity fracture identification result with the imaging logging data respectively to determine whether the identification rate of the resistivity fracture identification result for low-angle fractures is higher than a preset threshold; If it is determined that the identification rate of the resistivity fracture identification result for low-angle fractures is higher than the preset threshold, then the porosity fracture identification result corresponding to the well section where the calculated ratio of deep and shallow dual laterolog differences is less than zero is used as the low-angle fracture identification result. For the well sections that do not meet the condition of being less than zero, all the porosity fracture identification curves are set to zero, where the ratio of deep and shallow dual laterolog differences is calculated based on the resistivity data; If it is determined that the identification rate of the resistivity fracture identification result for low-angle fractures is not higher than the preset threshold, then the porosity fracture identification result of the well section where the ratio of deep and shallow dual laterolog differences is greater than or equal to zero is set to zero, and the porosity fracture identification result after the zeroing process is used as the low-angle fracture identification result.
2. The low-angle seam recognition method according to claim 1, wherein The resistivity data described above includes: shallow lateral resistivity, deep lateral resistivity, and the true formation resistivity value after invasion correction; The porosity data described above includes: neutron porosity, density porosity, and acoustic porosity.
3. The low-angle seam recognition method according to claim 2, wherein, The generating a resistivity fracture identification result based on the resistivity data includes: Generating a resistivity fracture identification result using the resistivity invasion correction difference ratio method based on the resistivity data; The generating a porosity fracture identification result based on the porosity data includes: Generating a porosity fracture identification result using the three-porosity identification method based on the porosity data.
4. The low-angle seam recognition method according to claim 3, wherein, The generating a resistivity fracture identification result based on the resistivity data includes: Generating a determined ratio of deep and shallow dual laterolog differences based on the shallow lateral resistivity value and the true formation resistivity value after invasion correction in the resistivity data and the following formula; Generating a resistivity fracture identification result based on the determined ratio of deep and shallow dual laterolog differences; Among them, RTC is the ratio of the deep and shallow dual laterolog differences, R lls is the shallow laterolog resistivity value, R t is the true formation resistivity value after invasion correction; Wherein, R t = 2.589R lld - 1.589R lls and R lld is the deep lateral resistivity value.
5. The low-angle seam recognition method according to claim 4, wherein The generating a porosity fracture identification result based on the porosity data includes: Determining the total porosity based on the neutron porosity and density porosity in the porosity data and the following formula; Determining the secondary porosity based on the determined total porosity, acoustic porosity, and the following formula; Generating a porosity fracture identification result based on the determined secondary porosity; Among them, Φ N is the neutron porosity, Φ D is the density porosity, θ S is the acoustic porosity, θ T is the total porosity; Among them, R p is the secondary porosity.
6. The low-angle seam recognition method according to claim 1, characterized in that, The preset threshold is 70%.
7. A low-angle seam recognition device, characterized in that, The device described above includes: A data acquisition module for obtaining resistivity data, porosity data, and imaging logging data of the target area; A resistivity fracture identification module for generating a resistivity fracture identification result based on the resistivity data; A porosity fracture identification module for generating a porosity fracture identification result based on the porosity data; A low-angle fracture identification module, configured to identify low-angle fractures based on the resistivity fracture identification result, the porosity fracture identification result, and the imaging logging data; The low-angle fracture identification module includes: A comparison unit, configured to compare the resistivity fracture identification result and the porosity fracture identification result with the imaging logging data respectively to determine whether the identification rate of the resistivity fracture identification result for low-angle fractures is higher than a preset threshold; An identification unit, configured to determine that if the identification rate of the resistivity fracture identification result for low-angle fractures is higher than the preset threshold, then use the porosity fracture identification result corresponding to the well section where the ratio of deep and shallow laterolog differences is less than zero as the low-angle fracture identification result. For the well sections that do not meet the condition of being less than zero, all the porosity fracture identification curves are set to zero, where the ratio of deep and shallow laterolog differences is calculated based on resistivity data; If it is determined that the identification rate of the resistivity fracture identification result for low-angle fractures is not higher than the preset threshold, then zero out the porosity fracture identification result of the well section where the ratio of deep and shallow laterolog differences is greater than or equal to zero, and use the zeroed porosity fracture identification result as the low-angle fracture identification result.
8. The low-angle fracture identification device according to claim 7, wherein The resistivity data includes: shallow lateral resistivity, deep lateral resistivity, and the true formation resistivity value after invasion correction; The porosity data includes: neutron porosity, density porosity, and acoustic porosity.
9. The low-angle seam recognition device according to claim 8, wherein The resistivity fracture identification module generates a resistivity fracture identification result based on the resistivity data, including: Generating a resistivity fracture identification result using the resistivity invasion correction difference ratio method based on the resistivity data; The porosity fracture identification module generates a porosity fracture identification result based on the porosity data, including: Generating a porosity fracture identification result using the three-porosity identification method based on the porosity data.
10. The low-angle seam recognition device according to claim 9, characterized in that, The resistivity fracture identification module includes: A difference ratio determination unit, configured to generate and determine the ratio of deep and shallow laterolog differences according to the shallow lateral resistivity value and the true formation resistivity value after invasion correction in the resistivity data and the following formula; A resistivity fracture identification unit, configured to generate a resistivity fracture identification result based on the determined ratio of deep and shallow laterolog differences; wherein, RTC is the ratio of the deep and shallow laterolog differences, R lls is the shallow laterolog resistivity value, and R t is the true formation resistivity value after invasion correction; Wherein, R t = 2.589R lld - 1.589R lls and R lld is the deep lateral resistivity value.
11. The low-angle seam recognition device according to claim 9, characterized in that, The porosity fracture identification module includes: A total porosity determination unit, configured to determine the total porosity according to the neutron porosity and density porosity in the porosity data and the following formula; Secondary porosity, configured to determine the secondary porosity according to the determined total porosity, acoustic porosity, and the following formula; A pore fracture identification unit, configured to generate a porosity fracture identification result based on the determined secondary porosity; Among them, θ N is the neutron porosity, θ D is the density porosity, Φ S is the acoustic porosity, Φ T is the total porosity; Among them, R p is the secondary porosity.
12. The low-angle seam recognition device according to claim 7, wherein The preset threshold is 70%.
13. A computer device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method according to any one of claims 1 to 6.
14. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program for executing the method according to any one of claims 1 to 6.