29 results about "Elastic wave velocity" patented technology

A method of determining rock matrix and fracture permeability to perform a reservoir fluid flow simulation. The method includes: obtaining a petrophysical characterization of a formation sample, measuring a first and second elastic wave velocity of the formation sample at a first and second frequency, calculating an average rock matrix density and an average porosity of the formation sample, determining a set of calculated elastic wave velocities for the first and second frequencies over a range of candidate permeabilities, and determining a rock matrix permeability and a fracture permeability based on the set of calculated elastic wave velocities at the first and second frequency respectively. The method further including performing a reservoir fluid flow simulation based, at least in part, on at least one of the rock matrix permeability and the fracture permeability.

The invention relates to a diversion tunnel complex environment concrete internal quality detection method based on elastic waves, and belongs to the technical field of concrete engineering nondestructive testing. The concrete lining thickness, the defect depth, the void position and the concrete elasticity modulus are calculated according to the wave velocity of the elastic wave, and the concrete strength can be converted through the concrete elasticity modulus. The method can detect and evaluate the concrete quality of the diversion tunnel lining project in a complex environment, provides targeted measures in the aspects of thicker tunnel lining, humid environment, noise interference and the like, and has the characteristics of convenience, high accuracy, rapidness and high efficiency, so that the concrete quality of the diversion tunnel lining project in the complex environment is scientifically and accurately detected and evaluated, and necessary prevention measures are taken for the tunnel project with hidden dangers. And disasters are avoided, so that safe production and operation of tunnel engineering are guaranteed.

The invention belongs to the technical field of landslide prediction, and relates to a shallow earthquake-stress joint inversion tunnel portal landslide prediction method, which comprises the following steps of: obtaining slope body first-arrival wave travel time data through a high-frequency seismic source array, constructing a geologic model containing a rock mass integrity coefficient matrix, identifying rock mass structure deterioration through elastic wave velocity distribution, and predicting the landslide at a tunnel portal through shallow earthquake-stress joint inversion. The method comprises the following steps: locking a potential unstable structure weak area in advance, then embedding a rock mass integrity coefficient matrix as priori knowledge into a stress field for solving, and by establishing a coupling relationship between seismic wave propagation and stress, inverting slope stress field distribution, capturing an early stress redistribution signal caused by structural degradation of a deep rock mass, and finally obtaining an initial stress redistribution signal of the deep rock mass. Early perception of mechanical instability precursor is realized; and then calculating strain energy density distribution of the potential slip surface based on the stress and strain data, and accurately judging the potential slip surface before the slip surface is not completely cut through by judging whether the strain energy density exceeds a threshold value, whether a continuous sheet is formed and whether the expansion rate exceeds the standard, thereby avoiding lagging caused by dependence on a displacement signal.

The invention discloses a TTI (Transmission Time Interval) medium rotating staggered grid finite difference three-dimensional forward modeling method, which mainly comprises the following steps of: establishing a TTI medium three-dimensional elastic wave velocity stress equation, determining a rotating staggered grid finite difference format for solving the equation, calculating the distribution of a three-dimensional velocity stress field, loading a carrier field to simulate a seismic source, and setting an absorption boundary. And finite difference numerical dispersion is eliminated. According to the method, TTI anisotropic medium three-dimensional seismic wave field simulation is realized, the problems of boundary reflection, numerical dispersion, poor calculation stability and the like existing in finite difference numerical simulation are solved, and the propagation law and kinematics characteristics of seismic waves in an anisotropic medium can be comprehensively known.

The invention discloses an acoustic monitoring method for a beam-direction stress state near a foundation surface of a concrete dam. The acoustic monitoring method comprises the following steps: S1, measuring a time difference between a sound wave transmitting signal and a sound wave receiving signal; s2, calculating the sound wave velocity of sound waves propagating along the upstream and downstream directions of the foundation surface of the dam body; s3, elastic wave velocity and viscous wave velocity are separated; s4, monitoring the sound wave velocity for a long time, and synchronously monitoring the reservoir water level to obtain a time-varying curve of the residual error; and S5, when the slope of the residual error changing curve changes suddenly along with the time, an alarm signal is sent out. According to the method, the elastic wave velocity is separated, and the influence of the strain rate effect is filtered out, so that the accurate measurement and monitoring of the beam-direction stress near the foundation surface of the concrete dam are realized.

The invention relates to the field of rock mechanics testing, in particular to a rock deformation and elastic wave velocity monitoring device and method under dynamic high temperature and high pressure, and the device comprises a clamp, a deformation monitoring assembly and an elastic wave velocity monitoring assembly. The deformation monitoring assembly is used for measuring thermal expansion deformation of the rock sample in the length, width and thickness directions at dynamic high temperature; the elastic wave velocity monitoring assembly is arranged on the rock sample and used for monitoring the elastic wave velocity of the rock sample. According to the monitoring device and method, at any moment, the thermal expansion of the clamp and the thermal expansion of the rod-shaped monitoring assembly offset each other, the thermal expansion deformation of the button and the thermal expansion deformation of the annular frame also offset each other, and the reading change of the displacement and deformation sensor directly and truly reflects the thermal expansion deformation of the rock sample; measurement errors caused by self deformation in the dynamic temperature rise process in the prior art are eliminated, and the accuracy of deformation data capture in the whole process from normal temperature to high temperature is ensured.

The present invention discloses a self-excited detection method for the reinforcement quality of a vibro-compacted gravel pile composite foundation. The pre-treated foundation is subjected to elastic wave velocity testing and test pile test velocity zoning, and at the same time, excess pore pressure monitoring is performed to obtain excess pore pressure; a detector and a collector are selected to establish a gravel pile construction excitation signal acquisition system, and a spatial array of detectors is arranged; elastic wave signals are collected in real time, filtered and processed, and signal identification is performed, the arrival time of the elastic wave signal is extracted, a travel time equation group is constructed, and the elastic wave velocity of the velocity partition is iteratively inverted; a corresponding relationship between the effective confining pressure, relative density, and elastic wave velocity of the soil is established, the relative density of the soil is input, and the reinforcement quality is judged by comparison. The present invention iteratively processes the travel time equation by velocity inversion and analysis to invert the elastic wave velocity and then interprets the velocity to determine the relative density of the soil, thereby realizing real-time detection of the reinforcement quality of the gravel pile composite foundation, and has the advantages of a large detection range, a wide detection time domain, and intelligent detection.

The application discloses a kind of embankment quality elastic wave tomography detection methods, comprising: synchronous acquisition compactor vibration source at source end reference signal and the ground surface response signal of measured region surface, establish the unified time base of both;Correlation function of source end reference signal and ground surface response signal is calculated, and the coarse estimation time is determined based on peak value.Extract feature frequency component from signal and obtain analytic signal phase information, calculate the phase difference of ground surface response signal at coarse estimation time relative to source end reference signal, and use the phase difference to determine fine time correction amount.Use fine time correction amount to correct coarse estimation time, obtain final elastic wave time and generate elastic wave velocity distribution image by inversion.The application breaks through the limitation of sampling rate on measurement accuracy by collecting real source end waveform and combining phase refinement technology, and improves the resolution and reliability of compaction quality detection.

The invention relates to a device and method for detecting the loose state of the internal structure of a high-speed rail supporting layer, the device comprises an elastic wave exciter arranged along one side line of the supporting layer and an elastic wave detector arranged along the other side line of the supporting layer, the elastic wave detector is in communication connection with a data acquisition instrument, the data acquisition instrument is installed on a mobile platform, and the mobile platform is connected with the elastic wave exciter. A computer connected with the data acquisition instrument is arranged on the mobile platform, and the elastic wave exciter is used for vertically knocking the side surface of a high-speed rail supporting layer so as to excite elastic waves; the elastic wave detector is used for detecting elastic wave speed response signals and transmitting the elastic wave speed response signals to the computer through the data acquisition instrument; and the computer is used for processing and analyzing the elastic wave speed response signal and outputting an internal structure state detection result of the supporting layer. Compared with the prior art, the method has the advantages that the propagation characteristics of elastic waves are utilized, the opposite penetrating waveform section is formed by changing the positions of the wave detector and the seismic source, and the loose state of the internal structure of the supporting layer can be detected and analyzed losslessly, rapidly and accurately.

To provide a speed estimation device, an evaluation system, and a speed estimation method that can improve the accuracy of estimating elastic wave velocity.SOLUTION: A speed estimation device of an embodiment has a spectrum information extraction unit, a peak frequency extraction unit, a range determination unit, a specification unit, and an elastic wave velocity estimation unit. The spectrum information extraction unit extracts, for each of different measurement points, spectrum information of one or more elastic waves generated from impact on the measurement points in a target member. The peak frequency extraction unit extracts one or more peak frequencies in every spectrum information extracted for each of the measurement points. The range determination unit determines a frequency range of the elastic wave on the basis of the one or more peak frequencies obtained for each of the measurement points. The specification unit specifies the peak frequency at each measurement point on the basis of the spectrum information and the frequency range. The elastic wave velocity estimation unit estimates the elastic wave velocity at each measurement point on the basis of the peak frequency and thickness information of the target member.SELECTED DRAWING: Figure 8

According to one embodiment, a speed estimation apparatus includes a spectrum information extractor, a peak frequency extractor, a range determiner, a specifier, and an elastic wave speed estimation unit. The spectrum information extractor extracts spectrum information of one or more elastic waves generated due to an impact at different measurement points of a target member for each measurement point. The peak frequency extractor extracts one or more peak frequencies in each of the pieces of spectrum information extracted for each of the measurement points. The range determiner determines a frequency range of the elastic waves based on the one or more peak frequencies. The specifier specifies a peak frequency at each measurement point based on the spectrum information and the frequency range. The elastic wave speed estimator configured to estimate an elastic wave speed at each measurement point based on the peak frequency and thickness information of the target member.

This application belongs to the field of landslide prediction technology, and relates to a method for predicting landslides at tunnel entrances using shallow seismic-stress joint inversion. This invention acquires first-arrival wave travel time data of the slope using a high-frequency seismic source array, constructs a geological model including a rock mass integrity coefficient matrix, identifies rock mass structural deterioration through elastic wave velocity distribution, and identifies potential unstable structural weak areas in advance. Then, the rock mass integrity coefficient matrix is ​​embedded as prior knowledge into the stress field solution. By establishing the coupling relationship between seismic wave propagation and stress, the stress field distribution of the slope is inverted, capturing early stress redistribution signals caused by structural deterioration in deep rock masses, achieving early detection of precursors to mechanical instability. Furthermore, based on stress and strain data, the strain energy density distribution of the potential slip surface is calculated. By judging whether the strain energy density exceeds a threshold, whether it is continuous, and whether the expansion rate exceeds the standard, the potential slip surface can be accurately identified before it is fully connected, avoiding the lag caused by relying on displacement signals.

The invention provides an ancient building exterior wall tile weathering degree detection and evaluation method, which comprises the following steps: acquiring a missing volume value, a surface hardness value and an elastic wave velocity value of an ancient building exterior wall tile needing to be evaluated, and calculating a volume missing rate, a surface hardness decrease rate and an elastic wave velocity decrease rate of the ancient building exterior wall tile; according to the relative importance degrees of the volume missing rate, the surface hardness decrease rate and the elastic wave velocity decrease rate in the ancient building weathering degree evaluation, determining the volume missing rate weight, the surface hardness decrease rate weight and the elastic wave velocity decrease rate weight by using an analytic hierarchy process, thereby calculating the weathering index of the ancient building exterior wall tile. Determining the weathering grade of the ancient building exterior wall brick according to the weathering index, and evaluating the weathering degree of the ancient building exterior wall brick by using the weathering grade. The method comprehensively considers three factors of the volume loss rate, the surface hardness decrease rate and the elastic wave velocity decrease rate of the external wall tile to scientifically and comprehensively evaluate the weathering degree of the external wall tile of the ancient building.

The present invention discloses a method and device for predicting elastic wave velocity in a strongly transversely isotropic medium, and relates to the technical field of geophysical oil and gas exploration. An elastic wave velocity model is constructed using the rock attribute information of the target section of the reservoir, and the elastic wave velocity is predicted based on the elastic wave velocity model. The present invention provides an elastic wave velocity model in which the elastic wave velocity uses the fracture weakness as an independent variable, which is convenient for determining the elastic wave velocity based on the fracture parameters and can be applied to strong anisotropy conditions. Moreover, in the case of strong anisotropy, the seismic elastic wave velocity prediction results provided by the present invention are more consistent with the precise values ​​and have a wider range of applications. In the process of establishing the elastic wave velocity model, the rock attribute information of the target section of the reservoir is fully utilized to accurately construct the corresponding rock physics model, which can accurately determine the corresponding reservoir rock and the fluid properties contained therein based on the seismic response characteristics.

The invention discloses an engineering geophysical prospecting data-based embankment safe and rapid grading method, which relates to the technical field of embankment safe and rapid grading methods, and comprises the steps of testing embankment soil state parameters, establishing a relation between geophysical prospecting data and soil state, carrying out embankment safe and rapid grading and the like. Through soil sample testing, the relation between the compactness, the water content and the like of a soil body and indexes such as dielectric constants, elastic wave velocities and electrical resistivity is established; through a forward modeling test or a physical test, relevance between multi-source geophysical prospecting parameters such as radar waves, elastic waves and resistivity and core factors such as embankment body structural integrity, embankment foundation stability and seepage conditions reflecting embankment structure safety is established, and a geophysical prospecting data analysis principle is formed; and adopting an analytic hierarchy process, and introducing a dike aging factor to construct a discrimination matrix to calculate a dike comprehensive score. According to the method, safe, rapid and dynamic grading of the dike can be realized, the dike risk level is pre-warned, and guidance is provided for rescue power investment.

The application discloses a karst collapse hidden danger area multi-method cooperative detection method and system, relates to the technical field of geographic detection, and comprises the following steps: arranging an observation network in a target hidden danger area, respectively implementing an electric excitation sequence and a mechanical vibration excitation sequence, collecting potential difference data and vibration waveform data, inverting the potential difference data and the vibration waveform data at each collection time, analyzing the phase relationship and the change form of the resistivity sequence and the elastic wave velocity sequence with respect to the corresponding excitation time sequence, and dividing the area into regions with different response characteristics according to the characteristic differences of the multi-physical field dynamic response modes in space. The application can extract the multi-physical field dynamic response mode representing the dynamics process such as groundwater seepage and particle migration, realize direct description of the dynamics process activity and evolution time sequence in the hidden danger area, and improve the ability of identifying the hidden danger development stage from the physical response.

The invention discloses a real-time monitoring and risk early warning method for corrosion damage of a deep soft rock roadway, and relates to the technical field of roadway monitoring. The method comprises the following steps: starting a data acquisition system, calculating a rock mass elastic wave velocity attenuation rate, a support rotation angle and a corrosion rate, and standardizing overpressure and impulse data to generate a real-time curve; real-time prediction of the risk level is realized based on the LSTM neural network; and selecting parameters sensitive to the corrosion holes, establishing a development index of the corrosion holes, and generating a corrosion damage report of each monitoring unit. According to the method, a relation model between monitoring data and the corrosion damage degree is established by fusing triple triggering conditions of LSTM neural network prediction, a multi-load coupling curve and a key parameter threshold value; the mechanical response and damage development process of the roadway under the corrosion effect is simulated, the future deformation and damage trend of the roadway is predicted, and the accuracy of an analysis result is improved.

Systems and methods for determining a connected porosity and a non-connected porosity in a fluid-saturated hydrocarbon reservoir are disclosed. The methods include obtaining at least one rock sample from the fluid-saturated hydrocarbon reservoir, determining, using a petrophysical sample analyzer, at least one petrophysical parameter of the rock sample, and measuring at least one of an elastic wave velocity and an elastic wave attenuation for each of a plurality of wave frequencies. The methods further include determining, using a computer processor, the connected porosity and the non-connected porosity of the rock sample using a dual-porosity single-permeability model based, at least in part, on at least one petrophysical parameter and at least one of the elastic wave velocity and the elastic wave attenuation.

The invention discloses an ultrasonic oblique incidence solid elastic wave speed measuring device and method.The device comprises a transducer, a three-axis positioning control system (4), an angle adjusting mechanism (6) and a host (1), the transducer is arranged in liquid, emits ultrasonic waves to a to-be-measured solid (10) at an adjustable incidence angle, collects refraction signals and uploads the refraction signals to the host (1) through a data collecting system (5), and the host (1) sends the refraction signals to the three-axis positioning control system (4); the host (1) is connected with the three-axis positioning control system (4), the three-axis positioning control system (4) is connected with the angle adjusting mechanism (6), and the angle of the transducer is adjusted through the angle adjusting mechanism (6). According to the invention, high-efficiency and non-destructive measurement of solid longitudinal wave and transverse wave velocity can be realized.

The invention discloses a method and a device for establishing a physical model of fracture-developed rock. The method comprises the following steps: determining an elastic parameter, a density and a flexibility parameter of a porous rock matrix according to mineral components, volume contents and pore parameters in combination with the density and the elastic parameter of the mineral components; according to the truncation Gaussian function corresponding to the azimuth angle of each group of cracks and a Bond transformation matrix used for converting the local coordinate system of the cracks into an observation coordinate system, in combination with the flexibility parameter of the porous rock matrix, determining the flexibility parameter and the elastic parameter of the dry rock skeleton; based on the density of the porous rock matrix and the elastic parameters and fluid attributes of the dry rock skeleton, determining the elastic parameters and density of the fluid saturated rock; and determining the elastic wave velocity based on the elastic parameter and density of the fluid saturated rock. According to the method, multi-group fracture distribution description and rock physical modeling can be realized.

To easily measure the tension force of an anchor.SOLUTION: An evaluation method includes the steps for: acquiring detection data, from a sensor detecting elastic waves, which indicates the detection result of an elastic wave at the time of being input to an anchor and the detection result of an elastic wave reflected within the anchor and input to the sensor; calculating the velocity of the elastic wave propagating through the anchor on the basis of the detection data; and calculating a tension force pe occurring in the anchor on the basis of the following equation (1), where Vp0 is a reference velocity of the elastic wave propagating through the anchor measured in advance, Vpa is a calculated elastic wave velocity, and θ is a correlation coefficient.SELECTED DRAWING: Figure 6

The invention belongs to the field of foundation treatment, and discloses a soil body compaction degree monitoring method and a related device, and the method comprises the steps: collecting a vibration wave signal of a vibration wheel of a road roller operating on the surface of a to-be-detected soil body; carrying out spatial domain filtering and frequency domain filtering processing to obtain a preprocessed vibration wave signal; based on the pre-processed vibration wave signal, combining a pre-constructed elastic wave velocity-density dynamic inversion model to obtain the total density of the to-be-measured soil body; the elastic wave velocity-density dynamic inversion model is deduced according to an expanded Biot fluid saturated porous medium wave equation; the expanded Biot fluid saturated porous medium wave equation is a Biot fluid saturated porous medium wave equation into which unsaturated soil phase constitutive and equivalent soil particle skeleton modulus are introduced; according to the total density of the to-be-measured soil body, the compaction degree of the to-be-measured soil body is obtained through calculation; according to the invention, high-precision monitoring of the compaction degree of the soil body is realized, continuous and real-time monitoring of the construction process of the adaptive road roller can be realized, and the monitoring efficiency and precision are both considered.

The invention discloses a multi-method cooperative detection method and system for a karst collapse hidden danger area, and relates to the technical field of geographic detection, and the method comprises the steps: laying an observation network in a target hidden danger area, respectively implementing an electric excitation sequence and a mechanical vibration excitation sequence, and collecting potential difference data and vibration waveform data; performing inversion on the potential difference data and the vibration waveform data at each acquisition moment, analyzing the phase relation and the change form of the resistivity sequence and the elastic wave speed sequence relative to the corresponding excitation time sequence, and dividing areas with different response characteristics according to the characteristic difference of the multi-physical field dynamic response mode in the space; according to the method, the multi-physical-field dynamic response mode representing the dynamic processes such as underground water seepage and particle migration can be extracted, direct depiction of the activity degree and the evolution time sequence of the dynamic processes in the hidden danger area is achieved, and the capacity of identifying the development stage of the hidden danger from physical response is improved.

The application relates to a CFST arch bridge tube concrete defect intelligent identification method based on an MCMC model, which comprises the following steps: detecting a CFST arch bridge rib (1) section by using an elastic wave, acquiring detection data, and the detection data comprising an elastic wave velocity; inputting the detected elastic wave velocity into a probability model for quantitative identification; the probability model comprising: judging the compactness of the detection data according to the probability model; when the detection data is located in a 95% confidence interval range of the probability model, the compactness of the section detection data is good; and when the detection data is located outside the 95% confidence interval of the probability model, the compactness of the section detection data is poor; the method has simple steps, is convenient to operate, and has good effects.

The invention provides a rock mass elastic wave velocity large-scale rapid prediction method based on a seismic wave method. The rock mass elastic wave velocity large-scale rapid prediction method comprises the steps that the real elastic wave velocity of the earth surface is acquired; acquiring a real elastic wave velocity in the tunnel; preprocessing the real elastic wave velocity data of the earth surface and the real elastic wave velocity data in the tunnel; constructing a rock mass elastic wave velocity prediction model based on the real elastic wave velocity of the earth surface and the real elastic wave velocity in the tunnel; and performing elastic wave velocity large-scale rapid prediction on the rock mass based on the rock mass elastic wave velocity prediction model. According to the rock mass elastic wave velocity large-scale rapid prediction method based on the seismic wave method provided by the invention, the problems of low efficiency and insufficient precision of obtaining the rock mass elastic wave velocity in the prior art can be solved.

The invention discloses an elastic wave tomography detection method for embankment quality, which comprises the following steps of: synchronously acquiring a source end reference signal at a vibration source of a road roller and an earth surface response signal on the surface of a detected area, and establishing a unified time reference of the source end reference signal and the earth surface response signal; and calculating a cross-correlation function of the source end reference signal and the surface response signal, and determining coarse estimation arrival time based on a peak value. And extracting a characteristic frequency component from the signal, obtaining and analyzing signal phase information, calculating the phase difference of the earth surface response signal relative to the source end reference signal at the roughly estimated time, and determining the fine time correction by using the phase difference. And correcting the roughly estimated arrival time by using the fine time correction to obtain final elastic wave arrival time, and performing inversion to generate an elastic wave velocity distribution image. According to the invention, by collecting the real source end waveform and combining the phase refinement technology, the limitation of the sampling rate on the measurement precision is broken through, and the resolution and reliability of compaction quality detection are improved.

The invention relates to the technical field of building anchorage devices, in particular to a nonlinear fluctuation analysis method for working stress of a prestressed anchorage device in an isotropic structure, which comprises the following steps of: 1, establishing an anchorage device model: establishing the prestressed anchorage device model of an isotropic material, and assuming that the anchorage device is subjected to unidirectional tensile stress in an isotropic medium; step 2, acoustic elasticity theory analysis: analyzing the relationship between stress and elastic wave velocity in the prestressed anchorage device based on a nonlinear fluctuation theory; supposing that the medium is a continuum, the small disturbance of the sound wave is superposed on the static large deformation of the medium; the one-dimensional acoustic-elastic relation in a one-way stress medium is deduced according to the nonlinear wave equation, the working stress of the pre-stressed anchor cable can be accurately monitored through a nonlinear wave analysis method, the too high or too low stress condition can be found and processed in time, structural faults are prevented, and the reliability of the pre-stressed anchor cable is improved. And the safety and reliability of the project are improved.

The application discloses a method for rapidly simulating three-dimensional scattered waves in a sound wave far detection transversely isotropic stratum, and belongs to the field of seismic wave forward simulation. The application is used to solve the technical problem that it is difficult to efficiently and universally simulate three-dimensional scattered waves in single well imaging of the transversely isotropic stratum in the prior art. The application comprises the following steps: obtaining parameters required for simulation; deriving an analytical solution of a sound source far field wave under a transversely isotropic medium according to a steepest descent integral method, and obtaining an incident wave field, wherein saddle points of quasi-SV waves and quasi-P waves are solved by using a Newton-Raphson method; establishing an elastic wave velocity-stress staggered grid, introducing the incident wave field into a calculation space boundary by using a total field / scattered field technology, and performing numerical simulation on a local region containing a scatterer to obtain a near field scattered wave field; calculating a virtual source wave field based on the analytical solution; and converting the near field scattered wave field and the virtual source wave field into a scattered wave field at a position of a well receiver based on a fluid-solid coupling reciprocity theorem. The application is suitable for a transversely isotropic medium with arbitrary parameters.

The invention provides a loose sandstone displacement process digital rock core construction and elastic wave velocity simulation method, and belongs to the technical field of oil and gas geophysical exploration, and the method comprises the following steps: constructing an initial digital rock core; continuously calculating and operating water / oil saturation values of all positions of the array in a traversal mode, so as to realize establishment of digital rock cores with different water / oil saturations; on the basis of digital rock cores with different water / oil saturations, continuously calculating and operating substance types at all positions of an array by adopting a traversal mode, and realizing construction of digital rock core framework weakening in the displacement process; initial conditions and boundary conditions are established, constitutive and wave equations of the digital rock core containing the multiphase fluid are adopted, and elastic wave velocity calculation of the digital rock core is achieved. According to the method, the digital rock core in the loose sandstone multiphase fluid displacement process is constructed, the elastic wave velocity of the loose sandstone multiphase fluid displacement process is simulated through the rotary staggered grid finite difference technology, and key technical support can be provided for quantitative parameter interpretation of remaining oil distribution.