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.展开更多
This paper addresses the Multi-Vehicle Routing Problem with Time Windows and Simultaneous Pickup and Delivery(MVRPTWSPD),aiming to optimize logistics distribution routes and minimize total costs.A vehicle routing opti...This paper addresses the Multi-Vehicle Routing Problem with Time Windows and Simultaneous Pickup and Delivery(MVRPTWSPD),aiming to optimize logistics distribution routes and minimize total costs.A vehicle routing optimization model is developed based on the operational requirements of the KS Logistics Center,focusing on minimizing vehicle dispatch,loading and unloading,operating,and time window penalty costs.The model incorporates constraints such as vehicle capacity,time windows,and travel distance,and is solved using a genetic algorithm to ensure optimal route planning.Through MATLAB simulations,34 customer points are analyzed,demonstrating that the simultaneous pickup and delivery model reduces total costs by 30.13%,increases vehicle loading rates by 20.04%,and decreases travel distance compared to delivery-only or pickup-only models.The results demonstrate the significant advantages of the simultaneous pickup and delivery mode in reducing logistics costs and improving vehicle utilization,offering valuable insights for enhancing the operational efficiency of the KS Logistics Center.展开更多
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.展开更多
In this paper, a hybrid genetic algorithm （GA） is proposed for the traveling salesman problem （TSP） with pickup and delivery （TSPPD）. In our algorithm, a novel pheromone-based crossover operator is advanced that...In this paper, a hybrid genetic algorithm （GA） is proposed for the traveling salesman problem （TSP） with pickup and delivery （TSPPD）. In our algorithm, a novel pheromone-based crossover operator is advanced that utilizes both local and global information to construct offspring. In addition, a local search procedure is integrated into the GA to accelerate convergence. The proposed GA has been tested on benchmark instances, and the computational results show that it gives better convergence than existing heuristics.展开更多
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.展开更多
The weapon transportation support scheduling problem on aircraft carrier deck is the key to restricting the sortie rate and combat capability of carrier-based aircraft.This paper studies the problem and presents a nov...The weapon transportation support scheduling problem on aircraft carrier deck is the key to restricting the sortie rate and combat capability of carrier-based aircraft.This paper studies the problem and presents a novel solution architecture.Taking the interference of the carrier-based aircraft deck layout on the weapon transportation route and precedence constraint into consideration,a mixed integer formulation is established to minimize the total objective,which is constituted of makespan,load variance and accumulative transfer time of support unit.Solution approach is developed for the model.Firstly,based on modeling the carrier aircraft parked on deck as convex obstacles,the path library of weapon transportation is constructed through visibility graph and Warshall-Floyd methods.We then propose a bi-population immune algorithm in which a population-based forward/backward scheduling technique,local search schemes and a chaotic catastrophe operator are embedded.Besides,the randomkey solution representation and serial scheduling generation scheme are adopted to conveniently obtain a better solution.The Taguchi method is additionally employed to determine key parameters of the algorithm.Finally,on a set of generated realistic instances,we demonstrate that the proposed algorithm outperforms all compared algorithms designed for similar optimization problems and can significantly improve the efficiency,and that the established model and the bi-population immune algorithm can effectively respond to the weapon support requirements of carrier-based aircraft under different sortie missions.展开更多
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.展开更多
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.展开更多
基金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.
文摘This paper addresses the Multi-Vehicle Routing Problem with Time Windows and Simultaneous Pickup and Delivery(MVRPTWSPD),aiming to optimize logistics distribution routes and minimize total costs.A vehicle routing optimization model is developed based on the operational requirements of the KS Logistics Center,focusing on minimizing vehicle dispatch,loading and unloading,operating,and time window penalty costs.The model incorporates constraints such as vehicle capacity,time windows,and travel distance,and is solved using a genetic algorithm to ensure optimal route planning.Through MATLAB simulations,34 customer points are analyzed,demonstrating that the simultaneous pickup and delivery model reduces total costs by 30.13%,increases vehicle loading rates by 20.04%,and decreases travel distance compared to delivery-only or pickup-only models.The results demonstrate the significant advantages of the simultaneous pickup and delivery mode in reducing logistics costs and improving vehicle utilization,offering valuable insights for enhancing the operational efficiency of the KS Logistics Center.
基金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.
文摘In this paper, a hybrid genetic algorithm （GA） is proposed for the traveling salesman problem （TSP） with pickup and delivery （TSPPD）. In our algorithm, a novel pheromone-based crossover operator is advanced that utilizes both local and global information to construct offspring. In addition, a local search procedure is integrated into the GA to accelerate convergence. The proposed GA has been tested on benchmark instances, and the computational results show that it gives better convergence than existing heuristics.
文摘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.
基金the financial support of the National Natural Science Foundation of China(No.52102453)。
文摘The weapon transportation support scheduling problem on aircraft carrier deck is the key to restricting the sortie rate and combat capability of carrier-based aircraft.This paper studies the problem and presents a novel solution architecture.Taking the interference of the carrier-based aircraft deck layout on the weapon transportation route and precedence constraint into consideration,a mixed integer formulation is established to minimize the total objective,which is constituted of makespan,load variance and accumulative transfer time of support unit.Solution approach is developed for the model.Firstly,based on modeling the carrier aircraft parked on deck as convex obstacles,the path library of weapon transportation is constructed through visibility graph and Warshall-Floyd methods.We then propose a bi-population immune algorithm in which a population-based forward/backward scheduling technique,local search schemes and a chaotic catastrophe operator are embedded.Besides,the randomkey solution representation and serial scheduling generation scheme are adopted to conveniently obtain a better solution.The Taguchi method is additionally employed to determine key parameters of the algorithm.Finally,on a set of generated realistic instances,we demonstrate that the proposed algorithm outperforms all compared algorithms designed for similar optimization problems and can significantly improve the efficiency,and that the established model and the bi-population immune algorithm can effectively respond to the weapon support requirements of carrier-based aircraft under different sortie missions.
基金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.
基金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.
