CN112987096A - Method for calculating sound velocity of high-argillaceous sandstone - Google Patents

Method for calculating sound velocity of high-argillaceous sandstone Download PDF

Info

Publication number
CN112987096A
CN112987096A CN202110289819.4A CN202110289819A CN112987096A CN 112987096 A CN112987096 A CN 112987096A CN 202110289819 A CN202110289819 A CN 202110289819A CN 112987096 A CN112987096 A CN 112987096A
Authority
CN
China
Prior art keywords
argillaceous
modulus
equivalent
sandstone
content
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110289819.4A
Other languages
Chinese (zh)
Other versions
CN112987096B (en
Inventor
韩学辉
江佳洋
韩振鹏
刘智桢
罗振
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dongying Yongsheng Energy Technology Co ltd
Karamay Onco Energy Technology Co ltd
China University of Petroleum East China
Original Assignee
Dongying Yongsheng Energy Technology Co ltd
Karamay Onco Energy Technology Co ltd
China University of Petroleum East China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dongying Yongsheng Energy Technology Co ltd, Karamay Onco Energy Technology Co ltd, China University of Petroleum East China filed Critical Dongying Yongsheng Energy Technology Co ltd
Priority to CN202110289819.4A priority Critical patent/CN112987096B/en
Publication of CN112987096A publication Critical patent/CN112987096A/en
Application granted granted Critical
Publication of CN112987096B publication Critical patent/CN112987096B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention discloses a method for calculating the sound velocity of high-mud sandstone, which comprises the following steps: s1, determining a mud distribution mode in the high-mud sandstone; s2 quantitatively distinguishing the content of the cemented shale; s3, calculating the content of the structural argillaceous substance and the content of the dispersed argillaceous substance; s4 calculating the equivalent skeleton modulus of quartz and cemented argillaceous based on the corrected CCT model; s5, replacing a part of quartz with the structural argillaceous substance as a skeleton by using a DEM model, and calculating the equivalent modulus of the structural argillaceous substance after replacing the skeleton; s6, filling the dispersed argillaceous substances into the pores based on the unconsolidated sandstone model, and calculating to obtain the equivalent elastic modulus of the high-argillaceous sandstone under the constraint of the structural argillaceous substances, the cemented argillaceous substances and the dispersed argillaceous substances; s7 calculating the velocity V of longitudinal wavep1Velocity V of sum transverse waves1Compared with the prior art, the method for calculating the sound velocity of the high-argillaceous sandstone has higher accuracy.

Description

Method for calculating sound velocity of high-argillaceous sandstone
Technical Field
The invention belongs to the technical field of petroleum and gas exploration, and relates to a method for calculating the sound velocity of high-argillaceous sandstone.
Background
Argillaceous sandstone reservoirs are one of the most important reservoirs. The sound velocity is one of the most important rock physical characteristic parameters of the argillaceous sandstone, can reflect the properties of a framework and pore fluid, and is an important physical property for evaluating parameters such as the porosity and the saturation of a reservoir by earthquake and well logging.
The high-mud sandstone contains mud in various distribution forms: cementing, dispersing, and structuring muds, which have different effects on the acoustic velocity properties of the muddy sandstone, need to be distinguished when studying their acoustic properties.
In the existing cemented sandstone model with the constraint of cemented shale content, the shale is assumed to exist only in a pore space, and is divided into two types of cemented shale and dispersed shale, wherein the cemented shale plays a role in cementing to strengthen a framework, the dispersed shale plays a role in filling to reduce the porosity, and the two types of shale can increase the sound velocity of the argillaceous sandstone.
A method for distinguishing between cemented and dispersed argillaceous bodies has been proposed by observing the contact relationship and relative distribution between argillaceous bodies and particles in cast body sheets, such as korean scholar et al: the continuously distributed argillaceous matter in contact with two or more particles is a cemented argillaceous matter; based on the criterion, the cemented shale content of the artificial argillaceous sandstone is estimated based on a pixel pickup method, and the cemented shale content is used as an input parameter of a cemented sandstone model to optimize a CCT model, compared with an original model, the sound speed error of the method is reduced by 20%, and the prediction accuracy is remarkably improved (Han Chi, Nie Jun, Guo Jun Xin and the like 2020. quantitative estimation of the contact cemented shale in the argillaceous sandstone and influence on sandstone elasticity. geophysical report 63 (4): 1654-. However, the above model ignores the influence of structural mud on the sound velocity, and cannot explain the phenomenon that the sound velocity shows a decreasing trend along with the increase of the content of the mud, so that the model is not suitable for calculating the sound velocity of the high-mud sandstone any more.
Therefore, it is necessary to establish a method for calculating the sound velocity of the high-shale-content sandstone under the dual constraints of the cemented shale content and the structural shale content so as to accurately forward predict the sound velocity of the high-shale-content sandstone and provide a basis for applying earthquake and logging sound velocity to invert parameters such as porosity and the like.
Disclosure of Invention
The invention aims to provide a method for calculating the sound velocity of high-argillaceous sandstone, which is characterized in that the argillaceous of the high-argillaceous sandstone is divided into three parts of cemented argillaceous, dispersed argillaceous and structural argillaceous, and the concepts of cemented argillaceous for cementing, dispersed argillaceous for filling and structural argillaceous for replacing a quartz framework are introduced into the calculation of the modulus of the high-argillaceous sandstone on the basis of an improved CCT theory, a differential equivalent medium model and an unconsolidated sandstone model, so that the method for calculating the sound velocity of the high-argillaceous sandstone in different argillaceous distribution modes is established.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method for calculating the sound velocity of high-argillaceous sandstone comprises the following steps:
s1, determining a mud distribution mode in the high-mud sandstone by using the sample slice identification result;
s2, quantitatively distinguishing the content V of the cemented shale by using a dyeing mark-pixel pickup methodcement clay；
S3, calculating the content V of the structural argillaceous substance by using a substance balance methodstructure clayAnd dispersed argillaceous content Vdispersed clay；
S4, calculating the equivalent volume modulus K of the quartz and the cemented argillaceous based on the corrected CCT modelcctAnd equivalent shear modulus Gcct；
S5, replacing partial quartz with the structural argillaceous substance as a skeleton by utilizing the DEM model, and calculating the equivalent volume modulus K of the structural argillaceous substance after replacing the skeletoncctdemAnd equivalent shear modulus Gcctdem；
S6, filling the dispersed argillaceous substances into the pores based on the unconsolidated sandstone model, and calculating to obtain the equivalent bulk elastic modulus K of the high-argillaceous sandstone under the constraint of the structural argillaceous substances, the cemented argillaceous substances and the dispersed argillaceous substanceseffAnd equivalent modulus of elasticity in shear Geff；
S7, calculating the longitudinal wave velocity V of the high-mud sandstone modulus obtained in the step S6p1Velocity V of sum transverse waves1。
Preferably, the argillaceous sandstone in the step S1 has an argillaceous content of 30-50%.
Preferably, the mud distribution manner in step S1 is cementing mud, dispersing mud and structural mud.
Preferably, the step S1 further includes measuring the porosity of the high argillaceous sandstone
Figure BDA0002976529960000036
Density rho, longitudinal wave velocity VpTransverse wave velocity VsAnd a mud content Vsh。
Preferably, the structural argillaceous content V in step S3structure clayThe calculation formula of (2) is as follows:
Vg+Vstructure clay+Φ0＝1 (1)
wherein, VgIs the content of quartz skeleton, Vstructure clayIs structured in terms of mud content, phi0Is the critical porosity.
Preferably, the dispersed argillaceous V of step S3dispersed clayThe calculation formula of (2) is as follows:
Vsh＝Vdispersed clay+Vstructure clay+Vcement clay (2)
wherein, VshIs the total argillaceous content.
Preferably, step S4 equivalent bulk modulus K of said quartz and cemented sludgecctAnd equivalent shear modulus GcctThe calculation formula of (2) is as follows:
Figure BDA0002976529960000031
Figure BDA0002976529960000032
wherein,
Figure BDA0002976529960000033
cement distribution mode 1:
Figure BDA0002976529960000034
cement distribution mode 2:
Figure BDA0002976529960000035
Figure BDA0002976529960000041
Figure BDA0002976529960000042
in the formula, alpha1The ratio of the radius of the cementing plane to the radius of the sandstone particles is the normalized cementing radius;
ε is the ratio of the center thickness of the cement to the radius of the particles, i.e., the normalized center thickness of the cement;
Φ0critical porosity for loose sandstone;
phi is the true porosity of the loose sandstone;
Gcand G is the shear modulus of the shale cement and the quartz particles respectively;
vcand v is the poisson's ratio of the cement and the quartz particles, respectively;
cement distribution mode 1 is a case where the cement is distributed only between particles; cement distribution pattern 2 is the case when the cement is uniformly distributed on the particle surface.
Preferably, the equivalent modulus of the structural argillaceous matrix in step S5 after replacing the skeleton is: equivalent bulk modulus KcctdemAnd equivalent shear modulus GcctdemThe calculation formula of (2) is as follows:
Figure BDA0002976529960000043
Figure BDA0002976529960000044
in the formula, K*(0)＝K1，μ*(0)＝μ1As an initial condition for the coupled differential equation;
K*(0)＝K1，μ*(0)＝μ1as an initial condition for the coupled differential equation, K1，μ1Bulk and shear moduli for the background medium (i.e., the first phase, equivalent bulk and equivalent shear moduli for quartz and cemented muds); k2，μ2The volume modulus and shear modulus (i.e., the second phase, bulk modulus and shear modulus of the structural argillaceous mass) of the gradually added inclusion object; p and Q are geometric factors used to characterize the geometry of the filler, according to the Minear study, the structural argillaceous aspect ratio, replacing quartz as skeleton, is 1.0, the superscript 2 of P and Q means that this geometric factor is aimed at having an equivalent modulus K*And mu*(i.e. K)cctdemAnd Gcctdem) The main phase background medium of (2).
Preferably, the modulus of the high argillaceous sandstone in step S6 is: equivalent bulk modulus of elasticity KeffAnd equivalent modulus of elasticity in shear GeffThe calculation formula of (2) is as follows:
Figure BDA0002976529960000051
Figure BDA0002976529960000052
in the formula, KeffAnd GeffEquivalent bulk modulus and equivalent shear modulus at a specific porosity, respectively;
Φ is the existing porosity;
Φ1subtracting the cementitious shale content from the critical porosity;
Kcctdemand GcctdemRespectively an equivalent bulk modulus and an equivalent shear modulus calculated by using a CCT model and a DEM model,
wherein the cement content is the content of cemented shale, i.e.
Figure BDA0002976529960000053
KsAnd GsThe equivalent bulk modulus and the equivalent shear modulus of quartz and structural argillas are respectively obtained through calculation by a Voigt-Reuss-Hill model.
The invention also provides application of the calculation method in forward modeling prediction of the acoustic velocity of the high-shale-content sandstone and inversion estimation of the porosity of the high-shale-content sandstone through earthquake and well logging.
The invention has the beneficial effects that:
the method comprises the steps of dividing the argillaceous sandstone into cemented argillaceous, dispersed argillaceous and structural argillaceous, and introducing concepts of cementing argillaceous for cementing, filling dispersed argillaceous for filling and replacing a quartz framework with structural argillaceous in the calculation of the modulus of the high argillaceous sandstone based on an improved CCT theory, a differential equivalent medium model and an unconsolidated sandstone model.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Fig. 2 is a slice experiment result of argillaceous sandstone.
Figure 3 is a graph of identifying cementitious shale content using IPP software.
Figure 4 is a model of basal cemented sandstone loose sand.
Figure 5 is a diagram of a process for solving the modulus of high argillaceous sandstones.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention are further described below with reference to the accompanying drawings.
As shown in fig. 1-5, the invention provides a method for calculating the sound velocity of high-argillaceous sandstone, which comprises the following steps:
s1, determining a argillaceous distribution mode in the high-argillaceous sandstone according to the sample slice identification result, wherein as shown in FIG. 2, when the argillaceous content is small, the argillaceous is mainly distributed among particles, and the cemented argillaceous is in a dominant effect; as the content of the sludge continues to rise, the sludge is mainly distributed in the pore space, and the dispersed sludge plays a leading role at the moment; when the mud content reaches a certain critical point, as the mud content continues to rise, the mud is mainly distributed in the quartz framework, and the mud in the structure plays a leading role at the moment;
s2, for the sandstone with high mud content, assuming the framework of the sandstone to be cemented is the one with the porosity of phi0Equivalent spherical particles with the coordination number C of 0.36 and the average coordination number C of 9 are closely and randomly arranged, and the content V of the structural argillaceous substance is calculated according to a substance balance equationstructure clayThe specific calculation formula is as follows:
Vg+Vstructure clay+Φ0＝1 (1)
wherein, VgIs the content of quartz skeleton, Vstructure clayIs structured in terms of mud content, phi0Critical porosity;
then, the same artificial argillaceous sandstone casting body slice image is displayedQuantitative discrimination of cementitious shale content V using a dye mark-pixel pickup methodcement clayAs shown in fig. 3, the thick black line segment is cemented mud;
thirdly, calculating the dispersed argillaceous V according to a material balance equationdispersed clayThe specific calculation formula is as follows:
Vsh＝Vdispersed clay+Vstructure clay+Vcement clay (2)
wherein, VshIs the total argillaceous content.
S3, for the argillaceous sandstone sample, firstly considering the cementation effect of quartz and cemented argillaceous, namely adding the cemented argillaceous into the quartz to reduce the porosity and increase the effective modulus of the particle aggregate (skeleton), so that the equivalent bulk modulus K of the quartz and the cemented argillaceous is calculated by utilizing a continuous cementation theory (CCT model)cctAnd equivalent shear modulus GcctThe specific calculation formula is shown as follows;
Figure BDA0002976529960000071
Figure BDA0002976529960000072
Figure BDA0002976529960000073
cement distribution mode 1:
Figure BDA0002976529960000074
cement distribution mode 2:
Figure BDA0002976529960000075
Figure BDA0002976529960000076
Figure BDA0002976529960000081
in the formula, alpha1The ratio of the radius of the cementing plane to the radius of the sandstone particles is the normalized cementing radius;
ε is the ratio of the center thickness of the cement to the radius of the particles, i.e., the normalized center thickness of the cement;
Φ0critical porosity for loose sandstone;
phi is the true porosity of the loose sandstone;
Gcand G is the shear modulus of the shale cement and the quartz particles respectively;
vcand v is the poisson's ratio of the cement and the quartz particles, respectively;
wherein, the cement distribution mode is shown in figure 4: cement distribution mode 1 is a case where the cement content is small and is distributed only between particles; the cement distribution mode 2 is a case where the cement content is large and the cement is uniformly distributed on the surface of the particles, and refers to cement distribution modes 1 and 2 in a substrate-type cemented loose sandstone acoustic velocity correction model based on a CCT model (korea scholar, university of china, proceedings of petroleum university (nature science edition), vol.37, No. 4, 2013);
s4, performing skeleton replacement by using a differential equivalent medium theory (DEM model), and calculating the equivalent modulus of the structure mud after the skeleton is replaced: equivalent bulk modulus KcctdemAnd equivalent shear modulus GcctdemThe specific calculation formula is as follows:
Figure BDA0002976529960000082
Figure BDA0002976529960000083
in the formula, K*(0)＝K1，μ*(0)＝μ1As an initial condition for the coupled differential equation;
K*(0)＝K1，μ*(0)＝μ1as an initial condition for the coupled differential equation, K1，μ1Bulk and shear moduli for the background medium (i.e., the first phase, equivalent bulk and equivalent shear moduli for quartz and cemented muds); k2，μ2The volume modulus and shear modulus (i.e., the second phase, bulk modulus and shear modulus of the structural argillaceous mass) of the gradually added inclusion object;
p and Q are geometric factors used to characterize the geometry of the filler, according to the Minear study, the structural argillaceous aspect ratio, replacing quartz as skeleton, is 1.0, the superscript 2 of P and Q means that this geometric factor is aimed at having an equivalent modulus K*And mu*(i.e. K)cctdemAnd Gcctdem) The main phase background medium of (2).
S5, filling the dispersed argillaceous substances into the pores based on the unconsolidated sandstone model, and calculating to obtain the equivalent bulk elastic modulus K of the high-argillaceous sandstone under the constraint of the structural argillaceous substances, the cemented argillaceous substances and the dispersed argillaceous substanceseffAnd equivalent modulus of elasticity in shear GeffThe calculation formula is as follows:
Figure BDA0002976529960000091
Figure BDA0002976529960000092
in the formula, KeffAnd GeffEquivalent bulk modulus and equivalent shear modulus at a specific porosity, respectively;
Φ is the existing porosity;
Φ1subtracting the cementitious shale content from the critical porosity;
Kcctdemand GcctdemRespectively an equivalent bulk modulus and an equivalent shear modulus calculated by using a CCT model and a DEM model,
wherein the cement content is the content of cemented shale, i.e.
Figure BDA0002976529960000093
KsAnd GsThe specific solving process diagram of the equivalent skeleton bulk modulus and the equivalent shear modulus of quartz and structural argillaceous sand respectively obtained by calculation through a Voigt-Reuss-Hill model is shown in FIG. 5.
S6, calculating the longitudinal wave velocity V by using the modulus of the high-mud sandstone obtained in the step S5p1Velocity V of sum transverse waves1。
Example 1
The method for calculating the sound velocity of the artificial high-mud sandstone comprises the following steps:
s1, selecting 4 artificial sandstone samples with high mud content distribution of 30-50%, determining the distribution form of mud in the samples to be three types of cemented mud, dispersed mud and structural mud according to the identification result of the casting slice (as shown in figure 2), and measuring the porosity of the samples
Figure BDA0002976529960000101
Density rho, longitudinal wave velocity VpTransverse wave velocity VsAnd a argillaceous content Vsh；
S2, observing a sample slice under a mirror to determine a argillaceous distribution mode, adding a contact argillaceous distribution line by using Image processing software, wherein a thick black line is cemented argillaceous distributed among particles to play a cementing role, picking up pixels of the cemented argillaceous by using Image-pro-plus software as shown in FIG. 3, and determining the content V of the cemented argillaceouscement clay；
S3, determining the structural argillaceous content V by using a material balance methodstructure clayAnd dispersed argillaceous content Vdispersed clay；
The content of the structural argillaceous component in step S3Vstructure clayThe calculation formula of (2) is as follows:
Vg+Vstructure clay+Φ0＝1 (1)
wherein, VgIs the content of quartz skeleton, Vstructure clayIs structured in terms of mud content, phi0Is the critical porosity.
Preferably, the dispersed argillaceous V of step S3dispersed clayThe calculation formula of (2) is as follows:
Vsh＝Vdispersed clay+Vstructure clay+Vcement clay (2)
wherein, VshThe total mud content is;
s4, calculating to obtain equivalent bulk modulus K of the quartz and the cemented shale based on the improved CCT modelcctAnd equivalent shear modulus Gcct；
Equivalent framework modulus of the quartz and cemented mud: equivalent bulk modulus KcctAnd equivalent shear modulus GcctThe calculation formula of (2) is as follows:
Figure BDA0002976529960000111
Figure BDA0002976529960000112
wherein,
Figure BDA0002976529960000113
cement distribution mode 1:
Figure BDA0002976529960000114
cement distribution mode 2:
Figure BDA0002976529960000115
Figure BDA0002976529960000116
Figure BDA0002976529960000117
in the formula, alpha1The ratio of the radius of the cementing plane to the radius of the sandstone particles is the normalized cementing radius;
ε is the ratio of the center thickness of the cement to the radius of the particles, i.e., the normalized center thickness of the cement;
Φ0critical porosity for loose sandstone;
phi is the true porosity of the loose sandstone;
Gcand G is the shear modulus of the shale cement and the quartz particles respectively;
vcand v is the poisson's ratio of the cement and the quartz particles, respectively;
wherein, the cement distribution mode is shown in figure 4: cement distribution mode 1 is a case where the cement content is small and is distributed only between particles; the cement distribution mode 2 is a case where the cement content is large and the cement is uniformly distributed on the particle surface;
s5, calculating the equivalent modulus of the cemented sandstone after the structural argillaceous body replaces part of quartz to serve as a framework by utilizing a DEM model: equivalent bulk modulus KcctdemAnd equivalent shear modulus Gcctdem；
Equivalent bulk modulus K of the structure argillaceous substituted skeleton in step S5cctdemAnd equivalent shear modulus GcctdemThe calculation formula of (2) is as follows:
Figure BDA0002976529960000121
Figure BDA0002976529960000122
in the formula, K*(0)＝K1，μ*(0)＝μ1As an initial condition for the coupled differential equation;
K*(0)＝K1，μ*(0)＝μ1as an initial condition for the coupled differential equation, K1，μ1Bulk and shear moduli for the background medium (i.e., the first phase, equivalent bulk and equivalent shear moduli for quartz and cemented muds); k2，μ2The volume modulus and shear modulus (i.e., the second phase, bulk modulus and shear modulus of the structural argillaceous mass) of the gradually added inclusion object;
p and Q are geometric factors used to characterize the geometry of the filler, according to the Minear study, the structural argillaceous aspect ratio, replacing quartz as skeleton, is 1.0, the superscript 2 of P and Q means that this geometric factor is aimed at having an equivalent modulus K*And mu*(i.e. K)cctdemAnd Gcctdem) The inclusion material 2 in the primary phase background medium;
s6, filling the dispersed argillaceous substances into the pores based on the unconsolidated sandstone model, and calculating to obtain the equivalent bulk elastic modulus K of the high-argillaceous sandstone under the constraint of the structural argillaceous substances, the cemented argillaceous substances and the dispersed argillaceous substanceseffAnd equivalent modulus of elasticity in shear Geff；
The equivalent bulk modulus K of the high argillaceous sandstone in the step S6effAnd equivalent modulus of elasticity in shear GeffThe calculation formula of (2) is as follows:
Figure BDA0002976529960000123
Figure BDA0002976529960000124
Figure BDA0002976529960000131
in the formula, KeffAnd GeffAre respectively specific toEquivalent bulk modulus and equivalent shear modulus of elasticity at porosity;
Φ is the existing porosity;
Φ1subtracting the cementitious shale content from the critical porosity;
Kcctdemand GcctdemRespectively an equivalent bulk modulus and an equivalent shear modulus calculated by using a CCT model and a DEM model,
wherein the cement content is the content of cemented shale, i.e.
Figure BDA0002976529960000134
KsAnd GsThe specific solving process diagram of the modulus of the high-mud sandstone is shown in figure 5, wherein the specific solving process diagram is respectively the equivalent skeleton bulk modulus and the equivalent shear modulus of quartz and structural mud obtained by calculation through a Voigt-Reuss-Hill model;
s7, calculating the longitudinal wave velocity V by using the modulus of the high-mud sandstone obtained in the step S6p1Velocity V of sum transverse waves1。
The calculation results are shown in table 1.
Table 1 example 1 calculation results
Figure BDA0002976529960000132
Note that the mud content V in the surfaceshPorosity of the porous material
Figure BDA0002976529960000133
The content of the cementing mud and the content of the dispersed mud are both decimal.
As shown in Table 1, the longitudinal and transverse wave velocities calculated by applying the new model are consistent with the experimental measurement results, wherein the relative prediction error of the longitudinal wave velocity is less than 5%, and the relative prediction error of the transverse wave velocity is less than 15%.
It should be emphasized that the embodiments described herein are illustrative rather than restrictive, and thus the present invention is not limited to the embodiments described in the detailed description, but also includes other embodiments that can be derived from the technical solutions of the present invention by those skilled in the art.

Claims (10)

1. A method for calculating the sound velocity of high-argillaceous sandstone comprises the following steps:
s1, determining a mud distribution mode in the high-mud sandstone by using the sample slice identification result;
s2, quantitatively distinguishing the content V of the cemented shale by using a dyeing mark-pixel pickup methodcementclay；
S3, calculating the content V of the structural argillaceous substance by using a substance balance methodstructureclayAnd dispersed argillaceous content Vdispersedclay；
S4, calculating the equivalent volume modulus K of the quartz and the cemented argillaceous based on the corrected CCT modelcctAnd equivalent shear modulus Gcct；
S5, replacing partial quartz with the structural argillaceous substance as a skeleton by utilizing the DEM model, and calculating the equivalent volume modulus K of the structural argillaceous substance after replacing the skeletoncctdemAnd equivalent shear modulus Gcctdem；
S6, filling the dispersed argillaceous substances into the pores based on the unconsolidated sandstone model, and calculating to obtain the equivalent bulk elastic modulus K of the high-argillaceous sandstone under the constraint of the structural argillaceous substances, the cemented argillaceous substances and the dispersed argillaceous substanceseffAnd equivalent modulus of elasticity in shear Geff；
S7, calculating the modulus of the high-mud sandstone obtained in the step S6, and calculating the longitudinal wave velocity Vp1Velocity V of sum transverse waves1。
2. The calculation method according to claim 1, wherein the argillaceous sandstone of high and medium argillaceous content in step S1 has an argillaceous content of 30-50%.
3. The calculation method according to claim 1, wherein the argillaceous distribution manner in step S1 is cemented argillaceous, dispersed argillaceous, and structured argillaceous.
4. According to the claimsThe calculation method of claim 1, wherein the step S1 further comprises measuring the porosity of the high argillaceous sandstone
Figure FDA0002976529950000011
Density rho, longitudinal wave velocity VpTransverse wave velocity VsAnd a mud content Vsh。
5. The calculation method according to claim 1, wherein the structural argillaceous content V in step S3structureclayThe calculation formula of (2) is as follows:
Vg+Vstructureclay+Φ0＝1 (1)
wherein, VgIs the content of quartz skeleton, VstructureclayIs structured in terms of mud content, phi0Is the critical porosity.
6. The calculation method according to claim 1, wherein the dispersed argillaceous V in step S3dispersedclayThe calculation formula of (2) is as follows:
Vsh＝Vdispersedclay+Vstructureclay+Vcementclay (2)
wherein, VshIs the total argillaceous content.
7. The method of claim 1, wherein the equivalent bulk modulus K of the quartz and cemented sludge in step S4cctAnd equivalent shear modulus GcctThe calculation formula of (2) is as follows:
Figure FDA0002976529950000021
Figure FDA0002976529950000022
wherein,
Figure FDA0002976529950000023
cement distribution mode 1:
Figure FDA0002976529950000024
cement distribution mode 2:
Figure FDA0002976529950000025
Figure FDA0002976529950000026
Figure FDA0002976529950000027
in the formula,
α1the ratio of the radius of the cementing plane to the radius of the sandstone particles;
ε is the ratio of the center thickness of the cement to the particle radius;
Φ0critical porosity for loose sandstone;
phi is the true porosity of the loose sandstone;
Gcand G is the shear modulus of the shale cement and the quartz particles respectively;
vcand v is the poisson's ratio of the cement and the quartz particles, respectively;
cement distribution mode 1 is a case where the cement is distributed only between particles;
cement distribution pattern 2 is the case when the cement is uniformly distributed on the particle surface.
8. The method according to claim 1, wherein the equivalent bulk modulus K of the structural argillaceous device after skeleton replacement in step S5cctdemAnd equivalent shear modulus GcctdemThe calculation formula of (2) is as follows:
Figure FDA0002976529950000031
Figure FDA0002976529950000032
in the formula,
y is the volume fraction of the structural argillaceous material;
K*(0)＝K1，μ*(0)＝μ1as an initial condition for the coupled differential equation, K1，μ1Bulk and shear moduli for background media; k2，μ2The bulk modulus and shear modulus of the gradually added inclusion bodies;
p and Q are geometric factors used to characterize the geometry of the filler, according to the Minear study, the structural argillaceous aspect ratio, replacing quartz as skeleton, is 1.0, the superscript 2 of P and Q means that this geometric factor is aimed at having an equivalent modulus K*And mu*The main phase background medium of (2).
9. The method of claim 1, wherein the equivalent bulk modulus K of the high argillaceous sandstone in step S6effAnd equivalent modulus of elasticity in shear GeffThe calculation formula of (2) is as follows:
Figure FDA0002976529950000041
Figure FDA0002976529950000042
in the formula,
Keffand GeffEquivalent bulk modulus and equivalent shear modulus at a specific porosity, respectively;
Φ is the existing porosity;
Φ1subtracting the cementitious shale content from the critical porosity;
Kcctdemand GcctdemRespectively an equivalent bulk modulus and an equivalent shear modulus calculated by using a CCT model and a DEM model,
wherein the content of the cementing material is the content of cementing mud,
Figure FDA0002976529950000043
Ksand GsRespectively calculating the equivalent skeleton bulk modulus and shear modulus of quartz and structural argillaceous substances obtained through a Voigt-reus-Hill model.
10. Use of the calculation method according to any one of claims 1 to 9 for forward prediction of the acoustic velocity of high-mudcontent sandstone and/or for seismic and well logging inversion estimation of the porosity of high-mudcontent sandstone.
CN202110289819.4A 2021-03-15 2021-03-15 Method for calculating sound velocity of high-argillaceous sandstone Active CN112987096B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110289819.4A CN112987096B (en) 2021-03-15 2021-03-15 Method for calculating sound velocity of high-argillaceous sandstone

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110289819.4A CN112987096B (en) 2021-03-15 2021-03-15 Method for calculating sound velocity of high-argillaceous sandstone

Publications (2)

Publication Number Publication Date
CN112987096A true CN112987096A (en) 2021-06-18
CN112987096B CN112987096B (en) 2021-11-26

Family

ID=76333758

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110289819.4A Active CN112987096B (en) 2021-03-15 2021-03-15 Method for calculating sound velocity of high-argillaceous sandstone

Country Status (1)

Country Link
CN (1) CN112987096B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114414027A (en) * 2021-12-16 2022-04-29 东营市永昇能源科技有限责任公司 Method for calculating sound velocity of supporting conglomerate and non-supporting conglomerate

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130229892A1 (en) * 2010-09-14 2013-09-05 Statoil Asa Method of predicting the pressure sensitivity of seismic velocity within reservoir rocks
CN104570064A (en) * 2013-10-09 2015-04-29 中国石油化工股份有限公司 Method for calculating shear wave velocity of sandstone formation
CN108399270A (en) * 2017-02-08 2018-08-14 中国石油化工股份有限公司 A kind of method of anisotropy shale ratio in determining shale formation
CN109884696A (en) * 2019-03-20 2019-06-14 中国石油化工股份有限公司 Sandstone Computational approach of S-wave velocity based on calcite cementation function influence
CN110133722A (en) * 2019-06-06 2019-08-16 克拉玛依市昂科能源科技有限公司 A kind of application velocity of sound qualitative recognition mud distribution type new method
CN112489735A (en) * 2020-11-30 2021-03-12 中国海洋石油集团有限公司 Carbonate cement content calculating method based on logging and linear particle swarm optimization

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130229892A1 (en) * 2010-09-14 2013-09-05 Statoil Asa Method of predicting the pressure sensitivity of seismic velocity within reservoir rocks
CN104570064A (en) * 2013-10-09 2015-04-29 中国石油化工股份有限公司 Method for calculating shear wave velocity of sandstone formation
CN108399270A (en) * 2017-02-08 2018-08-14 中国石油化工股份有限公司 A kind of method of anisotropy shale ratio in determining shale formation
CN109884696A (en) * 2019-03-20 2019-06-14 中国石油化工股份有限公司 Sandstone Computational approach of S-wave velocity based on calcite cementation function influence
CN110133722A (en) * 2019-06-06 2019-08-16 克拉玛依市昂科能源科技有限公司 A kind of application velocity of sound qualitative recognition mud distribution type new method
CN112489735A (en) * 2020-11-30 2021-03-12 中国海洋石油集团有限公司 Carbonate cement content calculating method based on logging and linear particle swarm optimization

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
韩学辉: "基于CCT模型的基底式胶结疏松砂岩声波速度修正模型", 《中国石油大学学报（自然科学版）》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114414027A (en) * 2021-12-16 2022-04-29 东营市永昇能源科技有限责任公司 Method for calculating sound velocity of supporting conglomerate and non-supporting conglomerate

Also Published As

Publication number Publication date
CN112987096B (en) 2021-11-26

Similar Documents

Publication Publication Date Title
CN104950331B (en) A kind of porosity of sand mud reservoir and the earthquake prediction method of shale content
CN108399270B (en) Method for determining anisotropic shale proportion in shale stratum
EA027440B1 (en) Method of predicting the pressure sensitivity of seismic velocity within reservoir rocks
Duan et al. Lithology identification and reservoir characteristics of the mixed siliciclastic-carbonate rocks of the lower third member of the Shahejie formation in the south of the Laizhouwan Sag, Bohai Bay Basin, China
Jin et al. Rock-physics-model-based pore type characterization and its implication for porosity and permeability qualification in a deeply-buried carbonate reservoir, Changxing formation, Lower Permian, Sichuan Bain, China
US7664602B2 (en) Method and apparatus for quantitative determination of sediment physical properties of seabed and layered subbottom using single-channel seismic-reflection profiling data
CN112987096B (en) Method for calculating sound velocity of high-argillaceous sandstone
CN109581490A (en) A kind of Wave velocity estimation method and system of compact reservoir
Ham et al. Grain-scale tensile and shear strengths of glass beads cemented by MICP
Nabawy et al. Seismic reservoir characterization of the syn-rift lower Miocene Rudeis Formation in the July oilfield, Gulf of Suez basin, Egypt: implication for reservoir quality assessment
Belozerov et al. Concept of technology for determining the permeability and porosity properties of terrigenous reservoirs on a digital rock sample model
Zhou et al. The quantitative characterization of stylolites in the limestone reservoirs of the Lower Triassic Feixianguan Formation, northeastern Sichuan Basin: Insights to the influence of pressure solution on the quality of carbonate reservoirs
Guo et al. An NMR-assisted laboratory investigation of coal fines migration in fracture proppants during single-phase water production
Wu et al. Pore structure and movable fluid characteristics of typical sedimentary lithofacies in a tight conglomerate reservoir, Mahu depression, northwest China
Yu et al. Stress sensitivity of elastic moduli in high-porosity cemented sandstone—Heuristic models and experimental data
CN112228050B (en) Method for quantitatively evaluating macroscopic heterogeneity of tight oil reservoir and application thereof
Davoodi et al. A fresh look at the Lucia classification using mud-and grain-dominated reservoirs of the Persian Gulf
CN109884696A (en) Sandstone Computational approach of S-wave velocity based on calcite cementation function influence
Vernik et al. On some controversial issues in rock physics
Zhang et al. Construction and validation of the upscaling model of organic-rich shale by considering water-sensitivity effects
CN112946739B (en) Deep carbonate reservoir seismic rock physical template construction method and reservoir parameter prediction method in fracture-erosion hole double-hole system
Benjumea et al. Petrophysical characterization of carbonates (SE of Spain) through full wave sonic data
CN114414027B (en) Method for calculating sound velocity of supporting conglomerate and non-supporting conglomerate
Li et al. Study on Anisotropic Petrophysical Modeling of Shale: A Case Study of Shale Oil in Qingshankou Formation in Sanzhao Sag, Songliao Basin, China
Ambarsari et al. Reservoir quality determination by using critical porosity and volume of clay: a case study in the Talang Akar Formation in the NW Java Basin, Indonesia

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant