• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A method to predict rockburst using temporal trend test and its application

    2024-03-25 11:06:06YrongXueZhenleiLiDzhoSongXueqiuHeHongleiWngChoZhouJinqingChenAlekseiSobolev

    Yrong Xue,Zhenlei Li,Dzho Song,*,Xueqiu He,b,Honglei Wng,Cho Zhou,Jinqing Chen,Aleksei Sobolev

    a School of Civil and Resources Engineering, University of Science and Technology Beijing, Beijing,100083, China

    b Zhong-an Academy of Safety Engineering, Beijing,100083, China

    c China Energy Group Xinjiang Energy Co., Ltd, Urumqi, 830027, China

    d Khabarovsk Federal Research Center of the Far Eastern Branch of the Russian Academy of Sciences (KhFRC FEB RAS), 51 Turgenev Street, Khabarovsk, 680000,Russia

    Keywords: Rockburst Microseismicity Early warning Mann-Kendall trend test Confusion matrix Multi-indices fusion

    ABSTRACT Rockbursts have become a significant hazard in underground mining,underscoring the need for a robust early warning model to ensure safety management.This study presents a novel approach for rockburst prediction,integrating the Mann-Kendall trend test (MKT) and multi-indices fusion to enable real-time and quantitative assessment of rockburst hazards.The methodology employed in this study involves the development of a comprehensive precursory index library for rockbursts.The MKT is then applied to analyze the real-time trend of each index,with adherence to rockburst characterization laws serving as the warning criterion.By employing a confusion matrix,the warning effectiveness of each index is assessed,enabling index preference determination.Ultimately,the integrated rockburst hazard index Q is derived through data fusion.The results demonstrate that the proposed model achieves a warning effectiveness of 0.563 for Q,surpassing the performance of any individual index.Moreover,the model’s adaptability and scalability are enhanced through periodic updates driven by actual field monitoring data,making it suitable for complex underground working environments.By providing an efficient and accurate basis for decision-making,the proposed model holds great potential for the prevention and control of rockbursts.It offers a valuable tool for enhancing safety measures in underground mining operations.

    1.Introduction

    Rockburst,which is considered one of the most severe dynamic disasters in coal mines,refers to the sudden and forceful ejection of coal and rock masses caused by the accumulation of elastic deformation energy internal.It often leads to casualties and substantial property damage(Cai et al.,2020;Xue et al.,2021;Zhu et al.,2016,2018).With the continuous expansion of coal mining into deeper areas,the coal seam structure of the stope and the surrounding rock storage conditions around the roadway have become increasingly complex.Consequently,the internal dynamic response characteristics of coal and rock masses have become more intricate(He et al.,2020).As a result,the frequency of rockburst hazards has risen sharply,posing significant threats to the safety of individuals working in mines and their properties.

    In an effort to prevent and control rockburst disasters,the microseismic (MS)monitoring system has been widely adopted in underground mines,proving to be a powerful tool for rockburst prediction.Lu et al.(2015)investigated the evolutionary patterns of multi-parameter precursory characteristics before and after rockburst events.Cao et al.(2016) conducted a qualitative analysis of the evolution of microseismicity leading up to a catastrophic rockburst.Their study concluded that abnormal clustering of seismic sources,abnormal variations in daily total energy release,and event counts could be considered precursors to rockburst incidents.Yu et al.(2016) indicated that the daily maximum MS energy could be used to estimate rockburst intensity.Tang et al.(2018)proposed that the spatial and temporal concentration of MS events,along with a decrease in theb-value,can be regarded as precursors to the instability of the surrounding rock.He et al.(2019)conducted a study on the variations in MS and acoustic emission monitoring data prior to rockburst events in steeply inclined and extremely thick coal seams.They established the energy deviation value and total deviation high-value indicators as precursory warnings for rockburst.Li et al.(2021)observed a decreasing trend in MS events prior to rockbursts,accompanied by a “quiet period,” while electromagnetic radiation intensity increased and reached a maximum before rockbursts.Zhang et al.(2021) proposed that peak particle velocity could serve as an early warning index for rockburst,emphasizing the need to consider this indicator alongside postprocessing indicators such as event energy.Additionally,various indicators such as fractal dimension,b-value,apparent cumulative apparent volume,the energy index,cumulative released energy,and theEs/Epvalue of MS events have been analyzed in relation to rockburst events.These analyses provide valuable insights for rockburst prediction (Feng et al.,2016;Mondal and Roy,2019;Xie and Pariseau,1993;Xue et al.,2020).

    However,due to the complexity of rockburst occurrence mechanisms,the use of different warning indicators can lead to varied warning results for the same event.This is because these indicators reflect the evolution of rockburst precursors based on different principles.Consequently,assessing the actual hazard state can pose challenges for the mine personnel.Consequently,there have been notable efforts to develop a multi-indices rockburst early warning model to improve the accuracy of predictions.Cai et al.(2014) developed a multidimensional index system that incorporates MS information.They combined the comprehensive anomaly index method with theR-value method to provide a quantitatively description of the real-time rockburst hazard status.Liu et al.(2016) introduced a methodology for dynamic risk assessment and management of rockbursts in drill and blast tunnels.The approach utilizes quantitative MS indices to evaluate the probability of rockburst occurrence.Cai et al.(2018) proposed a fuzzy comprehensive rockburst risk evaluation model.The model incorporates components such as the Gaussian shape membership function,the confusion matrix and the maximum membership degree principle.Dou et al.(2018) utilized theR-value scoring method to determine the weights and estimate the critical values of various MS indicators,such as bursting strain energy,time-spacemagnitude independent information,and time-space-magnitude compound information.The study demonstrated that their comprehensive index effectively quantifies the pre-warning of rockburst risk.Cao et al.(2020) developed a probabilistic forecasting methodology for rockburst hazard at Coal Mine Velenje.Their approach integrated data-driven techniques with a physicsbased framework,utilizing MS monitoring data in conjunction with concurrent face advance records.The study achieved satisfactory results,highlighting the effectiveness of their methodology in rockburst hazard assessment.In a study by Feng et al.(2015),a microseismicity-based rockburst warning method was proposed.This method utilized real-time microseismic data and a rockburst warning formula.The formula included a rockburst database,selection of typical rockburst cases,functional relationships between microseismicity and rockbursts,optimal weighting coefficients,and dynamic updating.The method was successfully applied to rockburst warning in deep tunnels at the Jinping II hydropower project.Yin et al.(2021) proposed an integrated CNN-Adam-BO algorithm based on microseismic monitoring data for real-time prediction of rockburst intensity.Cao et al.(2022) presented a knowledge and data fusion-driven deep neural network called FDNet for coal burst prediction.The FDNet used the existing mine seismic model to extract explicit features and employed deep learning to automatically extract implicit features from mine microseismic data.This approach provided new insights in coal burst prediction.In addition,decision tree(Wang,2021;Zhao et al.,2021),support vector machine (Ji et al.,2020;Jin et al.,2022),neural network (Feng et al.,2019),and other algorithms have also been used to establish new multi-indices fusion early warning models for rockburst.

    Although significant advancements have been made in rockburst early warning models,the currently available single-indicator or multi-indicator warning methods mostly rely on threshold values to determine whether an alarm is triggered.However,the occurrence of rockburst is a complex nonlinear process with mechanisms that are not fully understood.This lack of understanding has created a challenge in improving the accuracy of existing rockburst warning methods.On-site empirical evidence indicates that there are anomalous patterns preceding rockburst incidents.However,these patterns are often less considered in existing warning methods.Therefore,to enhance the accuracy of warnings,it is crucial to establish a rockburst early warning model that incorporates multidimensional warning indicators and captures real-time changes in these indicators.

    To address the aforementioned limitations,this paper proposes a multi-indices fusion rockburst early warning model based on the Mann-Kendall trend test method (MKT).The model aims to improve the rockburst early warning precision by incorporating multiple indicators and utilizing MKT for trend analysis.The model initiates by applying MKT to capture the temporal changes of multidimensional rockburst precursory indices.It evaluates whether the changing trends conform to the characterization law of rockburst precursors,forming the foundation for early warning.The model continuously updates and selects the indices by utilizing the confusion matrix to reevaluate the early warning effectiveness of each index.This process ensures that the most effective indices are chosen for optimal performance in predicting rockburst events.The effectiveness of early warnings serves as the data fusion weight,leading to a quantified comprehensive rockburst hazard evaluation index.The proposed model has demonstrated successful applications at Kuangou Coal Mine in China,providing a timely and accurate decision-making tool for the prevention and control of underground rockburst incidents.

    2.Principles and methodology

    Pu et al.(2019,2020) have found that rockburst is a complex nonlinear process;that is,there is a degree of uncertainty in the occurrence of precursor anomalies of rockbursts.This uncertainty arises from the incomplete understanding of the mechanism and the difficulty in specifying relevant measurement parameters during the rockburst evolution.Nevertheless,rockburst is fundamentally the result of failure stress in coal and rock masses,leading to the formation of large internal cracks and energy release.Extensive research has identified various precursory indices that reflect the process of fissures initiation,propagation,convergence,and connection into macro-fractures in coal and rock masses.The abnormal change trend of each precursory index before rockburst can be regarded as a danger signal and has achieved good application results in infield practice.However,the complex geological environment and mining conditions in each mine give rise to unique characteristics for the application of individual indices in different regions.To enhance the accuracy of rockburst prediction,this paper proposes a multi-indices fusion rockburst early warning model based on the MKT.

    2.1.Overview of the multi-indices fusion rockburst early warning model

    As illustrated in Fig.1,the proposed model begins by collecting real-time monitoring data from the online rockburst monitoring system and then uses a certain length of the time window and sliding step to calculate the multidimensional rockburst precursory indices,which reflect the rockburst development process information.The model utilizes MKT to determine the trend of each index,ensuring adherence to the characterization law of rockburst precursors for early warning.The effectiveness of each index is evaluated and ranked using the confusion matrix.Selected indices with high effectiveness are combined to form the integrated rockburst hazard indexQ.Qis categorized into four levels of rockburst risk:no risk,weak risk,medium risk,and high risk,corresponding to specific ranges.More details about the development of the early warning model can be found in the following section.

    Fig.1.Multi-indices fusion rockburst early warning model.

    2.2.Details of the multi-indices fusion rockburst early warning model

    2.2.1.Pre-processing of real-time monitoring data

    Currently,online monitoring systems for underground dynamic disasters,such as MS,electromagnetic radiation,and acoustic emission,are widely used.These systems capture real-time information on the energy release of coal-rock mass breakage,seismic source location,roof pressure,and other parameters,providing insights into the evolution of coal-rock dynamic hazards.However,the data obtained from these monitoring systems often show irregular time intervals,making it challenging to perform data mining and analysis.To address this challenge,it is necessary to standardize the data obtained from different sources by structuring them with consistent dimensional granularity and attributes.

    By dividing the original monitoring data into fixed-length time windows and calculating the feature of each sample within the windows,it is possible to transform irregular time series into regular time series and obtain the trends of the features.This process is illustrated in Fig.2.Firstly,a sliding time window of length Δtis defined,and the monitoring data time series is divided intonsets of data,each with a length Δtand corresponding to the end moments the time window.The data set for the moment Tiis denoted as Xi[x1,x2,x3,…,xk](k≤t,0

    Fig.2.Pre-processing of irregular monitoring data (Cai et al.,2018).

    2.2.2.Trend determination of precursory indices

    Extracting meaningful precursory indices from original underground monitoring data is an effective approach for disaster prevention and control in coal-rock dynamic environments.Various precursory indices have been widely utilized,including thebvalue(Gutenberg,1956),total fault areaA(t) (Lu et al.,2015),lack of seismicMm(Utsu et al.,1995),energy deviationDE(He et al.,2019),among others.These indices often exhibit a continuous trend of increase or decrease before the occurrence of disasters like rockbursts,enabling timely hazard level assessment.However,the range of each index varies significantly under complex mining conditions and geological environments.Consequently,the conventional approach of using a single critical value lacks scalability,as it varies with specific mine conditions.This limitation necessitates a more effective time series trend test method for practical application of precursory indices,as relying solely on human experience to identify abnormal trends hampers the utilization efficiency of these indices.

    We performed MKT to determine the monotonic trends in the time series of precursory indices.MKT is rank-based and nonparametric,which means it does not rely on specific data distribution and emphasizes the relative order of magnitude rather than the actual data values.Furthermore,MKT is capable of handling extreme values,making it highly suitable for analyzing underground mine monitoring data.The MKT can be defined as

    wherenpresents the length of dataset Xi[x1,x2,x3,…,xn],xkcorresponding to the rank for thekth observations(k=1,2,3,…,n-1),andxjrepresents the rank for thejth observations(j=k+1,2,3,…,n).

    According to Mann and Kendall (Mann,1945;Kendall,1948),when the value ofnis greater than or equal to 10,the statisticSfollows an approximate normal distribution with the following mean and variance:

    wherenrepresents the length of dataset Xi[x1,x2,x3,…,xn],grepresents the number of these equal trend values or groups,andtprepresents the number of data values in thepth group.

    The Kendall standardized test statisticsZis calculated as

    If the calculated value ofZis greater than zero and passes the significance test (with a significance level of α=0.05 andZ=±1.96),it indicates an increasing trend in the data.Conversely,if the calculated value ofZis less than zero,it suggests a decreasing trend in the data.However,ifZdoes not pass the significance test,it implies that the data does not exhibit an obvious change trend.This approach allows for the statistical analysis of trend behavior in the data,providing valuable insights into the temporal variations of rockburst-related indicators.

    2.2.3.Evaluation of early warning effectiveness and preference of indices

    Numerous studies have focused on rockburst precursory indices,offering valuable insights into the evolution of rockbursts from multiple perspectives such as time,space,and intensity.These indices exhibit different response characterization laws influenced by various contributing factors.Consequently,the evaluation and selection of precursory indices play a crucial role in determining early warning levels during field applications.To enhance decisionmaking efficiency for underground personnel,it is important to choose appropriate indices and assign them appropriate weights.In this study,the effectiveness of each precursory index is evaluated and optimized usingF_scorein Table 1.The matrix includes true positive (TP) and true negative (TN) values,representing correct warnings,and false positive (FP) and false negative (FN) values,representing false and missed warnings.Initially,actual monitoring data and records of high-energy tremors and rockburst events are used to calculate and rank theF_scorefor each index in the precursory database.Indices with highF_scoresare selected as the foundation for the model.The indices are periodically re-evaluated and selected to ensure their efficiency in assisting the model’s early warning decisions.Further details can be found in Fig.1.

    Table 1 Calculation of the confusion matrix and early warning effectiveness (Fawcett,2006).

    2.2.4.Multi-indices fusion and rockburst risk classification

    Integrating various multidimensional coal-rock dynamic disaster early warning indices that consider temporal,spatial,and intensity aspects is crucial for creating a unified and quantitative real-time warning system with consistent criteria and thresholds.By applying the comprehensive anomaly index method (Cai et al.,2014),the integrated rockburst hazard indexQis constructed,and its calculation method is shown in Eq.(6):

    wherenrepresents the total number of preferred precursory indices,Fkrepresents theF_score,which reflects the early warning effectiveness of thek-th indicator.Wk(±)represents the anomaly membership of thek-th index,and they range from 0 to 1.The calculation ofWk(±)can be performed using the following method.

    For the positive precursory index,using MKT to determine its trend in the previous period and get

    For the negative precursory index,the same reasoning yields

    Based on theoretical analysis and numerous field experiments,the coal-rock dynamic hazard level can be categorized into four levels (Dou and He,2007).These hazard levels align with the classifications specified in the“Rules for Prevention and Control of Coal Mine Rockburst,China.” The corresponding hazard levels are presented in Table 2.

    Table 2 Classification standard for coal-rock dynamic hazards.

    3.Case study

    3.1.Engineering description and rockburst contributing factor analysis

    Kuangou Coal Mine (KCM) is situated in Urumqi,Xinjiang Province,China.The mine has six mineable coal seams ranging from shallow to deep: B4-2 (mined out),B4-1 (mined out),B3,B2(under mining),B1,and B0.Currently,mining activities are mainly concentrated in the I010203 working face of the B2 coal seam,and the general layout and stratigraphic structure diagram of KCM are shown in Figs.3 and 4.In Fig.3,the black line represents the working face and roadway in the B2 seam,while the blue line represents the working face and roadway in the B4-1 seam.The I010203 working face has an inclination width of 192 m and a strike length of 1,469 m.It employs the fully mechanized caving mining method,with a coal cutting thickness of 3.2 m and a roof coal release thickness of 6.3 m.The average burial depth is approximately 350 m.The working face is situated on the west side of the safety coal block,with a 15 m section pillar reserved between its north side and the I010201 goaf of the same coal seam.Moreover,the I010403 and I010405 goafs are situated around 70 m above the I010203 working face in the B4-1 coal seam,forming a “knifehandle”boundary.The B1 coal seam,which is yet to be mined,lies approximately 30 m below the B2 coal seam.Within the working face,there is a significant fault labeled F2-1,which has an inclination of 60°and a fault displacement of 2.9-6.8 m.According to the test,the average duration of dynamic fractureDTof B2 coal seam at different locations is 250.33ms,the elastic strain energy indexWETis 3.43,the bursting energy indexKEis 2.79,the uniaxial compressive strengthRCis 26.34 MPa,according to the Methods for Test,Monitoring and Prevention of Rockburst(GB/T 25217.2-2010),B2 coal seam has weak bursting liability.In the same way,the roof of B2 coal seam has strong bursting liability and the floor has weak bursting liability.

    Fig.3.Geological condition of KCM and the layout of MS monitoring system installed until March 7,2018.

    Fig.4.Schematic diagram of stratigraphic structure.

    KCM is equipped with the ARAMIS M/E MS monitoring system.The MS sensors have a sampling frequency of 500 Hz and a sensitivity of 110 Vs/m ± 10%.These sensors are capable of monitoring low-frequency high-energy MS events with an energy threshold of 100 J and a frequency range of 0-150 Hz.The positioning accuracy of the sensors is±20 m in theXandYdirections and±50 m in theZdirection.Based on the geological conditions and characteristics of the surrounding rocks at location I010203,2 MS sensors were arranged in the haulage roadway of the working face with a spacing of 150 m,marked as S2 and T3;1 MS sensor was arranged in the craft lane,marked as T16;2 MS sensors were arranged in the ventilation roadway with a spacing of 150 m,marked as T10 and T11(see Fig.3).To ensure that the working face remains within the detection range of the MS sensors,a specific adjustment is made when the sensors are located less than 50 m away from the working face.In such cases,the sensors are shifted 300 m ahead of the working face.This adjustment ensures that the working face remains within the coverage area of the MS sensors for effective monitoring.(Please refer to Khan et al.,2022 for more details about the MS system).From February 1,2018,to January 31,2019,a total of 15 rockburst events have occurred at the I010203 working face during this period,and the location of each event source is shown in Fig.5,and the information is shown in Table 3,which shows that:

    Table 3 Overview of each rockburst event.

    Fig.5.Source location diagram of each rockburst events.

    (1) Rockburst events occur at varying focal locations,with the majority concentrated within the range of 0-365 m from the leading working face.The density of focal locations is high both before and after the fault location,with a significant number of events occurring in the haulage roadway and some in the surrounding rock near the ventilation roadway.The range of rockburst events typically falls within 100-550 m from the working face.This indicates that the highest rockburst risk is observed within the 100-550 m range from the cutting hole.

    (2) The majority of rockburst sources are located in the coal seam and roof in the middle and lower sections of the working face.They also occur in the roof and section coal pillar in the lateral goaf.These areas are identified as the main elastic energy concentration areas,and the event density of the roof in the vertical direction is significantly higher than that of the floor (11 out of 15 times).The breaking of the hard roof is the main inducing factor of the rockburst of KCM.

    The analysis results show that the hard roof breakage at I010203 is the main inducing factor of rockburst events,and rockbursts are often preceded by high-energy tremors(energy>106J),which are a direct cause of rockbursts as a dynamic load disturbance(Manouchehrian and Cai,2017;Wang et al.,2018).Based on this,the MS data collected by the MS sensors at the I010203 working face are used as the basis for extracting the precursor indices of coal-rock dynamic hazards to predict the occurrence of highenergy tremors.These indices enable the prediction of highenergy tremors and subsequently facilitate the monitoring and warning of rockburst hazards.

    3.2.Rockburst precursory indices

    The effective utilization of MKT for monitoring and early warning of rockbursts relies on the availability of accurate precursory indices that possess clear physical significance.Numerous researchers (Cai et al.,2014,2018;He et al.,2019,2021;Liu et al.,2019;Lu et al.,2015;Ma et al.,2019;Qin et al.,2019;Xu et al.,2017;Yu et al.,2017) have proposed multidimensional precursory indices based on the “temporal-spatial-intensity” framework and these indices can effectively capture the evolving patterns preceding rockburst events,and their practical application in mining has yielded positive results.Therefore,KCM adopts a library of 18 precursory indices to monitor and predict rockburst events,as shown in Fig.6.For detailed descriptions each warning index,please refer to Xue et al.(2023).

    Fig.6.The library of precursory indices.

    From February 1,2018,to January 31,2019,the I010203 working face experienced a series of rockburst events,making it an ideal case study site.During this period,a total of 30065 MS events(Except September 24 to October 19,2018,with no data recorded because of working face closure) and 15 rockburst events (see Table 3).A time window of 15-d and a sliding step of 1-d are employed to calculate of the actual values for each index,and its time-sequence evolution graphs are shown in Fig.7.

    Fig.7.Temporal evolution of MS precursory indices before and after rockburst events.

    Prior to the occurrence of a rockburst,each precursory index shows abnormal changes,characterized by a continuous increase or decrease.The values of these indices fluctuate within a range of high or low values until the rockburst event eventually takes place.For example,A(t)(see Fig.7l)demonstrates a rapid increase in the period leading up to a rockburst event,exhibiting a notable highvalue anomaly,indicating that the internal crack propagation degree of coal and rock mass increases rapidly,the number of large cracks increases,and the fracture of coal and rock mass intensifies,which induces a high-energy mine tremor,and its value falls back to the normal level after energy release;whenSD(see Fig.7i) increases rapidly in time series,signifies a rapid increase in both the frequency and energy level of MS events.This increase indicates a more concentrated spatial distribution of these events and a heightened degree of internal rupture within the coal and rock mass.Consequently,it is accompanied by the occurrence of rockbursts;b(see Fig.7p) decreases rapidly before the occurrence of a rockburst,and then fluctuates in the low value range,indicating that the level of stress concentration within the coal rock mass and the degree of elastic energy accumulation increase,leading to a rapid increase in the proportion of large rupture,and the possibility of a high-energy tremor increases,making it more likely to induce a rockburst.Each index in the rockburst precursory indices library reflects the precursor evolution law of rockburst from the multidimensional “temporal-spatial-intensity” and corresponds to the existence of certain characterization laws,which can monitor and warn the rockburst risk.

    3.3.Evaluation of the early warning effectiveness and selection of indices

    To monitor the rockburst hazard in the mine,the study employed the MKT with a 15-d time window and a 1-d sliding step.This analysis was conducted to evaluate the change trend of each precursory index,considering the rockburst characterization law outlined in Section 3.2 and the associated change characteristics.For ease of reference,the 15 rockburst events mentioned in Table 3 were assigned numbers from 1 to 15 in chronological order.The results of assessing the trend of the precursory indices within a 15-d period preceding each rockburst event are presented in Table 4,where “I" indicates that the index has an increasing trend,“D" indicates that the index has a decreasing trend,and"-"indicates that the index does not pass the significance test and without a significant change trend.

    Table 4 Trend determination results of each precursory index before rockbursts and its early warning effectiveness.

    The results depicted in Fig.7 and Table 4 demonstrate that the change trends of the precursory indices,as determined by MKT,largely align with the qualitative judgment results and adhere to the characterization law derived for the precursory indices.This indicates that the MKT method is effective in monitoring and providing warnings for rockburst hazards.However,it should be noted that a single precursory index has limitations in capturing the complete evolutionary patterns of rockburst precursors.Consequently,relying solely on a single index may not always yield accurate early warnings for every rockburst event.For example,Qtexhibits a declining trend prior to events 1-6,10,and 12-15,and it has good early warning effectiveness as a negative rockburst precursory index accurately warning 11 out of 15 rockburst events.Furthermore,the different dimensions of each early warning index may contribute to vary early warning results for the same event.For example,prior to the onset of rockburst 8,DEin the intensity dimension andDFin the temporal dimension show an increasing trend as positive indices warn the danger,while λ in the spatial dimension also shows an increasing trend as a negative index and does not successfully warn of the occurrence of rockburst 8.Therefore,it is necessary to fuse the indices in the precursory library,the advantages of different dimensions of indices complement each other,which can be more efficient and accurate monitoring and warning of rockburst hazards.

    The precursory index is deemed accurate if it exhibits an abnormal trend within 5 d before the occurrence of a high-energy tremor (with energy greater than 106J).Specifically,positive precursory indices should demonstrate an increasing trend,while negative precursory indices should show a decreasing trend.Any deviation from these criteria would be classified as a missed or false alarm.The choice of a 5-d period for early warning indices is in accord with the practical experience of KCM,and it can be adjusted for other mines based on their specific requirements.To evaluate the early warning effectiveness of each index,the confusion matrix(refer to Table 1) is utilized.From the confusion matrix,theRecall,precision,andF_scoreof each index can be calculated.Table 4 presents the results of the early warning effectiveness for each index.

    Based on Table 4,the precursory indices can be ranked in terms of their early warning effectiveness as follows:Emax>DF>Eavg>Esum>DE>Fsum>AC=P(b)>F>b>λ>Qt>ΔF>ΔE>A(t)>Mm=A(b)>SD.Notably,the indices in the intensity dimension generally exhibit higher early warning effectiveness compared to those in the temporal and spatial dimensions.These findings suggest that the precursory indices related to intensity in rockburst events carry a higher level of information entropy compared to the temporal and spatial indices.The energy characteristics of MS play a crucial role in capturing the evolutionary process of rockburst precursors.These intensity-based indices demonstrate a high level of effectiveness in reflecting the underlying dynamics of rockburst events.As a result,they hold significant importance and should receive increased attention in practical applications for improved rockburst prediction and early warning systems.

    To maximize the effectiveness of early warning,the model selected preferred indices from each dimension.Out of the precursory indices library,eight indices (Emax,DF,Eavg,Esum,DE,Fsum,AC,P(b))were chosen based on theirF_scoresexceeding 0.2.These selected indices accounted for 74.9% of the total warning effectiveness and covered a significant portion of the intensity dimensions of rockburst hazard precursory indices.This selection ensures a strong foundation for the overall warning capability of the model.Regular re-evaluation and optimization of the early warning performance of each index in the precursory indices library are important to maintain the scalability and robustness of the model.Ideally,this re-evaluation should occur at least once every month.By dynamically updating the model based on actual field data,it can effectively assist relevant personnel in making rockburst warning decisions.This approach ensures that the model remains adaptable to changing conditions and maintains its efficiency in real-world applications.

    3.4.Multi-indices fusion early warning of rockburst

    3.4.1.Application effect test

    TheF_scoreof each rockburst precursory index is utilized to calculate the weight assigned to each index.This approach ensures that indices with higher warning effectiveness receive higher weights,while indices with lower warning effectiveness are given lower weights.By combining these weights with Eq.(6),Qis calculated.The time-series ofQis depicted in Fig.8.Observing the curve,it is evident that the majority ofQvalues exceed 0.5 within a 5-d period prior to the occurrence of high-energy tremors.This corresponds to the medium or strong rockburst hazard category described in Table 2.The results indicate a strong correlation between the index and the occurrence of rockburst.Moreover,it demonstrates the capability of the index to provide effective monitoring and early warning for rockburst events.

    Fig.8.Warning results of Q.

    TheQdemonstrated strong performance in terms ofF_score,Recall,andPrecision,with values of 0.563,0.500,and 0.643,respectively.These metrics indicate thatQoutperformed individual warning indices,exhibiting higherF_scoreandPrecision.This suggests thatQhas improved warning effectiveness,as depicted in Fig.9.However,theRecallofQis lower than the 5 single precursory indicesDF(Recall=0.667),P(b)(Recall=0.611),AC(Recall=0.611),DE(Recall=0.611) andFsum(Recall=0.556).This means that the accuracy ofQ’s early warning for the actual occurrence of highenergy tremors in the original monitoring data is lower than these five single indices,which are more conservative in the early warning process compared to the single indices,and will only warn when the predicted outcome is a high probability of rockburst events,which is shown asPrecisionis higher than other single indices.

    Fig.9.Comparing the effectiveness of early warning indices.

    In general,Qoffers several advantages compared to single precursory indices.By combining multiple indices from different dimensions,Qincorporates the strengths of each index and establishes warning criteria among the indices.This approach eliminates the influence of singular values and improves the accuracy of the warning,resulting in higher integrated warning effectiveness compared to single indices.Furthermore,Qserves as a unified quantitative evaluation index for rockburst hazard,overcoming the issue of conflicting early warning results from multiple sources.It provides clear quantitative classification criteria,as shown in Table 2,which significantly reduces the management cost associated with underground rockburst prevention and control.Moreover,Qis data-driven,as its base indices are regularly updated.This ensures its adaptability and stability under complex working conditions.By assisting personnel in making efficient and accurate disaster prevention and control decisions,Qhelps avoid major casualties and property damage.

    3.4.2.Comparison with previous models

    Previously,rockburst early warnings at KCM relied on monitoring high-energy MS events,specifically those with energy exceeding 105J.According to the criterion,if an MS event at the working face had energy surpassing 105J,it indicated a dangerous state,whereas energy below this threshold indicated a safe condition.To assess rockburst warning using this approach,the daily maximum MS energy data from February 1,2018,to January 31,2019,was utilized,and the results are shown in Fig.10.

    Fig.10.Results of previous early warning method in KCM.

    The effectiveness of the previous MS warning method was evaluated based on issuing a hazard warning within the same day.If a warning was issued on the same day as the hazard occurred,it was considered as a correct report;otherwise,it was classified as a missed or false alarm.The early warning results of the previous method are presented in Table 5,indicating that the method had a high missing rate.In 13 out of 18 d when actual hazards occurred,no warning was issued.This high missing rate could lead to severe consequences as the necessary measures were not taken in a timely manner.After calculating the metrics,it was found that the previous model had a significantly lower early warning effectiveness,withPrecisionof 0.128,Recallof 0.278 andF_scoreof 0.175 compared to the multi-indices fusion rockburst early warning model proposed in this paper.

    Table 5 Early warning results of the previous model.

    The low effectiveness of the previous model can be attributed to several factors.Firstly,the previous model solely relied on the daily maximum energy (Emax) index in the intensity dimension of rockburst hazard precursory indices.This approach neglected the temporal and spatial dimensions,resulting in the loss of important information regarding the evolution of rockburst precursors.Consequently,the overall warning effectiveness was compromised.Additionally,the previous model utilized a single critical value for issuing warnings,which limited its ability to fully leverage the potential of the precursory indices.In contrast,the new model takes into account the change trend information of theEmaxindex,resulting in an early warning effectiveness of 0.529,making it the most effective single warning indicator.In comparison,the previous model generated several missed alarms.Consequently,the new model achieved an early warning effectiveness 3.2 times higher than that of the previous model.

    In summary,the new model effectively enhances the accuracy of rockburst prediction by utilizing multiple dimensions of precursory indices and leveraging their full potential.

    4.Discussion

    Given the unclear mechanism of rockburst hazards in mines operating under complex working conditions,accurately and efficiently predicting and controlling these hazards remains challenging.While significant progress has been made in monitoring and early warning systems,there is still a need for a highly scalable and accurate rockburst early warning model.To address this,the present study proposes an MKT-based multi-indices fusion early warning model for rockburst hazards (refer to Fig.1).This model leverages multidimensional warning indicators to capture the precursor trends that reflect the evolution of rockbursts.It incorporates the use of a confusion matrix to evaluate the effectiveness of individual warning indices and periodically optimizes the indices based on this evaluation.The integrated rockburst hazard indexQ(refer to Fig.8)demonstrates that it can reach a medium-or high-hazard warning level within 5 d prior to the occurrence of most high-energy tremors.The warning effectiveness of the model reaches 0.563,representing a significant improvement compared to relying on a single warning index(refer to Fig.9).The model has the capability to dynamically adjust the selection of indices based on real-time monitoring data in the field.This adaptability allows it to effectively respond to changing working conditions and ensures the production of stable warning results.The practicality of the model is strengthened by providing reliable warning outcomes for the prevention and control of rockbursts in mining operations.

    In contrast to previous studies that primarily relied on exceeding hazard thresholds or qualitative human judgment to trigger warnings,this study employs mathematical-statistical methods to capture the changing trends of multidimensional warning indicators in real-time.This approach enables a comprehensive representation of the overall evolution process,including fissure initiation,propagation,convergence,and the formation of macro-fractures in coal and rock masses that ultimately lead to rockburst occurrences.Warnings are only issued when the indicator’s change trend clearly aligns with the rockburst precursor response.Furthermore,the study incorporates the use of a confusion matrix to comprehensively evaluate the effectiveness of each individual precursory index.By considering the warning effectiveness of each index,the model can select the most appropriate indices that align with the actual conditions of the mine.This selection process enhances the overall warning performance of the model and ensures that low-effectiveness indices,which may not adequately represent the mine’s specific conditions,do not negatively impact the warning results.The periodic re-evaluation and optimization of indices based on data-driven self-feedback further maintains the model’s effectiveness and sustainability within the complex underground environment.Moreover,this approach significantly reduces human operation and management costs.In practical applications,it is important to note that for mines with newly installed online monitoring systems,equal weights can be initially assigned to the warning indices.The remaining calculations should remain unchanged until sufficient data becomes available.

    In mining practice,rockburst accidents can lead to significant losses.Therefore,it is crucial for an early warning model to minimize both missed alarms (where the model fails to predict a rockburst event that occurs) and false alarms (where the model issues a warning but no rockburst occurs).The costs associated with these two scenarios can vary significantly for the mine.The evaluation of the warning model’s effectiveness relies on the use of the confusion matrix,specifically focusing onPrecisionandRecall.Precisionmeasures the accuracy of the model’s warnings relative to the total number of warnings issued,whileRecallmeasures the accuracy of predicting actual rockburst occurrences relative to the total number of actual rockbursts.TheF_scoreis the harmonic mean ofPrecisionandRecall,providing a comprehensive assessment of the model’s performance.Maximizing bothPrecisionandRecallis crucial for achieving optimal performance.However,achieving high values for both measures simultaneously can be challenging in practice.Mines often prioritize models with higherRecallto minimize the risk of significant losses.To address this,different warning thresholds can be set in the model to determine the changes in indices.The results of applying different warning thresholds are presented in Fig.11.These thresholds allow for finetuning the model’s sensitivity and strike a balance between missed forecasts and false alarms,taking into account the specific requirements and risk tolerance of the mine.

    Fig.11.Comparing the effectiveness of warnings under different thresholds.

    It can be seen that as the warning threshold increases,the model’sPrecisionincreases whileRecalldecreases.The highestF_scoreof 0.643 is achieved at a warning threshold value of 0.6.However,it is important to note that at this threshold,theRecalldecreases to 0.5,indicating that the model issues warnings only when it is more confident in its predictions,resulting in an increased number of missed alarms.On the other hand,theF_scoreof 0.513 is obtained at a warning threshold value of 0.4,which is higher than the threshold of 0.5.However,theRecallat this threshold drops significantly to 0.476 compared to theRecallof 0.643 at the threshold of 0.5.Considering a comprehensive evaluation,a warning threshold value of 0.5 forQappears to strike a balance betweenPrecisionandRecall,achieving a reasonableF_scoreof 0.643 while maintaining aRecallof 0.5.This threshold value aligns with the concept of“producer equilibrium”for mining enterprises,taking into account the trade-off between missed alarms and false alarms.Nevertheless,in practical applications,the threshold value can be adjusted based on the mine’s own risk tolerance level to best align with its specific interests and requirements.

    In the next step,we plan to use a “cost-sensitive” approach to further improve the rockburst early warning model,i.e.to apply a larger penalty to the model when it generates“missed”alarms and a smaller penalty when it generates “false” alarms so that the model can be trained to better suit the actual warning needs in the field and contribute to creating a safer environment in the mining industry.

    5.Conclusions

    Rockburst frequently result in significant casualties and property damage.Therefore,accurate monitoring and early warning systems are crucial for disaster prevention.This paper proposes a multi-indices fusion rockburst early warning model utilizing the MKT to enhance the quantitative and precise monitoring and warning of rockburst disasters.The main findings of this study are as follows:

    (1) A rockburst early warning model based on MKT is developed,incorporating 18 rockburst precursory indices with clear physical significance.The model assesses the conformity of their temporal trends with the rockburst precursor characterization law to determine warnings.The effectiveness of each index is evaluated and ranked using the confusion matrix,and the indices with higher effectiveness are given more weight.The multi-indices fusion approach is employed using a comprehensive anomaly index method,resulting in the integrated rockburst hazard indexQ.The values ofQcorrespond to four levels of rockburst hazard: none,weak,medium,and high.

    (2) The field application results demonstrate thatQ,serving as a unified quantitative evaluation index,combines the strengths of individual early warning indicators.It achieves an early warning effectiveness of 0.563,surpassing that of any single early warning index.Moreover,the indices from different dimensions provide complementary early warning criteria,mitigating conflicts in the results from multiple sources.The model outperforms the previous mine’s early warning model by 3.2 times in effectiveness,effectively supporting mine personnel in making accurate and efficient disaster prevention and control decisions.

    (3) The model has the capability to undergo periodic updates through self-feedback using real-time monitoring data from the field.This allows for the selection of indices that are most suitable for the complex and dynamic underground working environment,resulting in higher warning effectiveness.As a result,the model reduces the costs associated with human operation and management.Additionally,the model can provide stable early warning results based on online realtime monitoring data and exhibits strong scalability,making it easily adaptable to other mine sites.Overall,this model represents an effective and innovative approach for monitoring and early warning of rockbursts.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (Grant Nos.52011530037 and 51904019)and the Fundamental Research Funds for the Central Universities and the Youth Teacher International Exchange&Growth Program(Grant No.QNXM20210004).We also greatly appreciate the assistance provided by Kuangou coal mine,China Energy Group Xinjiang Energy Co.,Ltd.

    免费一级毛片在线播放高清视频 | 欧美久久黑人一区二区| 麻豆国产av国片精品| 久久精品成人免费网站| 久久精品成人免费网站| 男人操女人黄网站| 一进一出抽搐动态| 国产麻豆69| 一进一出抽搐动态| 国产精品亚洲av一区麻豆| 欧美另类亚洲清纯唯美| 亚洲国产精品久久男人天堂| 久久国产精品人妻蜜桃| 夜夜躁狠狠躁天天躁| 麻豆一二三区av精品| 成人国产综合亚洲| 日韩欧美一区视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 非洲黑人性xxxx精品又粗又长| 久久久久精品国产欧美久久久| 美女国产高潮福利片在线看| 色综合亚洲欧美另类图片| 18禁美女被吸乳视频| 亚洲人成电影免费在线| 亚洲国产欧美网| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区激情短视频| 男人操女人黄网站| 亚洲av成人一区二区三| 99国产极品粉嫩在线观看| 国产精品九九99| 久久久久九九精品影院| 国产精品电影一区二区三区| 波多野结衣高清无吗| tocl精华| 久久人妻福利社区极品人妻图片| 亚洲狠狠婷婷综合久久图片| 久久人妻熟女aⅴ| 青草久久国产| 欧美最黄视频在线播放免费| 神马国产精品三级电影在线观看 | 亚洲精品一卡2卡三卡4卡5卡| 国产午夜福利久久久久久| 午夜日韩欧美国产| 国产在线观看jvid| 中文字幕久久专区| 国产精品99久久99久久久不卡| 久久国产乱子伦精品免费另类| 国产三级在线视频| 18禁国产床啪视频网站| 亚洲中文av在线| 国产精品久久久久久人妻精品电影| 久久中文字幕一级| 中国美女看黄片| 成人国语在线视频| 国产在线精品亚洲第一网站| 色老头精品视频在线观看| 精品福利观看| 美女免费视频网站| 丰满的人妻完整版| 久久久久久国产a免费观看| 免费搜索国产男女视频| 99在线视频只有这里精品首页| 老司机福利观看| 久久久久久久久免费视频了| 国产乱人伦免费视频| 国产三级在线视频| 国产三级黄色录像| 久久久久久久精品吃奶| 久久亚洲真实| 午夜精品在线福利| 性欧美人与动物交配| 亚洲 欧美一区二区三区| 国产视频一区二区在线看| 国产不卡一卡二| 亚洲精品粉嫩美女一区| 久久久久九九精品影院| 首页视频小说图片口味搜索| 精品久久久久久久久久免费视频| 一二三四在线观看免费中文在| 一区二区三区激情视频| 黑人操中国人逼视频| 久久婷婷成人综合色麻豆| 成人三级黄色视频| 日本精品一区二区三区蜜桃| 久久久国产精品麻豆| 满18在线观看网站| 国产色视频综合| 亚洲精品国产一区二区精华液| 欧美乱码精品一区二区三区| 久久国产精品影院| 亚洲 国产 在线| 午夜a级毛片| 黑丝袜美女国产一区| 国产精品久久久久久人妻精品电影| 99精品久久久久人妻精品| 午夜久久久在线观看| 99久久99久久久精品蜜桃| 女人爽到高潮嗷嗷叫在线视频| 精品国产超薄肉色丝袜足j| 国产欧美日韩一区二区三| 国产伦人伦偷精品视频| 免费不卡黄色视频| 在线观看舔阴道视频| 午夜福利高清视频| 亚洲成人国产一区在线观看| 亚洲九九香蕉| 后天国语完整版免费观看| 精品午夜福利视频在线观看一区| 黄色毛片三级朝国网站| 麻豆久久精品国产亚洲av| 在线永久观看黄色视频| 满18在线观看网站| 十分钟在线观看高清视频www| 一级毛片女人18水好多| 国产视频一区二区在线看| 国产野战对白在线观看| 亚洲成av片中文字幕在线观看| 精品免费久久久久久久清纯| 午夜福利在线观看吧| 91av网站免费观看| 无人区码免费观看不卡| 精品一区二区三区av网在线观看| 岛国在线观看网站| 欧美+亚洲+日韩+国产| 亚洲男人天堂网一区| 琪琪午夜伦伦电影理论片6080| 黄片小视频在线播放| 日韩欧美一区二区三区在线观看| 国产av一区二区精品久久| 日本免费a在线| 69精品国产乱码久久久| 淫妇啪啪啪对白视频| 精品福利观看| 精品一品国产午夜福利视频| 午夜福利免费观看在线| 欧美乱妇无乱码| 校园春色视频在线观看| 免费在线观看视频国产中文字幕亚洲| 变态另类成人亚洲欧美熟女 | 美女高潮到喷水免费观看| 啪啪无遮挡十八禁网站| 国产精品美女特级片免费视频播放器 | 美女免费视频网站| 亚洲自拍偷在线| 黄色视频不卡| 少妇裸体淫交视频免费看高清 | 亚洲最大成人中文| 天天一区二区日本电影三级 | 久久久久久久久久久久大奶| 后天国语完整版免费观看| 欧美人与性动交α欧美精品济南到| 久久久国产成人免费| 国产欧美日韩一区二区精品| tocl精华| 成人18禁在线播放| 无限看片的www在线观看| 色综合婷婷激情| 亚洲欧美激情在线| 妹子高潮喷水视频| 巨乳人妻的诱惑在线观看| 免费高清在线观看日韩| 国产欧美日韩综合在线一区二区| 国产日韩一区二区三区精品不卡| 99国产精品一区二区三区| 亚洲精品国产区一区二| 美女免费视频网站| 一二三四社区在线视频社区8| 欧美日韩乱码在线| 亚洲一卡2卡3卡4卡5卡精品中文| 一区二区三区精品91| 日本撒尿小便嘘嘘汇集6| 欧美大码av| 热re99久久国产66热| 中文亚洲av片在线观看爽| 麻豆久久精品国产亚洲av| 一级a爱视频在线免费观看| 日韩精品青青久久久久久| avwww免费| av欧美777| 午夜福利免费观看在线| 欧美日韩福利视频一区二区| 午夜福利成人在线免费观看| 日本在线视频免费播放| 两个人看的免费小视频| 久久青草综合色| 91在线观看av| 一区二区三区国产精品乱码| 91麻豆av在线| 日韩一卡2卡3卡4卡2021年| bbb黄色大片| tocl精华| 成人手机av| 好男人在线观看高清免费视频 | 亚洲成av人片免费观看| 女人精品久久久久毛片| 免费不卡黄色视频| 午夜福利成人在线免费观看| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲精品av在线| 天天一区二区日本电影三级 | 久久久精品欧美日韩精品| 叶爱在线成人免费视频播放| 满18在线观看网站| 麻豆av在线久日| 19禁男女啪啪无遮挡网站| 日韩欧美国产一区二区入口| 九色亚洲精品在线播放| 一本大道久久a久久精品| 亚洲电影在线观看av| 天天躁狠狠躁夜夜躁狠狠躁| 18禁国产床啪视频网站| 亚洲人成电影观看| av超薄肉色丝袜交足视频| 变态另类丝袜制服| 国产三级在线视频| 一边摸一边抽搐一进一出视频| svipshipincom国产片| 窝窝影院91人妻| 国产片内射在线| 一区在线观看完整版| 欧美在线一区亚洲| 色播亚洲综合网| 校园春色视频在线观看| 麻豆成人av在线观看| 91麻豆精品激情在线观看国产| 十八禁人妻一区二区| 看免费av毛片| 韩国精品一区二区三区| 一夜夜www| 国产亚洲欧美精品永久| 97人妻天天添夜夜摸| 黄色女人牲交| 大型av网站在线播放| 欧美一级毛片孕妇| 亚洲欧美日韩另类电影网站| 在线观看免费视频网站a站| av在线天堂中文字幕| 亚洲人成电影免费在线| 亚洲欧美精品综合久久99| 成人18禁高潮啪啪吃奶动态图| 制服诱惑二区| 中国美女看黄片| 久久久久久久久久久久大奶| 久久中文字幕人妻熟女| 欧美老熟妇乱子伦牲交| 国产av精品麻豆| 精品国产乱码久久久久久男人| 国产一区二区在线av高清观看| 欧美成狂野欧美在线观看| 最新美女视频免费是黄的| 中出人妻视频一区二区| 天堂动漫精品| 免费看美女性在线毛片视频| 两性夫妻黄色片| 久久香蕉激情| 久久青草综合色| 亚洲五月色婷婷综合| 免费高清视频大片| 国产区一区二久久| 99国产精品99久久久久| 国产精品99久久99久久久不卡| ponron亚洲| 精品国产美女av久久久久小说| 黄色女人牲交| 国产男靠女视频免费网站| 精品国产一区二区三区四区第35| 国产亚洲精品第一综合不卡| 亚洲国产精品999在线| 极品教师在线免费播放| 国产单亲对白刺激| 在线十欧美十亚洲十日本专区| 国产高清激情床上av| 高清毛片免费观看视频网站| av天堂久久9| 搡老熟女国产l中国老女人| 涩涩av久久男人的天堂| 丝袜美腿诱惑在线| 少妇的丰满在线观看| 欧美激情久久久久久爽电影 | 老司机午夜福利在线观看视频| 日本在线视频免费播放| 欧美成人免费av一区二区三区| 一级毛片高清免费大全| 久久人人97超碰香蕉20202| 亚洲熟妇中文字幕五十中出| 久久香蕉激情| 黄色a级毛片大全视频| 老司机靠b影院| 琪琪午夜伦伦电影理论片6080| 在线观看日韩欧美| 一级毛片精品| 18禁观看日本| АⅤ资源中文在线天堂| 亚洲一卡2卡3卡4卡5卡精品中文| 91大片在线观看| 亚洲中文字幕日韩| 国产精品一区二区三区四区久久 | 亚洲精品在线观看二区| 变态另类丝袜制服| 欧美人与性动交α欧美精品济南到| 日日夜夜操网爽| 亚洲专区中文字幕在线| videosex国产| 人人妻人人澡人人看| 成人精品一区二区免费| ponron亚洲| 嫩草影院精品99| 三级毛片av免费| 国产亚洲精品综合一区在线观看 | 久99久视频精品免费| 啦啦啦观看免费观看视频高清 | av超薄肉色丝袜交足视频| 成人18禁在线播放| 少妇 在线观看| 国产精品,欧美在线| 久久性视频一级片| 国内精品久久久久精免费| 老熟妇仑乱视频hdxx| 黄色视频不卡| 国内精品久久久久久久电影| 国产成+人综合+亚洲专区| 国产精品久久久av美女十八| 免费女性裸体啪啪无遮挡网站| 中出人妻视频一区二区| 手机成人av网站| 精品一品国产午夜福利视频| 高清黄色对白视频在线免费看| 亚洲avbb在线观看| 日本黄色视频三级网站网址| 国产一区二区三区综合在线观看| 日韩有码中文字幕| 一区二区三区国产精品乱码| 久久久国产欧美日韩av| 日韩av在线大香蕉| 国产激情久久老熟女| 久热爱精品视频在线9| 亚洲成av片中文字幕在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 看黄色毛片网站| 国产在线观看jvid| 99riav亚洲国产免费| 韩国av一区二区三区四区| av视频免费观看在线观看| 黄片大片在线免费观看| 免费在线观看日本一区| 国产精品影院久久| 满18在线观看网站| 午夜久久久久精精品| 嫩草影视91久久| 久久久久久久久久久久大奶| 窝窝影院91人妻| 亚洲自拍偷在线| 视频在线观看一区二区三区| 中文字幕最新亚洲高清| 999久久久精品免费观看国产| 亚洲av电影在线进入| 免费久久久久久久精品成人欧美视频| 精品国产亚洲在线| 露出奶头的视频| 日本三级黄在线观看| 老熟妇仑乱视频hdxx| 搡老妇女老女人老熟妇| 久久影院123| av天堂在线播放| 中亚洲国语对白在线视频| 可以在线观看毛片的网站| 日本欧美视频一区| 久久九九热精品免费| 成人永久免费在线观看视频| 亚洲成人国产一区在线观看| 国产精品亚洲一级av第二区| 啦啦啦免费观看视频1| 亚洲色图 男人天堂 中文字幕| 亚洲国产日韩欧美精品在线观看 | 曰老女人黄片| 国产精品二区激情视频| 99久久精品国产亚洲精品| 波多野结衣巨乳人妻| 身体一侧抽搐| 国内毛片毛片毛片毛片毛片| svipshipincom国产片| 男人舔女人的私密视频| 久久天堂一区二区三区四区| 中亚洲国语对白在线视频| 久久 成人 亚洲| 国语自产精品视频在线第100页| 国产高清激情床上av| 国产麻豆成人av免费视频| 曰老女人黄片| 亚洲精品粉嫩美女一区| 国产激情久久老熟女| 成人av一区二区三区在线看| cao死你这个sao货| 日韩大尺度精品在线看网址 | 国产精品亚洲av一区麻豆| 免费无遮挡裸体视频| 99热只有精品国产| 最新在线观看一区二区三区| 亚洲精品中文字幕在线视频| 一夜夜www| 亚洲色图综合在线观看| 国产又色又爽无遮挡免费看| 99在线视频只有这里精品首页| 曰老女人黄片| av有码第一页| 国产av一区在线观看免费| 人人妻人人澡欧美一区二区 | 久久国产精品影院| 午夜福利免费观看在线| 亚洲第一电影网av| 男女午夜视频在线观看| 亚洲色图 男人天堂 中文字幕| 欧美乱色亚洲激情| 午夜免费观看网址| 亚洲一卡2卡3卡4卡5卡精品中文| 91九色精品人成在线观看| 深夜精品福利| 色av中文字幕| 日韩精品青青久久久久久| 中国美女看黄片| 国产亚洲精品一区二区www| 岛国视频午夜一区免费看| 大码成人一级视频| 每晚都被弄得嗷嗷叫到高潮| 久久亚洲真实| 91九色精品人成在线观看| 天堂动漫精品| √禁漫天堂资源中文www| 涩涩av久久男人的天堂| 真人一进一出gif抽搐免费| 岛国在线观看网站| 国产真人三级小视频在线观看| 国产成人av教育| 日本黄色视频三级网站网址| 搡老妇女老女人老熟妇| 精品久久久久久久人妻蜜臀av | 久久精品国产亚洲av高清一级| 中文字幕高清在线视频| 欧美激情 高清一区二区三区| 精品日产1卡2卡| 欧美日本亚洲视频在线播放| 大香蕉久久成人网| 欧美不卡视频在线免费观看 | 侵犯人妻中文字幕一二三四区| 亚洲av电影在线进入| 黄色a级毛片大全视频| avwww免费| av片东京热男人的天堂| 亚洲一区二区三区色噜噜| 日本 欧美在线| 9191精品国产免费久久| 天堂√8在线中文| 欧洲精品卡2卡3卡4卡5卡区| 国产xxxxx性猛交| 日本撒尿小便嘘嘘汇集6| 欧美一区二区精品小视频在线| 又大又爽又粗| 老司机深夜福利视频在线观看| 黄色视频不卡| www.精华液| av中文乱码字幕在线| 女警被强在线播放| 国产亚洲精品久久久久久毛片| 1024香蕉在线观看| 色在线成人网| 免费在线观看黄色视频的| 曰老女人黄片| 日韩欧美免费精品| 我的亚洲天堂| 国产精品久久久久久精品电影 | 女人被狂操c到高潮| 久久人人97超碰香蕉20202| 亚洲成av人片免费观看| 久久国产精品影院| 啪啪无遮挡十八禁网站| 老司机深夜福利视频在线观看| 国产精品日韩av在线免费观看 | 少妇被粗大的猛进出69影院| 久久精品国产99精品国产亚洲性色 | 美女国产高潮福利片在线看| 黄色视频,在线免费观看| 国产1区2区3区精品| 国产精品一区二区三区四区久久 | 久99久视频精品免费| 亚洲精品中文字幕在线视频| 亚洲国产精品成人综合色| 老司机深夜福利视频在线观看| 天天添夜夜摸| 99久久99久久久精品蜜桃| 琪琪午夜伦伦电影理论片6080| 我的亚洲天堂| 国产欧美日韩一区二区三| 法律面前人人平等表现在哪些方面| 中文字幕人妻熟女乱码| 久久精品91蜜桃| 大香蕉久久成人网| 久久香蕉精品热| 黄频高清免费视频| 黑人操中国人逼视频| or卡值多少钱| 一级毛片高清免费大全| 免费高清视频大片| 人妻久久中文字幕网| 欧美av亚洲av综合av国产av| 午夜老司机福利片| 日本免费一区二区三区高清不卡 | 在线观看舔阴道视频| 日本欧美视频一区| 热99re8久久精品国产| 欧美日韩黄片免| 久久精品91无色码中文字幕| 国产熟女xx| 精品人妻在线不人妻| 日本一区二区免费在线视频| 少妇粗大呻吟视频| 18禁美女被吸乳视频| 久久国产精品人妻蜜桃| 欧美老熟妇乱子伦牲交| 国产人伦9x9x在线观看| 国产一区在线观看成人免费| 久久久久久亚洲精品国产蜜桃av| 亚洲精品一区av在线观看| 国内精品久久久久久久电影| 欧美+亚洲+日韩+国产| 18禁黄网站禁片午夜丰满| 18禁观看日本| 久9热在线精品视频| 悠悠久久av| 男人操女人黄网站| 欧美亚洲日本最大视频资源| 久久午夜亚洲精品久久| 国产一区二区三区视频了| 麻豆国产av国片精品| 精品电影一区二区在线| 在线观看日韩欧美| 午夜视频精品福利| 免费在线观看黄色视频的| 精品电影一区二区在线| 无人区码免费观看不卡| 嫩草影视91久久| 亚洲欧美激情综合另类| 91成人精品电影| 亚洲第一欧美日韩一区二区三区| 波多野结衣av一区二区av| 国产单亲对白刺激| 免费少妇av软件| 熟女少妇亚洲综合色aaa.| 精品一品国产午夜福利视频| 多毛熟女@视频| 久久天堂一区二区三区四区| 欧美久久黑人一区二区| 亚洲成av片中文字幕在线观看| 久久人妻av系列| 少妇被粗大的猛进出69影院| 在线播放国产精品三级| 男女做爰动态图高潮gif福利片 | 老熟妇乱子伦视频在线观看| 大型av网站在线播放| 手机成人av网站| 丝袜人妻中文字幕| 亚洲精品国产区一区二| 女警被强在线播放| 午夜免费成人在线视频| 中文字幕色久视频| 精品国产一区二区三区四区第35| 多毛熟女@视频| 久久中文看片网| 日日爽夜夜爽网站| 90打野战视频偷拍视频| 美国免费a级毛片| 满18在线观看网站| 少妇 在线观看| 国产色视频综合| 黄色 视频免费看| bbb黄色大片| 男女午夜视频在线观看| 老司机午夜福利在线观看视频| 男女下面进入的视频免费午夜 | 国产av又大| 成人国产一区最新在线观看| 国产亚洲欧美98| 精品久久久久久,| 午夜福利成人在线免费观看| 婷婷六月久久综合丁香| 久久国产精品人妻蜜桃| 国产一区二区在线av高清观看| 色播在线永久视频| 搡老妇女老女人老熟妇| 亚洲专区国产一区二区| 在线观看免费日韩欧美大片| 少妇粗大呻吟视频| e午夜精品久久久久久久| 波多野结衣av一区二区av| 国产激情久久老熟女| 岛国视频午夜一区免费看| 一级毛片精品| 黄色视频,在线免费观看| 久久精品亚洲熟妇少妇任你| 母亲3免费完整高清在线观看| 精品少妇一区二区三区视频日本电影| 给我免费播放毛片高清在线观看| 999久久久国产精品视频| 国产精品亚洲av一区麻豆| 日韩精品免费视频一区二区三区| 欧美大码av| 亚洲中文av在线| 精品久久久久久久人妻蜜臀av | 免费人成视频x8x8入口观看| 咕卡用的链子| 天天添夜夜摸| 99久久99久久久精品蜜桃| 久久久久久大精品| 亚洲色图综合在线观看| 亚洲av美国av|