A method and device for determining shale content based on seismic elastic parameters

Abstract

The application discloses a method and device for determining shale content based on seismic elastic parameters. The method comprises the following steps: determining the distance between the to-be-estimated point and a first template line in a two-dimensional coordinate system of longitudinal wave impedance and transverse wave impedance according to the longitudinal wave impedance and the transverse wave impedance of the to-be-estimated point; and determining the shale content of the to-be-estimated point by means of distance and shale content difference in a proportional manner according to the distance, the standard distance between the first template line and a second template line, and the difference between the shale content corresponding to the first template line and the shale content corresponding to the second template line, wherein the first template line and the second template line are two parallel straight lines in the coordinate system, the shale content of each point on the first template line is the same, and the shale content of each point on the second template line is the same. The method can quantitatively calculate the shale content based on seismic elastic parameter inversion data (longitudinal wave impedance and transverse wave impedance), has strong operability, high precision, and can be better applied to production.

Description

Technical Field

[0001] This invention relates to the fields of oil and gas geophysical exploration and well logging technology, and in particular to a method and apparatus for determining clay content based on seismic elastic parameters. Background Technology

[0002] Accurate prediction of clay content (Vsh) has a significant impact on the prediction of other reservoir parameters, drilling design, and wellbore stability evaluation. It can be said that accurate clay content prediction plays a crucial role in finding high-quality oil and gas reservoirs.

[0003] Currently, there are two main methods for quantitative prediction of clay content using seismic data: (1) Directly constructing the relationship between clay content and elastic curves such as P-wave impedance by fitting well logging data, and applying this relationship to the elastic data volume for the prediction of the entire clay content. (2) Constructing the relationship between clay content and elastic parameters by using inverse rock physics modeling, and applying this relationship to the prediction of clay content. Summary of the Invention

[0004] The inventors found that the above methods are difficult to apply to the actual production of clay content because: the first method is simple and easy to implement, but the linear or nonlinear relationship of regression is difficult to apply to seismic inversion data, because seismic inversion seismic data often only yields parameters such as P-wave impedance and S-wave impedance; the second method is theoretically advanced, but involves many calculation parameters, which leads to instability in the solution process and makes it difficult to operate in practice.

[0005] In order to at least partially solve the technical problems existing in the prior art, the inventors made this invention, which provides a method and apparatus for determining clay content based on seismic elastic parameters through specific implementation methods. It can quantitatively calculate clay content based on seismic elastic parameter inversion data (longitudinal wave impedance, transverse wave impedance), which is highly operable, has high accuracy, and can be better applied to production.

[0006] In a first aspect, embodiments of the present invention provide a method for determining clay content based on seismic elastic parameters, comprising:

[0007] Based on the longitudinal and transverse wave impedances of the point to be estimated, determine the distance between the point to be estimated and the first template line in the two-dimensional coordinate system of longitudinal and transverse wave impedances.

[0008] Based on the distance, the standard distance between the first template line and the second template line, and the difference between the mud content corresponding to the first template line and the mud content corresponding to the second template line, the mud content of the point to be estimated is determined by proportionally dividing the distance by the difference in mud content. The first template line and the second template line are two parallel straight lines in the coordinate system. The mud content is the same at all points on the first template line and the mud content is the same at all points on the second template line.

[0009] Secondly, embodiments of the present invention provide a device for determining clay content based on seismic elastic parameters, comprising:

[0010] The distance determination module is used to determine the distance between the point to be estimated and the first template line in the two-dimensional coordinate system of the longitudinal and transverse wave impedances, based on the longitudinal and transverse wave impedances of the point to be estimated.

[0011] The mud content determination module is used to determine the mud content of the point to be estimated by means of the distance, the standard distance between the first template line and the second template line, and the difference between the mud content corresponding to the first template line and the mud content corresponding to the second template line, in a proportional manner to the difference in mud content. The first template line and the second template line are two parallel straight lines in the coordinate system. The mud content of each point on the first template line is the same, and the mud content of each point on the second template line is the same.

[0012] Thirdly, embodiments of the present invention provide a computer storage medium storing computer-executable instructions, which, when executed by a processor, implement the above-described method for determining mud content based on seismic elastic parameters.

[0013] Fourthly, embodiments of this disclosure provide a server, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the above-described method for determining mud content based on seismic elastic parameters.

[0014] The beneficial effects of the above-described technical solutions provided in the embodiments of the present invention include at least the following:

[0015] The method for determining clay content based on seismic elastic parameters provided in this invention starts directly from seismic elastic parameter inversion data (P-wave impedance, S-wave impedance), pre-establishes a template line for interpreting seismic elastic parameters for clay content, and calculates the clay content by solving the distance between the point to be estimated and the template line. This method is highly operable, has high accuracy, and can be better applied to production.

[0016] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description, claims, and drawings.

[0017] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description

[0018] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used in conjunction with embodiments of the invention to explain the invention and do not constitute a limitation thereof. In the drawings:

[0019] Figure 1 This is a flowchart of the method for establishing template lines for interpreting seismic elastic parameters in Embodiment 1 of the present invention;

[0020] Figure 2 This is a schematic diagram of the template line established in Embodiment 1 of the present invention;

[0021] Figure 3 This is a flowchart of the method for determining clay content based on seismic elastic parameters in Embodiment 2 of the present invention;

[0022] Figure 4 This is a schematic diagram showing the calculation results of mud content in Embodiment 2 of the present invention;

[0023] Figure 5 This is a schematic diagram of the device for determining mud content based on seismic elastic parameters in an embodiment of the present invention. Detailed Implementation

[0024] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0025] It should be understood that the terminology used in this invention is merely for describing particular embodiments and is not intended to limit the invention. Furthermore, with respect to numerical ranges in this invention, it should be understood that each intermediate value between the upper and lower limits of the range is also specifically disclosed. Every smaller range between any stated value or intermediate value within a stated range, and any other stated value or intermediate value within said range, is also included in this invention. The upper and lower limits of these smaller ranges may be independently included or excluded from the range.

[0026] Unless otherwise stated, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains. While only preferred methods and materials have been described herein, any methods and materials similar or equivalent to those described herein may be used in the implementation or testing of this invention. All references to this specification are incorporated by way of citation to disclose and describe methods and / or materials associated with those references. In the event of any conflict with any incorporated reference, the content of this specification shall prevail.

[0027] Example 1

[0028] Embodiment 1 of the present invention provides a method for establishing template lines for the interpretation of seismic elastic parameters, the process of which is as follows: Figure 1 As shown, it includes the following steps:

[0029] Step S11: Based on the first sample set, establish a rock physics model that characterizes the relationship between clay content, porosity and P-wave and S-wave velocities through multiple linear regression.

[0030] Each sample in the first sample set contains the clay content, porosity, P-wave velocity, and S-wave velocity at the same location, obtained through rock physics experiments and / or well logging data.

[0031] Based on the first sample set, a rock physics model characterizing the relationship between clay content, porosity, and P-wave and S-wave velocities is established through multiple linear regression as follows:

[0032] V P =AB·φ-C·V sh (1)

[0033] V S =DE·φ-FV sh (2)

[0034] In equations (1) and (2), V P V is the longitudinal wave velocity. S V is the transverse wave velocity, φ is the porosity, and V is the transverse wave velocity. sh denoted as mud content, and A, B, C, D, E, and F are fitting constants.

[0035] Step S12: Based on the second sample set, the third sample set is obtained by using the conversion relationship between P-wave velocity and P-wave impedance and the rock physics model.

[0036] Each sample in the second sample set contains the clay content and porosity at the same location, while each sample in the third sample set contains the clay content, porosity, longitudinal wave impedance, and transverse wave impedance at the same location.

[0037] Using the rock physics model, the third sample set is obtained through the following equations (3)-(5):

[0038] I P =V P ·ρ (3)

[0039] I S =V S ·ρ (4)

[0040] ρ=ρ ma *(1-V sh -φ)+ρ sh *V sh +ρ f *φ (5)

[0041] In equations (3)-(5), ρ ma ρ sh and ρ f Let ρ be the density of pure sandstone, pure mudstone, and fluid, respectively, where ρ is the reservoir density, and I is the density of the fluid. S For transverse wave impedance, I P This is the longitudinal wave impedance.

[0042] There are relatively few sample points in the study area that can be directly obtained, which include clay content, porosity, P-wave impedance, and S-wave impedance. By using a rock physics model, sample points that include clay content, porosity, P-wave impedance, and S-wave impedance can be calculated from the sample points that include clay content and porosity, thus expanding the number of samples in the third sample set.

[0043] Step S13: Determine multiple mud content values, extract multiple sub-sample sets from the third sample set, each sub-sample set corresponds to a mud content value, and project each sub-sample set into a two-dimensional coordinate system of longitudinal and transverse wave impedance with longitudinal wave impedance as the horizontal axis and transverse wave impedance as the vertical axis to obtain multiple template lines.

[0044] Determining multiple clay content values ​​can be achieved by fixing the clay content between 0 and 1, and then determining multiple clay content values ​​at set intervals (which can be flexibly set according to research needs). The largest of the determined clay content values ​​is 1, and the smallest is 0. See also Figure 2 As shown in the figure, the mud content V represents the mud content V. sh The template lines corresponding to 0.00, 0.20, 0.40, 0.60, 0.80, and 1.00 (from bottom to top, the clay content corresponding to the template lines increases sequentially) are shown. It can be seen that the template lines are parallel to each other, and their form is: Y = G·X + P j That is, G·X-Y+P j =0, where X is the longitudinal wave impedance I P Y is the transverse wave impedance I S G is the slope, P j This represents the intercept of each template line.

[0045] The distance between any two template lines can be determined by the following formula (6):

[0046]

[0047] In equation (6), D T P is the standard distance between template line 11 and template line 22. 11 and P 22 These are the intercepts of template line 11 and template line 22, respectively.

[0048] Calculations revealed that the ratio between the distance between any two template lines and the absolute value of the difference in their corresponding clay content is a constant. Therefore, two template lines can be selected as standards, and the distance between them and the corresponding difference in clay content can be calculated in advance. Then, by using the distance between the point to be estimated and one of the template lines, the clay content of the point to be estimated can be determined.

[0049] Based on the understanding of the above patterns, in subsequent practical applications, it is unnecessary to create too many template lines; only two standard template lines are needed. That is:

[0050] Based on steps S11 and S12, a first character sample set with a mud content of a first set value and a second character sample set with a mud content of a second set value are extracted from the third sample set. The extracted first character sample set and second character sample set are projected into the coordinate system to obtain the first template line corresponding to the first sub-sample set and the second template line corresponding to the second sub-sample set.

[0051] For ease of calculation, template lines are typically established corresponding to clay contents of 0 and 1. Otherwise, it is necessary to determine the relationship between the clay content of the point to be estimated and the clay content corresponding to the standard template line. Furthermore, it is necessary to determine whether the point to be estimated is above or below the standard template line. If it is above, the clay content of the point to be estimated is greater than the clay content corresponding to the standard template line; if it is below, the clay content of the point to be estimated is less than the clay content corresponding to the standard template line. However, by defining the standard template line (i.e., the first template line) as the template line with a clay content of 0, the step of determining the relationship between the clay content of the point to be estimated and the clay content corresponding to the standard template line can be eliminated, because the clay content will definitely be greater than or equal to 0.

[0052] In this embodiment of the invention, "distance" refers to vertical distance. Furthermore, the distance between two template lines refers to the vertical distance between them; the distance from the point to be estimated to the template line refers to the length of the vertical line segment obtained by drawing a perpendicular line from the point to be estimated to the template line. This distance is not a traditional spatial or temporal distance, but rather a distance in an extended sense, specifically a distance within the longitudinal and transverse wave impedance domains.

[0053] Example 2

[0054] Embodiment 2 of the present invention provides a method for determining clay content based on seismic elastic parameters, the process of which is as follows: Figure 3 As shown, it includes the following steps:

[0055] Step S31: Determine the distance between the point to be estimated and the first template line in the two-dimensional coordinate system of the longitudinal and transverse wave impedances based on the longitudinal and transverse wave impedances of the point to be estimated.

[0056] The distance between the point to be estimated and the first template line in the two-dimensional coordinate system of the longitudinal and transverse wave impedances is determined using the following formula (7):

[0057]

[0058] In equation (7), D i Let I be the distance between the point to be estimated with index i and the first template line. pi and I si Let G be the longitudinal wave impedance and transverse wave impedance of the point to be estimated, respectively; G be the slope of the template line in the coordinate system (the slopes of all template lines are equal); and P1 be the intercept of the first template line.

[0059] Step S32: Based on the distance, the standard distance between the first template line and the second template line, and the difference between the mud content corresponding to the first template line and the mud content corresponding to the second template line, determine the mud content of the point to be estimated by proportionally dividing the distance by the difference in mud content.

[0060] The first template line and the second template line are two parallel straight lines in the coordinate system. The mud content is the same at all points on the first template line and the mud content is the same at all points on the second template line.

[0061] The standard distance between the first template line and the second template line is predetermined by the following formula (8):

[0062]

[0063] In equation (8), D s G is the standard distance between the first template line and the second template line, G is the slope of the template line in the coordinate system, and P1 and P2 are the intercepts of the first template line and the second template line, respectively.

[0064] When the clay content at each point on the first template line is 0 and the clay content at each point on the second template line is 1, the clay content V of the point to be estimated is determined by the following formula (9). shi :

[0065]

[0066] In equation (9), D i D is the distance between the point to be estimated with index i and the first template line. sThis is the standard distance between the first template line and the second template line.

[0067] Previous methods for quantitatively predicting clay content have either low accuracy or, at best, limited applicability. This embodiment presents a method for determining clay content based on seismic elastic parameters. Using calibrated well or core data, a first and second sample set are obtained. A rock physical model characterizing the relationship between clay content, porosity, and P-wave and S-wave velocities is established. Furthermore, a seismic elastic parameter interpretation template line is created for clay content interpretation. The clay content is calculated by determining the distance from the point to be estimated to the template line. This method combines experimental data to construct a rock physical model, and the quantitative calculation of clay content can be achieved using relatively simple mathematical relationships, making it more suitable for industrial application.

[0068] The clay content was calculated using the above-mentioned method for determining clay content based on seismic elastic parameters. The calculation results are shown in the figure below. Figure 4 As can be seen, the verification well and the measured mud content have a good match, indicating that it has high accuracy.

[0069] Based on the inventive concept of this invention, embodiments of this invention also provide a device for determining clay content based on seismic elastic parameters, the structure of which is as follows: Figure 5 As shown, it includes:

[0070] The distance determination module 51 is used to determine the distance between the point to be estimated and the first template line in the two-dimensional coordinate system of the longitudinal and transverse wave impedances based on the longitudinal and transverse wave impedances of the point to be estimated.

[0071] The mud content determination module 52 is used to determine the mud content of the point to be estimated by means of the distance, the standard distance between the first template line and the second template line, and the difference between the mud content corresponding to the first template line and the mud content corresponding to the second template line, in a proportional manner to the difference in mud content. The first template line and the second template line are two parallel straight lines in the coordinate system. The mud content of each point on the first template line is the same, and the mud content of each point on the second template line is the same.

[0072] In some embodiments, the distance determination module 51, in determining the distance between the point to be estimated and the first template line in the two-dimensional coordinate system of longitudinal and transverse wave impedances, is used for:

[0073] The distance between the point to be estimated and the first template line in the two-dimensional coordinate system of longitudinal and transverse wave impedance is determined using the following formula:

[0074]

[0075] Among them, D i Let I be the distance between the point to be estimated with index i and the first template line. pi and ISi Let G be the longitudinal wave impedance and transverse wave impedance of the point to be estimated, respectively; G be the slope of the template line in the coordinate system; and P1 be the intercept of the first template line.

[0076] In some embodiments, the mud content determination module 52 determines the mud content of the point to be estimated by means of a distance proportional to the difference in mud content, for the following purposes:

[0077] The mud content V of the point to be estimated is determined by the following formula. shi :

[0078]

[0079] In the formula, D i D is the distance between the point to be estimated with index i and the first template line. s This is the standard distance between the first template line and the second template line.

[0080] In some embodiments, the above apparatus further includes a template line creation module 53, configured to: pre-create the first template line and the second template line in the following manner:

[0081] Based on the first sample set, a rock physics model characterizing the relationship between clay content, porosity, and P-wave and S-wave velocities is established through multiple linear regression. Each sample in the first sample set contains clay content, porosity, P-wave velocity, and S-wave velocity at the same location. Based on the second sample set, a third sample set is obtained using the conversion relationship between P-wave and S-wave velocity and P-wave and S-wave impedance, and the rock physics model. Each sample in the second sample set contains clay content and porosity at the same location, and each sample in the third sample set contains clay content, porosity, P-wave impedance, and S-wave impedance at the same location. A first-character sample set with a clay content of a first set value and a second-character sample set with a clay content of a second set value are extracted from the third sample set. The extracted first-character sample set and second-character sample set are projected into the coordinate system to obtain the first template line corresponding to the first sub-sample set and the second template line corresponding to the second sub-sample set.

[0082] In some embodiments, the template line establishment module 53, which establishes a rock physics model characterizing the relationship between clay content, porosity, and P-wave and S-wave velocities based on a first sample set through multiple linear regression, is used for:

[0083] Based on the first sample set, a rock physics model characterizing the relationship between clay content, porosity, and P-wave and S-wave velocities is established through multiple linear regression as follows:

[0084] V P =AB·φ-C·V sh

[0085] VS =DE·φ-F·V sh

[0086] Among them, V P V is the longitudinal wave velocity. S V is the transverse wave velocity, φ is the porosity, and V is the transverse wave velocity. sh denoted as mud content, and A, B, C, D, E, and F are fitting constants.

[0087] In some embodiments, the template line establishment module 53, which utilizes the conversion relationship between P-wave and S-wave impedance and the rock physics model to obtain a third sample set, is used for:

[0088] Using the aforementioned rock physics model, the third sample set is obtained through the following equation:

[0089] I P =V P ·ρ

[0090] I S =V S ·ρ

[0091] ρ=ρ ma *(1-V sh -φ)+ρ sh *V sh +ρ f *φ

[0092] Where, ρ ma ρ sh and ρ f Let ρ be the density of pure sandstone, pure mudstone, and fluid, respectively, where ρ is the reservoir density, and I is the density of the fluid. S For transverse wave impedance, I P This is the longitudinal wave impedance.

[0093] Regarding the apparatus in the above embodiments, the specific manner in which each module performs its operation has been described in detail in the embodiments related to the method, and will not be elaborated upon here.

[0094] Based on the inventive concept of the present invention, embodiments of the present invention also provide a computer storage medium storing computer-executable instructions, which, when executed by a processor, implement the above-described method for determining mud content based on seismic elastic parameters.

[0095] Based on the inventive concept of this invention, this embodiment of the invention also provides a server, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the above-described method for determining mud content based on seismic elastic parameters.

[0096] Unless otherwise specifically stated, terms such as processing, calculation, operation, determination, display, etc., may refer to the actions and / or processes of one or more processing or computing systems or similar devices that represent the manipulation and conversion of data representing physical (e.g., electronic) quantities within the registers or memory of the processing system into other data similarly representing physical quantities within the memory, registers, or other such information storage, transmission, or display devices of the processing system. Information and signals can be represented using any of a variety of different techniques and methods. For example, data, instructions, commands, information, signals, bits, symbols, and chips mentioned throughout the above description can be represented by voltage, current, electromagnetic waves, magnetic fields or particles, light fields or particles, or any combination thereof.

[0097] It should be understood that the specific order or hierarchy of steps in the disclosed process is an example of an exemplary method. Based on design preferences, it should be understood that the specific order or hierarchy of steps in the process may be rearranged without departing from the scope of this disclosure. The appended method claims provide elements of various steps in an exemplary order and are not intended to limit the scope to the specific order or hierarchy described.

[0098] In the detailed description above, various features are combined together in a single embodiment to simplify this disclosure. This approach to disclosure should not be construed as reflecting an intention that embodiments of the claimed subject matter require more features than those stated in each claim. Rather, as reflected in the appended claims, the invention is presented with fewer features than all of the features in a single disclosed embodiment. Therefore, the appended claims are hereby clearly incorporated into the detailed description, wherein each claim stands alone as a preferred embodiment of the invention.

[0099] Those skilled in the art will also understand that the various illustrative logic blocks, modules, circuits, and algorithm steps described in conjunction with the embodiments herein can be implemented as electronic hardware, computer software, or a combination thereof. To clearly illustrate the interchangeability between hardware and software, the various illustrative components, blocks, modules, circuits, and steps described above are generally described in terms of their functionality. Whether such functionality is implemented as hardware or software depends on the specific application and the design constraints imposed on the overall system. Those skilled in the art can implement the described functionality in alternative ways for each specific application; however, such implementation decisions should not be construed as departing from the scope of this disclosure.

[0100] The steps of the methods or algorithms described in conjunction with the embodiments herein can be directly embodied in hardware, software modules executed by a processor, or a combination thereof. The software modules can reside in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disks, removable disks, CD-ROMs, or any other form of storage medium well known in the art. An exemplary storage medium is connected to the processor, enabling the processor to read information from and write information to the storage medium. Of course, the storage medium can also be a component of the processor. The processor and storage medium can reside in an ASIC. The ASIC can reside in a user terminal. Alternatively, the processor and storage medium can exist as discrete components in the user terminal.

[0101] For software implementation, the techniques described in this application can be implemented using modules (e.g., procedures, functions, etc.) that perform the functions described in this application. This software code can be stored in memory units and executed by a processor. The memory units can be implemented within the processor or outside the processor; in the latter case, they are communicatively coupled to the processor via various means, as is well known in the art.

[0102] The foregoing description includes examples of one or more embodiments. It is certainly impossible to describe all possible combinations of components or methods in order to describe the above embodiments, but those skilled in the art will recognize that further combinations and arrangements of the various embodiments are possible. Therefore, the embodiments described herein are intended to cover all such changes, modifications, and variations that fall within the scope of the appended claims. Furthermore, the term “comprising” as used in the specification or claims is interpreted in a manner similar to the term “including,” as it is understood when used as a conjunction in the claims. Additionally, the use of any term “or” in the specification of the claims is intended to mean “non-exclusive or.” The terms “first” and “second” are used for descriptive purposes only and should not be construed as indicating or implying relative importance.

Claims

1. A method for determining clay content based on seismic elastic parameters, characterized in that, include: Based on the longitudinal and transverse wave impedances of the point to be estimated, determine the distance between the point to be estimated and the first template line in the two-dimensional coordinate system of longitudinal and transverse wave impedances. Based on the distance, the standard distance between the first template line and the second template line, and the difference between the mud content corresponding to the first template line and the mud content corresponding to the second template line, the mud content of the point to be estimated is determined by proportionally dividing the distance by the difference in mud content. The first template line and the second template line are two parallel straight lines in the coordinate system. The mud content of each point on the first template line is the same, and the mud content of each point on the second template line is the same. The first template line and the second template line are pre-established in the following manner: Based on the first sample set, a rock physics model characterizing the relationship between clay content, porosity, and P-wave and S-wave velocities is established through multiple linear regression. Each sample in the first sample set contains clay content, porosity, P-wave velocity, and S-wave velocity at the same location. Based on the second sample set, a third sample set is obtained using the conversion relationship between P-wave and S-wave velocity and P-wave and S-wave impedance, and the rock physics model. Each sample in the second sample set contains clay content and porosity at the same location, and each sample in the third sample set contains clay content, porosity, P-wave impedance, and S-wave impedance at the same location. A first-character sample set with a clay content of a first set value and a second-character sample set with a clay content of a second set value are extracted from the third sample set. The extracted first-character sample set and second-character sample set are projected into the coordinate system to obtain the first template line corresponding to the first sub-sample set and the second template line corresponding to the second sub-sample set. The rock physics model is as follows: （1）； （2）； In equations (1) and (2), For the longitudinal wave velocity, For transverse wave velocity, Porosity The mud content is represented by A, B, C, D, E, and F, which are fitting constants. The third sample set is: （3）； （4）； （5）； In equations (3)-(5), , and The densities of pure sandstone, pure mudstone, and fluid are given respectively. For reservoir density, For transverse wave impedance, This is the longitudinal wave impedance.

2. The method according to claim 1, characterized in that, The coordinate system has longitudinal wave impedance as the horizontal axis and transverse wave impedance as the vertical axis.

3. The method according to claim 2, characterized in that, Determining the distance between the point to be estimated and the first template line in the two-dimensional coordinate system of longitudinal and transverse wave impedance includes: The distance between the point to be estimated and the first template line in the two-dimensional coordinate system of longitudinal and transverse wave impedance is determined using the following formula (6): = （6）； In equation (6), Let i be the distance between the point to be estimated (i) and the first template line. and These are the longitudinal wave impedance and transverse wave impedance of the point to be estimated, respectively. The slope of the template line in the coordinate system is given. The intercept of the first template line is denoted as .

4. The method according to claim 2, characterized in that, The standard distance between the first template line and the second template line is predetermined by the following formula (7): （7）； In equation (7), This is the standard distance between the first template line and the second template line. The slope of the template line in the coordinate system is given. and These are the intercepts of the first template line and the second template line, respectively.

5. The method according to claim 1, characterized in that, The mud content at each point on the first template line is 0, and the mud content at each point on the second template line is 1.

6. The method according to claim 5, characterized in that, The method of determining the mud content of the point to be estimated by proportionally varying the distance to the difference in mud content includes: The mud content of the point to be estimated is determined by the following formula (8). : （8）； In equation (8), Let i be the distance between the point to be estimated (i) and the first template line. This is the standard distance between the first template line and the second template line.

7. A device for determining clay content based on seismic elastic parameters, characterized in that, The apparatus is used to perform the method for determining clay content based on seismic elastic parameters as described in claim 1, and the apparatus comprises: The distance determination module is used to determine the distance between the point to be estimated and the first template line in the two-dimensional coordinate system of the longitudinal and transverse wave impedances, based on the longitudinal and transverse wave impedances of the point to be estimated. The mud content determination module is used to determine the mud content of the point to be estimated by means of the distance, the standard distance between the first template line and the second template line, and the difference between the mud content corresponding to the first template line and the mud content corresponding to the second template line, in a proportional manner to the difference in mud content. The first template line and the second template line are two parallel straight lines in the coordinate system. The mud content of each point on the first template line is the same, and the mud content of each point on the second template line is the same.

8. A computer storage medium, characterized in that, The computer storage medium stores computer-executable instructions, which, when executed by a processor, implement the method for determining mud content based on seismic elastic parameters as described in any one of claims 1 to 6.

9. A server, characterized in that, include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the program, implements the method for determining clay content based on seismic elastic parameters as described in any one of claims 1 to 6.

Images