A P-wave velocity model was built in the central southern of the Tanlu Fault based on double-difference tomography.The results suggest the presence of a low-velocity anomaly extending from the surface to a depth of 25...A P-wave velocity model was built in the central southern of the Tanlu Fault based on double-difference tomography.The results suggest the presence of a low-velocity anomaly extending from the surface to a depth of 25 km around the Tanlu and Feixi Faults,representing fault-related fluids caused by partial melting.The relocated earthquakes indicate a significant concentration of seismic activity above 20 km around the Tanlu and Feixi Faults,suggesting that prominent fault systems possibly serve as conduits for the upward migration of deep minerals.The proposed geodynamic model,supported by geological and geophysical data,suggests that the migration of deep mineralized materials extends along the Tanlu Fault.The obtained results serve as a crucial foundation for elucidating the intricate process of mineralization in the central southern segment of the Tanlu Fault,thereby enhancing comprehension regarding the interaction among ore body formation,fault fluids,localized melting,and seismic activity.展开更多
The integration of IoT and Deep Learning(DL)has significantly advanced real-time health monitoring and predictive maintenance in prognostic and health management(PHM).Electrocardiograms(ECGs)are widely used for cardio...The integration of IoT and Deep Learning(DL)has significantly advanced real-time health monitoring and predictive maintenance in prognostic and health management(PHM).Electrocardiograms(ECGs)are widely used for cardiovascular disease(CVD)diagnosis,but fluctuating signal patterns make classification challenging.Computer-assisted automated diagnostic tools that enhance ECG signal categorization using sophisticated algorithms and machine learning are helping healthcare practitioners manage greater patient populations.With this motivation,the study proposes a DL framework leveraging the PTB-XL ECG dataset to improve CVD diagnosis.Deep Transfer Learning(DTL)techniques extract features,followed by feature fusion to eliminate redundancy and retain the most informative features.Utilizing the African Vulture Optimization Algorithm(AVOA)for feature selection is more effective than the standard methods,as it offers an ideal balance between exploration and exploitation that results in an optimal set of features,improving classification performance while reducing redundancy.Various machine learning classifiers,including Support Vector Machine(SVM),eXtreme Gradient Boosting(XGBoost),Adaptive Boosting(AdaBoost),and Extreme Learning Machine(ELM),are used for further classification.Additionally,an ensemble model is developed to further improve accuracy.Experimental results demonstrate that the proposed model achieves the highest accuracy of 96.31%,highlighting its effectiveness in enhancing CVD diagnosis.展开更多
We experimentally studied the effect of crack aperture on P-wave velocity, amplitude, anisotropy and dispersion. Experimental models were constructed based on Hudson＇s theory. Six crack models were embedded with equa...We experimentally studied the effect of crack aperture on P-wave velocity, amplitude, anisotropy and dispersion. Experimental models were constructed based on Hudson＇s theory. Six crack models were embedded with equal-radius penny-shaped crack inclusions in each layer. The P-wave velocity and amplitude were measured parallel and perpendicular to the layers of cracks at frequencies of 0.1 MHz to 1 MHz. The experiments show that as the crack aperture increases from 0.l mm to 0.34 mm, the amplitude of the P-waves parallel to the crack layers decreases linearly with increasing frequency and the P-wave velocity dispersion varies from 1.5% to 2.1%, whereas the amplitude of the P-wave perpendicular to the crack layers decreases quadratically with increasing frequency and the velocity dispersion varies from 1.9% to 4.7%. The variation in the velocity dispersion parallel and perpendicular to the cracks intensifies the anisotropy dispersion of the P-waves in the crack models （6.7% to 83%）. The P-wave dispersion strongly depends on the scattering characteristics of the crack apertures.展开更多
Finite-frequency travel time tomography is a newly developing method.The main procedure in this new method is to compute the traveltime sensitive kernel.The travel time of the same scatterer needs to be used for compu...Finite-frequency travel time tomography is a newly developing method.The main procedure in this new method is to compute the traveltime sensitive kernel.The travel time of the same scatterer needs to be used for computing the traveltime sensitive kernel many times.It is a time-consuming task.It is easy and fast to get the travel time from analytic equations in a simple model such as a homogenous or linear velocity media.However,most of the earth models are layered.It is cumbersome to get the travel time from analytic equations.In order to enhance the computation efficiency,we used the table look-up method to compute the finite-frequency travel time sensitive kernel for P-waves in a layered structure model.We chose the AK135 earth model for the velocity model.The table look-up method saved about 50% of the computation time.We enhanced the computation speed by using the table lookup method in the same velocity model,which was very useful for enhancing the computation efficiency for the finite-frequency travel time tomography.展开更多
The recently built China Digital Seismic Network consists of the China National Digital Seismic Network （CNDSN）, 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including...The recently built China Digital Seismic Network consists of the China National Digital Seismic Network （CNDSN）, 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850＋ broadband stations. It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth＇s deep interior besides its routine task of seismic monitoring. Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves. Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation. We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station. An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station. The method provides robust estimates including a measurement error calculated from background noise levels. We found that about one third of the stations have some sort of problems, including misorientation of the two horizontal components, mislabeling and polarity reversal in one or more components. These problems need to be taken into account for any rotation based seismic studies.展开更多
Engineering rock mass classification,based on empirical relations between rock mass parameters and engineering applications,is commonly used in rock engineering and forms the basis for designing rock structures.The ba...Engineering rock mass classification,based on empirical relations between rock mass parameters and engineering applications,is commonly used in rock engineering and forms the basis for designing rock structures.The basic data required may be obtained from visual observation and laboratory or field tests.However,owing to the discontinuous and variable nature of rock masses,it is difficult for rock engineers to directly obtain the specific design parameters needed.As an alternative,the use of geophysical methods in geomechanics such as seismography may largely address this problem.In this study,25 seismic profiles with the total length of 543 m have been scanned to determine the geomechanical properties of the rock mass in blocks Ⅰ,Ⅲ and Ⅳ-2 of the Choghart iron mine.Moreover,rock joint measurements and sampling for laboratory tests were conducted.The results show that the rock mass rating（RMR） and Q values have a close relation with P-wave velocity parameters,including P-wave velocity in field(V;).P-wave velocity in the laboratory(V;) and the ratio of V;V;(i.e.K;= V;/V;.However,Q value,totally,has greater correlation coefficient and less error than the RMR,In addition,rock mass parameters including rock quality designation（RQD）,uniaxial compressive strength（UCS）,joint roughness coefficient（JRC） and Schmidt number（RN） show close relationship with P-wave velocity.An equation based on these parameters was obtained to estimate the P-wave velocity in the rock mass with a correlation coefficient of 91%.The velocities in two orthogonal directions and the results of joint study show that the wave velocity anisotropy in rock mass may be used as an efficient tool to assess the strong and weak directions in rock mass.展开更多
Fast and accurate P-wave arrival picking significantly affects the performance of earthquake early warning(EEW)systems.Automated P-wave picking algorithms used in EEW have encountered problems of falsely picking up no...Fast and accurate P-wave arrival picking significantly affects the performance of earthquake early warning(EEW)systems.Automated P-wave picking algorithms used in EEW have encountered problems of falsely picking up noise,missing P-waves and inaccurate P-wave arrival estimation.To address these issues,an automatic algorithm based on the convolution neural network(DPick)was developed,and trained with a moderate number of data sets of 17,717 accelerograms.Compared to the widely used approach of the short-term average/long-term average of signal characteristic function(STA/LTA),DPick is 1.6 times less likely to detect noise as a P-wave,and 76 times less likely to miss P-waves.In terms of estimating P-wave arrival time,when the detection task is completed within 1 s,DPick′s detection occurrence is 7.4 times that of STA/LTA in the 0.05 s error band,and 1.6 times when the error band is 0.10 s.This verified that the proposed method has the potential for wide applications in EEW.展开更多
It is generally accepted that the uniaxial compressive strength(UCS)and P-wave velocity of rocks tend to decrease simultaneously with increasing temperature.However,based on a great number of statistical data and syst...It is generally accepted that the uniaxial compressive strength(UCS)and P-wave velocity of rocks tend to decrease simultaneously with increasing temperature.However,based on a great number of statistical data and systematic analysis of the microstructure variation of rocks with temperature rising and corresponding propagation mechanism of elastic wave,the results show that(1)There are three different trends for the changes of UCS and P-wave velocity of sandstone when heated from room temperature(20C or 25C)to 800C:(i)Both the UCS and P-wave velocity decrease simultaneously;(ii)The UCS increases initially and then decreases,while the P-wave velocity decreases continuously;and(iii)The UCS increases initially and then fluctuates,while the P-wave velocity continuously decreases.(2)The UCS changes at room temperaturee400C,400Ce600C,and 600Ce800C are mainly attributed to the discrepancy of microstructure characteristics and quartz content,the transformation plasticity of clay minerals,and the balance between the thermal cementation and thermal damage,respectively.(3)The inconsistency in the trends of UCS and P-wave velocity changes is caused by the change of quartz content,phase transition of water and certain minerals.展开更多
If the thickness of coal seams and the lithology of both roofs and floors of coal seams have not changed at all or only a little, then it is thought that the elastic anisotropy of coal seams depends mainly on fracture...If the thickness of coal seams and the lithology of both roofs and floors of coal seams have not changed at all or only a little, then it is thought that the elastic anisotropy of coal seams depends mainly on fractures and obeys the horizontally symmetric model of an azimuth anisotropy. For a fixed offset, the amplitude A of the reflection P-wave and the cosine of 2φ has an approximately linear relation, (φ is the source-detector azimuth with respect to the fracture strike. Based on this relationship, many things can be done, such as the extraction of macro bins, the correction of residual normal moveout, the formation of azimuth gather, the transformation and normalization of azimuth gathers and the extraction of reflection wave amplitudes of coal seams. The least squares method was used to inverse theoretically the direction and density of fractures of coal seams. The result is in good agreement with the regional geological structure, indicating that the azimuth anisotropic analysis of the P-wave is feasible in evaluating the density and direction of fractures in coal seams.展开更多
Model-driven and data-driven inversions are two prominent methods for obtaining P-wave impedance,which is significant in reservoir description and identification.Based on proper initial models,most model-driven method...Model-driven and data-driven inversions are two prominent methods for obtaining P-wave impedance,which is significant in reservoir description and identification.Based on proper initial models,most model-driven methods primarily use the limited frequency bandwidth information of seismic data and can invert P-wave impedance with high accuracy,but not high resolution.Conventional data-driven methods mainly employ the information from well-log data and can provide high-accuracy and highresolution P-wave impedance owing to the superior nonlinear curve fitting capacity of neural networks.However,these methods require a significant number of training samples,which are frequently insufficient.To obtain P-wave impedance with both high accuracy and high resolution,we propose a model-data-driven inversion method using Res Nets and the normalized zero-lag cross-correlation objective function which is effective for avoiding local minima and suppressing random noise.By using initial models and training samples,the proposed model-data-driven method can invert P-wave impedance with satisfactory accuracy and resolution.Tests on synthetic and field data demonstrate the proposed method’s efficacy and practicability.展开更多
Fluid and effective fracture identification in reservoirs is a crucial part of reservoir prediction.The frequency-dependent AVO inversion algorithms have proven to be effective for identifying fluid through its disper...Fluid and effective fracture identification in reservoirs is a crucial part of reservoir prediction.The frequency-dependent AVO inversion algorithms have proven to be effective for identifying fluid through its dispersion property.However,the conventional frequency-dependent AVO inversion algorithms based on Smith&Gidlow and Aki&Richards approximations do not consider the acquisition azimuth of seismic data and neglect the effect of seismic anisotropic dispersion in the actual medium.The aligned fractures in the subsurface medium induce anisotropy.The seismic anisotropy should be considered while accounting for the seismic dispersion properties through fluid-saturated fractured reservoirs.Anisotropy in such reservoirs is frequency-related due to wave-induced fluid-flow(WIFF)between interconnected fractures and pores.It can be used to identify fluid and effective fractures(fluid-saturated)by using azimuthal seismic data via anisotropic dispersion properties.In this paper,based on Rüger’s equation,we derived an analytical expression in the frequency domain for the frequencydependent AVOAz inversion in terms of fracture orientation,dispersion gradient of isotropic background rock,anisotropic dispersion gradient,and the dispersion at a normal incident angle.The frequency-dependent AVOAz equation utilizes azimuthal seismic data and considers the effect of both isotropic and anisotropic dispersion.Reassigned Gabor Transform(RGT)is used to achieve highresolution frequency division data.We then propose the frequency-dependent AVOAz inversion method to identify fluid and characterize effective fractures in fractured porous reservoirs.Through application to high-qualified seismic data of dolomite and carbonate reservoirs,the results show that the method is useful for identifying fluid and effective fractures in fluid-saturated fractured rocks.展开更多
Crack closure is one of the reasons inducing changes of P-wave velocity of rocks under compression.In this context,a method is proposed to investigate the relationships among P-wave velocity,pre-existing cracks,and co...Crack closure is one of the reasons inducing changes of P-wave velocity of rocks under compression.In this context,a method is proposed to investigate the relationships among P-wave velocity,pre-existing cracks,and confining pressure based on the discrete element method(DEM).Pre-existing open cracks inside the rocks are generated by the initial gap of the flat-joint model.The validity of the method is evaluated by comparing the P-wave velocity tested on a sandstone specimen with numerical result.As the crack size is determined by the diameter of particles,the effects of three factors,i.e.number,aspect ratio,and orientation of cracks on the P-wave velocity are discussed.The results show that P-wave velocity is controlled by the(i.e.number) of open micro-cracks,while the closure pressure is determined by the aspect ratio of crack.The reason accounting for the anisotropy of P-wave velocity is the difference in crack number in measurement paths.Both of the number and aspect ratio of cracks can affect the responses of P-wave velocity to the applied confining pressure.Under confining pressure,the number of open cracks inside rocks will dominate the lowest P-wave velocity,and the P-wave velocity of the rock containing narrower cracks is more sensitive to the confining pressure.In this sense,crack density is difficult to be back-calculated merely by P-wave velocity.The proposed method offers a means to analyze the effect of pre-existing cracks on P-wave velocity.展开更多
Laboratory tests were performed on Toyoura sand specimens to investigate the relationship between degree of saturation Sr, B-value and P-wave velocity Vp. Different types of pore water （de-aired water or tap water） ...Laboratory tests were performed on Toyoura sand specimens to investigate the relationship between degree of saturation Sr, B-value and P-wave velocity Vp. Different types of pore water （de-aired water or tap water） and pore gas （air or CO2） as well as different magnitudes of back pressure were used to achieve different Sr （or B-value）. The measured relationship between B-value and Vp was not consistent with the theoretical prediction. The measurement shows that the Vp value in the specimen flushed with de-aired water is independent of B-value （or St） and is always around the one in fully saturated condition. However, the Vp value in the specimen flushed with tap water increases with B-value, but the shape of the relationship between Vp and B-value is quite different from the theoretical prediction. The possible explanation for the discrepancy between laboratory measurement and theoretical prediction lies in that the air exists in the water as air bubbles and therefore the pore fluid （air-water mixture） is heterogeneous instead of homogenous assumed in the theoretical prediction.展开更多
As P-wave velocity is sensitive to the variations in coal reservoir parameters,it is possible to monitor the injected CO_(2)through P-wave velocity during CO_(2)sequestration in coal.However,the effects of CO_(2)on th...As P-wave velocity is sensitive to the variations in coal reservoir parameters,it is possible to monitor the injected CO_(2)through P-wave velocity during CO_(2)sequestration in coal.However,the effects of CO_(2)on the coal P-wave velocity under triaxial stress are not clearly discerned.In the present study,different boundary conditions and gases were utilised to investigate the factors affecting the P-wave velocity after the interaction of coal with CO_(2).Experiments with helium indicated that the pore pressure primarily affected the P-wave velocity by altering the effective stress.Experiments with CH4 and CO_(2)indicated that matrix swelling induced-cleats porosity decline significantly promoted P-wave velocity.Moreover,CO_(2)caused a wider scale and severe weakening of coal matrix than CH4,thereby significantly decreasing the P-wave velocity,and the decline in P-wave velocity increases with vitrinite content.Furthermore,experiments under different boundary conditions showed that with the boundary condition having more constraints,the decrement of pore pressure on P-wave velocity is more weaken,whereas the improvement of matrix swelling on P-wave velocity is more evident.This study contributes to understanding the mechanism of effect of CO_(2)on P-wave velocity under triaxial stress condition and provides guidance for monitoring CO_(2)sequestration in coal.展开更多
基金supported by the National Natural Science Foundation of China(Nos.42574119,42274083,41974049)partly supported by the Urban Geological Survey Project of Linyi,Shandong Province,China(No.SDGP371300202102000468).
文摘A P-wave velocity model was built in the central southern of the Tanlu Fault based on double-difference tomography.The results suggest the presence of a low-velocity anomaly extending from the surface to a depth of 25 km around the Tanlu and Feixi Faults,representing fault-related fluids caused by partial melting.The relocated earthquakes indicate a significant concentration of seismic activity above 20 km around the Tanlu and Feixi Faults,suggesting that prominent fault systems possibly serve as conduits for the upward migration of deep minerals.The proposed geodynamic model,supported by geological and geophysical data,suggests that the migration of deep mineralized materials extends along the Tanlu Fault.The obtained results serve as a crucial foundation for elucidating the intricate process of mineralization in the central southern segment of the Tanlu Fault,thereby enhancing comprehension regarding the interaction among ore body formation,fault fluids,localized melting,and seismic activity.
基金funded by Researchers Supporting ProjectNumber(RSPD2025R947),King Saud University,Riyadh,Saudi Arabia.
文摘The integration of IoT and Deep Learning(DL)has significantly advanced real-time health monitoring and predictive maintenance in prognostic and health management(PHM).Electrocardiograms(ECGs)are widely used for cardiovascular disease(CVD)diagnosis,but fluctuating signal patterns make classification challenging.Computer-assisted automated diagnostic tools that enhance ECG signal categorization using sophisticated algorithms and machine learning are helping healthcare practitioners manage greater patient populations.With this motivation,the study proposes a DL framework leveraging the PTB-XL ECG dataset to improve CVD diagnosis.Deep Transfer Learning(DTL)techniques extract features,followed by feature fusion to eliminate redundancy and retain the most informative features.Utilizing the African Vulture Optimization Algorithm(AVOA)for feature selection is more effective than the standard methods,as it offers an ideal balance between exploration and exploitation that results in an optimal set of features,improving classification performance while reducing redundancy.Various machine learning classifiers,including Support Vector Machine(SVM),eXtreme Gradient Boosting(XGBoost),Adaptive Boosting(AdaBoost),and Extreme Learning Machine(ELM),are used for further classification.Additionally,an ensemble model is developed to further improve accuracy.Experimental results demonstrate that the proposed model achieves the highest accuracy of 96.31%,highlighting its effectiveness in enhancing CVD diagnosis.
基金supported by the Major National Project Program (No.2011ZX05007-006)
文摘We experimentally studied the effect of crack aperture on P-wave velocity, amplitude, anisotropy and dispersion. Experimental models were constructed based on Hudson＇s theory. Six crack models were embedded with equal-radius penny-shaped crack inclusions in each layer. The P-wave velocity and amplitude were measured parallel and perpendicular to the layers of cracks at frequencies of 0.1 MHz to 1 MHz. The experiments show that as the crack aperture increases from 0.l mm to 0.34 mm, the amplitude of the P-waves parallel to the crack layers decreases linearly with increasing frequency and the P-wave velocity dispersion varies from 1.5% to 2.1%, whereas the amplitude of the P-wave perpendicular to the crack layers decreases quadratically with increasing frequency and the velocity dispersion varies from 1.9% to 4.7%. The variation in the velocity dispersion parallel and perpendicular to the cracks intensifies the anisotropy dispersion of the P-waves in the crack models （6.7% to 83%）. The P-wave dispersion strongly depends on the scattering characteristics of the crack apertures.
基金supported by the National Natural Science Foundation of China (Grant No. 90814013)
文摘Finite-frequency travel time tomography is a newly developing method.The main procedure in this new method is to compute the traveltime sensitive kernel.The travel time of the same scatterer needs to be used for computing the traveltime sensitive kernel many times.It is a time-consuming task.It is easy and fast to get the travel time from analytic equations in a simple model such as a homogenous or linear velocity media.However,most of the earth models are layered.It is cumbersome to get the travel time from analytic equations.In order to enhance the computation efficiency,we used the table look-up method to compute the finite-frequency travel time sensitive kernel for P-waves in a layered structure model.We chose the AK135 earth model for the velocity model.The table look-up method saved about 50% of the computation time.We enhanced the computation speed by using the table lookup method in the same velocity model,which was very useful for enhancing the computation efficiency for the finite-frequency travel time tomography.
基金supported by NSF grant EAR-063566(F.N.)National Natural Science Foundation of China grant 40774042(J.L.)
文摘The recently built China Digital Seismic Network consists of the China National Digital Seismic Network （CNDSN）, 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850＋ broadband stations. It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth＇s deep interior besides its routine task of seismic monitoring. Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves. Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation. We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station. An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station. The method provides robust estimates including a measurement error calculated from background noise levels. We found that about one third of the stations have some sort of problems, including misorientation of the two horizontal components, mislabeling and polarity reversal in one or more components. These problems need to be taken into account for any rotation based seismic studies.
文摘Engineering rock mass classification,based on empirical relations between rock mass parameters and engineering applications,is commonly used in rock engineering and forms the basis for designing rock structures.The basic data required may be obtained from visual observation and laboratory or field tests.However,owing to the discontinuous and variable nature of rock masses,it is difficult for rock engineers to directly obtain the specific design parameters needed.As an alternative,the use of geophysical methods in geomechanics such as seismography may largely address this problem.In this study,25 seismic profiles with the total length of 543 m have been scanned to determine the geomechanical properties of the rock mass in blocks Ⅰ,Ⅲ and Ⅳ-2 of the Choghart iron mine.Moreover,rock joint measurements and sampling for laboratory tests were conducted.The results show that the rock mass rating（RMR） and Q values have a close relation with P-wave velocity parameters,including P-wave velocity in field(V;).P-wave velocity in the laboratory(V;) and the ratio of V;V;(i.e.K;= V;/V;.However,Q value,totally,has greater correlation coefficient and less error than the RMR,In addition,rock mass parameters including rock quality designation（RQD）,uniaxial compressive strength（UCS）,joint roughness coefficient（JRC） and Schmidt number（RN） show close relationship with P-wave velocity.An equation based on these parameters was obtained to estimate the P-wave velocity in the rock mass with a correlation coefficient of 91%.The velocities in two orthogonal directions and the results of joint study show that the wave velocity anisotropy in rock mass may be used as an efficient tool to assess the strong and weak directions in rock mass.
基金National Natural Science Foundation of China under Grant Nos.51968016 and 5197083806the Guangxi Innovation Driven Development Project(Science and Technology Major Project,Grant No.Guike AA18118008).
文摘Fast and accurate P-wave arrival picking significantly affects the performance of earthquake early warning(EEW)systems.Automated P-wave picking algorithms used in EEW have encountered problems of falsely picking up noise,missing P-waves and inaccurate P-wave arrival estimation.To address these issues,an automatic algorithm based on the convolution neural network(DPick)was developed,and trained with a moderate number of data sets of 17,717 accelerograms.Compared to the widely used approach of the short-term average/long-term average of signal characteristic function(STA/LTA),DPick is 1.6 times less likely to detect noise as a P-wave,and 76 times less likely to miss P-waves.In terms of estimating P-wave arrival time,when the detection task is completed within 1 s,DPick′s detection occurrence is 7.4 times that of STA/LTA in the 0.05 s error band,and 1.6 times when the error band is 0.10 s.This verified that the proposed method has the potential for wide applications in EEW.
基金This work was supported by the National Natural Science Foundation of China(Grant No.41772333)the program of State Key Laboratory of Frozen Soil Engineering(Grant No.SKLFSE201713)the Shaanxi Province New-Star Talents Promotion Project of Science and Technology(Grant No.2019KJXX-049).
文摘It is generally accepted that the uniaxial compressive strength(UCS)and P-wave velocity of rocks tend to decrease simultaneously with increasing temperature.However,based on a great number of statistical data and systematic analysis of the microstructure variation of rocks with temperature rising and corresponding propagation mechanism of elastic wave,the results show that(1)There are three different trends for the changes of UCS and P-wave velocity of sandstone when heated from room temperature(20C or 25C)to 800C:(i)Both the UCS and P-wave velocity decrease simultaneously;(ii)The UCS increases initially and then decreases,while the P-wave velocity decreases continuously;and(iii)The UCS increases initially and then fluctuates,while the P-wave velocity continuously decreases.(2)The UCS changes at room temperaturee400C,400Ce600C,and 600Ce800C are mainly attributed to the discrepancy of microstructure characteristics and quartz content,the transformation plasticity of clay minerals,and the balance between the thermal cementation and thermal damage,respectively.(3)The inconsistency in the trends of UCS and P-wave velocity changes is caused by the change of quartz content,phase transition of water and certain minerals.
基金Projects 40574058 supported by the National Natural Science Foundation of China, 2005cb221500 the National Key Basic Research and Development(973) Program of China and 03(2007) the Scientific and Technological Project about Geology and Mineral Resources of Henan Land Resources Department
文摘If the thickness of coal seams and the lithology of both roofs and floors of coal seams have not changed at all or only a little, then it is thought that the elastic anisotropy of coal seams depends mainly on fractures and obeys the horizontally symmetric model of an azimuth anisotropy. For a fixed offset, the amplitude A of the reflection P-wave and the cosine of 2φ has an approximately linear relation, (φ is the source-detector azimuth with respect to the fracture strike. Based on this relationship, many things can be done, such as the extraction of macro bins, the correction of residual normal moveout, the formation of azimuth gather, the transformation and normalization of azimuth gathers and the extraction of reflection wave amplitudes of coal seams. The least squares method was used to inverse theoretically the direction and density of fractures of coal seams. The result is in good agreement with the regional geological structure, indicating that the azimuth anisotropic analysis of the P-wave is feasible in evaluating the density and direction of fractures in coal seams.
基金financially supported by the Important National Science&Technology Specific Project of China(Grant No.2017ZX05018-005)
文摘Model-driven and data-driven inversions are two prominent methods for obtaining P-wave impedance,which is significant in reservoir description and identification.Based on proper initial models,most model-driven methods primarily use the limited frequency bandwidth information of seismic data and can invert P-wave impedance with high accuracy,but not high resolution.Conventional data-driven methods mainly employ the information from well-log data and can provide high-accuracy and highresolution P-wave impedance owing to the superior nonlinear curve fitting capacity of neural networks.However,these methods require a significant number of training samples,which are frequently insufficient.To obtain P-wave impedance with both high accuracy and high resolution,we propose a model-data-driven inversion method using Res Nets and the normalized zero-lag cross-correlation objective function which is effective for avoiding local minima and suppressing random noise.By using initial models and training samples,the proposed model-data-driven method can invert P-wave impedance with satisfactory accuracy and resolution.Tests on synthetic and field data demonstrate the proposed method’s efficacy and practicability.
基金supported by the National Major Science and Technology Project of China(2016ZX05004003)the National Natural Science Foundation of China(41574103,41974120,U20B2015)Open Fund of State Key Laboratory of Coal Resources and Safe Mining(Grant No.SKLCRSM19KFA08)。
文摘Fluid and effective fracture identification in reservoirs is a crucial part of reservoir prediction.The frequency-dependent AVO inversion algorithms have proven to be effective for identifying fluid through its dispersion property.However,the conventional frequency-dependent AVO inversion algorithms based on Smith&Gidlow and Aki&Richards approximations do not consider the acquisition azimuth of seismic data and neglect the effect of seismic anisotropic dispersion in the actual medium.The aligned fractures in the subsurface medium induce anisotropy.The seismic anisotropy should be considered while accounting for the seismic dispersion properties through fluid-saturated fractured reservoirs.Anisotropy in such reservoirs is frequency-related due to wave-induced fluid-flow(WIFF)between interconnected fractures and pores.It can be used to identify fluid and effective fractures(fluid-saturated)by using azimuthal seismic data via anisotropic dispersion properties.In this paper,based on Rüger’s equation,we derived an analytical expression in the frequency domain for the frequencydependent AVOAz inversion in terms of fracture orientation,dispersion gradient of isotropic background rock,anisotropic dispersion gradient,and the dispersion at a normal incident angle.The frequency-dependent AVOAz equation utilizes azimuthal seismic data and considers the effect of both isotropic and anisotropic dispersion.Reassigned Gabor Transform(RGT)is used to achieve highresolution frequency division data.We then propose the frequency-dependent AVOAz inversion method to identify fluid and characterize effective fractures in fractured porous reservoirs.Through application to high-qualified seismic data of dolomite and carbonate reservoirs,the results show that the method is useful for identifying fluid and effective fractures in fluid-saturated fractured rocks.
基金partially supported by the National Natural Science Foundation of China(Grant Nos.41872210 and 51674238)the Open Research Fund of the State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences(Grant No.Z017006)。
文摘Crack closure is one of the reasons inducing changes of P-wave velocity of rocks under compression.In this context,a method is proposed to investigate the relationships among P-wave velocity,pre-existing cracks,and confining pressure based on the discrete element method(DEM).Pre-existing open cracks inside the rocks are generated by the initial gap of the flat-joint model.The validity of the method is evaluated by comparing the P-wave velocity tested on a sandstone specimen with numerical result.As the crack size is determined by the diameter of particles,the effects of three factors,i.e.number,aspect ratio,and orientation of cracks on the P-wave velocity are discussed.The results show that P-wave velocity is controlled by the(i.e.number) of open micro-cracks,while the closure pressure is determined by the aspect ratio of crack.The reason accounting for the anisotropy of P-wave velocity is the difference in crack number in measurement paths.Both of the number and aspect ratio of cracks can affect the responses of P-wave velocity to the applied confining pressure.Under confining pressure,the number of open cracks inside rocks will dominate the lowest P-wave velocity,and the P-wave velocity of the rock containing narrower cracks is more sensitive to the confining pressure.In this sense,crack density is difficult to be back-calculated merely by P-wave velocity.The proposed method offers a means to analyze the effect of pre-existing cracks on P-wave velocity.
基金Foundation item： Project（2012CB719803） supported by the National Basic Research Program of China Project（201011159098） supported by the Seed Funding for Basic Research Scheme from The University of Hong Kong, China
文摘Laboratory tests were performed on Toyoura sand specimens to investigate the relationship between degree of saturation Sr, B-value and P-wave velocity Vp. Different types of pore water （de-aired water or tap water） and pore gas （air or CO2） as well as different magnitudes of back pressure were used to achieve different Sr （or B-value）. The measured relationship between B-value and Vp was not consistent with the theoretical prediction. The measurement shows that the Vp value in the specimen flushed with de-aired water is independent of B-value （or St） and is always around the one in fully saturated condition. However, the Vp value in the specimen flushed with tap water increases with B-value, but the shape of the relationship between Vp and B-value is quite different from the theoretical prediction. The possible explanation for the discrepancy between laboratory measurement and theoretical prediction lies in that the air exists in the water as air bubbles and therefore the pore fluid （air-water mixture） is heterogeneous instead of homogenous assumed in the theoretical prediction.
基金supported by the National Natural Science Foundation of China(No.51974304)the Natural Science Foundation of Hebei Province(No.E2020402075)+2 种基金the 2nd Xplorer Prize sponsored by the Tencent Foundationthe Program for Changjiang Scholars and Innovative Research Team in University(No.IRT 17R103)the Qinglan Project of Jiangsu Province.
文摘As P-wave velocity is sensitive to the variations in coal reservoir parameters,it is possible to monitor the injected CO_(2)through P-wave velocity during CO_(2)sequestration in coal.However,the effects of CO_(2)on the coal P-wave velocity under triaxial stress are not clearly discerned.In the present study,different boundary conditions and gases were utilised to investigate the factors affecting the P-wave velocity after the interaction of coal with CO_(2).Experiments with helium indicated that the pore pressure primarily affected the P-wave velocity by altering the effective stress.Experiments with CH4 and CO_(2)indicated that matrix swelling induced-cleats porosity decline significantly promoted P-wave velocity.Moreover,CO_(2)caused a wider scale and severe weakening of coal matrix than CH4,thereby significantly decreasing the P-wave velocity,and the decline in P-wave velocity increases with vitrinite content.Furthermore,experiments under different boundary conditions showed that with the boundary condition having more constraints,the decrement of pore pressure on P-wave velocity is more weaken,whereas the improvement of matrix swelling on P-wave velocity is more evident.This study contributes to understanding the mechanism of effect of CO_(2)on P-wave velocity under triaxial stress condition and provides guidance for monitoring CO_(2)sequestration in coal.
