Relative sea level change analysis method, device, electronic device and storage medium

Abstract

The present invention relates to a relative sea level change analysis method, device, electronic device, and storage medium. The method comprises: obtaining well logging data of microresistivity imaging of sedimentary rock formations, preprocessing the well logging data to obtain a dynamic well logging image; obtaining grayscale values ​​of the dynamic well logging image, and calculating a well logging lithology curve of the microresistivity imaging based on the grayscale values; calculating high-order derivatives of the well logging lithology curve based on the well logging lithology curve, and determining the sedimentary rock formation depth and the cyclic layer interface depth based on the high-order derivatives; calculating the cyclic layer thickness and the average cyclic layer thickness based on the cyclic layer interface depth; calculating the cyclic layer offset thickness and the cyclic layer cumulative offset thickness based on the cyclic layer thickness and the cyclic layer average thickness; and analyzing the relative sea level change of the sedimentary rock formation based on the trend of the cumulative offset thickness of the cyclic layer. The present invention improves the accuracy and effectiveness of research on relative sea level change of sedimentary rock formations.

Description

Technical Field

[0001] The present invention relates to the technical field of marine geology and sequence stratigraphy research, and in particular to a relative sea level change analysis method, device, electronic equipment and storage medium. Background Art

[0002] Marine sedimentary rock formations record the changes and evolutionary characteristics of the marine environment during their deposition. Using sedimentary rock data to reconstruct relative sea level changes has always been a hot topic and a challenge in marine geology and sequence stratigraphy. Relative sea level change research primarily relies on data such as seismic data, field outcrops, drill cores, paleontology, and conventional well logging. However, seismic data has low resolution; outcrop rock data is difficult to obtain due to poor field conditions; drilling and coring are expensive and difficult to cover the entire target layer; paleontological data relies on extensive analytical testing, which is time-consuming and labor-intensive; and conventional well logging data are one-dimensional curves, which can be multi-solution prone. These data deficiencies significantly impact the precision and accuracy of relative sea level change results.

[0003] Therefore, how to provide a new relative sea level change analysis method, device, electronic equipment and storage medium to improve the efficiency and accuracy of relative sea level change research is a technical problem that needs to be solved urgently. Summary of the Invention

[0004] In view of this, it is necessary to provide a relative sea level change analysis method, device, electronic device and storage medium to solve the technical problems of high cost, low efficiency and low accuracy in the existing technology of relative sea level change research.

[0005] In order to solve the above problems, the present invention provides a relative sea level change analysis method, which, in one possible implementation, includes:

[0006] Acquiring logging data of microresistivity imaging of sedimentary rock formations, and preprocessing the logging data to obtain dynamic logging images of microresistivity imaging;

[0007] Acquiring the grayscale value of the logging dynamic image, and calculating the logging lithology curve of the microresistivity imaging based on the grayscale value;

[0008] Calculating high-order derivatives of the microresistivity imaging logging lithology curve, and determining the depth of the sedimentary rock layer and the depth of the sedimentary rock cycle layer interface based on the high-order derivatives;

[0009] Calculating the thickness of the cyclic sublayer and the average thickness of the cyclic sublayer within the depth range of the sedimentary rock layer based on the depth of the cyclic sublayer interface of the sedimentary rock;

[0010] Based on the cyclic sublayer thickness and the cyclic sublayer average thickness, calculating the cyclic sublayer offset thickness and the cyclic sublayer cumulative offset thickness within the depth range of the sedimentary rock layer;

[0011] The relative sea level changes of sedimentary rock strata are analyzed based on the changing trend of the cumulative offset thickness of the cyclic sub-layers.

[0012] In a possible implementation, the microresistivity imaging logging lithology curve includes:

[0013] in, is the logging lithology curve value of the microresistivity imaging corresponding to the current depth point, is the number of grayscale values ​​of the dynamic image of microresistivity imaging logging at the current depth point, The current depth point microresistivity imaging logging dynamic image Gray value of pixels, 0≤ ≤255, dimensionless.

[0014] In a possible implementation, the high-order derivatives of the microresistivity imaging logging lithology curve include the first-order derivative, the second-order derivative and the third-order derivative of the microresistivity imaging logging lithology curve;

[0015] The calculating of the high-order derivative of the microresistivity imaging logging lithology curve based on the logging lithology curve comprises:

[0016] The first derivative of the microresistivity imaging logging lithology curve is calculated based on the logging lithology curve:

[0017]

[0018] in, is the first derivative of the logging lithology curve of the microresistivity imaging corresponding to the current depth point, is the logging lithology curve value of the microresistivity imaging corresponding to the previous depth point, is the logging lithology curve value of the microresistivity imaging corresponding to the current depth point, The depth interval of the lithologic curve of microresistivity imaging logging;

[0019] The second-order derivative of the microresistivity imaging logging lithology curve is calculated based on the first-order derivative:

[0020]

[0021] in, is the second-order derivative of the microresistivity imaging logging lithology curve corresponding to the current depth point, is the first derivative of the microresistivity imaging logging lithology curve corresponding to the previous depth point;

[0022] The third-order derivative of the microresistivity imaging logging lithology curve is calculated based on the second-order derivative:

[0023]

[0024] in, is the third-order derivative of the microresistivity imaging logging lithology curve corresponding to the current depth point, It is the second derivative of the microresistivity imaging logging lithology curve corresponding to the previous depth point.

[0025] In a possible implementation, determining the sedimentary rock layer depth and the sedimentary rock cycle sub-layer interface depth based on the high-order derivatives includes:

[0026] The first depth value at which the second-order derivative of the microresistivity imaging logging lithology curve changes from negative to positive and the third-order derivative is not equal to 0 is determined as the depth of the bottom of the sedimentary rock layer, and the depth of the bottom interface of the sedimentary rock cycle layer is determined based on the depth of the bottom of the sedimentary rock layer;

[0027] The second depth value at which the second-order derivative of the microresistivity imaging logging lithology curve changes from positive to negative and the third-order derivative is not equal to 0 is determined as the top depth of the sedimentary rock layer; and based on the top depth of the sedimentary rock layer, the depth of the top interface of the sedimentary rock cycle layer is determined.

[0028] In a possible implementation, the calculating of the thickness of the cyclic sublayer and the average thickness of the cyclic sublayer within the depth range of the sedimentary rock layer based on the depth of the cyclic sublayer interface of the sedimentary rock includes:

[0029]

[0030]

[0031] in, is the sequence number of the cyclic sublayer from bottom to top, For the The depth of the bottom interface of the cyclic layer is For the The depth of the top interface of the cyclic layer, For the The thickness of the cyclic layer, is the total number of cyclic layers within the depth range of the sedimentary rock formation, is the average thickness of the cyclic sublayer.

[0032] In a possible implementation, the calculation of the cyclic sublayer offset thickness and the cyclic sublayer cumulative offset thickness within the depth range of the sedimentary rock layer based on the cyclic sublayer thickness and the cyclic sublayer average thickness includes:

[0033]

[0034]

[0035] in, For the The offset thickness of the small layer of each cycle is For the The cumulative offset thickness of the small layer in each cycle.

[0036] In a possible implementation, analyzing the relative sea level change of the sedimentary rock formation according to the change trend of the cumulative offset thickness of the cyclic sublayer includes:

[0037] If the cumulative offset thickness of the cyclic sublayer increases gradually from bottom to top, the relative sea level change of the sedimentary rock strata is in an upward state;

[0038] If the cumulative offset thickness of the cyclic sublayer gradually decreases from bottom to top, the relative sea level change of the sedimentary rock bottom layer is in a downward state.

[0039] In a second aspect, the present invention further provides a relative sea level change analysis device, comprising:

[0040] A preprocessing module is used to obtain logging data of microresistivity imaging of sedimentary rock formations, preprocess the logging data, and obtain a dynamic logging image of microresistivity imaging;

[0041] a curve calculation module, configured to obtain the grayscale value of the well logging dynamic image and calculate the well logging lithology curve of the microresistivity imaging based on the grayscale value;

[0042] a depth calculation module, configured to calculate the higher-order derivatives of the microresistivity imaging logging lithology curve, and determine the depth of the sedimentary rock layer and the depth of the sedimentary rock cycle sub-layer interface based on the higher-order derivatives;

[0043] A thickness calculation module is used to calculate the thickness of the cyclic sub-layer and the average thickness of the cyclic sub-layer within the depth range of the sedimentary rock layer based on the depth of the cyclic sub-layer interface of the sedimentary rock;

[0044] An offset thickness calculation module is used to calculate the offset thickness of the cyclic sub-layer and the cumulative offset thickness of the cyclic sub-layer within the depth range of the sedimentary rock layer based on the cyclic sub-layer thickness and the average thickness of the cyclic sub-layer;

[0045] The analysis module is used to analyze the relative sea level changes of sedimentary rock strata based on the changing trend of the cumulative offset thickness of the cyclic sub-layers.

[0046] In a third aspect, the present invention also provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps in the above-mentioned relative sea level change analysis method when executing the computer program.

[0047] In a fourth aspect, the present invention further provides a computer storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps in the relative sea level change analysis method as described above are implemented.

[0048] The present invention obtains logging data of microresistivity imaging of sedimentary rock formations, pre-processes the logging data to form a dynamic image of microresistivity imaging logging, calculates the grayscale average of the dynamic image of microresistivity imaging logging, forms a high-resolution microresistivity imaging logging lithology curve, calculates high-order derivatives of the microresistivity imaging logging lithology curve, determines the depth of the rock formation, automatically picks up the depth of the cyclical sub-layer interface, calculates the thickness of the cyclical sub-layer within the depth range of the sedimentary rock formation, the average thickness of the cyclical sub-layer, the offset thickness of the cyclical sub-layer and the cumulative offset thickness of the cyclical sub-layer, and finally analyzes the relative sea level change of the sedimentary rock formation according to the vertical change trend of the cumulative offset thickness of the cyclical sub-layer from bottom to top, thereby improving the accuracy and effectiveness of the research on the relative sea level change of the sedimentary rock formation. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative work.

[0050] Figure 1 A method flow chart of an embodiment of a relative sea level change analysis method provided by the present invention;

[0051] Figure 2 A schematic diagram of microresistivity imaging logging data of a portion of a sedimentary rock formation provided by one embodiment of the present invention;

[0052] Figure 3 A schematic diagram of a microresistivity imaging logging dynamic image of a portion of a well section in a sedimentary rock formation, a microresistivity imaging logging lithology curve, and the first, second, and third derivatives of the lithology curve, provided in one embodiment of the present invention;

[0053] Figure 4 A schematic diagram of the depth of rock layers and the depth of cyclic sub-layer interfaces in a sedimentary rock formation provided by one embodiment of the present invention;

[0054] Figure 5 A schematic diagram of the research results of the cyclical sub-layer offset thickness, the cyclical sub-layer cumulative offset thickness, and the relative sea level change within the depth range of a sedimentary rock layer provided by one embodiment of the present invention;

[0055] Figure 6A schematic structural diagram of an embodiment of a relative sea level change analysis device provided by the present invention;

[0056] Figure 7 A schematic structural diagram of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0057] The preferred embodiments of the present invention will be described in detail below in conjunction with the accompanying drawings, wherein the accompanying drawings constitute a part of this application and are used together with the embodiments of the present invention to illustrate the principles of the present invention, and are not used to limit the scope of the present invention.

[0058] The terms "first," "second," and so on, used in the embodiments of the present invention are for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Therefore, technical features designated as "first" or "second" may explicitly or implicitly include at least one such feature.

[0059] References herein to "embodiments" mean that a particular feature, structure, or characteristic described in connection with the embodiments may be included in at least one embodiment of the present invention. The appearance of this phrase in various places in the specification does not necessarily refer to the same embodiment, nor does it constitute a separate or alternative embodiment that is mutually exclusive of other embodiments. It is understood, both explicitly and implicitly, by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0060] A specific embodiment of the present invention discloses a relative sea level change analysis method, comprising:

[0061] Step 101: acquiring logging data of microresistivity imaging of sedimentary rock formations, preprocessing the logging data, and obtaining a dynamic logging image of microresistivity imaging;

[0062] Step 102: obtaining the grayscale value of the logging dynamic image, and calculating the logging lithology curve of the microresistivity imaging based on the grayscale value;

[0063] Step 103: Calculate the high-order derivatives of the microresistivity imaging logging lithology curve, and determine the depth of the sedimentary rock layer and the depth of the sedimentary rock cycle interface based on the high-order derivatives;

[0064] Step 104: Calculating the thickness of the cyclic sub-layer and the average thickness of the cyclic sub-layer within the depth range of the sedimentary rock layer based on the depth of the cyclic sub-layer interface of the sedimentary rock;

[0065] Step 105: Calculate the cyclic sub-layer offset thickness and the cyclic sub-layer cumulative offset thickness within the depth range of the sedimentary rock layer based on the cyclic sub-layer thickness and the cyclic sub-layer average thickness;

[0066] Step 106: Analyze the relative sea level change of the sedimentary rock strata based on the change trend of the cumulative offset thickness of the cyclic sub-layer.

[0067] Among them, the lithology of sedimentary rock formations includes mudstone, sandstone, dolomite and limestone. The raw data of microresistivity imaging logging are preprocessed, which usually includes removing outliers in the logging data and enhancing the dynamic logging images. The raw data of microresistivity imaging logging of the entire well section are obtained from the sedimentary rock formation. For some data, please refer to Figure 2 .

[0068] The present invention obtains logging data of microresistivity imaging of sedimentary rock formations, pre-processes the logging data to form a dynamic image of microresistivity imaging logging, calculates the grayscale average of the dynamic image of microresistivity imaging logging, forms a high-resolution microresistivity imaging logging lithology curve, calculates high-order derivatives of the microresistivity imaging logging lithology curve, determines the depth of the rock formation, automatically picks up the depth of the cyclical sub-layer interface, calculates the thickness of the cyclical sub-layer within the depth range of the sedimentary rock formation, the average thickness of the cyclical sub-layer, the offset thickness of the cyclical sub-layer and the cumulative offset thickness of the cyclical sub-layer, and finally analyzes the relative sea level change of the sedimentary rock formation according to the vertical change trend of the cumulative offset thickness of the cyclical sub-layer from bottom to top, thereby improving the accuracy and effectiveness of the research on the relative sea level change of the sedimentary rock formation.

[0069] In one embodiment of the present invention, the microresistivity imaging logging lithology curve includes: It can be understood that based on the dynamic image of micro-resistivity imaging logging, the grayscale average value of the dynamic image of micro-resistivity imaging logging is calculated to form a high-resolution micro-resistivity imaging logging lithology curve. is the logging lithology curve value of the microresistivity imaging corresponding to the current depth point, is the number of grayscale values ​​of the dynamic image of microresistivity imaging logging at the current depth point, The current depth point microresistivity imaging logging dynamic image Gray value of pixels, 0≤ ≤255, dimensionless.

[0070] In one embodiment of the present invention, the high-order derivatives of the microresistivity imaging logging lithology curve include the first-order derivative of the microresistivity imaging logging lithology curve, the second-order derivative of the microresistivity imaging logging lithology curve, and the third-order derivative of the microresistivity imaging logging lithology curve;

[0071] Calculate the high-order derivatives of microresistivity imaging logging lithology curves based on logging lithology curves, including:

[0072] Calculate the first derivative of the microresistivity imaging logging lithology curve based on the logging lithology curve:

[0073]

[0074] in, is the first derivative of the logging lithology curve of the microresistivity imaging corresponding to the current depth point, is the logging lithology curve value of the microresistivity imaging corresponding to the previous depth point, is the logging lithology curve value of the microresistivity imaging corresponding to the current depth point, The depth interval of the lithologic curve of microresistivity imaging logging;

[0075] Calculate the second-order derivative of the microresistivity imaging logging lithology curve based on the first-order derivative:

[0076]

[0077] in, is the second-order derivative of the microresistivity imaging logging lithology curve corresponding to the current depth point, is the first derivative of the microresistivity imaging logging lithology curve corresponding to the previous depth point;

[0078] Calculate the third-order derivative of the microresistivity imaging logging lithology curve based on the second-order derivative:

[0079]

[0080] in, is the third-order derivative of the microresistivity imaging logging lithology curve corresponding to the current depth point, It is the second derivative of the microresistivity imaging logging lithology curve corresponding to the previous depth point.

[0081] Determine the depth of sedimentary rock layers and the depth of sedimentary rock cycle interface based on high-order derivatives, including:

[0082] The first depth value at which the second-order derivative of the microresistivity imaging logging lithology curve changes from negative to positive and the third-order derivative is not equal to 0 is determined as the depth of the bottom of the sedimentary rock layer, and based on the depth of the bottom of the sedimentary rock layer, the depth of the bottom interface of the sedimentary rock cycle layer is determined;

[0083] The second depth value where the second-order derivative of the microresistivity imaging logging lithology curve changes from positive to negative and the third-order derivative is not equal to 0 is determined as the top depth of the sedimentary rock layer; and based on the top depth of the sedimentary rock layer, the depth of the top interface of the sedimentary rock cycle layer is determined.

[0084] It can be understood that by calculating the first derivative, second derivative and third derivative of the microresistivity imaging logging lithology curve. Figure 3 As shown, Figure 3An embodiment of the present invention provides a dynamic image of microresistivity imaging logging of a part of a well section of a sedimentary rock formation, a microresistivity imaging logging lithologic curve, and a schematic diagram of the first-order derivative, second-order derivative, and third-order derivative of the lithologic curve. The depth of the rock formation is determined based on the microresistivity imaging logging lithologic curve. If the depth value of the second-order derivative of the lithologic curve changes from negative to positive and the third-order derivative is not equal to 0, it is the depth of the bottom of the rock formation; if the depth value of the second-order derivative of the lithologic curve changes from positive to negative and the third-order derivative is not equal to 0, it is the depth of the top of the rock formation. The depth of the cyclic sub-layer interface is automatically picked up. The depth of the bottom of the sandstone, dolomite or limestone rock layer is the bottom interface of the cyclic sub-layer; the depth of the top of the mudstone rock layer overlying the sandstone, dolomite or limestone is the depth of the top interface of the cyclic sub-layer. As Figure 4 As shown, Figure 4 A schematic diagram of the rock layer depth and cyclic layer interface depth of a partial well section of a sedimentary rock formation provided by one embodiment of the present invention.

[0085] In one embodiment of the present invention, the thickness of the cyclic sub-layer and the average thickness of the cyclic sub-layer within the depth range of the sedimentary rock layer are calculated based on the depth of the cyclic sub-layer interface of the sedimentary rock, including:

[0086]

[0087]

[0088] in, is the sequence number of the cyclic sublayer from bottom to top, For the The depth of the bottom interface of the cyclic layer is For the The depth of the top interface of the cyclic layer, For the The thickness of the cyclic layer, is the total number of cyclic layers within the depth range of the sedimentary rock formation, is the average thickness of the cyclic sublayer.

[0089] It can be understood that, based on the depth of the cyclic layer interface, the thickness of the cyclic layer within the depth range of the sedimentary rock layer can be calculated, and the average thickness of the cyclic layer can be further calculated.

[0090] In one embodiment of the present invention, based on the cyclic layer thickness and the cyclic layer average thickness, calculating the cyclic layer offset thickness and the cyclic layer cumulative offset thickness within the depth range of the sedimentary rock layer includes:

[0091]

[0092]

[0093] in, For the The offset thickness of the small layer of each cycle is For the The cumulative offset thickness of the small layer in each cycle is Dimensionless.

[0094] According to the variation trend of the cumulative offset thickness of the cyclic sub-layer, the relative sea level change of the sedimentary rock strata is analyzed, including:

[0095] If the cumulative offset thickness of the cyclic sublayer increases gradually from bottom to top, the relative sea level change of the sedimentary rock strata is in an upward state;

[0096] If the cumulative offset thickness of the cyclic sublayer gradually decreases from bottom to top, the relative sea level change of the sedimentary rock bottom layer is in a downward state.

[0097] like Figure 5 As shown, Figure 5 A schematic diagram of the research results of the cyclical sub-layer offset thickness, cyclical sub-layer cumulative offset thickness and relative sea level change within the depth range of a sedimentary rock layer provided by one embodiment of the present invention. Figure 5 It can be seen from the figure that according to the vertical variation trend of the cumulative offset thickness of the cyclic small layer from bottom to top, the relative sea level change of the sedimentary rock strata is analyzed. The cumulative offset thickness of the cyclic small layer gradually increases from bottom to top, which means the relative sea level rises, and the cumulative offset thickness of the cyclic small layer gradually decreases from bottom to top, which means the relative sea level falls.

[0098] More specifically, in one embodiment of the present invention, referring to Figure 4 , a total of 517 cyclic layers were picked up, with a layer density of 8.86 layers / m, a minimum thickness of 2.42 cm, a maximum thickness of 74.38 cm, and an average thickness of 11.28 cm. According to the trend of the cumulative offset thickness of the cyclic layers from bottom to top, Figure 5 , a total of 5 periods of relative sea level rise and fall characteristics were identified.

[0099] In order to better implement the relative sea level change analysis method in the embodiment of the present invention, based on the relative sea level change analysis method, correspondingly, Figure 6 As shown, the embodiment of the present invention further provides a relative sea level change analysis device 600, which includes:

[0100] A preprocessing module 601 is used to obtain logging data of microresistivity imaging of sedimentary rock formations, preprocess the logging data, and obtain dynamic logging images of microresistivity imaging;

[0101] The curve calculation module 602 is used to obtain the gray value of the logging dynamic image and calculate the logging lithology curve of the microresistivity imaging based on the gray value;

[0102] Depth calculation module 603, used to calculate the high-order derivatives of the microresistivity imaging logging lithology curve, and determine the depth of the sedimentary rock layer and the depth of the sedimentary rock cycle layer interface based on the high-order derivatives;

[0103] Thickness calculation module 604, used to calculate the thickness of the cyclic sub-layer and the average thickness of the cyclic sub-layer within the depth range of the sedimentary rock layer based on the depth of the cyclic sub-layer interface of the sedimentary rock;

[0104] The offset thickness calculation module 605 is used to calculate the offset thickness of the cyclic sub-layer and the cumulative offset thickness of the cyclic sub-layer within the depth range of the sedimentary rock layer based on the cyclic sub-layer thickness and the average thickness of the cyclic sub-layer;

[0105] The analysis module 606 is used to analyze the relative sea level change of the sedimentary rock strata according to the change trend of the cumulative offset thickness of the cyclic sub-layer.

[0106] The relative sea level change analysis device 600 provided in the above embodiment can implement the technical solution described in the above relative sea level change analysis method embodiment. The specific implementation principles of the above modules or units can refer to the corresponding contents in the above relative sea level change analysis method embodiment, which will not be repeated here.

[0107] like Figure 7 As shown, the present invention also provides an electronic device 700. The electronic device 700 includes a processor 701, a memory 702 and a display 703. Figure 7 Only some of the components of the electronic device 700 are shown, but it should be understood that it is not required to implement all of the shown components, and more or fewer components may be implemented instead.

[0108] In some embodiments, the processor 701 may be a central processing unit (CPU), a microprocessor, or other data processing chip, configured to execute program codes or process data stored in the memory 702, such as the relative sea level change analysis method of the present invention.

[0109] In some embodiments, processor 701 may be a single server or a server group. The server group may be centralized or distributed. In some embodiments, processor 701 may be local or remote. In some embodiments, processor 701 may be implemented on a cloud platform. In one embodiment, the cloud platform may include a private cloud, a public cloud, a hybrid cloud, a community cloud, a distributed cloud, an internal cloud, multiple clouds, or any combination thereof.

[0110] In some embodiments, the memory 702 may be an internal storage unit of the electronic device 700, such as a hard disk or memory of the electronic device 700. In other embodiments, the memory 702 may also be an external storage device of the electronic device 700, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. equipped on the electronic device 700.

[0111] Furthermore, the memory 702 may include both an internal storage unit of the electronic device 700 and an external storage device. The memory 702 is used to store application software installed in the electronic device 700 and various data.

[0112] In some embodiments, display 703 can be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, or an OLED (Organic Light-Emitting Diode) touchscreen. Display 703 is used to display information about electronic device 700 and to display a visual user interface. Components 701-703 of electronic device 700 communicate with each other via a system bus.

[0113] In some embodiments, when the processor 701 executes the relative sea level change analysis program in the memory 702, the following steps may be implemented:

[0114] Acquire logging data of microresistivity imaging of sedimentary rock formations, pre-process the logging data, and obtain dynamic logging images of microresistivity imaging;

[0115] Obtaining the grayscale value of the logging dynamic image and calculating the logging lithology curve of the microresistivity imaging based on the grayscale value;

[0116] Calculate the high-order derivatives of the microresistivity imaging logging lithology curve and determine the depth of the sedimentary rock layer and the depth of the sedimentary rock cycle interface based on the high-order derivatives;

[0117] Based on the depth of the sedimentary rock cyclic layer interface, the thickness of the cyclic layer and the average thickness of the cyclic layer within the depth range of the sedimentary rock layer are calculated;

[0118] Based on the thickness of the cyclic layer and the average thickness of the cyclic layer, the offset thickness of the cyclic layer and the cumulative offset thickness of the cyclic layer within the depth range of the sedimentary rock layer are calculated;

[0119] The relative sea level changes of sedimentary rock strata are analyzed based on the changing trend of the cumulative offset thickness of the cyclic sub-layers.

[0120] It should be understood that, when the processor 701 executes the relative sea level change analysis program in the memory 702 , in addition to the above functions, it can also implement other functions. For details, please refer to the description of the corresponding method embodiment above.

[0121] Furthermore, the embodiment of the present invention does not specifically limit the type of the electronic device 700 mentioned. The electronic device 700 may be a portable electronic device such as a mobile phone, a tablet computer, a personal digital assistant (PDA), a wearable device, a laptop computer, or the like. Exemplary embodiments of portable electronic devices include, but are not limited to, portable electronic devices equipped with IOS, Android, Microsoft, or other operating systems. The above-mentioned portable electronic devices may also be other portable electronic devices, such as a laptop computer with a touch-sensitive surface (e.g., a touch panel). It should also be understood that in some other embodiments of the present invention, the electronic device 700 may not be a portable electronic device, but a desktop computer with a touch-sensitive surface (e.g., a touch panel).

[0122] Accordingly, an embodiment of the present application also provides a computer-readable storage medium, which is used to store computer-readable programs or instructions. When the program or instructions are executed by a processor, it can implement the steps or functions in the relative sea level change analysis method provided in the above-mentioned method embodiments.

[0123] Those skilled in the art will appreciate that all or part of the process steps of the above-described embodiments can be implemented by instructing related hardware through a computer program, and the program can be stored in a computer-readable storage medium, such as a magnetic disk, an optical disk, a read-only memory, or a random access memory.

[0124] The above description is only a preferred specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any changes or substitutions that can be easily thought of by any technician familiar with this technical field within the technical scope disclosed by the present invention should be covered by the scope of protection of the present invention.

Claims

1. A relative sea level change analysis method, characterized in that: include: Acquiring logging data of microresistivity imaging of sedimentary rock formations, and preprocessing the logging data to obtain dynamic logging images of microresistivity imaging; Acquiring a grayscale value of the dynamic logging image, and calculating a microresistivity imaging logging lithology curve based on the grayscale value; Calculating high-order derivatives of the microresistivity imaging logging lithology curve, and determining the depth of the sedimentary rock layer and the depth of the sedimentary rock cycle layer interface based on the high-order derivatives; Calculating the thickness of the cyclic sublayer and the average thickness of the cyclic sublayer within the depth range of the sedimentary rock layer based on the depth of the cyclic sublayer interface of the sedimentary rock; Based on the cyclic sublayer thickness and the cyclic sublayer average thickness, calculating the cyclic sublayer offset thickness and the cyclic sublayer cumulative offset thickness within the depth range of the sedimentary rock layer; The relative sea level changes of sedimentary rock strata are analyzed based on the changing trend of the cumulative offset thickness of the cyclic sub-layers.

2. The relative sea level change analysis method according to claim 1, characterized in that: The microresistivity imaging logging lithology curve includes: , in, is the logging lithology curve value of the microresistivity imaging corresponding to the current depth point, is the number of grayscale values ​​of the dynamic image of microresistivity imaging logging at the current depth point, The current depth point microresistivity imaging logging dynamic image Gray value of pixels, 0≤ ≤255, dimensionless.

3. The relative sea level change analysis method according to claim 1, characterized in that: The high-order derivatives of the microresistivity imaging logging lithology curve include the first-order derivative of the microresistivity imaging logging lithology curve, the second-order derivative of the microresistivity imaging logging lithology curve and the third-order derivative of the microresistivity imaging logging lithology curve; The calculating of the high-order derivative of the microresistivity imaging logging lithology curve based on the logging lithology curve comprises: The first derivative of the microresistivity imaging logging lithology curve is calculated based on the logging lithology curve: in, is the first derivative of the logging lithology curve of the microresistivity imaging corresponding to the current depth point, is the logging lithology curve value of the microresistivity imaging corresponding to the previous depth point, is the logging lithology curve value of the microresistivity imaging corresponding to the current depth point, The depth interval of the lithologic curve of microresistivity imaging logging; The second-order derivative of the microresistivity imaging logging lithology curve is calculated based on the first-order derivative: in, is the second-order derivative of the microresistivity imaging logging lithology curve corresponding to the current depth point, is the first derivative of the microresistivity imaging logging lithology curve corresponding to the previous depth point; The third-order derivative of the microresistivity imaging logging lithology curve is calculated based on the second-order derivative: in, is the third-order derivative of the microresistivity imaging logging lithology curve corresponding to the current depth point, It is the second derivative of the microresistivity imaging logging lithology curve corresponding to the previous depth point.

4. The relative sea level change analysis method according to claim 3, characterized in that: The method of determining the depth of the sedimentary rock layer and the depth of the sedimentary rock cycle interface based on the high-order derivatives includes: The first depth value at which the second-order derivative of the microresistivity imaging logging lithology curve changes from negative to positive and the third-order derivative is not equal to 0 is determined as the depth of the bottom of the sedimentary rock layer, and the depth of the bottom interface of the sedimentary rock cycle layer is determined based on the depth of the bottom of the sedimentary rock layer; The second depth value at which the second-order derivative of the microresistivity imaging logging lithology curve changes from positive to negative and the third-order derivative is not equal to 0 is determined as the top depth of the sedimentary rock layer; and based on the top depth of the sedimentary rock layer, the depth of the top interface of the sedimentary rock cycle layer is determined.

5. The relative sea level change analysis method according to claim 4, characterized in that: The step of calculating the thickness of the cyclic sub-layer and the average thickness of the cyclic sub-layer within the depth range of the sedimentary rock layer based on the depth of the cyclic sub-layer interface of the sedimentary rock comprises: in, is the sequence number of the cyclic sublayer from bottom to top, For the The depth of the bottom interface of the cyclic layer is For the The depth of the top interface of the cyclic layer, For the The thickness of the cyclic layer, is the total number of cyclic layers within the depth range of the sedimentary rock formation, is the average thickness of the cyclic sublayer.

6. The relative sea level change analysis method according to claim 5, characterized in that: The step of calculating the offset thickness of the cyclic sublayer and the cumulative offset thickness of the cyclic sublayer within the depth range of the sedimentary rock layer based on the cyclic sublayer thickness and the average cyclic sublayer thickness includes: in, For the The offset thickness of the small layer of each cycle is For the The cumulative offset thickness of the small layer in each cycle.

7. The relative sea level change analysis method according to claim 6, characterized in that: The analysis of the relative sea level change of the sedimentary rock strata based on the change trend of the cumulative offset thickness of the cyclic sub-layers includes: If the cumulative offset thickness of the cyclic sublayer increases gradually from bottom to top, the relative sea level change of the sedimentary rock strata is in an upward state; If the cumulative offset thickness of the cyclic sublayer gradually decreases from bottom to top, the relative sea level change of the sedimentary rock bottom layer is in a downward state.

8. A relative sea level change analysis device, characterized in that: include: A preprocessing module is used to obtain logging data of microresistivity imaging of sedimentary rock formations, preprocess the logging data, and obtain a dynamic logging image of microresistivity imaging; a curve calculation module, configured to obtain the grayscale value of the well logging dynamic image and calculate the well logging lithology curve of the microresistivity imaging based on the grayscale value; a depth calculation module, configured to calculate the higher-order derivatives of the microresistivity imaging logging lithology curve, and determine the depth of the sedimentary rock layer and the depth of the sedimentary rock cycle sub-layer interface based on the higher-order derivatives; A thickness calculation module is used to calculate the thickness of the cyclic sub-layer and the average thickness of the cyclic sub-layer within the depth range of the sedimentary rock layer based on the depth of the cyclic sub-layer interface of the sedimentary rock; An offset thickness calculation module is used to calculate the offset thickness of the cyclic sub-layer and the cumulative offset thickness of the cyclic sub-layer within the depth range of the sedimentary rock layer based on the cyclic sub-layer thickness and the average thickness of the cyclic sub-layer; The analysis module is used to analyze the relative sea level changes of sedimentary rock strata based on the changing trend of the cumulative offset thickness of the cyclic sub-layers.

9. An electronic device, characterized in that: comprising a memory and a processor, wherein, The memory is used to store programs; The processor is coupled to the memory and is configured to execute the program stored in the memory to implement the steps of the relative sea level change analysis method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that Used to store computer-readable programs or instructions, which, when executed by a processor, can implement the steps of the relative sea level change analysis method described in any one of claims 1 to 7.

Images