62 results about "Bulk modulus" patented technology

The invention relates to the technical field of oil and gas geophysical exploration, and particularly discloses a tight sandstone gas reservoir fluid mobility gradient attribute prediction method and system.The method comprises the steps that a reflection coefficient expression containing a fluid mobility item is constructed, and a formula for calculating fluid mobility through frequency dispersion attributes is deduced; constructing a fluid mobility gradient attribute related to the frequency and the incident angle of the seismic wave; the method comprises the following steps: establishing a double-pore rock physical model, generating a synthetic seismic record, applying frequency variation AVO inversion to the synthetic seismic record, and extracting a bulk modulus dispersion attribute and a shear modulus dispersion attribute; and calculating the fluid fluidity changing along with the frequency and the incident angle, and calculating the gradient attribute of the fluid fluidity based on a gradient attribute calculation formula and the fluid fluidity. According to the method, seismic wave frequency dispersion and attenuation response caused by permeability and gas saturation are fully utilized, the seismic gradient attribute of the fluid fluidity and the calculation method are constructed, and high-precision prediction of the fluid fluidity by applying seismic signal elasticity and frequency dispersion rock physical information is realized.

This invention discloses a quantitative evaluation method for the homogeneity of lightweight aggregate concrete based on rapid ultrasonic testing, belonging to the field of engineering materials testing technology. The method includes: S1, obtaining the shear wave velocity and longitudinal wave velocity at multiple test points of a series of lightweight aggregate concrete reference specimens; S2, calculating the shear modulus and bulk modulus at each test point of each reference specimen; S3, calculating the elastic modulus at each test point of each reference specimen; S4, constructing a relationship between the lightweight aggregate volume content, shear wave velocity, and longitudinal wave velocity based on the elastic modulus obtained in step S3; S5, obtaining the shear wave velocity and longitudinal wave velocity at each test point of the specimen to be tested, and calculating the lightweight aggregate volume content at each test point of the specimen to be tested based on the relationship; S6, evaluating the homogeneity of the specimen to be tested based on the lightweight aggregate volume content calculated in step S5. This invention solves the technical problem of insufficient accuracy in ultrasonic testing results for the homogeneity of lightweight aggregate concrete.

The present application relates to oil and gas geophysical exploration technology field, specifically disclose a kind of tight sandstone gas reservoir fluid mobility gradient attribute prediction method and system, the method includes constructing fluid mobility term's reflection coefficient expression, and deduce the formula for calculating fluid mobility using dispersion attribute, construct and frequency and seismic wave incidence angle related fluid mobility gradient attribute;Establish double-pore petrophysical model, generate synthetic seismogram, apply frequency-variable AVO inversion to the synthetic seismogram, extract bulk modulus dispersion attribute and shear modulus dispersion attribute;Fluid mobility changing with frequency and incidence angle is calculated, based on gradient attribute calculation formula and fluid mobility, fluid mobility gradient attribute is calculated.The present application makes full use of the response of the seismic wave dispersion and attenuation caused by permeability and gas saturation, constructs the seismic gradient attribute and calculation method of fluid mobility, realizes the high-precision prediction of fluid mobility by applying the elastic and dispersion petrophysical information of seismic signal.

The application provides a method for determining the longitudinal wave velocity after fluid replacement in an anisotropic medium and related equipment. The method comprises the following steps: determining an average bulk modulus based on a matrix bulk modulus; determining a first bulk modulus when the medium contains fluid based on first property data; determining a second bulk modulus when the medium does not contain fluid based on the first property data, the matrix bulk modulus, the average bulk modulus and the first bulk modulus; determining a second vertical stiffness coefficient when the medium does not contain fluid based on the first property data, the average bulk modulus, the first bulk modulus and the second bulk modulus; obtaining second property data of each predetermined replacement fluid; and determining the longitudinal wave velocity after the fluid in the medium is replaced by the predetermined replacement fluid based on the first property data, the average bulk modulus, the second bulk modulus, the second vertical stiffness coefficient and the second property data. The technical problem that the longitudinal wave velocity determined after fluid replacement in an anisotropic medium is not accurate in the prior art is solved.

The present application belongs to the field of material testing, and particularly relates to a viscoelastic material bulk modulus measuring device and method based on an electric measuring method, wherein the measuring device comprises a pressure tank and a viscoelastic material testing support; the viscoelastic material testing support comprises a fixing frame and a plurality of sets of deformation testing clamps; the deformation testing clamp comprises a fixed elastic sheet, a testing elastic sheet and a clamp block, the fixed elastic sheet is provided with a hollow slot in the middle, the testing elastic sheet is arranged through the hollow slot, and the fixed elastic sheet and the testing elastic sheet are perpendicular to each other, a strain gauge is pasted on the testing elastic sheet, adjacent one end of the fixed elastic sheet and the testing elastic sheet is fixedly arranged on the fixing frame, and adjacent the other end is fixedly arranged on the clamp block, and the clamp block is provided with an abutting rod perpendicular to the arrangement direction of the testing elastic sheet. The present application realizes direct and high signal-to-noise ratio measurement of the bulk strain of the viscoelastic material under a hydrostatic pressure environment by precisely isolating non-measurement direction interference through the double elastic sheet structure of the fixed elastic sheet and the testing elastic sheet of the deformation testing clamp.

The present application relates to the geological exploration technical field, more particularly, to a deep water super shallow layer natural gas reservoir rock physics shear wave velocity prediction method, comprising: calculating equivalent bulk modulus, shear modulus and density of a plurality of mineral mixed matrix according to shale content; then calculating skeleton bulk modulus and shear modulus; then calculating saturated rock bulk modulus and shear modulus through a model, and calculating saturated rock density through rock matrix density; distributing longitudinal and transverse velocities of suspended state and bearing state, and calculating equivalent longitudinal and transverse wave velocities; using equivalent longitudinal wave velocity to construct an inversion objective function with longitudinal wave velocity as a known condition; comparing actual longitudinal wave velocity with calculated longitudinal wave velocity; realizing nonlinear inversion of transverse wave velocity according to an algorithm, and obtaining transverse wave velocity of a to-be-inverted work area. The present application simultaneously considers rock physics models of particle suspended state and bearing state, and realizes transverse wave velocity prediction and deep water super shallow layer natural gas reservoir transverse wave velocity direct inversion.

The invention provides a seismic wave frequency dispersion and attenuation prediction method considering a cross-scale attenuation mechanism, and the method comprises the steps: calculating a longitudinal wave speed according to the physical parameters of a basic rock, and carrying out the reverse deduction of a first dry rock volume modulus through the longitudinal wave speed; directly calculating to obtain a second dry rock bulk modulus according to the physical parameters of the basic rock; establishing a dry rock volume modulus model of a cross-scale effect, and fusing the obtained first dry rock volume modulus and the second dry rock volume modulus in a weighted summation mode to obtain a fused dry rock volume modulus; and substituting the bulk modulus of the dry rock into a Biot equation, obtaining a relationship between frequency and wave number according to a plane wave analysis principle, and calculating longitudinal wave velocity dispersion and attenuation according to a definition formula. According to the seismic wave dispersion and attenuation prediction method considering the cross-scale attenuation mechanism, the problem of more complex multi-fluid and multi-pore type seismic wave dispersion and attenuation prediction is solved.

The invention provides a modeling method and a saturation prediction method of a rock physical model containing non-uniformly distributed hydrate sediments. Comprising the steps of obtaining model parameters, calculating the porosity of a mineral skeleton based on the model parameters, calculating the bulk modulus and the shear modulus of a dry hydrate skeleton, calculating the bulk modulus and the shear modulus of a dry mineral skeleton, and calculating the bulk modulus and the shear modulus of a dry composite skeleton composed of the dry mineral skeleton and the dry hydrate skeleton. And calculating the rock bulk modulus and shear modulus when the fluid is saturated, calculating the rock density value when the fluid is saturated, and calculating the stratum longitudinal wave velocity and the stratum transverse wave velocity. According to the method, on the basis of the forming mechanism and distribution characteristics of the non-uniformly distributed hydrates, the non-uniformity of the hydrates and mineral particles is considered in the modeling process, the hydrates and the mineral particles are regarded as two independent rock frameworks, and the isotropic rock physical model suitable for the non-uniformly distributed hydrates is constructed.

The invention relates to the technical field of oil and gas exploration, and discloses a hole structure parameter inversion method and device, electronic equipment and a medium. The method comprises the following steps: acquiring key rock physical parameters; processing logging data of the typical well, and calculating a hole structure parameter curve of the typical well; mathematical optimization is carried out, and coefficients of hole structure parameter constraint conditions are determined; pre-stack elastic parameter inversion is used to acquire a longitudinal wave impedance data body and a transverse wave impedance data body of the work area to be inverted; based on hole structure parameter constraint conditions, establishing a nonlinear mapping relation of longitudinal wave impedance and transverse wave impedance about the porosity, and solving a porosity inversion result by using a golden section algorithm; and calculating inversion results of the volume modulus hole structure parameters and the shear modulus hole structure parameters. According to the technical scheme, density is not used, inversion of hole structure parameters is achieved only through longitudinal wave impedance and transverse wave impedance, and input data are more stable; the use conditions of hole structure parameter inversion are broadened, and the stability and precision of hole structure parameter inversion are improved.

This invention belongs to the field of advanced metallic materials manufacturing, and relates to the design and preparation method of a high-plasticity Mg-Sn-Al magnesium alloy. First, a candidate system of Mg-xSn-yAl alloys is established. First-principles methods based on density functional theory are used to perform VCA modeling on each candidate composition and calculate elastic constants. Combined with the Voigt-Reuss-Hill model, parameters such as bulk modulus B, shear modulus G, Poisson's ratio, Pugh ratio, and Cauchy pressure are obtained. Using a Poisson's ratio greater than 0.36, a Pugh ratio greater than 3.4, and a Cauchy pressure greater than 18 GPa as screening criteria, the composition range of the high-plasticity alloy is determined to be Mg-(1-3.2)Sn-(0.7-0.8)Al (wt.%). Subsequently, alloy plates are prepared through melting and casting, hot extrusion, staged rolling, intermediate annealing, and pulsed current post-treatment, and tensile property testing and microstructure characterization are performed. The results show that the selected alloys have superior plasticity, with the Mg-3.2Sn-0.79Al alloy achieving an elongation of 15.6%. This invention enables rapid screening and preparation of high-ductility magnesium alloys, which can reduce trial-and-error costs and improve alloy design efficiency.

The application belongs to the field of rock physics and geophysical exploration, and relates to a simulation method for wettability affecting dispersion and attenuation of partially saturated sandstone, which comprises the following steps: S1. establishing a relationship between capillary pressure and water saturation for wettability correction; S2. calculating equivalent pore throat radius and effective permeability according to the relationship between capillary pressure and water saturation; S3. calculating characteristic frequency and complex bulk modulus under different wettability conditions based on a partially saturated patchy medium model; S4. simulating and comparing the characteristics of longitudinal wave velocity and attenuation under different wettability and water saturation conditions, analyzing the modulation law of frequency-dependent acoustic response caused by wettability change, and realizing quantitative characterization of the wettability effect. The step-by-step modeling method proposed in the application, combined with the wettability correction parameter represented by the contact angle, can systematically calculate the continuous coupling process from capillary pressure to equivalent pore throat structure and equivalent permeability to acoustic response, and realize quantitative analysis of the wettability change on the dispersion and attenuation law.

The invention relates to the technical field of petroleum and natural gas hydrate exploration, and discloses a method and device for identifying coexistence and occurrence types of hydrate and free gas and a medium. The method comprises the following steps: establishing a temperature-depth relational expression of stratums below the seabed; determining the depths of the bottom boundaries of the stable zones of the type I and type II hydrates; calculating a bulk modulus ksat, a shear modulus Gsat, a longitudinal wave velocity Vp, a transverse wave velocity Vs and a saturated water formation resistivity R0; comparing logging data; determining a hydrate-containing layer; determining a hydrate and free gas coexistence layer; and determining a pore filling type hydrate and free gas coexistence layer and a fracture filling type hydrate and free gas coexistence layer. According to the technical scheme, on the basis of intersection analysis of the porosity and the resistivity difference and conjoint analysis of various logging data, the hydrate and free gas coexisting stratum and the occurrence form are rapidly and effectively recognized, and the optimal target and the reservoir characteristic difference can be rapidly determined on a hydrate drilling site.

The embodiment of the application relates to the technical field of seismic exploration, and discloses a brittle prediction method, device, equipment and medium for a rock reservoir, the prediction method comprising the following steps: based on well logging data of a target rock reservoir, a target layer section of the target rock reservoir is constructed according to rock components and rock physical characteristic parameters, so as to construct a rock physical model for representing the relationship between the volume modulus, shear modulus, rock density and the content of main brittle minerals of the rock; a model of a brittle indicator factor is constructed according to the rock physical model and the change of the content of the main brittle minerals, the brittle indicator factor being used for representing the relative brittle degree of the rock mass; and the brittle indicator factor is converted into actual brittleness for representing the actual brittle degree of the rock mass. The brittle prediction method disclosed by the application fully utilizes the rock physical characteristic analysis of the target rock reservoir to analyze the influence of the brittle minerals on the brittle prediction, and effectively improves the rationality and accuracy of the rock brittle prediction.

This application relates to a method for calculating the internal pressure and output force of a hydraulic muscle using multi-factor coupling, comprising the following steps: obtaining the initial geometric parameters and rubber material parameters of the hydraulic muscle; constructing a multi-factor coupled output model based on the initial geometric parameters and rubber material parameters; and obtaining the input parameters of the hydraulic muscle, including the internal pressure data of the hydraulic muscle at the previous moment, the hydraulic oil input flow rate data at the current moment, the length data at the current moment, and the equivalent bulk modulus of elasticity data. This method, by constructing a multi-factor coupled output model, achieves high-precision calculation of the internal pressure and output force of the hydraulic muscle, significantly improving the physical realism of the simulation results. Furthermore, the initial geometric parameters and rubber material parameters required for model construction can be obtained by consulting product manuals or through simple measurements, solving the core pain points of existing data-driven methods, such as poor versatility and high calibration costs.

This invention discloses a method, apparatus, electronic device, and program product for predicting transverse and longitudinal wave velocities in rocks. The method includes: obtaining the porosity, density, and volume contents of quartz, clay, carbonate rocks, and pyrite in a marine-continental transitional shale rock; determining the equivalent matrix elastic modulus based on the porosity and the volume contents of quartz, clay, carbonate rocks, and pyrite; determining the volume modulus and shear modulus of the organic matter framework based on the volume modulus and shear modulus of the main phase material and the volume contents, volume modulus, and shear modulus of inclusions; obtaining the volume modulus and shear modulus of the shale matrix minerals based on the equivalent matrix elastic modulus, the volume modulus and shear modulus of the organic matter framework; and determining the volume modulus and shear modulus of the shale under saturated fluid conditions by mixing fluids, shale matrix minerals, and hard pores and fractures. This method can accurately predict transverse and longitudinal wave velocities.

The present application relates to a kind of double-scale full-band petrophysical modeling method, which is according to the inversion of target area rock P-wave and S-wave velocity, obtains the distribution of rock soft porosity and soft pore aspect ratio under different pressure;According to the distribution of rock soft porosity and soft pore aspect ratio, the frequency-dependent rock wet skeleton considering all aspect ratio soft pore jet flow effect is calculated using jet flow model;The frequency-dependent rock wet skeleton is brought into the frequency-dependent bulk modulus model of patch-saturated rock, and the petrophysical model considering micro and macro scale wave-induced flow effect is obtained.The present application comprehensively considers the influence of micro jet flow effect and mesoscale patch saturation, simulates the dispersion attenuation phenomenon caused by complex pore structure and complex fluid distribution of actual rock, and is more in line with the pore structure and fluid distribution in actual reservoir rock.

The invention provides a method, equipment and medium for estimating the acoustic velocity and density of multi-component mixed gas, and the method comprises the steps: introducing a temperature dependence function according to a Wilke gas mixing rule, and obtaining an absorption coefficient and a rejection coefficient of each component gas; defining the temperature dependence function to obtain an absorption coefficient and a rejection coefficient of the mixed gas; based on a temperature dependence function, improving a Peng-Robinson state equation to obtain the density and bulk modulus of the multi-component mixed gas, and then calculating the acoustic velocity of the multi-component mixed gas; and according to a Lohrenz-Bray-Clark model, the viscosity of the multi-component mixed gas is calculated. According to the method, a Peng-Robinson state equation is improved by introducing a temperature dependence function and a Wilke gas mixing rule so as to calculate the density and acoustic velocity of the multi-component mixed gas under different temperature-pressure conditions, and the viscosity of the multi-component mixed gas is obtained by combining and utilizing a Lohrenz-Bray-Clark model. And the accuracy of estimating the density, the acoustic velocity and the viscosity of the multi-component mixed gas is improved.

The application provides a fluid bulk modulus parameter extraction method, a fluid identification method, equipment and a medium, wherein the parameter extraction method comprises the following steps: dividing a seismic gather into at least three angle intervals, selecting a typical angle for each interval, and obtaining a corresponding partial stack for each interval; obtaining an elastic impedance curve and a curve of K f , f m and p for each typical angle from well logging data; performing elastic impedance inversion on the elastic impedance curve and the partial stack to obtain elastic impedance data of the typical angle; extracting wellside elastic impedance data from the corresponding elastic impedance data body by well logging; constructing a calculation relationship by using the curve of K f , f m and p and the elastic impedance curve of each typical angle, and bringing the curve of K f , f m and p and the wellside elastic impedance data into the calculation relationship to extract fluid bulk modulus. The technical scheme of the application can improve the resolution of the elastic impedance data body, reduce fluid identification artifacts, and improve the drilling rate of fluid inversion.

The present application belongs to the technical field of rock physics, and is especially a longitudinal wave velocity calculation method under a gas-water system. The method comprises the following steps: S1: obtaining longitudinal and transverse wave velocities of a sample under normal pressure and under a limit high pressure through experiments; S2: calculating a skeleton bulk modulus and a clay water relaxation time according to a sample XRD mineral analysis result; S3: converting the wave velocities into moduli; S4: calculating saturated volcanic lava longitudinal wave velocities considering frequency dispersion according to a Gurevich formula; S5: measuring a nuclear magnetic T2 spectrum of the sample under a saturated state; S6: calculating bulk moduli under different saturation degrees; and S7: converting to obtain longitudinal wave velocities under different saturation degrees. The present application explains the reason why the volcanic lava saturation degree-longitudinal wave velocity rule presents a "convex" shape from a mechanism, and effectively calculates volcanic lava longitudinal wave velocities under different saturation degrees.

The application discloses a three-parameter-based dry rock volume modulus calculation method, device, equipment and medium. The method comprises the following steps: determining a first physical characteristic parameter of a target dry rock; and determining a first dry rock volume modulus by using a relational model according to the first physical characteristic parameter of the target dry rock. The method can solve the problem that the adaptability of the existing dry rock volume modulus calculation method is poor, so that the seismic wave velocity calculation formed on the basis is not accurate enough, and the calculation precision of the dry rock volume modulus is improved.

An acoustic device for achieving a zero static bulk modulus, and a method. The acoustic device may comprise: an acoustic cavity, having a compressible fluid medium therein; a composite plate, arranged to vibratably close an opening on at least one side of the acoustic cavity; and an instability introducing member, coupled to the composite plate, and arranged to provide negative stiffness to the composite plate during vibration of the composite plate, so as to counteract positive stiffness of both the acoustic cavity and the composite plate, thereby achieving a critical state of an acoustic unit, the critical state indicating that a static bulk modulus of the acoustic unit approaches or equals zero.

The invention discloses a high-entropy rare earth aluminate based on light / heavy rare earth combination, the chemical formula of the high-entropy rare earth aluminate meets RE4Al2O9, RE is a rare earth element combination and is composed of a light rare earth element and three heavy rare earth elements, the number of the rare earth elements is four, and the substance amount ratio of the light rare earth element to the three heavy rare earth elements is 1: 1: 1: 1; the material is of a single-phase structure and has high-temperature stability. The preparation method comprises the following steps: firstly, carrying out solid-phase sintering, and then carrying out ball milling and sieving. When the material is used as a thermal barrier coating material, hot pressed sintering is carried out, and the Young modulus of 180-210 GPa, the shear modulus of 70-90 GPa, the bulk modulus of 130-150 GPa, the Vickers hardness of 9.1-9.2 GPa, the bending strength of 190-220 MPa and the fracture toughness of 2.9-3.2 MPa.m < 1 / 2 > are obtained. The thermal conductivity is 1.61-1.87 W.m <-1 >. K <-1 > in the temperature range of 25-1000 DEG C, and the thermal conductivity is 2.07-2.22 W.m <-1 >. K <-1 > at the high temperature of 1200 DEG C.

This invention provides a method for calculating the bulk modulus of elasticity of an ultra-high pressure common rail system. The method includes: constructing an ultra-high pressure common rail system test system; selecting measuring points X and Y; recording the pressure signal curves of measuring points X and Y during injection under fixed camshaft speed, injector pulse width, and high-pressure fuel line length; calculating the ultra-high pressure wave propagation velocity u = L / t under each operating condition; calculating the pressure wave propagation velocity at different fuel pressures and temperatures; calculating the fuel density corresponding to the pressure wave propagation velocity under the corresponding pressure and temperature conditions; obtaining the formula for calculating the bulk modulus of elasticity of the ultra-high pressure common rail system under different fuel pressure and temperature conditions; and verifying the calculation results. The method for calculating the bulk modulus of elasticity of an ultra-high pressure common rail system described in this invention can accurately calculate the bulk modulus of elasticity of the system under different fuel pressure and temperature conditions.

The invention belongs to the field of material testing, and particularly relates to a viscoelastic material bulk modulus measuring device and method based on an electrical measurement method, and the measuring device comprises a pressure tank and a viscoelastic material testing support; the viscous-elastic material testing bracket comprises a fixing frame and a plurality of groups of deformation testing clamps; the deformation test clamp comprises a fixed elastic piece, a test elastic piece and a clamp block, an empty groove is formed in the middle of the fixed elastic piece, the test elastic piece penetrates through the empty groove, the fixed elastic piece and the test elastic piece are perpendicular to each other, a strain gauge is pasted on the test elastic piece, the adjacent ends of the fixed elastic piece and the test elastic piece are fixedly arranged on the fixed frame, and the clamp block is fixed to the other adjacent ends. The clamp block is provided with an abutting rod which is perpendicular to the arrangement direction of the test elastic sheets. Interference in the non-measurement direction is accurately isolated through the double-elastic-sheet structure of the fixing elastic sheet and the testing elastic sheet of the deformation testing clamp, and direct and high-signal-to-noise-ratio measurement of the volume strain of the viscous-elastic material in the hydrostatic pressure environment is achieved.

The application discloses a kind of ultrasonic method-based light aggregate concrete compressive strength nondestructive testing methods, belong to engineering material testing field, this method includes S1, the ultrasonic transverse wave velocity and longitudinal wave velocity of each light aggregate concrete reference piece are detected;S2, the equivalent bulk modulus of each light aggregate concrete reference piece is calculated;S3, correction KT model is established;S4, the mapping relationship equation set of light aggregate concrete reference piece and correction KT model is established, and correction KT model is optimized adjustment;S5, the ultrasonic transverse wave velocity and longitudinal wave velocity of light aggregate concrete test piece, equivalent bulk modulus are obtained, and the volume fraction of light aggregate concrete test piece is calculated;S6, according to the volume fraction of light aggregate concrete test piece and the mapping relationship of compressive strength of mixing amount parameter, the compressive strength of light aggregate concrete test piece is calculated.The application solves the technical problem that the error of light aggregate concrete compressive strength nondestructive testing is higher.

The invention discloses a metamaterial with an evanescent shear modulus and a topological optimization design method of the metamaterial, and belongs to the technical field of aerospace structures and advanced material design. The method comprises the following steps: establishing a design domain of a periodic unit cell and discretizing; constructing a topological optimization mathematical model taking a relaxation function as an optimization target, and converting the problem of maximizing the ratio of the bulk modulus to the shear modulus, which is difficult to process directly, into a weighted optimization problem which can be solved steadily; predicting equivalent performance in combination with an energy-based homogenization method, and performing iterative solution by adopting an improved moving asymptote method and a novel floating projection function; and through error verification and parameterization reconstruction, outputting a geometric model with a clear boundary. According to the method, novel metamaterials such as two-dimensional quasi-single-mode / dual-mode and three-dimensional quasi-three-mode metamaterials with extremely high B / G ratio can be designed, seamless connection from optimization design to additive manufacturing is guaranteed, and an innovative material design solution is provided for high-performance vibration isolation and acoustic regulation and control components in the aerospace field.

The application discloses a kind of suitable for deep brackish aquifer seismic reflection wave response high fidelity simulation calculation method, the method includes S10: according to brackish aquifer development characteristics, establish brackish aquifer pore fluid rock physics model, calculate the viscosity, density, longitudinal wave velocity and bulk modulus of brackish aquifer pore fluid;S20: according to brackish aquifer development characteristics, establish brackish aquifer multi-scale rock physics model, obtain brackish aquifer equivalent stiffness matrix;S30: establish geology-seismic forward modeling;S40: utilize propagation matrix method to calculate brackish aquifer reflection and transmission coefficient matrix;S50: calculate brackish aquifer reflection wave seismic response.The application combines seismic rock physics modeling and seismic forward simulation technology, considers the viscoelasticity characteristics and thickness of brackish aquifer, is more in line with actual geological conditions, low in cost, high in accuracy, can quickly and accurately calculate deep brackish aquifer seismic reflection wave response, provides theoretical support for the evaluation and optimization of carbon storage space.

The invention provides a shale oil-containing inelastic reservoir seismic frequency-dependent reflectance calculation method, which comprises the following steps of: 1, inputting volume fractions and elastic parameters of different mineral components in a shale oil-containing inelastic reservoir rock physical model, and performing parameter conversion; 2, calculating the elastic modulus of the dry rock skeleton; 3, calculating the bulk modulus of the multiphase fluid mixture; step 4, performing equivalent filler complex modulus calculation; 5, calculating a non-elastic modulus and an oiliness identification factor of the oil-containing shale oil reservoir; and step 6, calculating speed dispersion and attenuation caused by wave-induced fluid flow. According to the method for calculating the seismic frequency-dependent reflectance of the shale oil-containing inelastic reservoir, the correspondence between the shale oil reservoir oiliness and the seismic frequency-dependent reflectance is disclosed, the rock physical characteristics and the response mechanism of the shale oil reservoir are clarified, and the method is of great significance to shale oil reservoir prediction.

This invention discloses a method and apparatus for shale rock physical modeling. The method includes: determining a maturity index curve based on kerogen maturity correlation curves; determining the density, bulk modulus, and shear modulus curves of kerogen based on discrete data of kerogen modulus and the maturity index curve; determining the porosity curves of organic and inorganic matter based on the adsorbed gas and free gas content curves and the maturity index and density curves of kerogen; determining the dry rock density, P-wave and S-wave velocities, and shear modulus curves based on the density, bulk modulus, and shear modulus curves of kerogen, clay minerals, and other minerals, as well as the total porosity curve, and the porosity curves of organic and inorganic matter; and determining a shale rock physical model based on the dry rock density, P-wave and S-wave velocities, and shear modulus curves, and the bulk modulus and density curves of the mixed fluid. This method enables the reasonable establishment of a shale rock physical model based on the changes in organic matter during the maturation process.

The present application relates to coal reservoir stress prediction technology field, specifically to a kind of coal reservoir stress prediction method.The method includes according to the depth distribution of drilling filtrate invasion formation and cleat direction, configure array sonic logging source distance;After source distance configuration, multiple logging is collected under each source distance P-wave velocity and S-wave velocity, and formation bulk density, cleat porosity is collected;Longitudinal wave velocity ratio is calculated by source distance, and far well boundary is selected, further identify far well P-wave velocity, far well S-wave velocity, near well P-wave velocity and near well S-wave velocity;Inversion dry skeleton bulk modulus and dry skeleton shear modulus;Near well wave velocity is replaced by fluid to obtain corrected P-wave velocity, corrected S-wave velocity and original gas-water two-phase state formation bulk density;Elastic parameter is calculated and substituted into geostress calculation model, and geostress prediction result is obtained.The present application effectively reduces the near well acoustic distortion caused by drilling fluid filtrate, and improves the accuracy of coal reservoir stress prediction.