A review of monitoring,calculation,and simulation methods for ground subsidence induced by coal mining

Subsidence data acquisition methods are crucial to mining subsidence research and an essential component of achieving the goal of environmentally friendly coal mining.The origin and history of the existing methods of field monitoring,calcula-tion,and simulation were introduced.It summarized and analyzed the main applications,flaws and solutions,and improve-ments of these methods.Based on this analysis,the future developing directions of subsidence data acquisition methods were prospected and suggested.The subsidence monitoring methods have evolved from conventional ground monitoring to combined methods involving ground-based,space-based,and air-based measurements.While the conventional methods are mature in technology and reliable in accuracy,emerging remote sensing technologies have obvious advantages in terms of reducing field workload and increasing data coverage.However,these remote sensing methods require further technological development to be more suitable for monitoring mining subsidence.The existing subsidence calculation methods have been applied to various geological and mining conditions,and many improvements have already been made.In the future,more attention should be paid to unifying the studies of calculation methods and mechanical principles.The simulation methods are quite dependent on the similarity of the model to the site conditions and are generally used as an auxiliary data source for subsidence studies.The cross-disciplinary studies between subsidence data acquisition methods and other technologies should be given serious consideration,as they can be expected to lead to breakthroughs in areas such as theories,devices,software,and other aspects.

Equivalent Method of Integrated Power Generation System of Wind, Photovoltaic and Energy Storage in Power Flow Calculation and Transient Simulation

针对工程实际开展风光储联合发电系统在潮流计算和机电暂态仿真中的等值方法研究，旨在为大容量风光储联合发电系统的并网仿真分析奠定基础。将潮流计算的等值分为单元机组和集电系统2部分来研究。单元机组等值采用根据不同控制模式选取不同节点类型的方法，针对集电系统等值提出基于损耗不变原则的方法。等值模型和详细模型的算例结果表明，潮流计算等值方法具有较好的精度。在机电暂态仿真动态等值中，基于实际工程计算的最严重工况分析原则，提出运行在满出力点的单机“倍乘”等值模型，为工程计算中的风光储联合发电系统动态等值提供了一种解决方案。

Active learning accelerated Monte-Carlo simulation based on the modified K-nearest neighbors algorithm and its application to reliability estimations

This paper proposes an active learning accelerated Monte-Carlo simulation method based on the modified K-nearest neighbors algorithm.The core idea of the proposed method is to judge whether or not the output of a random input point can be postulated through a classifier implemented through the modified K-nearest neighbors algorithm.Compared to other active learning methods resorting to experimental designs,the proposed method is characterized by employing Monte-Carlo simulation for sampling inputs and saving a large portion of the actual evaluations of outputs through an accurate classification,which is applicable for most structural reliability estimation problems.Moreover,the validity,efficiency,and accuracy of the proposed method are demonstrated numerically.In addition,the optimal value of K that maximizes the computational efficiency is studied.Finally,the proposed method is applied to the reliability estimation of the carbon fiber reinforced silicon carbide composite specimens subjected to random displacements,which further validates its practicability.

Simulation and Techno-Economic Performance of a Novel Charge Calculation and Melt Optimization Planning Model for Steel Making

Process algorithm, numerical model and techno-economic assessment of charge calculation and furnace bath optimization for target alloy for induction furnace-based steelmaking is presented in this study. The developed algorithm combines the make-to-order (MTO) and charge optimization planning (COP) of the steel melting shop in the production of target steel composition. Using a system-level approach, the unit operations involved in the melting process were analyzed with the purpose of initial charge calculation, prevailing alloy charge prediction and optimizing the sequence of melt chemistry modification. The model performance was established using real-time production data from a cast iron-based foundry with a 1- and 2-ton induction furnace capacity and a medium carbon-based foundry with a 10- and 15-ton induction furnace capacity. A simulation engine (CastMELT) was developed in Java IDE with a MySQL database for continuous interaction with changing process parameters to run the model for validation. The comparison between the model prediction and production results was analyzed for charge prediction, melt modification and ferroalloy optimization and possible cost savings. The model performance for elemental charge prediction and calculation purpose with respect to the charge input (at overall scrap meltdown) gave R-squared, Standard Error, Pearson correlation and Significance value of (0.934, 0.06, 0.97, 0.0003) for Carbon prediction, (0.962, 0.06, 0.98, 0.00009) for Silicon prediction, (0.999, 0.048, 0.999, 9E -11) for Manganese Prediction, and (0.997, 0.076, 0.999, 6E -7) for Chromium prediction respectively. Correlation analysis for melt modification (after charging of ferroalloy) using the model for after-alloying spark analysis compared with the target chemistry is at 99.82%. The results validate the suitability of the developed model as a functional system of induction furnace melting for combined charge calculation and melt optimization Techno-economic evaluation results showed that 0.98% - 0.25% ferroalloy saving per ton of melt is possible using the model. This brings about an annual production cost savings of 100,000 $/y in foundry A (medium carbon steel) and 20,000 $/y in foundry B (cast iron) on the use of different ferroalloy materials.

Calculation Method of Passive Residual Heat Removal Heat Exchanger and Numerical Simulation

The tube inside and outside heat transfer mechanism of Passive Residual Heat Removal Heat Exchanger (PRHR HX) was analyzed. The calculation method of this special heat exchanger under natural convection condition in In-containment Refueling Water Storage Tank (IRWST) was carried out. The single-tube coupling model three-dimensional natural circulation in the IRWST was simulated numerically using Fluent. The heat transfer and flow characteristics of the fluid in IRWST were obtained. The comparison of the results between theoretical arithmetic and numerical simulation showed that the theoretical calculation method is suitable for the heat transfer calculation of PRHR HX.

Phase equilibrium of Cd_(1-x)Zn_xS alloys studied by first-principles calculations and Monte Carlo simulations

The first-principles calculations based on density functional theory combined with cluster expansion techniques and Monte Carlo （MC） simulations were used to study the phase diagrams of both wurtzite （WZ） and zinc-blende （ZB） Cdl_xZnxS alloys. All formation energies are positive for WZ and ZB Cdl-xZnxS alloys, which means that the Cdl-xZnxS alloys are unstable and have a tendency to phase separation. For WZ and ZB Cdl_xZnxS alloys, the consolute temperatures are 655 K and 604 K, respectively, and they both have an asymmetric miscibility gap. We obtained the spatial distributions of Cd and Zn atoms in WZ and ZB Cd0.sZn0.sS alloys at different temperatures by MC simulations. We found that both WZ and ZB phases of Cdo.sZn0.sS alloy exhibit phase segregation of Cd and Zn atoms at low temperature, which is consistent with the phase diagrams.

An enhanced image binarization method incorporating with Monte-Carlo simulation

We proposed an enhanced image binarization method.The proposed solution incorporates Monte-Carlo simulation into the local thresholding method to address the essential issues with respect to complex background,spatially-changed illumination,and uncertainties of block size in traditional method.The proposed method first partitions the image into square blocks that reflect local characteristics of the image.After image partitioning,each block is binarized using Otsu’s thresholding method.To minimize the influence of the block size and the boundary effect,we incorporate Monte-Carlo simulation into the binarization algorithm.Iterative calculation with varying block sizes during Monte-Carlo simulation generates a probability map,which illustrates the probability of each pixel classified as foreground.By setting a probability threshold,and separating foreground and background of the source image,the final binary image can be obtained.The described method has been tested by benchmark tests.Results demonstrate that the proposed method performs well in dealing with the complex background and illumination condition.

Analysis of typhoon wind hazard in Shenzhen City by Monte-Carlo Simulation

As one of the most serious natural disasters,many typhoons affect southeastern China every year.Taking Shenzhen,a coastal city in southeast China as an example,we employed a Monte-Carlo simulation to generate a large number of virtual typhoons for wind hazard analysis.By analyzing 67-year historical typhoons data from 1949 to 2015 using the Best Track Dataset for Tropical Cyclones over the Western North Pacific recorded by the Shanghai Typhoon Institute,China Meteorological Administration(CMASTI),typhoon characteristic parameters were extracted and optimal statistical distributions established for the parameters in relation to Shenzhen.We employed the Monte-Carlo method to sample each distribution to generate the characteristic parameters of virtual typhoons.In addition,the Yah Meng(YM)wind field model was introduced,and the sensitivity of the YM model to several parameters discussed.Using the YM wind field model,extreme wind speeds were extracted from the virtual typhoons.The extreme wind speeds for different return periods were predicted and compared with the current structural code to provide improved wind load information for wind-resistant structural design.

STUDY ON MESO-MECHANICAL SIMULATION OF GROUTING FOR UNFAVORABLE GEOLOGICAL CONDITIONS

The numerical simulation program of PFC2D(Particle Flow Code in 2 Dimension)particle flow based on the flow-solid coupling principle and,on its built-in FISHTANK function library and FISH language,defines the flow equation and pressure equation of fluid domain respectively,and carries out numerical simulation calculations on the diffusion process and,on the morphology and particle displacement of slurry during the slurry injection process.By adjusting the parameters of hist,n_bond,s_bond and measure in the PFC command flow,the tracking of granular body displacement changes is achieved,and the mesoscopic mechanism such as the diffusion law of soil slurry at different depths and the change of formation porosity is revealed.The numerical calculations show that:the grouting pressure has a significant effect on the alteration and destruction of the formation structure,and the fracturing effect becomes gradually worse with increasing adhesive strength,while the porosity increases significantly with increasing grouting pressure.Based on the elastic-plastic theory of the Mohr-Colomb criterion to theoretically derive the stress field of the soil around the borehole,it is pointed out that the mechanical mechanism of annular tension and radial compression is the fundamental reason for the appearance of fracturing grouting action mode.The increase of slurry viscosity is beneficial to improve the grouting effect of fracturing-compacting grouting,while the increase of friction coefficient has little effect on the grouting effect.The comparative analysis of the laboratory tests shows that the PFC2D simulation of the grouting process is feasible.

Temperature Field in Laser Line Scanning Thermography: Analytical Calculation and Experiment

The temperature field in laser line scanning thermography is investigated comprehensively in this work,including analytical calculation and experiment.Firstly,the principle of laser line scanning thermography is analyzed.On this basis,a physical laser line scanning model is proposed.Afterwards,based on Fourier transform(FT)and segregation variablemethod(SVM),the heat conduction differential equation in laser line scanning thermography is derived in detail.The temperature field of the composite-based coatings model with defects is simulated numerically.The results show that the laser line scanning thermography can effectively detect the defects in the model.The correctness of the analytical calculation is verified by comparing the surface temperature distribution obtained by analytical calculation and numerical simulation.Additionally,an experiment is carried out and the changeable surface temperature obtained by analytical calculation is compared with the experimental results.It shows that the error of maximum temperature obtained by analytical calculation and by experiment is 8%with high accuracy,which proves the correctness of analytical calculation and enriches the theoretical foundation of laser line scanning thermography.

Monte-Carlo Simulation of Surface Crack Growth Rate for Offshore Structural Steel E36-Z35

The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.

Principle of Line Configuration and Monte-Carlo Simulation for Shared Multi-Channel System

Based on the steady-state solution of finite-state birth and death process, the principle of line configuration for shared multi-channel system is analyzed. Call congestion ratio equation and channel utilization ratio equation are deduced, and visualized data analysis is presented. The analy-sis indicates that, calculated with the proposed equations, the overestimate for call congestion ratio and channel utilization ratio can be rectified, and thereby the cost of channels can be saved by 2000 in a small system.With MATLAB programming, line configuration methods are provided. In order to generally and intuitively show the dynamic running of the system, and to analyze,promote and improve it, the system is simulated using M/M/n/n/m queuing model and Monte-Carlo method. In addition, the simulation validates the correctness of the theoretical analysis and optimizing configuration method.

Simulating Calculation of Frictional Force between Mold and Billet

Assumption about the inner surface of mold is made according to the forming mechanisms of frictional force. The frictional force between billet and mold can be analysed by calculating the thermal shrinkage and the shearing strength of the meniscus.The model is used to calculate the maximum drawing speed at different conditions. The results are very satisfactory.

Study on the Relationship between Lamellar Spacing and Growth Rate inthe Regular Eutectic Growthby Monte-Carlo Simulation

A modified Monte-Carlo(MC) method to simulate the regular growth of binary eutectic alloys is presented. It is found that the growth rate has a linear dependence on the chemical potential difference between the solid and liquid; the relation between the lamellar spacing A and growth rate R accords well with the prediction of Jackson-Hunt(JH) theory unless the growth rate is very low.

MONTE-CARLO SIMULATION FOR ATOMIC DEPOSITION OF AMORPHOUS ELECTROLESS Ni_(80)P_(20) COATING

Atomic growth process and structure of Amorphous Electroless Coating have been studied, using Monte-Carlo simulation method. The simulation results of amorphous Ni80P20 coating show that PDFs are in accordance with practical values. The migrations of adatoms in coating's growth are different from that of solidification of amorphous materials. In some cases, the migrated adatoms in the process of growth of amorphous coating are not enough to occupy all vacancies and traps, so the amorphous coating is micro-porous. The immovable probability k and the largest migration distance of adatoms, which lie on the electroless bath components, affect the PDF, volume density and microporosity remarkably.

Monte-Carlo Simulation of Response Functions for Natural Gamma-Rays in LaBr_3 Detector System with Complex Borehole Configurations

Usually, there are several methods, e.g. experiment, interpolation experiment-based, analytic function, and Monte-Carlo simulation, to calculate the response functions in LaBr3（Ce） detectors. In logging applications, the experiment-based methods cannot be adopted because of their limitations. Analytic function has the advantage of fast calculating speed, but it is very difficult to take into account many effects that occur in practical applications. On the contrary, Monte-Carlo simulation can deal with physical and geometric configurations very tactfully. It has a distinct advantage for calculating the functions with complex configurations in borehole. A new application of LaBr3（Ce） detector is in natural gamma-rays borehole spectrometer for uranium well logging. Calculation of response functions must consider a series of physical and geometric factors under complex logging conditions, including earth formations and its relevant parameters, different energies, material and thickness of the casings, the fluid between the two tubes, and relative position of the LaBr3（Ce） crystal to steel ingot at the front of logging tube. The present work establishes Monte-Carlo simulation models for the above-mentioned situations, and then performs calculations for main gamma-rays from natural radio-elements series. The response functions can offer experimental directions for the design of borehole detection system, and provide technique basis and basic data for spectral analysis of natural gamma-rays, and for sonrceless calibration in uranium quantitative interpretation.

Computer Simulating Calculation on the Microstructure Evolutions during Hot Strip Rolling of C-Mn Steels

The program to predict the microstructure evolutions during hot strip rolling of C-M n steels has been developed in this paper, BV using this program, the microstructure changes with the processing parameters were analysed in detail. showing not only a good agreement of prediction with the measured values, but also entirely possibility to optimize hot strip rolling precess by computer simulation

Monte Carlo simulation based thinning and calculating method for noninvasive blood glucose sensing by near-infrared spectroscopy

Previous results show that the floating reference theory(FRT)is an effective tool to reduce the infuence of interference factors on noninvasive blood glucose sensing by near infrared spectros-copy(NIRS).It is the key to measure the floating reference point(FRP)precisely for the application of FRT.Monte Carlo(MC)simulation has been introduced to quantitatively in-vestigate the effects of positioning errors and light source drifts on measuring FRP.In this article,thinning and calculating method(TCM)is proposed to quantify the positioning error.Mean-while,the normalization process(NP)is developed to significantly reduce the error induced by light source drift.The results according to TCM show that 7 purm deviations in positioning can generate about 10.63%relative error in FRP.It is more noticeable that 1%fluctuation in light source intensity may lead to 12.21%relative errors.Gratifyingly,the proposed NP model can effectively reduce the error caused by light source drift.Therefore,the measurement system for FRPs must meet that the positioning error is less than 7 purm,and the light source drift is kept within 1%.Furthermore,an improvement for measurement system is proposed in order to take advantage of the NP model.

Monte-Carlo simulation of a stochastic differential equation

For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo（MC） techniques are considered to be more effective than other algorithms, such as finite element method or finite difference method. The inhomogeneity of diffusion coefficient strongly limits the use of different numerical techniques. For better convergence, methods with higher orders have been kept forward to allow MC codes with large step size. The main focus of this work is to look for operators that can produce converging results for large step sizes. As a first step, our comparative analysis has been applied to a general stochastic problem.Subsequently, our formulization is applied to the problem of pitch angle scattering resulting from Coulomb collisions of charge particles in the toroidal devices.