Accurate and rapid determination of source locations is of great significance for surface microseismic monitoring.Traditional methods,such as diffraction stacking,are time-consuming and challenging for real-time monit...Accurate and rapid determination of source locations is of great significance for surface microseismic monitoring.Traditional methods,such as diffraction stacking,are time-consuming and challenging for real-time monitoring.In this study,we propose an approach to locate microseismic events using a deep learning algorithm with surface data.A fully convolutional network is designed to predict source locations.The input data is the waveform of a microseismic event,and the output consists of three 1D Gaussian distributions representing the probability distribution of the source location in the x,y,and z dimensions.The theoretical dataset is generated to train the model,and several data augmentation methods are applied to reduce discrepancies between the theoretical and field data.After applying the trained model to field data,the results demonstrate that our method is fast and achieves comparable location accuracy to the traditional diffraction stacking location method,making it promising for real-time microseismic monitoring.展开更多
A landslide always results from a progressive process of slope deformation. In recent years, an increasing number of slope instabilities have occurred with regard to human engineering activities such as hydropower or ...A landslide always results from a progressive process of slope deformation. In recent years, an increasing number of slope instabilities have occurred with regard to human engineering activities such as hydropower or traffic construction in mountainous area, which cause even greater casualties and economic loss compared with the natural hazards. The development of such earth surface process may hold long period with mechanisms still not fully understood. Using monitoring technology is an effective and intuitive approach to assist analyzing the slope deformation process and their driving factors. This study presents an engineering slope excavated during the construction of Changheba Hydropower Station, which is located in the upper reaches of Dadu River, Sichuan Province, southwest China. The engineering slope experienced and featured a five-year continuous deformation which caused continuous high risks to the engineering activities. We conducted in-depth analysis for such a long-term deformation process based on ground and subsurface monitoring data, collected successive data with a series of monitoring equipment such as automated total station, borehole inclinometers and other auxiliary apparatus, and identified the deformation process based on the comprehensive analysis of monitoring data as well as field investigation. After analyzing the effects of engineering activities and natural factors on the continuous deformation, we found that the overburden strata provided deformable mass while the excavation-produced steep terrain initiated the slope deformation in limit equilibrium state over a long period of time;afterwards, the intense rainwater accelerated slope deformation in the rainy season.展开更多
Microseismic (MS) source location plays an important role in MS monitoring. This paper proposes a MS source location method based on particle swarm optimization (PSO) and multi-sensor arrays, where a free weight joint...Microseismic (MS) source location plays an important role in MS monitoring. This paper proposes a MS source location method based on particle swarm optimization (PSO) and multi-sensor arrays, where a free weight joints the P-wave first arrival data. This method adaptively adjusts the preference for “superior” arrays and leverages “inferior” arrays to escape local optima, thereby improving the location accuracy. The effectiveness and stability of this method were validated through synthetic tests, pencil-lead break (PLB) experiments, and mining engineering applications. Specifically, for synthetic tests with 1 μs Gaussian noise and 100 μs large noise in rock samples, the location error of the multi-sensor arrays jointed location method is only 0.30 cm, which improves location accuracy by 97.51% compared to that using a single sensor array. The average location error of PLB events on three surfaces of a rock sample is reduced by 48.95%, 26.40%, and 55.84%, respectively. For mine blast event tests, the average location error of the dual sensor arrays jointed method is 62.74 m, 54.32% and 14.29% lower than that using only sensor arrays 1 and 2, respectively. In summary, the proposed multi-sensor arrays jointed location method demonstrates good noise resistance, stability, and accuracy, providing a compelling new solution for MS location in relevant mining scenarios.展开更多
A case study of seismic interferometry applied to a small microseismic monitoring network is here presented.The main objectives of this study are(i)to quantify the lateral variability of shear-wave ve-locities in the ...A case study of seismic interferometry applied to a small microseismic monitoring network is here presented.The main objectives of this study are(i)to quantify the lateral variability of shear-wave ve-locities in the studied area,and(ii)to investigate the bias produced by noise directionality and non-stationarity in the velocity estimate.Despite the limited number of stations and the short-period char-acter of the seismic sensors,the empirical Green's functions were retrieved for all station pairs using two years of passive data.Both group and phase velocities were derived,the former using the widespread frequency-time analysis,the latter through the analysis of the real part of the cross-spectra.The main advantage of combining these two methods is a more accurate identification of higher modes,resulting in a reduction of ambiguity during picking and data interpretation.Surface wave tomography was run to obtain the spatial distribution of group and phase velocities for the same wavelengths.The low standard deviation of the results suggests that the sparse character of the network does not limit the applicability of the method,for this specific case.The obtained maps highlight the presence of a lower velocity area that extends from the centre of the network towards southeast.Group and phase velocity dispersion curves have been jointly inverted to retrieve as many shear-wave velocity profiles as selected station pairs.While the average model can be used for a more accurate location of the local natural seismicity,the associated standard deviations give us an indication of the lateral heterogeneity of seismic velocities as a function of depth.Finally,the same velocity analysis was repeated for different time windows in order to quantify the error associated to variations in the noise field.Errors as large as 4%have been found,related to the unfavorable orientation of the receiver pairs with respect to strongly directional noise sources,and to the very short time widows.It was shown that using a one-year time window these errors arereduced to 0.3%.展开更多
Detection of structural changes from an opera- tional process is a major goal in machine condition moni- toring. Existing methods for this purpose are mainly based on retrospective analysis, resulting in a large detec...Detection of structural changes from an opera- tional process is a major goal in machine condition moni- toring. Existing methods for this purpose are mainly based on retrospective analysis, resulting in a large detection delay that limits their usages in real applications. This paper presents a new adaptive real-time change detection algorithm, an extension of the recent research by combin- ing with an incremental sliding-window strategy, to handle the multi-change detection in long-term monitoring of machine operations. In particular, in the framework, Hil- bert space embedding of distribution is used to map the original data into the Re-producing Kernel Hilbert Space （RK_HS） for change detection; then, a new adaptive threshold strategy can be developed when making change decision, in which a global factor （used to control the coarse-to-fine level of detection） is introduced to replace the fixed value of threshold. Through experiments on a range of real testing data which was collected from an experimental rotating machinery system, the excellent detection performances of the algorithm for engineering applications were demonstrated. Compared with state-of- the-art methods, the proposed algorithm can be more suitable for long-term machinery condition monitoring without any manual re-calibration, thus is promising in modern industries.展开更多
In this work,we report long-term trends in the abundance and breeding performance of Adélie penguins(Pygoscelis adeliae)nesting in three Antarctic colonies(i.e.,at Martin Point,South Orkneys Islands;Stranger Poin...In this work,we report long-term trends in the abundance and breeding performance of Adélie penguins(Pygoscelis adeliae)nesting in three Antarctic colonies(i.e.,at Martin Point,South Orkneys Islands;Stranger Point/Cabo Funes,South Shetland Islands;and Esperanza/Hope Bay in the Antarctic Peninsula)from 1995/96 to 2022/23.Using yearly count data of breeding groups selected,we observed a decline in the number of breeding pairs and chicks in crèche at all colonies studied.However,the magnitude of change was higher at Stranger Point than that in the remaining colonies.Moreover,the index of breeding success,which was calculated as the ratio of chicks in crèche to breeding pairs,exhibited no apparent trend throughout the study period.However,it displayed greater variability at Martin Point compared to the other two colonies under investigation.Although the number of chicks in crèche of Adélie penguins showed a declining pattern,the average breeding performance was similar to that reported in gentoo penguin colonies,specifically,those undergoing a population increase(even in sympatric colonies facing similar local conditions).Consequently,it is plausible to assume a reduction of the over-winter survival as a likely cause of the declining trend observed,at least in the Stranger Point and Esperanza colonies.However,we cannot rule out local effects during the breeding season affecting the Adélie population of Martin Point.展开更多
[Objective] This study aimed to investigate the effects of long-term differ- ent fertilization in three types of soils on wheat yield and soil nutrient variation in Shandong Province. [Method] A 30-year located experi...[Objective] This study aimed to investigate the effects of long-term differ- ent fertilization in three types of soils on wheat yield and soil nutrient variation in Shandong Province. [Method] A 30-year located experiment in Jinan of Shandong Province was selected and the results of soil nutrient and crop yield in 1984, 1987, 1988, 1989, 1990, 2001, 2005, 2006, 2007 and 2010 were measured and collected. In this study, five treatments： CK, NP, NK, PK and NPK of the located experiment were selected. [Result] The three types of soils in wheat yields decreased signifi- cantly in the first several years and in 2006. Wheat yields of the treatments with P fertilizers were obviously higher than those without P fertilizers; it was shown that phosphorus is the primary nutritional factor for high-yielding of wheat. The highest yield is from cinnamon soil, followed by that from brown soil, and the lowest pro- duction is from fluvo-aquic soil. Under the same fertilization, the influence of other factors on wheat yield of brown soil is the smallest, while cinnamon soil is vulnera- ble to the influence of external conditions, resulting in larger fluctuation of annual wheat yield. The alkali-hydro nitrogen contents of three kinds of soils first de- creased, then raised, and at last reduced apparently. Since 2007, the change of al- kali-hydro nitrogen content appeared rebounded. The available P contents of no- phosphorus treatments decreased over time while those of the treatments with P fertilizers increased at first, then decreased, and after that kept relatively stable. The available K contents of no K treatments decreased slowly. The treatments of PK and NK had higher available K content than NPK treatment. [Conclusion] Thus, it is an effective fertilization measure to improve the wheat yield by supplying reasonable phosphate fertilizer and nitrogen fertilizer and making up potassium fertilizer.展开更多
To reveal variation of organic matter content in red paddy soil, Through a 27 years-located fertilization experiment in red paddy soil, the content and composi- tion of organic matter in paddy soil were studied in thi...To reveal variation of organic matter content in red paddy soil, Through a 27 years-located fertilization experiment in red paddy soil, the content and composi- tion of organic matter in paddy soil were studied in this paper. The result showed that： the dynamics of soil organic matter of the different fertilization treatments showed significant differences, in the premise of equal nutrient （nitrogen and phos- phorus and potassium）, combining application of organic-inorganic was benefited for the accumulation of organic matter in paddy soil than without fertilizer treatment or chemical fertilizer treatment; the dynamics of soil humic acid, HA and FA of chemi- cal fertilizer only and Combining application of organic-inorganic treatments had basi- cally the same trend, But the contents of humic acid, HA and FA of combining ap- plication of organic-inorganic treatments had been higher than that without fertilizer and chemical fertilizer treatment. Moreover combining application of organic-inorganic treatments was benefited for improving the contents of HA and FA, but decreasing HA/FA ratio as an extension of time. Combining application of organic-inorganic was benefited for improving the contents ol readily oxidizable organic matter. And the contents of soil organic matter in long-term experiment and the contents of readily oxidizable O.M were highly significant positive correlation and the contents of soil total nitrogen, avail nitrogen and potassium were significant positive correlation; the contents of soil readily oxidizable O.M and the contents of soil total nitrogen, avail- able P and rice yield were significant positive correlation. Thereinto, the correlation （r=0.818 1 ） between the rice yield and soil readily oxidizable O.M was higher than the correlation （r=0.802 0） between the rice yield and soil organic matter. It showed the soil readily oxidized organic matter had greater contribution to the rice yield.展开更多
Passive microseismic monitoring(PMM)serves as a fundamental technology for assessing hydraulic fracturing(HF)effectiveness,with a key focus on accurate and efficient phase detection/arrival picking and source location...Passive microseismic monitoring(PMM)serves as a fundamental technology for assessing hydraulic fracturing(HF)effectiveness,with a key focus on accurate and efficient phase detection/arrival picking and source location.In PMM data processing,the data-driven paradigm(deep learning based)outperforms the model-driven paradigm in characteristic extraction but lacks quality control and uncertainty quantification.Monte Carlo Dropout,a Bayesian uncertainty quantification technique,performs stochastic neuron deactivation through multiple forward propagation samplings.Therefore,this study proposes a deep learning neural network incorporating uncertainty quantification with manual quality control integration,establishing an optimized workflow spanning automated phase detection to robust source location.The methodology implementation comprises two principal components:(1)The MDNet employing Monte Carlo Dropout strategy enabling simultaneous phase detection/arrival picking and unce rtainty estimation;(2)an integrated hybrid-driven workflow with a traveltime-based inve rsion method for source location.Validation with field data demonstrates that MD-Net achieves superior performance under low signal-to-noise ratio conditions,maintaining detection accuracy exceeding 99%for both P-and S-waves.The phase arrival picking precision shows significant improvement,with a 40%reduction in standard deviation compared to the baseline model(P-S time difference decreasing from12.0 ms to 7.1 ms),while providing quantifiable uncertainty metrics for manual calibration.Source location results further reveal that our hybrid-driven workflow produces more physically plausible event distributions,with 100%of microseismic eve nts clustering along the primary fracture expanding direction.This performance surpasses traditional cross-correlation methods and single/multi-trace data-driven me thods in spatial rationality.This study establishes an inte rpretable,high-pre cision automated framework for HF-PMM applications,demonstrating potential for extension to diverse geological settings and monitoring configurations.展开更多
Knowledge of the locations of seismic sources is critical for microseismic monitoring. Time-window-based elastic wave interferometric imaging and weighted- elastic-wave （WEW） interferometric imaging are proposed and...Knowledge of the locations of seismic sources is critical for microseismic monitoring. Time-window-based elastic wave interferometric imaging and weighted- elastic-wave （WEW） interferometric imaging are proposed and used to locate modeled microseismic sources. The proposed method improves the precision and eliminates artifacts in location profiles. Numerical experiments based on a horizontally layered isotropic medium have shown that the method offers the following advantages： It can deal with Iow-SNR microseismic data with velocity perturbations as well as relatively sparse receivers and still maintain relatively high precision despite the errors in the velocity model. Furthermore, it is more efficient than conventional traveltime inversion methods because interferometric imaging does not require traveltime picking. Numerical results using a 2D fault model have also suggested that the weighted-elastic-wave interferometric imaging can locate multiple sources with higher location precision than the time-reverse imaging method.展开更多
Downhole microseismic data has the significant advantages of high signal-to-noise ratio and well-developed P and S waves and the core component of microseismic monitoring is microseismic event location associated with...Downhole microseismic data has the significant advantages of high signal-to-noise ratio and well-developed P and S waves and the core component of microseismic monitoring is microseismic event location associated with hydraulic fracturing in a relatively high confidence level and accuracy.In this study,we present a multidimensional DIRECT inversion method for microseismic locations and applicability tests over modeling data based on a downhole microseismic monitoring system.Synthetic tests inidcate that the objective function of locations can be defined as a multi-dimensional matrix space by employing the global optimization DIRECT algorithm,because it can be run without the initial value and objective function derivation,and the discretely scattered objective points lead to an expeditious contraction of objective functions in each dimension.This study shows that the DIRECT algorithm can be extensively applied in real downhole microseismic monitoring data from hydraulic fracturing completions.Therefore,the methodology,based on a multidimensional DIRECT algorithm,can provide significant high accuracy and convergent efficiency as well as robust computation for interpretable spatiotemporal microseismic evolution,which is more suitable for real-time processing of a large amount of downhole microseismic monitoring data.展开更多
Micro-seismic phenomena, occurring when rock masses are subjected to forces and failures, allow the determination of their unstable states and failure zones by analyzing micro-seismic signals. We first present the pri...Micro-seismic phenomena, occurring when rock masses are subjected to forces and failures, allow the determination of their unstable states and failure zones by analyzing micro-seismic signals. We first present the principles of micro-seismic monitoring and location, as well as an underground explosion-proof micro-seismic monitoring system. Given a practical engineering application, we describe the application of micro-seismic monitoring technology in determining the height of a "two-zone" overburden, i.e., a caving zone and a fracture zone, the width of a coal-pillar section and the depth of failure of a floor. The workfaces monitored accomplished safe and highly efficient mining based on our micro-seismic monitoring results and provide direct proof of the reliability and validity of micro-seismic monitoring technology.展开更多
Locating seismic events is a central task for earthquake monitoring.Compared to arrival-based location methods,waveformbased location methods do not require picking phase arrivals and are more suitable for locating se...Locating seismic events is a central task for earthquake monitoring.Compared to arrival-based location methods,waveformbased location methods do not require picking phase arrivals and are more suitable for locating seismic events with noisy waveforms.Among waveform-based location methods,one approach is to stack different attributes of P and S waveforms around arrival times corresponding to potential event locations and origin times,and the maximum stacking values are assumed to indicate the correct event location and origin time.In this study,to obtain a high-resolution location image,we improve the waveform-based location method by applying a hybrid multiplicative imaging condition to characteristic functions of seismic waveforms.In our new stacking method,stations are divided into groups;characteristic functions of seismic waveforms recorded at stations in the same group are summed,and then multiplied among groups.We find that this approach can largely eliminate the cumulative effects of noise in the summation process and thus improve the resolution of location images.We test the new method and compare it to three other stacking methods,using both synthetic and real datasets that are related to induced seismicity occurring in petroleum/gas production.The test results confirm that the new stacking method can provide higher-resolution location images than those derived from currently used methods.展开更多
The earthquake real-time monitoring system of the Chinese National Digital Seismic Network has been in operation since"the Ninth Five-year Plan"period,and the stability of the system has been well tested.In ...The earthquake real-time monitoring system of the Chinese National Digital Seismic Network has been in operation since"the Ninth Five-year Plan"period,and the stability of the system has been well tested.In recent years,with the continuous improvement of monitoring technology and increase of public demands,the original real-time monitoring system needs to be upgraded and improved in terms of timeliness,stability,accuracy and ease of operation.Therefore,by accessing a total of more than 1,000 seismic stations,reducing the seismic trigger threshold of the monitoring system,eliminating the false trigger stations and optimizing the seismic waveform display interface,the current earthquake monitoring demands can be satisfied on the basis of ensuring the stable operation of the system.展开更多
Under two rock strata combination conditions, over 10,000 microseismic events were received with microseismic location monitoring technology which possessed by the author＇s studying team, used in fully mechanized coa...Under two rock strata combination conditions, over 10,000 microseismic events were received with microseismic location monitoring technology which possessed by the author＇s studying team, used in fully mechanized coal face of Huafeng Mine of Xinwen Coal Mining Group Co., Shandong Province. On the basis of the achievement of the location results, the conclusions were drawn as follows： On the basis of the achievement of 3D strata fracturing situation and the section plane of microseimic events in different areas, the relationship between spatial structure of overlying strata and mining pressure field was found, and we might describe distribution range of dynamic pressure of advance pressure and lateral stress around long face, and range of structure ad-tivation. Quantitative guidance to prevent dynamic disasters was provided. The practice in coal mine got a effective results. According to the FLAG3D soft numerical simulation of diameter drilling hole （the diameter is 300 mm） to relieve pressure in specified geological condition in Huafeng Mine, the right distance of two dirlls is 2.5 m and the right depth is 12 m. The research pro-vided basic guiding and practical experiences for the underground microseismic monitoring and disaster prevention in side slopes or tunnels engineering.展开更多
To trace the potential hazards of open-pit slope in Longshou mine,global positioning system（GPS） is applied to monitoring ground movement and deformation induced by transition from open-pit to underground mining.Thr...To trace the potential hazards of open-pit slope in Longshou mine,global positioning system（GPS） is applied to monitoring ground movement and deformation induced by transition from open-pit to underground mining.Through long-term monitoring from 2003 to 2008,huge amounts of data were acquired.Monitoring results show that large-scale ground movement and deformation have occurred in mining area,and the movement area is ellipse-shaped.The displacement boundary of settlement trough is 2.0 km long along the exploratory line,and 1.5 km long along the strike of ore body.GPS monitoring results basically agree with the practical deformation state of open-pit slope.It is indicated that the long-term GPS monitoring is an effective way to understand the mechanism of ground movement and deformation in mine area. 更多展开更多
This study was conducted to investigate the effects of long-term located fertilization on soil phosphorus,the changes of soil available phosphorus(OlsenP),the evolution of soil total phosphorus(TP)and the ratio change...This study was conducted to investigate the effects of long-term located fertilization on soil phosphorus,the changes of soil available phosphorus(OlsenP),the evolution of soil total phosphorus(TP)and the ratio change of Olsen-P to TP(PAC)by 33-year fertilization experiments in winter wheat-summer maize rotation system in Shandong fluvo-aquic soil.Eight treatments were designed as no fertilization(CK),nitrogen fertilizer(N),nitrogen and phosphate fertilizer(NP),nitrogen and potassium fertilizer(NK),phosphate and potassium fertilizer(PK),nitrogen-phosphate-potassium fertilizer(NPK),reduced NPK fertilizer(N(15)PK),and increased NPK fertilizer(N_(25)PK).Meanwhile,eight organic fertilizer-added treatments were designed based on the application of inorganic fertilizer the same as the above ones.The results showed that TP,Olsen-P and PAC of treatments added with organic fertilizer were higher than those without organic fertilizer,and those of the treatments applied with phosphate fertilizer were higher than those of no phosphate fertilizer.With the increase of years,soil P pool decreased due to crop absorption,nutrient loss and morphological transformation and other causes under the treatments of without and only phosphate fertilizer,while remained stable under the treatments added with organic fertilizer.The PAC values were generally lower in fluvo-aquic soil,and it could be improved by the application of organic fertilizer.On the whole,the application of chemical phosphate fertilizer combined with organic fertilizer could improve the phosphorus content in soil and ensure the supply of phosphorus nutrition.This study would provide scientific basis for fertilization management and soil fertility in fluvo-aquic soil.展开更多
A systematic study concerning the effects of a long-term stationary fertilization on content and property of soil humus in fluvo-aquic soil sampled from Malan Farm, Xinji City, Hebei, and arid red soil and paddy red s...A systematic study concerning the effects of a long-term stationary fertilization on content and property of soil humus in fluvo-aquic soil sampled from Malan Farm, Xinji City, Hebei, and arid red soil and paddy red soil sampled from the Institute of Red Soil, Jinxian County, Jiangxi was conducted. The results showed that long-term fertilization had effects not only on the content and composition of soil humus, but also on the physico-chemical property of humus. With applying organic manure or combined application of organic manure and chemical fertilizer, E4 and E6 values of humic acid decreased in fluvo-aquic soil and arid red soil, but increased in paddy red soil. In paddy red soil, E4 and E6 values of humic acid increased also with a single application of chemical fertilizer, but E4 and E6 values had less change of humic acid in fluvo-aquic soil and arid red soil. The effects on the visible spectroscopic property of f ulvic acid were different from that of humic acid. Long-term application of organic manure or combined application of organic manure and chemical fertilizer could increase E4 and E6 values of f ulvic acid in three types of soil. Single application of chemical fertilizer had less effect on the E4 and E6. Long-term fertilization could also influence the ultraviolet spectroscopic property of humus. With a single application of organic manure or combined application of organic manure and chemical fertilizer, the ultraviolet absorbance of humic acid and f ulvic acid increased in the three types of soil. But this effect was obvious only in short wave length, and the effect could decrease if the wave length increased. With a single application of chemical fertilizer the ultraviolet absorbance of fulvic acid could increase, but it of humic acid increased only in fluvo-aquic soil. Long-term application of organic manure or combined application of organic manure and chemical fertilizer could increase the content of total acidic groups, carboxy groups and phenolic hydroxy groups of humic acid and fuvic acid in the three types of soil. Single application of chemical fertilizer had less effect on the content of total acidic groups, carboxy groups and phenolic hydroxy groups of humic acid and fuvic acid in the three types of soil.展开更多
This study was aimed at assessing the locational compliance status of petroleum handling facilities in the Niger Delta to the specifications of the Nigerian Upstream Petroleum Regulatory Commission (NUPRC) locational ...This study was aimed at assessing the locational compliance status of petroleum handling facilities in the Niger Delta to the specifications of the Nigerian Upstream Petroleum Regulatory Commission (NUPRC) locational standards. A cross-sectional research design was employed using a Standard checklist of seventeen compliance specifications of NUPRC. A walk-through survey was carried out in 118 identified facilities in three locations: Eket in Akwa-Ibom State, Port Harcourt in Rivers State and Warri in Delta State, respectively. The data obtained were analyzed using inferential and descriptive statistics, Kendall Tau-B and Principal Component analyses. The results indicate that the Petroleum Product facilities complied with 6 (35.29%) out of 17 specifications while 11 (64.70%) specifications were violated. Locations compliance of stations in the Niger Delta region is generally poor, but Stations in Port Harcourt recorded a higher compliance rate compared to Eket and Warri. Results of Kendall’s tau-b and Principal Component analysis indicated positive association between all the land space locational compliances. This study attributes the poor compliance rate to rapid urbanization, overpopulation, proliferation of filling stations and poor monitoring by regulatory agencies. It is recommended that, NUPRC should improve its monitoring and enforce regulatory operational specifications.展开更多
基金supported by National Natural Science Foundation of China Grant(No.42004040,42474092,U2239204,and 42304145)Natural Science Foundation of Jiangxi Province Grant(20242BAB25190 and 20232BAB213077).
文摘Accurate and rapid determination of source locations is of great significance for surface microseismic monitoring.Traditional methods,such as diffraction stacking,are time-consuming and challenging for real-time monitoring.In this study,we propose an approach to locate microseismic events using a deep learning algorithm with surface data.A fully convolutional network is designed to predict source locations.The input data is the waveform of a microseismic event,and the output consists of three 1D Gaussian distributions representing the probability distribution of the source location in the x,y,and z dimensions.The theoretical dataset is generated to train the model,and several data augmentation methods are applied to reduce discrepancies between the theoretical and field data.After applying the trained model to field data,the results demonstrate that our method is fast and achieves comparable location accuracy to the traditional diffraction stacking location method,making it promising for real-time microseismic monitoring.
基金funded by the Second Tibetan Plateau Scientific Expedition and Research Program(STEP)(2019QZKK0904)National Natural Science Foundation of China(42077266,41825018,42090051,41941018,41902289)Strategic Priority Research Program of the Chinese Academy of Sciences(XDA23090402)。
文摘A landslide always results from a progressive process of slope deformation. In recent years, an increasing number of slope instabilities have occurred with regard to human engineering activities such as hydropower or traffic construction in mountainous area, which cause even greater casualties and economic loss compared with the natural hazards. The development of such earth surface process may hold long period with mechanisms still not fully understood. Using monitoring technology is an effective and intuitive approach to assist analyzing the slope deformation process and their driving factors. This study presents an engineering slope excavated during the construction of Changheba Hydropower Station, which is located in the upper reaches of Dadu River, Sichuan Province, southwest China. The engineering slope experienced and featured a five-year continuous deformation which caused continuous high risks to the engineering activities. We conducted in-depth analysis for such a long-term deformation process based on ground and subsurface monitoring data, collected successive data with a series of monitoring equipment such as automated total station, borehole inclinometers and other auxiliary apparatus, and identified the deformation process based on the comprehensive analysis of monitoring data as well as field investigation. After analyzing the effects of engineering activities and natural factors on the continuous deformation, we found that the overburden strata provided deformable mass while the excavation-produced steep terrain initiated the slope deformation in limit equilibrium state over a long period of time;afterwards, the intense rainwater accelerated slope deformation in the rainy season.
基金Project(SICGM2023301) supported by the State Key Laboratory of Strata Intelligent Control and Green Mining Co-founded by Shandong Province and the Ministry of Science and Technology,ChinaProject(SMDPC202202) supported by the Key Laboratory of Mining Disaster Prevention and Control,ChinaProject(U21A2030) supported by the National Natural Science Foundation of China。
文摘Microseismic (MS) source location plays an important role in MS monitoring. This paper proposes a MS source location method based on particle swarm optimization (PSO) and multi-sensor arrays, where a free weight joints the P-wave first arrival data. This method adaptively adjusts the preference for “superior” arrays and leverages “inferior” arrays to escape local optima, thereby improving the location accuracy. The effectiveness and stability of this method were validated through synthetic tests, pencil-lead break (PLB) experiments, and mining engineering applications. Specifically, for synthetic tests with 1 μs Gaussian noise and 100 μs large noise in rock samples, the location error of the multi-sensor arrays jointed location method is only 0.30 cm, which improves location accuracy by 97.51% compared to that using a single sensor array. The average location error of PLB events on three surfaces of a rock sample is reduced by 48.95%, 26.40%, and 55.84%, respectively. For mine blast event tests, the average location error of the dual sensor arrays jointed method is 62.74 m, 54.32% and 14.29% lower than that using only sensor arrays 1 and 2, respectively. In summary, the proposed multi-sensor arrays jointed location method demonstrates good noise resistance, stability, and accuracy, providing a compelling new solution for MS location in relevant mining scenarios.
基金This study is part of the project 2021RUAPON-REACT EU-Finanziamento PON“Ricerca e Innovazione”20142020,grant n.19-G-12543-2,funded by the Italian Ministry of University and Research(MUR)This study was developed in the frame of“The Geosciences for Sustainable Development”project(Budget Minis-tero dell'Universita e della Ricerca-Dipartimenti di Eccellenza 2023-2027,code n.C93C23002690001).
文摘A case study of seismic interferometry applied to a small microseismic monitoring network is here presented.The main objectives of this study are(i)to quantify the lateral variability of shear-wave ve-locities in the studied area,and(ii)to investigate the bias produced by noise directionality and non-stationarity in the velocity estimate.Despite the limited number of stations and the short-period char-acter of the seismic sensors,the empirical Green's functions were retrieved for all station pairs using two years of passive data.Both group and phase velocities were derived,the former using the widespread frequency-time analysis,the latter through the analysis of the real part of the cross-spectra.The main advantage of combining these two methods is a more accurate identification of higher modes,resulting in a reduction of ambiguity during picking and data interpretation.Surface wave tomography was run to obtain the spatial distribution of group and phase velocities for the same wavelengths.The low standard deviation of the results suggests that the sparse character of the network does not limit the applicability of the method,for this specific case.The obtained maps highlight the presence of a lower velocity area that extends from the centre of the network towards southeast.Group and phase velocity dispersion curves have been jointly inverted to retrieve as many shear-wave velocity profiles as selected station pairs.While the average model can be used for a more accurate location of the local natural seismicity,the associated standard deviations give us an indication of the lateral heterogeneity of seismic velocities as a function of depth.Finally,the same velocity analysis was repeated for different time windows in order to quantify the error associated to variations in the noise field.Errors as large as 4%have been found,related to the unfavorable orientation of the receiver pairs with respect to strongly directional noise sources,and to the very short time widows.It was shown that using a one-year time window these errors arereduced to 0.3%.
基金Supported by National Natural Science Foundation of China(Grant Nos.61403232,61327003)Shandong Provincial Natural Science Foundation of China(Grant No.ZR2014FQ025)Young Scholars Program of Shandong University,China(YSPSDU,2015WLJH30)
文摘Detection of structural changes from an opera- tional process is a major goal in machine condition moni- toring. Existing methods for this purpose are mainly based on retrospective analysis, resulting in a large detection delay that limits their usages in real applications. This paper presents a new adaptive real-time change detection algorithm, an extension of the recent research by combin- ing with an incremental sliding-window strategy, to handle the multi-change detection in long-term monitoring of machine operations. In particular, in the framework, Hil- bert space embedding of distribution is used to map the original data into the Re-producing Kernel Hilbert Space （RK_HS） for change detection; then, a new adaptive threshold strategy can be developed when making change decision, in which a global factor （used to control the coarse-to-fine level of detection） is introduced to replace the fixed value of threshold. Through experiments on a range of real testing data which was collected from an experimental rotating machinery system, the excellent detection performances of the algorithm for engineering applications were demonstrated. Compared with state-of- the-art methods, the proposed algorithm can be more suitable for long-term machinery condition monitoring without any manual re-calibration, thus is promising in modern industries.
基金Nacional de Promoción Científica y Tecnológica(Grant:PICTO 2010-0111)the Instituto Antártico Argentino-Dirección Nacional del Antártico(PINST-05)provided financial and logistical support.
文摘In this work,we report long-term trends in the abundance and breeding performance of Adélie penguins(Pygoscelis adeliae)nesting in three Antarctic colonies(i.e.,at Martin Point,South Orkneys Islands;Stranger Point/Cabo Funes,South Shetland Islands;and Esperanza/Hope Bay in the Antarctic Peninsula)from 1995/96 to 2022/23.Using yearly count data of breeding groups selected,we observed a decline in the number of breeding pairs and chicks in crèche at all colonies studied.However,the magnitude of change was higher at Stranger Point than that in the remaining colonies.Moreover,the index of breeding success,which was calculated as the ratio of chicks in crèche to breeding pairs,exhibited no apparent trend throughout the study period.However,it displayed greater variability at Martin Point compared to the other two colonies under investigation.Although the number of chicks in crèche of Adélie penguins showed a declining pattern,the average breeding performance was similar to that reported in gentoo penguin colonies,specifically,those undergoing a population increase(even in sympatric colonies facing similar local conditions).Consequently,it is plausible to assume a reduction of the over-winter survival as a likely cause of the declining trend observed,at least in the Stranger Point and Esperanza colonies.However,we cannot rule out local effects during the breeding season affecting the Adélie population of Martin Point.
基金Supported by the Special Fund for Agro-scientific Research in the Public Interest of China(201203030,201203050)Special Fund for "Taishan Scholar" Construction Engineering "Agricultural Nonpoint Source Pollution Prevention and Control"~~
文摘[Objective] This study aimed to investigate the effects of long-term differ- ent fertilization in three types of soils on wheat yield and soil nutrient variation in Shandong Province. [Method] A 30-year located experiment in Jinan of Shandong Province was selected and the results of soil nutrient and crop yield in 1984, 1987, 1988, 1989, 1990, 2001, 2005, 2006, 2007 and 2010 were measured and collected. In this study, five treatments： CK, NP, NK, PK and NPK of the located experiment were selected. [Result] The three types of soils in wheat yields decreased signifi- cantly in the first several years and in 2006. Wheat yields of the treatments with P fertilizers were obviously higher than those without P fertilizers; it was shown that phosphorus is the primary nutritional factor for high-yielding of wheat. The highest yield is from cinnamon soil, followed by that from brown soil, and the lowest pro- duction is from fluvo-aquic soil. Under the same fertilization, the influence of other factors on wheat yield of brown soil is the smallest, while cinnamon soil is vulnera- ble to the influence of external conditions, resulting in larger fluctuation of annual wheat yield. The alkali-hydro nitrogen contents of three kinds of soils first de- creased, then raised, and at last reduced apparently. Since 2007, the change of al- kali-hydro nitrogen content appeared rebounded. The available P contents of no- phosphorus treatments decreased over time while those of the treatments with P fertilizers increased at first, then decreased, and after that kept relatively stable. The available K contents of no K treatments decreased slowly. The treatments of PK and NK had higher available K content than NPK treatment. [Conclusion] Thus, it is an effective fertilization measure to improve the wheat yield by supplying reasonable phosphate fertilizer and nitrogen fertilizer and making up potassium fertilizer.
基金Supported by National Key Technology R&D Program(2006BAD02A04)~~
文摘To reveal variation of organic matter content in red paddy soil, Through a 27 years-located fertilization experiment in red paddy soil, the content and composi- tion of organic matter in paddy soil were studied in this paper. The result showed that： the dynamics of soil organic matter of the different fertilization treatments showed significant differences, in the premise of equal nutrient （nitrogen and phos- phorus and potassium）, combining application of organic-inorganic was benefited for the accumulation of organic matter in paddy soil than without fertilizer treatment or chemical fertilizer treatment; the dynamics of soil humic acid, HA and FA of chemi- cal fertilizer only and Combining application of organic-inorganic treatments had basi- cally the same trend, But the contents of humic acid, HA and FA of combining ap- plication of organic-inorganic treatments had been higher than that without fertilizer and chemical fertilizer treatment. Moreover combining application of organic-inorganic treatments was benefited for improving the contents of HA and FA, but decreasing HA/FA ratio as an extension of time. Combining application of organic-inorganic was benefited for improving the contents ol readily oxidizable organic matter. And the contents of soil organic matter in long-term experiment and the contents of readily oxidizable O.M were highly significant positive correlation and the contents of soil total nitrogen, avail nitrogen and potassium were significant positive correlation; the contents of soil readily oxidizable O.M and the contents of soil total nitrogen, avail- able P and rice yield were significant positive correlation. Thereinto, the correlation （r=0.818 1 ） between the rice yield and soil readily oxidizable O.M was higher than the correlation （r=0.802 0） between the rice yield and soil organic matter. It showed the soil readily oxidized organic matter had greater contribution to the rice yield.
基金funded by the Deep Earth Probe and Mineral Resources Exploration-National Science and Technology Major Project(Grant No.2024ZD1002503)。
文摘Passive microseismic monitoring(PMM)serves as a fundamental technology for assessing hydraulic fracturing(HF)effectiveness,with a key focus on accurate and efficient phase detection/arrival picking and source location.In PMM data processing,the data-driven paradigm(deep learning based)outperforms the model-driven paradigm in characteristic extraction but lacks quality control and uncertainty quantification.Monte Carlo Dropout,a Bayesian uncertainty quantification technique,performs stochastic neuron deactivation through multiple forward propagation samplings.Therefore,this study proposes a deep learning neural network incorporating uncertainty quantification with manual quality control integration,establishing an optimized workflow spanning automated phase detection to robust source location.The methodology implementation comprises two principal components:(1)The MDNet employing Monte Carlo Dropout strategy enabling simultaneous phase detection/arrival picking and unce rtainty estimation;(2)an integrated hybrid-driven workflow with a traveltime-based inve rsion method for source location.Validation with field data demonstrates that MD-Net achieves superior performance under low signal-to-noise ratio conditions,maintaining detection accuracy exceeding 99%for both P-and S-waves.The phase arrival picking precision shows significant improvement,with a 40%reduction in standard deviation compared to the baseline model(P-S time difference decreasing from12.0 ms to 7.1 ms),while providing quantifiable uncertainty metrics for manual calibration.Source location results further reveal that our hybrid-driven workflow produces more physically plausible event distributions,with 100%of microseismic eve nts clustering along the primary fracture expanding direction.This performance surpasses traditional cross-correlation methods and single/multi-trace data-driven me thods in spatial rationality.This study establishes an inte rpretable,high-pre cision automated framework for HF-PMM applications,demonstrating potential for extension to diverse geological settings and monitoring configurations.
基金supported by the R&D of Key Instruments and Technologies for Deep Resources Prospecting(No.ZDYZ2012-1)National Natural Science Foundation of China(No.11374322)
文摘Knowledge of the locations of seismic sources is critical for microseismic monitoring. Time-window-based elastic wave interferometric imaging and weighted- elastic-wave （WEW） interferometric imaging are proposed and used to locate modeled microseismic sources. The proposed method improves the precision and eliminates artifacts in location profiles. Numerical experiments based on a horizontally layered isotropic medium have shown that the method offers the following advantages： It can deal with Iow-SNR microseismic data with velocity perturbations as well as relatively sparse receivers and still maintain relatively high precision despite the errors in the velocity model. Furthermore, it is more efficient than conventional traveltime inversion methods because interferometric imaging does not require traveltime picking. Numerical results using a 2D fault model have also suggested that the weighted-elastic-wave interferometric imaging can locate multiple sources with higher location precision than the time-reverse imaging method.
基金financially supported by the National Natural Science Foundation of China (Grant No. 41807296 and No. 41802006)Natural science found for universities of Anhui province (Grant No. KJ2017A036)
文摘Downhole microseismic data has the significant advantages of high signal-to-noise ratio and well-developed P and S waves and the core component of microseismic monitoring is microseismic event location associated with hydraulic fracturing in a relatively high confidence level and accuracy.In this study,we present a multidimensional DIRECT inversion method for microseismic locations and applicability tests over modeling data based on a downhole microseismic monitoring system.Synthetic tests inidcate that the objective function of locations can be defined as a multi-dimensional matrix space by employing the global optimization DIRECT algorithm,because it can be run without the initial value and objective function derivation,and the discretely scattered objective points lead to an expeditious contraction of objective functions in each dimension.This study shows that the DIRECT algorithm can be extensively applied in real downhole microseismic monitoring data from hydraulic fracturing completions.Therefore,the methodology,based on a multidimensional DIRECT algorithm,can provide significant high accuracy and convergent efficiency as well as robust computation for interpretable spatiotemporal microseismic evolution,which is more suitable for real-time processing of a large amount of downhole microseismic monitoring data.
基金supported by the National Basic Research Program of China (No. 2010CB226805)the National Natural Science Foundation of China (No. 50874103)+1 种基金the State Laboratory Fund (No. SKLGDUEK0905)the Talent Introduction Fund of Anhui University of Science & Technology
文摘Micro-seismic phenomena, occurring when rock masses are subjected to forces and failures, allow the determination of their unstable states and failure zones by analyzing micro-seismic signals. We first present the principles of micro-seismic monitoring and location, as well as an underground explosion-proof micro-seismic monitoring system. Given a practical engineering application, we describe the application of micro-seismic monitoring technology in determining the height of a "two-zone" overburden, i.e., a caving zone and a fracture zone, the width of a coal-pillar section and the depth of failure of a floor. The workfaces monitored accomplished safe and highly efficient mining based on our micro-seismic monitoring results and provide direct proof of the reliability and validity of micro-seismic monitoring technology.
基金supported by National Key R&D Program of China under grant 2018YFC1504102the National Natural Science Foundation of China under grant 41961134001。
文摘Locating seismic events is a central task for earthquake monitoring.Compared to arrival-based location methods,waveformbased location methods do not require picking phase arrivals and are more suitable for locating seismic events with noisy waveforms.Among waveform-based location methods,one approach is to stack different attributes of P and S waveforms around arrival times corresponding to potential event locations and origin times,and the maximum stacking values are assumed to indicate the correct event location and origin time.In this study,to obtain a high-resolution location image,we improve the waveform-based location method by applying a hybrid multiplicative imaging condition to characteristic functions of seismic waveforms.In our new stacking method,stations are divided into groups;characteristic functions of seismic waveforms recorded at stations in the same group are summed,and then multiplied among groups.We find that this approach can largely eliminate the cumulative effects of noise in the summation process and thus improve the resolution of location images.We test the new method and compare it to three other stacking methods,using both synthetic and real datasets that are related to induced seismicity occurring in petroleum/gas production.The test results confirm that the new stacking method can provide higher-resolution location images than those derived from currently used methods.
基金the China Earthquake Network Center Seismic Network Department Daily Operation and Maintenance Funding Support(1950411001)
文摘The earthquake real-time monitoring system of the Chinese National Digital Seismic Network has been in operation since"the Ninth Five-year Plan"period,and the stability of the system has been well tested.In recent years,with the continuous improvement of monitoring technology and increase of public demands,the original real-time monitoring system needs to be upgraded and improved in terms of timeliness,stability,accuracy and ease of operation.Therefore,by accessing a total of more than 1,000 seismic stations,reducing the seismic trigger threshold of the monitoring system,eliminating the false trigger stations and optimizing the seismic waveform display interface,the current earthquake monitoring demands can be satisfied on the basis of ensuring the stable operation of the system.
基金Supported by the National Natural Science Foundation of China(40674017,50534080)the Open Fund of State Key Laboratory of Coal Resources and Safe Mining(2007-04)+1 种基金the Open Doctor Innovation Fund of Shandong Province(200703020)the Chinese Postdoctoral Science Fundation(20080440304)
文摘Under two rock strata combination conditions, over 10,000 microseismic events were received with microseismic location monitoring technology which possessed by the author＇s studying team, used in fully mechanized coal face of Huafeng Mine of Xinwen Coal Mining Group Co., Shandong Province. On the basis of the achievement of the location results, the conclusions were drawn as follows： On the basis of the achievement of 3D strata fracturing situation and the section plane of microseimic events in different areas, the relationship between spatial structure of overlying strata and mining pressure field was found, and we might describe distribution range of dynamic pressure of advance pressure and lateral stress around long face, and range of structure ad-tivation. Quantitative guidance to prevent dynamic disasters was provided. The practice in coal mine got a effective results. According to the FLAG3D soft numerical simulation of diameter drilling hole （the diameter is 300 mm） to relieve pressure in specified geological condition in Huafeng Mine, the right distance of two dirlls is 2.5 m and the right depth is 12 m. The research pro-vided basic guiding and practical experiences for the underground microseismic monitoring and disaster prevention in side slopes or tunnels engineering.
基金Supported by the National Natural Science Foundation of China (40972197,41002107, 41030750)the Program of Knowledge Innovation of the Chinese Academy of Sciences (KZCX2-YW-Q03-02)
文摘To trace the potential hazards of open-pit slope in Longshou mine,global positioning system（GPS） is applied to monitoring ground movement and deformation induced by transition from open-pit to underground mining.Through long-term monitoring from 2003 to 2008,huge amounts of data were acquired.Monitoring results show that large-scale ground movement and deformation have occurred in mining area,and the movement area is ellipse-shaped.The displacement boundary of settlement trough is 2.0 km long along the exploratory line,and 1.5 km long along the strike of ore body.GPS monitoring results basically agree with the practical deformation state of open-pit slope.It is indicated that the long-term GPS monitoring is an effective way to understand the mechanism of ground movement and deformation in mine area. 更多
基金Supported by Sub-project of National Key Research and Development Program of China in the 13thFive-Year Plan of China(2016YFD0300804-5)Special Fund for Agro-scientific Research in the Public Interest(201503106)Special Construction Project of"Overseas Taishan Scholar"
文摘This study was conducted to investigate the effects of long-term located fertilization on soil phosphorus,the changes of soil available phosphorus(OlsenP),the evolution of soil total phosphorus(TP)and the ratio change of Olsen-P to TP(PAC)by 33-year fertilization experiments in winter wheat-summer maize rotation system in Shandong fluvo-aquic soil.Eight treatments were designed as no fertilization(CK),nitrogen fertilizer(N),nitrogen and phosphate fertilizer(NP),nitrogen and potassium fertilizer(NK),phosphate and potassium fertilizer(PK),nitrogen-phosphate-potassium fertilizer(NPK),reduced NPK fertilizer(N(15)PK),and increased NPK fertilizer(N_(25)PK).Meanwhile,eight organic fertilizer-added treatments were designed based on the application of inorganic fertilizer the same as the above ones.The results showed that TP,Olsen-P and PAC of treatments added with organic fertilizer were higher than those without organic fertilizer,and those of the treatments applied with phosphate fertilizer were higher than those of no phosphate fertilizer.With the increase of years,soil P pool decreased due to crop absorption,nutrient loss and morphological transformation and other causes under the treatments of without and only phosphate fertilizer,while remained stable under the treatments added with organic fertilizer.The PAC values were generally lower in fluvo-aquic soil,and it could be improved by the application of organic fertilizer.On the whole,the application of chemical phosphate fertilizer combined with organic fertilizer could improve the phosphorus content in soil and ensure the supply of phosphorus nutrition.This study would provide scientific basis for fertilization management and soil fertility in fluvo-aquic soil.
文摘A systematic study concerning the effects of a long-term stationary fertilization on content and property of soil humus in fluvo-aquic soil sampled from Malan Farm, Xinji City, Hebei, and arid red soil and paddy red soil sampled from the Institute of Red Soil, Jinxian County, Jiangxi was conducted. The results showed that long-term fertilization had effects not only on the content and composition of soil humus, but also on the physico-chemical property of humus. With applying organic manure or combined application of organic manure and chemical fertilizer, E4 and E6 values of humic acid decreased in fluvo-aquic soil and arid red soil, but increased in paddy red soil. In paddy red soil, E4 and E6 values of humic acid increased also with a single application of chemical fertilizer, but E4 and E6 values had less change of humic acid in fluvo-aquic soil and arid red soil. The effects on the visible spectroscopic property of f ulvic acid were different from that of humic acid. Long-term application of organic manure or combined application of organic manure and chemical fertilizer could increase E4 and E6 values of f ulvic acid in three types of soil. Single application of chemical fertilizer had less effect on the E4 and E6. Long-term fertilization could also influence the ultraviolet spectroscopic property of humus. With a single application of organic manure or combined application of organic manure and chemical fertilizer, the ultraviolet absorbance of humic acid and f ulvic acid increased in the three types of soil. But this effect was obvious only in short wave length, and the effect could decrease if the wave length increased. With a single application of chemical fertilizer the ultraviolet absorbance of fulvic acid could increase, but it of humic acid increased only in fluvo-aquic soil. Long-term application of organic manure or combined application of organic manure and chemical fertilizer could increase the content of total acidic groups, carboxy groups and phenolic hydroxy groups of humic acid and fuvic acid in the three types of soil. Single application of chemical fertilizer had less effect on the content of total acidic groups, carboxy groups and phenolic hydroxy groups of humic acid and fuvic acid in the three types of soil.
文摘This study was aimed at assessing the locational compliance status of petroleum handling facilities in the Niger Delta to the specifications of the Nigerian Upstream Petroleum Regulatory Commission (NUPRC) locational standards. A cross-sectional research design was employed using a Standard checklist of seventeen compliance specifications of NUPRC. A walk-through survey was carried out in 118 identified facilities in three locations: Eket in Akwa-Ibom State, Port Harcourt in Rivers State and Warri in Delta State, respectively. The data obtained were analyzed using inferential and descriptive statistics, Kendall Tau-B and Principal Component analyses. The results indicate that the Petroleum Product facilities complied with 6 (35.29%) out of 17 specifications while 11 (64.70%) specifications were violated. Locations compliance of stations in the Niger Delta region is generally poor, but Stations in Port Harcourt recorded a higher compliance rate compared to Eket and Warri. Results of Kendall’s tau-b and Principal Component analysis indicated positive association between all the land space locational compliances. This study attributes the poor compliance rate to rapid urbanization, overpopulation, proliferation of filling stations and poor monitoring by regulatory agencies. It is recommended that, NUPRC should improve its monitoring and enforce regulatory operational specifications.
