• <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.

    青青草视频在线视频观看| 在线观看国产h片| 18禁裸乳无遮挡动漫免费视频| 在线观看三级黄色| 我要看日韩黄色一级片| 丰满少妇做爰视频| 美女cb高潮喷水在线观看| 97精品久久久久久久久久精品| 纯流量卡能插随身wifi吗| 中文字幕人妻熟人妻熟丝袜美| 黄色配什么色好看| 日本av免费视频播放| 伦理电影免费视频| 秋霞在线观看毛片| 欧美变态另类bdsm刘玥| 国产精品国产三级专区第一集| 有码 亚洲区| 亚洲欧美清纯卡通| 国产在线一区二区三区精| 国产欧美亚洲国产| 日韩中字成人| 王馨瑶露胸无遮挡在线观看| 久久久a久久爽久久v久久| 夫妻性生交免费视频一级片| 我的女老师完整版在线观看| 热re99久久国产66热| 亚洲精品中文字幕在线视频 | 成人亚洲欧美一区二区av| 26uuu在线亚洲综合色| 高清欧美精品videossex| 国产乱来视频区| 久久毛片免费看一区二区三区| 男的添女的下面高潮视频| 亚洲av福利一区| 一本大道久久a久久精品| 中文字幕亚洲精品专区| 精品一区二区三卡| 国产深夜福利视频在线观看| 久热久热在线精品观看| 亚洲精品456在线播放app| 国产欧美日韩一区二区三区在线 | 久久久午夜欧美精品| √禁漫天堂资源中文www| 亚洲国产色片| 欧美高清成人免费视频www| 精品人妻偷拍中文字幕| 成人国产av品久久久| freevideosex欧美| 午夜免费鲁丝| 久久精品国产亚洲av天美| 看十八女毛片水多多多| 日本wwww免费看| 午夜福利视频精品| 国产深夜福利视频在线观看| 欧美xxⅹ黑人| 伦理电影免费视频| 久热久热在线精品观看| 国产白丝娇喘喷水9色精品| 久久精品国产自在天天线| 一区二区三区四区激情视频| 亚洲精品一区蜜桃| 国产熟女欧美一区二区| 一边亲一边摸免费视频| 亚洲人与动物交配视频| 午夜福利,免费看| 国产精品成人在线| 天堂8中文在线网| 韩国av在线不卡| 亚洲欧美成人精品一区二区| 欧美日韩精品成人综合77777| 啦啦啦中文免费视频观看日本| 国产一级毛片在线| 人人妻人人澡人人爽人人夜夜| 人人妻人人澡人人看| 一二三四中文在线观看免费高清| 精品亚洲成国产av| 男女啪啪激烈高潮av片| 亚洲情色 制服丝袜| 热re99久久精品国产66热6| 欧美激情国产日韩精品一区| 国产欧美日韩综合在线一区二区 | 久久亚洲国产成人精品v| 日本与韩国留学比较| 日本av手机在线免费观看| 欧美成人精品欧美一级黄| 日韩,欧美,国产一区二区三区| 3wmmmm亚洲av在线观看| 中文字幕精品免费在线观看视频 | 久久久久久久国产电影| 欧美三级亚洲精品| 丰满迷人的少妇在线观看| 热re99久久精品国产66热6| 精品少妇黑人巨大在线播放| 久久久久久久久大av| 美女福利国产在线| 国产美女午夜福利| 香蕉精品网在线| 免费黄网站久久成人精品| 在线观看免费视频网站a站| 美女大奶头黄色视频| 国产亚洲最大av| 韩国高清视频一区二区三区| 夫妻性生交免费视频一级片| 国国产精品蜜臀av免费| 嘟嘟电影网在线观看| 国产成人精品福利久久| 亚洲av男天堂| 51国产日韩欧美| 欧美老熟妇乱子伦牲交| 观看av在线不卡| 国产精品久久久久久av不卡| 久久精品熟女亚洲av麻豆精品| 天天躁夜夜躁狠狠久久av| 久久久久人妻精品一区果冻| 女的被弄到高潮叫床怎么办| 久久国产精品大桥未久av | 亚洲精品一二三| 深夜a级毛片| 在线观看免费高清a一片| 亚洲不卡免费看| 精品国产露脸久久av麻豆| 欧美精品国产亚洲| 久热久热在线精品观看| 日本黄大片高清| 成人特级av手机在线观看| 另类精品久久| 欧美xxⅹ黑人| 久久精品久久精品一区二区三区| 国产免费又黄又爽又色| av免费观看日本| 亚洲高清免费不卡视频| 人人妻人人看人人澡| 国产一区二区三区av在线| 免费观看无遮挡的男女| 日韩视频在线欧美| 国产 一区精品| 国产精品久久久久久av不卡| 国产午夜精品一二区理论片| 99热国产这里只有精品6| 丰满少妇做爰视频| 亚洲精华国产精华液的使用体验| 菩萨蛮人人尽说江南好唐韦庄| 大香蕉97超碰在线| 亚洲在久久综合| 搡女人真爽免费视频火全软件| 亚洲四区av| 久久毛片免费看一区二区三区| 日本色播在线视频| 久久国内精品自在自线图片| 蜜桃在线观看..| 欧美日韩精品成人综合77777| 插阴视频在线观看视频| 久久这里有精品视频免费| 国产黄频视频在线观看| 午夜福利影视在线免费观看| 少妇的逼水好多| 婷婷色综合大香蕉| 人妻少妇偷人精品九色| 久久久国产精品麻豆| av又黄又爽大尺度在线免费看| 秋霞在线观看毛片| 美女大奶头黄色视频| 啦啦啦在线观看免费高清www| 99热6这里只有精品| 熟女av电影| 亚洲一区二区三区欧美精品| 一区二区av电影网| 欧美精品国产亚洲| 丰满少妇做爰视频| 国产在线一区二区三区精| 欧美精品一区二区大全| 色视频在线一区二区三区| 国产乱人偷精品视频| 午夜福利网站1000一区二区三区| 久久综合国产亚洲精品| 国产免费又黄又爽又色| 国产精品无大码| 午夜老司机福利剧场| 久久这里有精品视频免费| 精华霜和精华液先用哪个| 能在线免费看毛片的网站| 久热这里只有精品99| 亚洲美女黄色视频免费看| 成人无遮挡网站| 国产av一区二区精品久久| 3wmmmm亚洲av在线观看| 午夜激情福利司机影院| 人体艺术视频欧美日本| 久久人人爽av亚洲精品天堂| 一个人免费看片子| 在线观看免费日韩欧美大片 | 乱人伦中国视频| 91久久精品电影网| av卡一久久| 秋霞伦理黄片| 色吧在线观看| 特大巨黑吊av在线直播| 天天躁夜夜躁狠狠久久av| 一二三四中文在线观看免费高清| 99久久精品一区二区三区| 看非洲黑人一级黄片| 男人和女人高潮做爰伦理| 国产毛片在线视频| 一级a做视频免费观看| 特大巨黑吊av在线直播| 男人舔奶头视频| 热re99久久国产66热| 国产伦精品一区二区三区四那| 婷婷色综合大香蕉| 韩国av在线不卡| 大陆偷拍与自拍| 亚洲欧洲精品一区二区精品久久久 | 国产综合精华液| 老女人水多毛片| 国产精品偷伦视频观看了| 亚洲综合精品二区| 日本黄色日本黄色录像| 丰满乱子伦码专区| 日本91视频免费播放| 亚洲精品,欧美精品| 国产免费一区二区三区四区乱码| 亚洲,欧美,日韩| 亚洲精品乱久久久久久| 久久久久久人妻| 久久免费观看电影| 亚洲va在线va天堂va国产| 亚洲精品一二三| 国产精品国产三级国产av玫瑰| 建设人人有责人人尽责人人享有的| 看十八女毛片水多多多| 亚洲情色 制服丝袜| 国产一区有黄有色的免费视频| 亚洲精品日韩av片在线观看| 乱系列少妇在线播放| 成人亚洲精品一区在线观看| 18禁在线播放成人免费| 日韩一区二区视频免费看| 99热国产这里只有精品6| 国内精品宾馆在线| 午夜影院在线不卡| 久久久久久久久久人人人人人人| 五月伊人婷婷丁香| 一级毛片久久久久久久久女| a 毛片基地| 欧美一级a爱片免费观看看| 永久免费av网站大全| 欧美精品国产亚洲| 乱系列少妇在线播放| 国产精品不卡视频一区二区| 亚洲真实伦在线观看| 人人澡人人妻人| 欧美日韩亚洲高清精品| 五月开心婷婷网| www.色视频.com| av不卡在线播放| 美女中出高潮动态图| 日韩av不卡免费在线播放| 91久久精品国产一区二区三区| 亚洲国产成人一精品久久久| 国产精品成人在线| 狂野欧美白嫩少妇大欣赏| 国产熟女欧美一区二区| 好男人视频免费观看在线| 人人妻人人爽人人添夜夜欢视频 | 久久久久久久久久久丰满| 只有这里有精品99| 亚洲欧美日韩另类电影网站| 22中文网久久字幕| 日韩亚洲欧美综合| 色网站视频免费| 久久99热这里只频精品6学生| 日韩一区二区视频免费看| 女人精品久久久久毛片| 乱人伦中国视频| 欧美日韩在线观看h| 亚洲国产色片| 在线看a的网站| 啦啦啦中文免费视频观看日本| 少妇的逼水好多| 国产精品人妻久久久影院| 纯流量卡能插随身wifi吗| 汤姆久久久久久久影院中文字幕| 国产美女午夜福利| 伊人亚洲综合成人网| 亚洲欧美清纯卡通| 大陆偷拍与自拍| 91在线精品国自产拍蜜月| 久久毛片免费看一区二区三区| 日韩中文字幕视频在线看片| 新久久久久国产一级毛片| 欧美日韩精品成人综合77777| 亚洲欧美日韩另类电影网站| 一级,二级,三级黄色视频| 99久久精品国产国产毛片| 国产在线免费精品| av国产精品久久久久影院| 少妇的逼水好多| 一级片'在线观看视频| 涩涩av久久男人的天堂| 九九久久精品国产亚洲av麻豆| 亚洲精品中文字幕在线视频 | 国产精品嫩草影院av在线观看| 色94色欧美一区二区| 国产精品女同一区二区软件| 中文资源天堂在线| 搡老乐熟女国产| 大又大粗又爽又黄少妇毛片口| 色吧在线观看| av在线老鸭窝| 亚洲精品乱久久久久久| tube8黄色片| 国产伦精品一区二区三区视频9| 男人狂女人下面高潮的视频| 国产精品偷伦视频观看了| 欧美xxxx性猛交bbbb| 成年美女黄网站色视频大全免费 | 久久国产精品男人的天堂亚洲 | 色视频www国产| 在线精品无人区一区二区三| 国产欧美日韩综合在线一区二区 | 永久免费av网站大全| 午夜福利影视在线免费观看| 一级毛片黄色毛片免费观看视频| 夜夜骑夜夜射夜夜干| av天堂中文字幕网| 好男人视频免费观看在线| 精品酒店卫生间| 国产黄片视频在线免费观看| 成年人免费黄色播放视频 | a级一级毛片免费在线观看| 亚洲自偷自拍三级| 中文字幕久久专区| 国产欧美另类精品又又久久亚洲欧美| 99热这里只有是精品50| 99热网站在线观看| 精品亚洲成国产av| 日韩中文字幕视频在线看片| 欧美一级a爱片免费观看看| 久久久久精品久久久久真实原创| 中文字幕制服av| 亚洲美女黄色视频免费看| 日韩一区二区三区影片| av播播在线观看一区| 麻豆精品久久久久久蜜桃| 亚洲精品亚洲一区二区| 久久精品国产亚洲网站| freevideosex欧美| 性色av一级| 亚洲精品,欧美精品| 人妻系列 视频| 男人舔奶头视频| 狂野欧美白嫩少妇大欣赏| 美女大奶头黄色视频| 麻豆成人午夜福利视频| 久久久久久久久久久久大奶| 亚洲人与动物交配视频| 亚洲一区二区三区欧美精品| 国产精品国产三级国产专区5o| 国产黄色免费在线视频| 国产亚洲一区二区精品| 欧美日韩在线观看h| 在线观看www视频免费| 蜜臀久久99精品久久宅男| 亚洲欧美一区二区三区国产| 蜜臀久久99精品久久宅男| 日韩一区二区三区影片| 国产真实伦视频高清在线观看| 天堂俺去俺来也www色官网| 最新的欧美精品一区二区| 秋霞在线观看毛片| 如日韩欧美国产精品一区二区三区 | 六月丁香七月| 久久久久久久久久久免费av| 性色avwww在线观看| 大陆偷拍与自拍| 国语对白做爰xxxⅹ性视频网站| av天堂中文字幕网| 99久久精品国产国产毛片| 日本免费在线观看一区| 女性被躁到高潮视频| 永久免费av网站大全| 亚洲欧美中文字幕日韩二区| 亚洲一区二区三区欧美精品| 一级黄片播放器| 久久99精品国语久久久| 国产精品久久久久久精品古装| 亚洲欧美精品专区久久| 特大巨黑吊av在线直播| 少妇人妻 视频| 自线自在国产av| 亚洲怡红院男人天堂| 亚洲av在线观看美女高潮| 久久久久久伊人网av| 涩涩av久久男人的天堂| 一级片'在线观看视频| h视频一区二区三区| 色哟哟·www| 少妇的逼水好多| 青春草亚洲视频在线观看| 欧美+日韩+精品| 午夜免费鲁丝| 久久女婷五月综合色啪小说| 亚洲综合精品二区| 中文天堂在线官网| 麻豆乱淫一区二区| 国产熟女午夜一区二区三区 | 亚洲性久久影院| 日韩中文字幕视频在线看片| 韩国高清视频一区二区三区| 九九在线视频观看精品| 亚洲欧洲精品一区二区精品久久久 | 丝袜脚勾引网站| 天天躁夜夜躁狠狠久久av| 99热6这里只有精品| 国产午夜精品久久久久久一区二区三区| 美女视频免费永久观看网站| 亚洲伊人久久精品综合| 亚洲三级黄色毛片| 噜噜噜噜噜久久久久久91| 在线精品无人区一区二区三| 亚洲高清免费不卡视频| 日本黄大片高清| 最近2019中文字幕mv第一页| 少妇人妻一区二区三区视频| 亚洲av在线观看美女高潮| 18禁裸乳无遮挡动漫免费视频| 日韩电影二区| 18禁在线播放成人免费| 国产白丝娇喘喷水9色精品| av国产久精品久网站免费入址| 夜夜看夜夜爽夜夜摸| 男人爽女人下面视频在线观看| 日韩人妻高清精品专区| 美女视频免费永久观看网站| 国产成人91sexporn| 在线观看一区二区三区激情| 狂野欧美激情性xxxx在线观看| 久久av网站| 久久6这里有精品| 99精国产麻豆久久婷婷| 丝袜脚勾引网站| 自拍偷自拍亚洲精品老妇| 纯流量卡能插随身wifi吗| 99热网站在线观看| 久久这里有精品视频免费| 久久综合国产亚洲精品| av福利片在线观看| 人妻少妇偷人精品九色| 国产一区二区三区av在线| 亚洲精品乱码久久久v下载方式| 亚州av有码| 日韩一区二区视频免费看| 9色porny在线观看| 黄色日韩在线| 久久久久精品性色| 国产一区二区三区av在线| 亚洲图色成人| 2018国产大陆天天弄谢| 亚洲久久久国产精品| 精品卡一卡二卡四卡免费| 国产精品一区二区性色av| 亚洲经典国产精华液单| 成人无遮挡网站| 91久久精品电影网| 91精品国产九色| av免费在线看不卡| 亚洲人成网站在线观看播放| 一二三四中文在线观看免费高清| 大话2 男鬼变身卡| 热re99久久精品国产66热6| 国产精品秋霞免费鲁丝片| 国产成人精品一,二区| av一本久久久久| 久久鲁丝午夜福利片| 日韩免费高清中文字幕av| 日韩av不卡免费在线播放| 极品教师在线视频| av又黄又爽大尺度在线免费看| 18+在线观看网站| 91精品国产国语对白视频| 日本黄大片高清| 人妻系列 视频| 亚洲欧美成人综合另类久久久| 黄色欧美视频在线观看| 久久人人爽av亚洲精品天堂| 精品国产露脸久久av麻豆| 97精品久久久久久久久久精品| 国产一区二区在线观看日韩| 亚洲av成人精品一区久久| 青春草亚洲视频在线观看| 黄色日韩在线| 最近的中文字幕免费完整| 少妇丰满av| 欧美97在线视频| 久久综合国产亚洲精品| 高清午夜精品一区二区三区| 免费在线观看成人毛片| 日本vs欧美在线观看视频 | 一二三四中文在线观看免费高清| 久久99一区二区三区| 韩国高清视频一区二区三区| 女人精品久久久久毛片| 欧美老熟妇乱子伦牲交| 丰满少妇做爰视频| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久精品古装| 国产黄频视频在线观看| 国产精品欧美亚洲77777| 成人美女网站在线观看视频| 一级毛片 在线播放| 老熟女久久久| 日韩精品有码人妻一区| 亚洲欧美一区二区三区国产| 久久久久久久久久久久大奶| 51国产日韩欧美| 免费播放大片免费观看视频在线观看| 久久婷婷青草| 91久久精品国产一区二区成人| 我的老师免费观看完整版| 啦啦啦啦在线视频资源| 丰满人妻一区二区三区视频av| 国产精品人妻久久久久久| 亚洲高清免费不卡视频| 午夜日本视频在线| 亚洲精品一区蜜桃| 国产高清国产精品国产三级| 欧美区成人在线视频| 日韩一区二区三区影片| 我的女老师完整版在线观看| 久久av网站| 美女cb高潮喷水在线观看| 亚洲欧美精品专区久久| 亚洲av成人精品一区久久| 视频中文字幕在线观看| 美女脱内裤让男人舔精品视频| 人妻系列 视频| 一级毛片 在线播放| 亚洲无线观看免费| 欧美 亚洲 国产 日韩一| 日本av免费视频播放| 国产日韩欧美在线精品| 日本爱情动作片www.在线观看| 国产精品一区二区性色av| 国产精品蜜桃在线观看| 国产精品免费大片| 黄色日韩在线| 国产视频内射| 亚州av有码| 伊人久久精品亚洲午夜| 久久久久精品性色| 草草在线视频免费看| 少妇被粗大的猛进出69影院 | 桃花免费在线播放| 日日摸夜夜添夜夜添av毛片| 我要看黄色一级片免费的| 国产真实伦视频高清在线观看| 免费少妇av软件| 久久97久久精品| 久久久久久久久久久久大奶| 国产在线一区二区三区精| 亚洲av电影在线观看一区二区三区| 精品国产露脸久久av麻豆| 热re99久久精品国产66热6| 在线观看三级黄色| 少妇高潮的动态图| 国产成人aa在线观看| 亚洲美女视频黄频| 国产精品人妻久久久久久| 亚洲欧美清纯卡通| 国产精品无大码| 成人综合一区亚洲| 啦啦啦视频在线资源免费观看| 日韩一本色道免费dvd| 极品少妇高潮喷水抽搐| 人人妻人人添人人爽欧美一区卜| 边亲边吃奶的免费视频| 国产在线一区二区三区精| 18禁在线播放成人免费| 高清黄色对白视频在线免费看 | 97超视频在线观看视频| 日日摸夜夜添夜夜添av毛片| 在现免费观看毛片| 亚洲av成人精品一二三区| 欧美老熟妇乱子伦牲交| 老司机亚洲免费影院| 又大又黄又爽视频免费| 午夜免费男女啪啪视频观看| 亚洲va在线va天堂va国产| 亚洲精品日韩在线中文字幕| 国产色爽女视频免费观看| 高清不卡的av网站| 亚洲精品乱码久久久久久按摩| 久久精品国产鲁丝片午夜精品| 青青草视频在线视频观看| 国产女主播在线喷水免费视频网站| 黑人巨大精品欧美一区二区蜜桃 | 九色成人免费人妻av| 夫妻性生交免费视频一级片| 丰满饥渴人妻一区二区三| 日韩三级伦理在线观看| 国产欧美亚洲国产| 国产成人aa在线观看| 色吧在线观看| 欧美日本中文国产一区发布| 国产中年淑女户外野战色| 亚洲国产精品专区欧美| 欧美老熟妇乱子伦牲交| 中文乱码字字幕精品一区二区三区| 亚洲美女视频黄频| 色5月婷婷丁香| 人妻少妇偷人精品九色| 成人国产av品久久久| 精品酒店卫生间| 三级国产精品片| 国产国拍精品亚洲av在线观看|