Method for obtaining conversion relationship between dynamic and static elastic parameters

Abstract

A method for obtaining a conversion relationship between dynamic and static elastic parameters includes: Step S1, acquiring horizontal cores at different depths of the destination formation; Step S2, measuring the dynamic elastic parameters of the horizontal core under different pressures; Step S3, measuring the static elastic parameters of the horizontal core under different pressures; Step S4, measuring the clay content of the horizontal core; Step S5 establishing a function relationship of the ratio between the dynamic and static elastic parameters with the formation pressure and clay content; and completing the conversion between the dynamic and static elastic parameters. The technical solution provided by the present invention takes full account of the influence of the formation stress and the clay content on the conversion rule of dynamic and static elastic parameters and is of great significance for improving the logging evaluation accuracy of rock mechanical parameters.

Claims

1. A method for obtaining a conversion relationship between dynamic and static elastic parameters, comprising the following steps: Step S1, acquiring horizontal cores at different depths of the destination formation; Step S2, measuring with a device the dynamic elastic parameters of the horizontal cores under different effective formation stresses; Step S3, measuring with a device the static elastic parameters of the horizontal cores under different effective formation stresses; Step S4, measuring with a device the clay contents of the horizontal cores; Step S5, establishing a function relationship of the ratio between the dynamic and static elastic parameters with the effective formation stress and clay content; and completing the conversion between the dynamic and static elastic parameters; and Step S6, using the function relationship to determine, with a control unit, rock mechanical parameters of the core.

2. The method according to claim 1, wherein in step S1, the acquiring horizontal cores at different depths in the destination formation includes the following steps: Step S101, m pieces of cores are obtained by drilling in the horizontal direction at the same formation depth D1, denoted as D.sub.11, D.sub.12, D.sub.13, . . . , and D.sub.1m, respectively; Step S102, each of the m pieces of cores drilled at the formation depth D1 is divided into two smaller pieces, denoted as D.sub.11A, D.sub.11B, D.sub.12A, D.sub.12B, D.sub.13A, D.sub.13B, . . . , D.sub.1mA, and D.sub.1mB, respectively; Step S103: following the operation of Step S101 and Step S102, m pieces of cores at the formation depth Dn are obtained by drilling, with each of the m pieces of cores at Dn being divided into two smaller pieces, denoted as D.sub.n1A, D.sub.n1B, D.sub.n2A, D.sub.n2B, D.sub.n3A, D.sub.n3B, . . . , D.sub.nmA, and D.sub.nmB, respectively.

3. The method according to claim 2, wherein in step S101, the core is at least 8 cm in length and 2.5-5 cm in diameter.

4. The method according to claim 2, wherein in step S102, when each core is divided into two smaller pieces, each core of the two smaller pieces is at least 4 cm in length and 2.5-5 cm in diameter.

5. The method according to claim 2, wherein in step S2, the measuring the dynamic elastic parameters of the horizontal cores under different effective formation stresses includes the following steps: Step S201, based on the pressure of the destination formation, m effective formation stress values are set in an ascending order, denoted as P1, P2, P3, . . . , and Pm, respectively; wherein the effective formation stress interval [P1, Pm] includes the pressure of the entire destination formation; Step S202, under the effective formation stress of P1, the obtained horizontal cores D.sub.11A, D.sub.21A, D.sub.31A, D.sub.n1A are measured and the dynamic elastic parameters thereof are calculated and respectively denoted as Ed.sub.11, Ed.sub.21, Ed.sub.31, . . . , and Ed.sub.n1; following the above operation, the obtained horizontal cores D.sub.1mA, D.sub.2mA, D.sub.3mA, . . . , D.sub.nmA are measured and the dynamic elastic parameters thereof are calculated and respectively denoted as Ed.sub.1m, Ed.sub.2m, Ed.sub.3m, . . . , and Ed.sub.nm.

6. The method according to claim 5, wherein in Step 202, the measured items include a density, a longitudinal wave velocity and a shear wave velocity.

7. The method according to claim 5, wherein in Step S3, the measuring the static elastic parameters of the horizontal core under different effective formation stresses includes the following steps: Step S301: based on the pressure of the destination formation, m effective formation stress values are set in an ascending order, denoted as P1, P2, P3, . . . , Pm respectively; wherein the effective formation stress interval [P1, Pm] includes the effective formation stress of the entire destination formation; Step S302, under the effective formation stress of P1, the obtained horizontal cores D.sub.11B, D.sub.21B, D.sub.31B, D.sub.n1B are measured and the static elastic parameters thereof are calculated and respectively denoted as Es.sub.11, Es.sub.21, Es.sub.31, . . . , and Es.sub.n1; following the above operation, under the effective formation stress of Pm, the obtained horizontal cores D.sub.1mB, D.sub.2mB, D.sub.3mB, D.sub.nmB are measured and the static elastic parameters thereof are calculated and respectively denoted as Es.sub.1m, Es.sub.2m, Es.sub.3m, . . . , and Es.sub.nm.

8. The method according to claim 7, wherein in Step S302, the measured items include stress and strain.

9. The method according to claim 8, wherein, in Step S4, the measuring the clay content of the horizontal core includes the following steps: Step S401, the horizontal cores D.sub.11B, D.sub.12B, D.sub.13B, . . . , and D.sub.1mB which have been measured for the static elastic parameters are subjected to X-ray diffraction measurement to obtain the clay content Vcl1 at the formation depth D1; Step S402, following the operation of Step S401, the horizontal cores D.sub.n1B, D.sub.n2B, D.sub.n3B, . . . , and D.sub.nmB which have been measured for the static elastic parameters are subjected to X-ray diffraction measurement to obtain the clay content Vcln at the formation depth Dn.

10. The method according to claim 8, wherein, in Step S5, the establishing a function relationship of the ratio between the dynamic and static elastic parameters with the effective formation stress and clay content includes the following steps: Step S501: the ratios between the dynamic and static elastic parameters of the core sample at the formation depth D1 under the m effective formation stress values among P1-Pm are calculated, denoted as E.sub.d11/E.sub.s11, E.sub.d12/E.sub.s12, E.sub.d13/E.sub.s13, . . . , and E.sub.d1m/E.sub.s1m, respectively; a relationship of the ratio between the dynamic and static elastic parameters at the formation depth of D1 with the effective formation stress is formulated; following the above operation, the ratios between the dynamic and static elastic parameters of the core sample at the formation depth Dn under the m effective formation stress values among P1-Pm are calculated, denoted as E.sub.dn1/E.sub.sn1, E.sub.dn2/E.sub.sn2, E.sub.dn3/E.sub.sn3, . . . , and E.sub.dnm/E.sub.snm, respectively; a relationship of the ratio between the dynamic and static elastic parameters at the formation depth of Dn with the effective formation stress is formulated; Step S502, when the effective formation stress value is P1, the ratios between the dynamic and static elastic parameters of the core sample at the n formation depths among D1-Dn are calculated, denoted as E.sub.d11/E.sub.s11, E.sub.d21/E.sub.s21, E.sub.d31/E.sub.s31, . . . , and E.sub.dn1/E.sub.sn1, respectively; according to the obtained clay contents Vcl1, Vcl2, Vcl3, . . . , and Vcln at the n formation depths among D1-Dn, a relationship of the ratio between the dynamic and static elastic parameters with the clay content under the effective formation stress value of P1 is formulated; following the above operation, when the effective formation stress value is Pm, the ratios between the dynamic and static elastic parameters of the core sample at the n formation depths among D1-Dn are calculated, denoted as E.sub.d1m/E.sub.s1m, E.sub.d2m/E.sub.s2m, E.sub.d3m/E.sub.s3m, . . . , and E.sub.dnm/E.sub.snm, respectively; according to the obtained clay contents Vcl1, Vcl2, Vcl3, . . . , and Vcln at the n formation depths among D1-Dn, a relationship of the ratio between the dynamic and static elastic parameters with the clay content under the effective formation stress value of Pm is formulated; Step S503: according to the relationships obtained in step S501 and step S502, a function relationship of the ratio between the dynamic and static elastic parameters with the effective formation stress and clay content E.sub.d/E.sub.s=f(P,Vcl) is established.

11. The method according to claim 5, wherein in Step S302, the measured items include stress and strain.

12. The method according to claim 5, wherein, in Step S4, the measuring the clay content of the horizontal core includes the following steps: Step S401, the horizontal cores D.sub.11B, D.sub.12B, D.sub.13B, . . . , and D.sub.1mB which have been measured for the static elastic parameters are subjected to X-ray diffraction measurement to obtain the clay content Vcl1 at the formation depth D1; Step S402, following the operation of Step S401, the horizontal cores D.sub.n1B, D.sub.n2B, D.sub.n3B, . . . , and D.sub.nmB which have been measured for the static elastic parameters are subjected to X-ray diffraction measurement to obtain the clay content Vcln at the formation depth Dn.

13. The method according to claim 5, wherein, in Step S5, the establishing a function relationship of the ratio between the dynamic and static elastic parameters with the effective formation stress and clay content includes the following steps: Step S501: the ratios between the dynamic and static elastic parameters of the core sample at the formation depth D1 under the m effective formation stress values among P1-Pm are calculated, denoted as E.sub.d11/E.sub.s11, E.sub.d12/E.sub.s12, E.sub.d13/E.sub.s13, . . . , and E.sub.d1m/E.sub.s1m, respectively; a relationship of the ratio between the dynamic and static elastic parameters at the formation depth of D1 with the effective formation stress is formulated; following the above operation, the ratios between the dynamic and static elastic parameters of the core sample at the formation depth Dn under the m effective formation stress values among P1-Pm are calculated, denoted as E.sub.dn1/E.sub.sn1, E.sub.dn2/E.sub.sn2, E.sub.dn3/E.sub.sn3, . . . , E.sub.dnm/E.sub.snm, respectively; a relationship of the ratio between the dynamic and static elastic parameters at the formation depth of Dn with the effective formation stress is formulated; Step S502, when the effective formation stress value is P1, the ratios between the dynamic and static elastic parameters of the core sample at the n formation depths among D1-Dn are calculated, denoted as E.sub.d11/E.sub.s11, E.sub.d21/E.sub.s21, E.sub.d31/E.sub.s31, . . . , and E.sub.dn1/E.sub.sn1, respectively; according to the obtained clay contents Vcl1, Vcl2, Vcl3, . . . , and Vcln at the n formation depths among D1-Dn, a relationship of the ratio between the dynamic and static elastic parameters with the clay content under the effective formation stress value of P1 is formulated; following the above operation, when the effective formation stress value is Pm, the ratios between the dynamic and static elastic parameters of the core sample at the n formation depths among D1-Dn are calculated, denoted as E.sub.d1m/E.sub.s1m, E.sub.d2m/E.sub.s2m, E.sub.d3m/E.sub.s3m, . . . , and E.sub.dnm/E.sub.snm, respectively; according to the obtained clay contents Vcl1, Vcl2, Vcl3, . . . , and Vcln at the n formation depths among D1-Dn, a relationship of the ratio between the dynamic and static elastic parameters with the clay content under the effective formation stress value of Pm is formulated; Step S503: according to the relationships obtained in Step S501 and Step S502, a function relationship of the ratio between the dynamic and static elastic parameters with the effective formation stress and clay content E.sub.d/E.sub.s=f(P,Vcl) is established.

14. The method according to claim 2, wherein in Step S3, the measuring the static elastic parameters of the horizontal cores under different effective formation stresses includes the following steps: Step S301: based on the pressure of the destination formation, m effective formation stress values are set in an ascending order, denoted as P1, P2, P3, . . . , and Pm respectively; wherein the effective formation stress interval [P1, Pm] includes the effective formation stress of the entire destination formation; Step S302, under the effective formation stress of P1, the obtained horizontal cores D.sub.11B, D.sub.21B, D.sub.31B, . . . D.sub.n1B are measured and the static elastic parameters thereof are calculated and respectively denoted as Es.sub.11, Es.sub.21, Es.sub.31, . . . , and Es.sub.n1; following the above operation, under the effective formation stress of Pm, the obtained horizontal cores D.sub.1mB, D.sub.2mB, D.sub.3mB, D.sub.nmB are measured and the static elastic parameters thereof are calculated and respectively denoted as Es.sub.1m, Es.sub.2m, Es.sub.3m, . . . , and Es.sub.nm.

15. The method according to claim 14, wherein in Step S302, the measured items include stress and strain.

16. The method according to claim 14, wherein, in Step S4, the measuring the clay content of the horizontal core includes the following steps: Step S401, the horizontal cores D.sub.11B, D.sub.12B, D.sub.13B, . . . , and D.sub.1mB which have been measured for the static elastic parameters are subjected to X-ray diffraction measurement to obtain the clay content Vcl1 at the formation depth D1; Step S402, following the operation of Step S401, the horizontal cores D.sub.n1B, D.sub.n2B, D.sub.n3B, . . . , and D.sub.nmB which have been measured for the static elastic parameters are subjected to X-ray diffraction measurement to obtain the clay content Vcln at the formation depth Dn.

17. The method according to claim 14, wherein, in Step S5, the establishing a function relationship of the ratio between the dynamic and static elastic parameters with the effective formation stress and clay content includes the following steps: Step S501: the ratio between the dynamic and static elastic parameters of the core sample at the formation depth D1 under the m effective formation stress values among P1-Pm is calculated, denoted as E.sub.d11/E.sub.s11, E.sub.d12/E.sub.s12, E.sub.d13/E.sub.s13, . . . , and E.sub.d1m/E.sub.s1m, respectively; a relationship of the ratio between the dynamic and static elastic parameters at the formation depth of D1 with the effective formation stress is formulated; Following the above operation, the ratio between the dynamic and static elastic parameters of the core sample at the formation depth Dn under the m effective formation stress values among P1-Pm is calculated, denoted as E.sub.dn1/E.sub.sn1, E.sub.dn2/E.sub.sn2, E.sub.dn3/E.sub.sn3, . . . , E.sub.dnm/E.sub.snm, respectively; a relationship of the ratio between the dynamic and static elastic parameters at the formation depth of Dn with the effective formation stress is formulated; Step S502, when the effective formation stress value is P1, the ratio between the dynamic and static elastic parameters of the core sample at n formation depths in D1-Dn is calculated, denoted as, E.sub.d11/E.sub.s11, E.sub.d21/E.sub.s21, E.sub.d31/E.sub.s31, . . . , and E.sub.dn1/E.sub.sn1 respectively; according to the obtained clay contents Vcl1, Vcl2, Vcl3, . . . , and Vcln of the n formation depths in D1-Dn, a relationship of the ratio between the dynamic and static elastic parameters with the clay content at the effective formation stress value of P1 is formulated; Following the above operation, when the effective formation stress value is Pm, the ratio between the dynamic and static elastic parameters of the core sample at n formation depths in D1-Dn is calculated, denoted as, E.sub.d1m/E.sub.s1m, E.sub.d2m/E.sub.s2m, E.sub.d3m/E.sub.s3m, . . . , and E.sub.dnm/E.sub.snm, respectively; according to the obtained clay contents Vcl1, Vcl2, Vcl3, . . . , and Vcln of the n formation depths in D1-Dn, a relationship of the ratio between the dynamic and static elastic parameters with the clay content at the effective formation stress value of Pm is formulated; Step S503: according to the relationship obtained in step S501 and step S502, a function relationship of the ratio between the dynamic and static elastic parameters with the effective formation stress and clay content E.sub.d/E.sub.s=f(P,Vcl) is established.

18. The method according to claim 2, wherein in Step S4, the measuring the clay contents of the horizontal cores includes the following steps: Step S401, the horizontal cores D.sub.11B, D.sub.12B, D.sub.13B, . . . , and D.sub.1mB which have been measured for the static elastic parameters are subjected to X-ray diffraction measurement to obtain the clay content Vcl1 at the formation depth D1; Step S402, following the operation of Step S401, the horizontal cores D.sub.n1B, D.sub.n2B, D.sub.n3B, . . . , and D.sub.nmB which have been measured for the static elastic parameters are subjected to X-ray diffraction measurement to obtain the clay content Vcln at the formation depth Dn.

19. The method according to claim 18, wherein, in Step S5, the establishing a function relationship of the ratio between the dynamic and static elastic parameters with the effective formation stress and clay content includes the following steps: Step S501: the ratios between the dynamic and static elastic parameters of the core sample at the formation depth D1 under the m effective formation stress values among P1-Pm are calculated, denoted as E.sub.d11/E.sub.s11, E.sub.d12/E.sub.s12, E.sub.d13/E.sub.s13, . . . , and E.sub.d1m/E.sub.s1m, respectively; a relationship of the ratio between the dynamic and static elastic parameters at the formation depth of D1 with the effective formation stress is formulated; following the above operation, the ratios between the dynamic and static elastic parameters of the core sample at the formation depth Dn under the m effective formation stress values among P1-Pm are calculated, denoted as E.sub.dn1/E.sub.sn1, E.sub.dn2/E.sub.sn2, E.sub.dn3/E.sub.sn3, . . . , E.sub.dnm/E.sub.snm, respectively; a relationship of the ratio between the dynamic and static elastic parameters at the formation depth of Dn with the effective formation stress is formulated; Step S502, when the effective formation stress value is P1, the ratios between the dynamic and static elastic parameters of the core sample at the n formation depths among D1-Dn are calculated, denoted as E.sub.d11/E.sub.s11, E.sub.d21/E.sub.s21, E.sub.d31/E.sub.s31, . . . , and E.sub.dn1/E.sub.sn1, respectively; according to the obtained clay contents Vcl1, Vcl2, Vcl3, . . . , and Vcln at the n formation depths among D1-Dn, a relationship of the ratio between the dynamic and static elastic parameters with the clay content under the effective formation stress value of P1 is formulated; following the above operation, when the effective formation stress value is Pm, the ratios between the dynamic and static elastic parameters of the core sample at the n formation depths among D1-Dn are calculated, denoted as E.sub.d1m/E.sub.s1m, E.sub.d2m/E.sub.s2m, E.sub.d3m/E.sub.s3m, . . . , and E.sub.dnm/E.sub.snm, respectively; according to the obtained clay contents Vcl1, Vcl2, Vcl3, . . . , and Vcln at the n formation depths among D1-Dn, a relationship of the ratio between the dynamic and static elastic parameters with the clay content under the effective formation stress value of Pm is formulated; Step S503: according to the relationships obtained in step S501 and step S502, a function relationship of the ratio between the dynamic and static elastic parameters with the effective formation stress and clay content E.sub.d/E.sub.s=f(P,Vcl) is established.

20. The method according to claim 2, wherein, in Step S5, the establishing a function relationship of the ratio between the dynamic and static elastic parameters with the effective formation stress and clay content includes the following steps: Step S501: the ratios between the dynamic and static elastic parameters of the core sample at the formation depth D1 under the m effective formation stress values among P1-Pm are calculated, denoted as E.sub.d11/E.sub.s11, E.sub.d12/E.sub.s12, E.sub.d13/E.sub.s13, . . . , and E.sub.d1m/E.sub.s1m, respectively; a relationship of the ratio between the dynamic and static elastic parameters at the formation depth of D1 with the effective formation stress is formulated; following the above operation, the ratios between the dynamic and static elastic parameters of the core sample at the formation depth Dn under the m effective formation stress values among P1-Pm are calculated, denoted as E.sub.dn1/E.sub.sn1, E.sub.dn2/E.sub.sn2, E.sub.dn3/E.sub.sn3, . . . , and E.sub.dnm/E.sub.snm, respectively; a relationship of the ratio between the dynamic and static elastic parameters at the formation depth of Dn with the effective formation stress is formulated; Step S502, when the effective formation stress value is P1, the ratios between the dynamic and static elastic parameters of the core sample at the n formation depths among D1-Dn are calculated, denoted as E.sub.d11/E.sub.s11, E.sub.d21/E.sub.s21, E.sub.d31/E.sub.s31, . . . , and E.sub.dn1/E.sub.sn1, respectively; according to the obtained clay contents Vcl1, Vcl2, Vcl3, . . . , and Vcln at the n formation depths among D1-Dn, a relationship of the ratio between the dynamic and static elastic parameters with the clay content under the effective formation stress value of P1 is formulated; following the above operation, when the effective formation stress value is Pm, the ratios between the dynamic and static elastic parameters of the core sample at the n formation depths among D1-Dn are calculated, denoted as E.sub.d1m/E.sub.s1m, E.sub.d2m/E.sub.s2m, E.sub.d3m/E.sub.s3m, . . . , and E.sub.dnm/E.sub.snm, respectively; according to the obtained clay contents Vcl1, Vcl2, Vcl3, . . . , and Vcln at the n formation depths among D1-Dn, a relationship of the ratio between the dynamic and static elastic parameters with the clay content under the effective formation stress value of Pm is formulated; Step S503: according to the relationships obtained in Step S501 and Step S502, a function relationship of the ratio between the dynamic and static elastic parameters with the effective formation stress and clay content E.sub.d/E.sub.s=f(P,Vcl) is established.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a schematic flow chart of a method for converting dynamic and static elastic parameters according to an example of the present invention;

(2) FIG. 2 is a schematic structural view of an apparatus for measuring the longitudinal and transverse wave velocities of a core according to an example of the present invention;

(3) FIG. 3 is a graph of the relationship of the ratio between the dynamic and static elastic parameters with formation pressure according to an example of the present invention;

(4) FIG. 4 is a graph of the relationship of the ratio between the dynamic and static elastic parameters with the clay content when the formation pressure is 20 Mpa according to an example of the present invention;

(5) FIG. 5 is a graph of the relationship of the ratio between the dynamic and static elastic parameters with the clay content when the formation pressure is 40 Mpa according to an example of the present invention;

(6) FIG. 6 is a comparison diagram of the practical effects according to an example of the present invention;

DESCRIPTION OF THE KEY REFERENCE NUMBERS IN THE DRAWINGS

(7) 1. displacement meter; 2: first probe; 3: second probe; 4. measuring tube; 5. pore fluid control unit; 6: temperature control unit; 7: confining pressure control unit; 8: data acquisition and analysis unit.

DETAILED DESCRIPTION OF THE INVENTION

(8) For better understanding of the technical features, objects and advantageous effects of the present invention, the technical solutions of the present invention are now described in details below, which cannot be construed as limiting to the scope of the embodiments of the present invention.

EXAMPLES

(9) This example provides a method for obtaining a conversion relationship between dynamic and static elastic parameters (the flow of which is shown in FIG. 1), which includes the following steps:

(10) Step S1, acquiring horizontal cores at different depths of the destination formation:

(11) Horizontal plunger-like cores were obtained by drilling at different depths in the exploration area, with a set consisted of m=4 horizontal small cores obtained by drilling at the same formation depth; A total of n=8 sets of cores were obtained, denoted as D.sub.11, D.sub.12, D.sub.13, D.sub.14; D.sub.21, D.sub.22, D.sub.23, D.sub.24; . . . ; and D.sub.81, D.sub.82, D.sub.83, D.sub.84; wherein, each horizontal core was 2.5 cm in diameter and 8-10 cm in length.

(12) The above cores were cut, respectively, with each divided into two smaller pieces and the cutting surface parallel to the bottom surface of the plunger-like core, and the location of cutting point must ensure that each of the two pieces after cutting was 4-5 cm in length; the cores satisfying such conditions after cutting were respectively denoted as D.sub.11A, D.sub.11B, D.sub.12A, D.sub.12B, . . . , D.sub.14A, D.sub.14B; D.sub.21A, D.sub.21B, D.sub.22A, D.sub.22B, . . . , D.sub.24A, D.sub.24B; . . . , and D.sub.81A, D.sub.81B, D.sub.82A, D.sub.82B, . . . , D.sub.84A, D.sub.84B; wherein, D.sub.11A, D.sub.11B indicated that after D.sub.11 was cut into two smaller pieces, one of them was denoted as D.sub.11A, and the other D.sub.11B, respectively; and so on.

(13) Step S2, measuring the dynamic elastic parameters of the horizontal cores under different pressures:

(14) Based on the formation pressure, four formation pressure values were set in an ascending order, denoted as P1, P2, P3, and P4 (unit: MPa), wherein, P1=20 MPa, P2=30 MPa, P3=40 MPa, and P4=60 MPa; the pressure range of 20 MPa-60 MPa encompassed the pressure across the entire exploration area.

(15) Under the condition of a pressure of P1=20 MPa, the horizontal cores D.sub.11A, D.sub.21A, D.sub.31A, . . . , and D.sub.81A obtained in step S1 were subjected to density and longitudinal and shear wave velocity measurements under saturated brine conditions (for the longitudinal and transverse wave velocity measurements, reference was made to SY/T 6351-2012).

(16) The device for measuring the longitudinal and shear wave velocities of the core was not particularly limited and could be a conventional measuring device in the art, such as the device shown in FIG. 2 (which was one of the conventional measurement devices in the art). The device included a displacement meter 1, a first probe 2, a second probe 3, a measuring tube 4, a pore fluid control unit 5, a temperature control unit 6, a confining pressure control unit 7, and a data acquisition and analysis unit 8. During the measurement, the measure conditions were regulated by the pore fluid control unit 5, the temperature control unit 6, and the confining pressure control unit 7. In this device, the displacement meter 1 could perform deformation measurement; the probes at both ends of the sample could allow the fluid passing through on one hand and allow energy being transmitted and received on the other hand, while the data resulted from the measurements could be sent to the data acquisition and analysis unit 8. This measuring device was only an exemplary device of the embodiment of the present invention. Other measuring devices having the same function or capable of fulfilling the same purpose could also be used for the measurement of the longitudinal and shear wave velocities of cores.

(17) The measurement results are shown in Table 1, with their dynamic Young's moduli Ed.sub.11, Ed.sub.21, Ed.sub.31, . . . , and Ed.sub.81 calculated; wherein, the dynamic Young's moduli were calculated with an equation as shown in Equation 1:

(18) $\begin{array}{cc}E=\frac{\rho V_s^2\left(3V_p^2-4V_s^2\right)}{V_p^2-V_2^2}& \mathrm{Equation}1\end{array}$

(19) In Equation 1, E is dynamic Young's modulus, KPa; ρ is density, g/cm.sup.3, V.sub.p is longitudinal wave velocity, m/s, V.sub.s is shear wave velocity, m/s,

(20) TABLE-US-00001 TABLE 1 Experimental shear Dynamic confining Longitudinal wave Young's Sample pressure Density wave velocity velocity modulus No. (MPa) g/cm.sup.3 (m/s) (m/s) (GPa) D.sub.11A 20 2.65 5663 3426 75.35 D.sub.12A 30 2.65 5704 3456 76.59 D.sub.13A 40 2.65 5728 3472 77.29 D.sub.14A 60 2.65 5742 3484 77.76 D.sub.21A 20 2.64 5641 3400 74.14 D.sub.22A 30 2.64 5674 3418 74.96 D.sub.23A 40 2.64 5688 3430 75.43 D.sub.24A 60 2.64 5697 3437 75.72 D.sub.31A 20 2.61 5186 2905 56.01 D.sub.32A 30 2.61 5243 2954 57.74 D.sub.33A 40 2.61 5272 2980 58.66 D.sub.34A 60 2.61 5309 3009 59.71 D.sub.41A 20 2.56 5138 2914 54.91 D.sub.42A 30 2.56 5225 2977 57.16 D.sub.43A 40 2.56 5279 3014 58.53 D.sub.44A 60 2.56 5314 3041 59.49 D.sub.51A 20 2.57 5045 2885 53.78 D.sub.52A 30 2.57 5087 2912 54.75 D.sub.53A 40 2.57 5110 2926 55.27 D.sub.54A 60 2.57 5126 2934 55.59 D.sub.61A 20 2.54 5048 2970 55.36 D.sub.62A 30 2.54 5164 3051 58.24 D.sub.63A 40 2.54 5224 3092 59.75 D.sub.64A 60 2.54 5273 3125 60.98 D.sub.71A 20 2.54 5030 2977 55.39 D.sub.72A 30 2.54 5108 3042 57.59 D.sub.73A 40 2.54 5163 3086 59.12 D.sub.74A 60 2.54 5201 3111 60.05 D.sub.81A 20 2.48 4764 2653 44.53 D.sub.82A 30 2.48 4841 2718 46.53 D.sub.83A 40 2.48 4895 2757 47.79 D.sub.84A 60 2.48 4935 2782 48.63

(21) Then, under the condition of the pressure of P2=30 MPa, the horizontal cores D.sub.12A, D.sub.22A, D.sub.32A, . . . , and D.sub.82A obtained in Step S1 were subjected to density and longitudinal and shear wave velocity measurements under saturated brine conditions, and their dynamic Young's moduli Ed.sub.12, Ed.sub.22, Ed.sub.32, . . . , and Ed.sub.82 were calculated.

(22) Similarly, under the condition of the pressure of P4=60 MPa, the horizontal cores D.sub.14A, D.sub.24A, D.sub.34A, . . . , and D.sub.84A obtained in the step S1 were subjected to density and longitudinal and shear wave velocity measurements under saturated brine conditions, and their dynamic Young's moduli Ed.sub.14, Ed.sub.24, Ed.sub.34, . . . , and Ed.sub.84 were calculated.

(23) Step S3, measuring the static elastic parameters of the horizontal cores under different pressures:

(24) Under the condition of the pressure of P1=20 MPa, the horizontal cores D.sub.11B, D.sub.21B, D.sub.31B, . . . , and D.sub.81B obtained in Step S1 were subjected to stress and strain measurements under saturated saline conditions (for the stress and strain measurements, reference was made to GBT23561.9-2009), and their static Young's moduli were calculated according to the measurement results (as shown in Table 2); wherein the static Young's moduli were calculated with an equation as shown in Equation 2:

(25) $\begin{array}{cc}E=\frac{{\sigma }_{\mathrm{ab}}}{{\mathrm{.Math.}}_{\mathrm{ab}}}& \mathrm{Equation}2\end{array}$

(26) In Equation 2:

(27) σ.sub.ab—the stress difference between the end point and starting point of the straight line segment of the curve of stress vs. axial strain, in megapascals (MPa);

(28) ε.sub.ab—the strain difference between the end point and starting point of the straight line segment of the curve of stress vs. axial strain, in percentage.

(29) Similarly, under the conditions of pressure P2=30, 40 and 60 MPa, the horizontal cores D.sub.12B, D.sub.22B, D.sub.32B, . . . , D.sub.82B, . . . , D.sub.14B, D.sub.24B, D.sub.34B, . . . , D.sub.84B were subjected to stress and strain measurements under saturated saline conditions, and their static Young's moduli were calculated.

(30) TABLE-US-00002 TABLE 2 Experimental Static confining Stress Strain Young's Sample pressure difference difference modulus No. (MPa) σ.sub.ab (MPa) ε.sub.ab (%) (GPa) D.sub.11B 20 106.1708 0.2095 50.6 D.sub.12B 30 104.3521 0.1904 54.9 D.sub.13B 40 157.5117 0.2544 61.8 D.sub.14B 60 141.2113 0.2132 66.2 D.sub.21B 20 134.5789 0.2666 50.5 D.sub.22B 30 119.5818 0.2121 56.3 D.sub.23B 40 163.8272 0.2579 63.6 D.sub.24B 60 227.4188 0.3252 70 D.sub.31B 20 99.0613 0.2811 35.3 D.sub.32B 30 80.4307 0.2162 37.2 D.sub.33B 40 112.8493 0.291 38.8 D.sub.34B 60 107.6294 0.263 40.9 D.sub.41B 20 63.0874 0.1874 33.7 D.sub.42B 30 68.5192 0.1944 35.3 D.sub.43B 40 83.4911 0.2185 38.2 D.sub.44B 60 97.2092 0.2513 38.7 D.sub.51B 20 91.0993 0.275 33.1 D.sub.52B 30 94.6616 0.279 33.9 D.sub.53B 40 121.111 0.347 34.9 D.sub.54B 60 128.4008 0.358 35.9 D.sub.61b 20 91.6754 0.2513 36.5 D.sub.62B 30 88.6596 0.2226 39.8 D.sub.63B 40 97.4503 0.2286 42.6 D.sub.64B 60 101.1177 0.2189 46.1 D.sub.71B 20 76.6523 0.2076 36.9 D.sub.72B 30 96.5619 0.2431 39.7 D.sub.73B 40 107.8026 0.238 45.3 D.sub.74B 60 108.8003 0.2155 50.5 D.sub.81B 20 51.78 0.1859 27.9 D.sub.82B 30 90.3 0.2986 30.2 D.sub.83B 40 141.08 0.441 32 D.sub.84B 60 187.13 0.6021 31.1

(31) Step S4, measuring the clay contents of the horizontal cores;

(32) The divided samples of the core D.sub.11B, D.sub.12B, . . . , D.sub.14B after the stress and strain measurement were collected and subjected to X-ray diffraction measurement, and the obtained clay content values represented the clay contents at the formation depth D1;

(33) Similarly, the divided samples of the horizontal cores at other formation depths after the stress and strain measurement were respectively collected and subjected to X-ray diffraction measurement, and the obtained clay content values represented the clay contents at the corresponding formation depths, as shown in Table 3.

(34) TABLE-US-00003 TABLE 3 Sample No. Total clay content D1 3.5 D2 2 D3 19.8 D4 20.8 D5 16.7 D6 12.5 D7 11.8 D8 14.2

(35) Step S5, establishing a function relationship of the ratio between the dynamic and static elastic parameters with the formation pressure and clay contents:

(36) On the basis of the measured results at different formation depths under different pressure values, a relationship of the ratio between the dynamic and static elastic parameters (the ratio between the dynamic and static elastic parameters here referred to the ratio between the dynamic and static Young's modulus) with the formation pressure was established (as shown in FIG. 3). Meanwhile, a relationship of the ratio between the dynamic and static elastic parameters under a fixed pressure with the clay content was established in connection with the results of clay content measurements at various formation depths (as shown in FIG. 4 and FIG. 5 which demonstrated the relationship of the ratio between the dynamic and static elastic parameters and the clay content when the formation pressure was 20 MPa and 40 MPa respectively).

(37) According to the relationship of the ratio between the dynamic and static elastic parameters with the formation pressure, together with the relationship thereof with the clay content under different formation pressures, a function relationship of the ratio between the dynamic and static elastic parameters with the clay content and formation pressure was established, as shown in Equation 2:

(38) $\begin{array}{cc}\frac{E_d}{E_s}=\left(A+B.Math.V_{\mathrm{clay}}\right).Math.P_{\mathrm{eff}}^{C.Math.V_{\mathrm{clay}}+D}& \mathrm{Equation}2\end{array}$

(39) In Equation 2, E.sub.d is dynamic elastic parameter (also referred to as dynamic Young's modulus), E.sub.s is static elastic parameter (also referred to as static Young's modulus), V.sub.clay is clay content, P.sub.eff is formation pressure, and A, B, C, and D are empirical coefficients obtained from experimental results. In this example, A=3.5, B=−8.14, C=1.35, and D=−0.3.

(40) FIG. 6 is a comparison diagram of the effect of application in a well in the exploration area. Here, a total of 11 pieces of information was included, with Esta in the 9.sup.th representing the static Young's modulus calculated by the method provided in the example of the present invention, Esta2 representing the static Young's modulus calculated by a traditional dynamic and static linear regression method, and Esta_Core representing the static Young's modulus obtained from core measurements. The comparison results suggested that the results calculated by the method provided in the example of the present invention were closest to the results from the core measurements. The average relative error at three data points was 18%, as oppose to the average relative error from the conventional method of up to 45%. The statistical results from seven data points in three wells in the area demonstrated that the relative error of the method provided by the present invention was 13% in comparison with the results from the core measurements, as oppose to the relative error from the conventional method of up to 28%.

(41) The foregoing describes only the preferred embodiments of the present invention, but is not intended to limit the present invention. Various modifications and changes can be made to the embodiments of the present invention for those skilled in the art. Any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention are intended to be included in the protection scope of the present invention.