Search Documents

By R. R. Song, Q. Li, H. Xing, X. C. Wei

Geomechanical modelling of fault stability is an important part of CO2 capture and storage (CCS) research program to ensure its viability and long-term safety. Storage of CO2 in deep subsurface reservoirs requires estimates of sustainable fluid pressures that will not induce fracturing or fault reactivation which could cause fault slip. More importantly, fault slip gives rise to seismic waves, which inflict severe damage to the cap rock of the reservoirs and may result in leakage pathways for the buoyant plume and thus compromise the integrity of seal formations. In this study, we generated a numerical model tuned towards conditions usually encountered in the Shenhua CCS demonstration project in Ordos basin. In particular, with the purpose of evaluating the sensitivity of the parameters as well as their relationships with the moment magnitude of the seismicity, a variance-based sensitivity analysis was adopted for importance ranking of the uncertain input parameters of the given model. However, such flow models could involve one or several spatial inputs, for instance, fault dip, distance between fault and the injection well, friction angle of fault, permeability of fault, Young's Modulus of fault, injecting pore pressure and the stress field. Studying the sensitivity to each of these spatial inputs motivated the present work. The results led to the following conclusions: fault depth and fault dip significantly affect the behavior of the fault and the fault slip becomes violent with the fault depth increasing and fault dip decreasing. Furthermore, the results indicate that slip time increases with the fault dip decreasing and decreases with permeability of fault increasing. In addition, the calculation of seismic moment and moment magnitude of an earthquake demonstrates that it rises significantly with the decrease of fault dip and the increase of fault depth, and it decreases with the rising stress field. The results show that variations of the only controllable parameter injection pressure do not show great influence to the fault behavior compared with other geological factors, which imply site selection is extremely important for the CCS projects.

Jan 1, 2015

By R. P. Townsend

It is a widely accepted "fact" that a Koepe friction hoist is generally the best solution for production hoisting when compared to a double drum hoists. This would appear to be true if one looks at the initial capital cost of just the hoist, and also if one compares the size of the motors selected for the duty. But is this the whole story? What are the longer term costs involved in the hoist's operation, and what are the relative efficiency's of the two designs? Perhaps the Koepe hoist design enjoys a better reputation than it deserves.

Jan 1, 2011

By Wael Abdellah

Mine haulage drifts are the primary access to the mining blocks of an orebody in a multilevel mining system of a tabular ore deposit. The stability and functionality of haulage drifts are thus crucial to the success of a mining operation. Drift instability could lead to serious consequences such as injuries, production delays and higher operational cost. In this paper, the performance of the haulage drift stability is evaluated on the basis of the primary rock support system comprising 1.8 m resin grouted rebars in the drift walls and 2.1 m long bolt in the back. Two failure criteria namely the extent of Mohr-Coulomb yield zones around the drift, and linear elastic brittle shear failure are adopted and compared, with respect to lower level and same-level mining and filling sequence in the vicinity of the haulage drift. The case study is one of the #1 Shear East zone of the Garson Mine in Sudbury, Ontario. Random Monte Carlo (RMC) technique, which has been previously established to assess the haulage drift stability, is employed in this investigation. The RMC technique is used in conjunction with finite difference modelling software FLAC for random assignment of model input parameters in the FLAC grid. Deterministic model results reveal the drift behaviour in terms of deformation (convergence) and mining-induced stress distribution. Comparison between the two different criteria is carried out to determine which method is most suitable in evaluating unsatisfactory drift performance. The results are presented in terms of probability of instability and categorized with respect to failure condition and mining step. It is shown that the brittle shear failure condition based on linear elastic response calls for enhanced support system at a later mining step than what Mohr- Coulomb yielding condition does. The paper presents a detailed discussion of the stochastic analysis results.

Jun 1, 2012

By N. Erarslan, M. Ghamgosar

The size of rock fragments depends on the nature of the loading applied the strength of the minerals and the texture of the rock. Research in the literature demonstrates that there is a distinct difference in the rock fracturing process for static and cyclic loading. The tensile strength of brittle rocks is also one of the most significant parameters influencing the rock fracturing process and the propagation of micro-fractures on the crushing, blasting and cutting of rock. The Brazilian indirect tensile test, one of the methods approved by the International Society for Rock Mechanics (ISRM), has been used to determine rock fracture behaviour under different cyclic loading modes. In this study, the dimensions of Brisbane Tuff specimens have been selected according to the method outlined by the ISRM. The main aim of this study is to evaluate the effect of cyclic loading amplitude and frequency on the Fracture Process Zone (FPZ) at the micro-, meso- and macro-scales. The post-failure behaviour of the Brisbane Tuff specimens was obtained under cyclic loading with frequencies varied from 1 to 10 Hz. The amplitude was varied from 10% of Ultimate Tensile Strength (UTS) up to 80% UTS. Laboratory investigations of rock fragment size revealed that different loading frequencies and amplitudes produce different FPZ dimensions. Multivariate Analysis of Variance (MANOVA) results demonstrated that the FPZ dimension is relatively sensitive to the frequency change. Fragment size analysis showed that the FPZ is a consequence of the initiation, propagation and coalescence of micro-fractures and depends on the nature of the loading. Laboratory evidence showed that the FPZ dimension, for the cyclic tests performed on Brisbane tuff, is ranged from 0.07 to 20 mm. The achievement of this study highlights a key finding in the understanding of fracture patterns and propagation models in the FPZ under cyclic loading, which is the dominate mechanical action in the oscillator disc cutter (ODC).

Jan 1, 2015

By W. A. Gow, V. F. Harrison, H. H. Mccreedy

A semi-continuous bench-scale, countercurrent-decantation system, in which the ore was moved countercurrently to the flow of leaching solution, was used to conduct acid-leaching tests on uranium ore with the aid of bacteria. This method, which involved a low degree of agitation, was successfully applied to extract uranium in practical amounts and at practical rates. The acidity of the leaching solution was the main variable investigated in the bench-scale tests done at the Mines Branch, Ottawa. Over 90 11er cent of the uranium was extracted in a six-stage system. The results suggest that this technique could be ap-plied in the utilization of bacterial oxidation and leaching as a method for treating Elliot Lake uranium ore.

Jan 1, 1969

The dynamic breakage process of rocks by the impact of a hydraulic hammer is more complex than one normally assumes. It is even more difficult when one has to conduct the breakage under water. In such cases one does not know exactly what the surface topography and composition are at the bottom of the sea or river and, hence, it is more difficult to select the best position of the hammer before starting the hammering. Typically, an industrial hammer has a diameter of about 250 mm, a length of about 1.25 to 1.75 m and a mass of several tons. The hammer frequency can be up to 50 impacts per minute. The general principles of rock hammering can be easily summarized, but this knowledge is not sufficient to adapt adequately the method for less common rock conditions (e.g. less brittle rock, presence of vertical and subvertical fractures, etc.). That is the reason why we initiated a research project, whereby a scale model has been developed to conduct experiments in the laboratory and whereby discrete element simulations are being conducted. Both aim at a better understanding of the basic mechanisms of rock hammering. The scaled experimental model allows varying the impact energy and the static weight independently. Within the experiment, the hammer is in contact with the rock specimen under a certain confining pressure and with a certain contact pressure due to the static weight. A weight falls on the hammer, creating the dynamic action or impact; the height over which the weight falls can be varied. The penetration is recorded after each impact. After a certain penetration the test is stopped and the internal crack patterns are described. Some simulations have been conducted, allowing the modeling of individual fracture initiation and growth during the penetration of the hammer by either a constant displacement rate (using the UDEC-code) or by introducing the dynamic process (using the PFC-code). Some of the observed features in the laboratory experiments can be visualized in the simulations. However, further research is needed to allow the simulations being an adequate tool to well understand the dynamic process in laboratory experiments, as well as the in-situ behavior.

Aug 1, 2013

By R. H. Schonewille

One of the many challenges in developing new metallurgical processes is choosing the best option for validating process design parameters and establishing the process operating window. Historically this has been achieved by constructing pilot or demonstration plants of considerable size and complexity, to enable integrated operation of all key process units for a significant duration. While this would appear to be the best option for minimizing technical risk, it is often impractical from a perspective of cost, schedule, resources or bulk sample acquisition. Such an undertaking can even provide a false sense of security, as the optics of running a complete scaled-down version of a commercial plant suggest a high level of process understanding, regardless of the challenges encountered during operation. Xstrata Nickel has undertaken numerous process development and improvement projects over the past 15 years to support green field project development and continuous improvement in existing operations. The approaches used have included; computer simulation, cold flow modeling, construction of pilot test units in third party facilities, commercial plant referencing and the more conventional construction of integrated test plants. Regardless of the strategy employed, the keys to success have included; managing stakeholder expectations, involving the commercial plant design engineers in the planning stage, anticipating subsequent data and sample requirements and most importantly, constructing and operating complex demonstration plants only as a last resort. This paper provides examples of various process development approaches used by Xstrata Nickel and discusses the rationale specific to each project.

Jan 1, 2011

By J. E. Dutrizac

Canada is currently the world's largest supplier of zinc concentrates and is also a major producer of zinc metal. This dominant position in the world's zinc market is a fairly recent circumstance and, in fact, the country's zinc industry was founded only in the early part of the twentieth century. By contrast, the American zinc industry was well established by the mid-nineteenth century, largely based on technology originally imported from the still older European practice. Initial attempts to create a viable Canadian zinc smelter were based on the traditional practice and were unsuccessful; success came only with new electrolytic technology, much of which was developed domestically. The purpose of this paper is to outline the early attempts to produce metallic zinc in Canada and to trace the gradual development of these first efforts into a viable zinc industry based on modern science and technology.

Jan 1, 1999

By O. Hall

STOPING at Noranda" is a sequel to "Mining at Noranda". "Mining at Noranda", presented in 1937 (1), sketches Noranda's history, metallurgical problems, exploration, mining methods, mining details, pillar extraction, and costs. "Stoping at Noranda" deals in greater detail with our stoping, the use of sub-level drifts instead of benches, and the drilling out of the stopes with jointed rock drill steel and diamond drills. Sun-LEVEL DRIFT STOPING: Sub-level drift stoping is almost universal throughout the mine and is the standard sroping method in both sulphide ore and rhyolite flux ore. The features that have led to the adoption of this standard are: (1) Continuous symmetrical drilling in sub-level drifts with no delay for blasting. (2) Large-tonnage blasts Saturday night after the drillers have left. (3) Comparative freedom from gas and smoke during the working days. ( 4) The daily safety achieved by having all stope drilling done in drifts. (5) The safety and security achieved by drilling out stopes, floors, and pillars before there is any weighting or cracking.

Jan 1, 1939

By Dennis Z. F. Mraz

"Because rock salt and pot ash formations contain particles capable of elastic deformation, the material has usually been compared with various rheological models containing elastic elements. Extensive laboratory testing of salt and potash samples resulted in various theoretical conclusions based on conditions quite different from those in a deep mine. It is felt that sufficient data have now been collected from underground measurements and observations to demonstrate the inelastic behaviour of salt masses in a deep mine. The conditions in deep salt strata can be defined as near hydrostatic, because the acting pressures are very high compared to the negligible elastic limit of the salt. The pressure field and creep around the single opening is a function of its shape, size and age. The same factors, together with extraction, determine the behaviour of the mining a rea. The application of the correct analytical method in mine design can significantly improve the excavation techniques and predetermine the long-term behaviour of openings."

Jan 1, 1973

By Wm. D. Jr. Sheldon

THIS paper presents practical data to aid in the selection of a mine ventilating fan, and outlines considerations to be made to ensure satisfactory operation and optimum efficiency of the fan. Much has been published on the general subject of mine ventilation, but the writer has found few specific data on factors and considerations which govern the selection of a fan-data which should be available to the operator if his fan is to give satisfactory service within the range of the high-efficiency portion of its performance curve and deliver the specified air capacity when operating against the anticipated range of resistance of the mine. The total quantity of fresh air furnished by the ventilation system must be adequate to provide, and to ensure at all times, safe and healthful working conditions in the mine. It is assumed that the mine operator is able to determine the air quantities required for the various parts of the mine.

Jan 1, 1941

By K. Forsberg, A. Rasmuson, R. Rodriguez Varela, J. Martinez De La Cruz

"This study analyzes the separation of rare earth elements by means of a supported liquid membrane (SLM) with the final objective of recovering REE from an apatite concentrate (Alemrajabi, Forsberg, & Rasmuson, 2015; Mohammadi, Forsberg, Martinez de la Cruz, & Rasmuson, 2015). In this study, the feed phase consists of an aqueous phase containing REE, the membrane phase contains the organic solvent and carrier (D2EHPA and EHEHPA), and the stripping phase is a hydrochloric acid solution. An evaluation of mass fluxes, separation efficiency, and membrane performance is conducted. In order to analyze the suitability of liquid membranes in the extraction of REE, a deeper insight into the mechanisms causing instability is needed. In this study, surface tension is considered as the key factor of Marangoni instabilities and spontaneous emulsification.INTRODUCTIONThe demand of REE and their compounds has increased considerably due to their special physicochemical features that make them crucial elements for numerous high-tech applications, including wind turbines, hybrid cars, mobile phones, and defense technologies. They are also extensively utilized in traditional industries such as nuclear energy, petroleum, metallurgy, agriculture, and textiles (Xie, Zhang, Dreisinger, & Dyole, 2014; Zhang, Wu, Bian, & Zeng, 2014).Solvent extraction is widely used for the separation of REE (Xie et al., 2014; Zhao, Huo, Pan, Li, & Liu, 2014). Di-(2-ethylhexyl)-phosphoric acid (D2EHPA) and 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (EHEHPA) are two of the most effective extractants for the separation of rare earth elements among organophosphorus extractants (Gupta & Krishnamurthy, 2005). However, the difficulty in stripping loaded D2EHPA extractant, especially for heavy REE, has limited its utilization. Yttrium can be stripped at lower acidities in EHEHPA systems compared to D2EHPA systems (Gupta et al., 2005), but the extraction efficiency is lower. Conventional solvent extraction by mixer-settlers is restricted by its large solvent inventory and loading limits, and suffers from phase disengagement, emulsification, loss of solvent, and flooding problems (Swain, Jeong, Lee, & Lee, 2007). In this context, supported liquid membranes (SLMs) present a potential alternative due to their advantages, being low solvent inventory, process simplicity, low energy and operational costs, great enrichment ratio, and combination of extraction and stripping into one single step. A SLM is a system in which two aqueous phases are separated by an organic liquid membrane suspended in the pores of a porous support by capillary forces. One of the aqueous phases is the feed solution with the solute and the other is the stripping or receiving solution. The solute diffuses through the liquid membrane and then is extracted into the stripping solution (Kislik, 2010)."

Jan 1, 2016

By W. J. Adams

It is recognized that within life cycle assessment that the eco-toxicity potentials (TPs) derived for metals require additional research to improve their reliability. Current models used to develop TPs do not include important fate process for metals such as speciation and bioavailability. Toxicity potentials are standard values used in LCA to enable a comparison of potential toxic impacts between substances. Huijbrechts (2001) calculated toxicity potentials for 181 substances utilizing the USES model for six impact categories, i.e. freshwater aquatic ecotoxicity, marine aquatic ecotoxicity, freshwater sediment ecotoxicity, marine sediment ecotoxicity, terrestrial ecotoxicity and human toxicity (not considered here). Marine ecotox potential for metals is typically high using this approach, whereas most open ocean metal concentrations are very low. A sensitivity analysis of the USES-LCA model was undertaken using copper as an example. Results indicate: fate processes on the global scale are driving the marine water and sediment TPs for copper most; the second most important fate parameter is sedimentation from the water and burial in deep sediment layers; inclusion of speciation and bio-availability processes would be very beneficial; more accurate estimates of input parameters should be used, and the PNECwater and PNECsediment and all partition coefficients for metals should be updated in the model.

Jan 1, 2007

By M. Bustin, A. Hosseinian, G. Preisig

Hydraulic fracturing and hydraulic shearing techniques are essential for the development of unconventional oil & gas and geothermal resources. Such techniques will also play a key role in the development of CO2 sequestration, and management of high stresses encountered with deep mining. However, the optimal design of hydraulic injections is often limited by continuum treatments of the rock mass without consideration of the influence of the natural fractures present. These may interact with the hydraulic fracture to promote offsets or arrest, or may cause the natural fractures to undergo shear and dilation. Given the limited access to the rock afforded by deep boreholes common to oil & gas and geothermal projects, advanced numerical modeling is critical for understanding key mechanisms and interactions. In this work, we first present a coupled hydromechanical distinct-element model focusing on the propagation of a hydraulic fracture. Normal and shear displacements arising from the generated hydraulic fractures are subsequently integrated into a fracture flow model for detailed analysis of specific flow behaviour and fracture permeability. For this, a scanned 3-D profile of a natural fracture is used to correlate surface roughness to dilation as a function of shearing. Fracture flow is computed at different time steps corresponding to different fracture propagation states and geometries. Modeling results indicate that the opening of a tensile hydraulic fracture is accompanied by millimetre-scale shear displacements associated with slip and wedging between adjacent blocks. This may leads to a mismatch between the upper and lower walls of the fracture if asperities are irregular and strong enough. However, fluid-induced tensile fractures have low Joint Roughness Coefficients (JRC), and associated millimetre-scale shear displacements are too small to ensure a significant offset of asperities that prevent a fully-elastic closure when the hydraulic fracture is depressurized. Injection of proppant is thus critical to keeping hydraulic fracture apertures open. Future work is planned that will see the models calibrated against field data including microseismic event locations, magnitudes and source mechanisms, tiltmeter deformations, stress change, and pore pressures.

Jan 1, 2015

By Michael Li

Paste backfill as a means of supporting underground mining is gaining popularity. Despite this trend, engineering and design of mine paste backfill systems are still largely based on empirical methods rather than established scientific principles. Because of a lack of sound understanding of mineral paste flow, some degree of trial and error is often involved during commissioning and operating of new paste fill systems. New operations often have to be built with flexibilities and contingencies, escalating the overall project costs. This paper reviews the applicability of some scientific and engineering principles to mineral paste flows. Recent experimental data are used to illustrate important points. The phenomenon of mineral paste flow is examined under the context of improving the efficiency and reliability of underground delivery. Gaps in our understanding of mineral paste rheology are pointed out and from those some future research directions naturally emerge. Status and problems of the current engineering and design practice are also briefly discussed.

Nov 1, 2002

By Y. W. Pan, Y. R. Liu

The disturbed natural rock mass tends to equilibrium state accompanying with non-equilibrium evolution process which has significant influence on the long-term stability of engineering structure. The deformation development in the non-equilibrium evolution is determined by two intrinsic processes, trend to equilibrium and damage development, which can be illustrated by evolution of excavated tunnel. Thus the creep and damage must be taken into account in study of non-equilibrium evolution process. Within Rice’s irreversible thermodynamics with internal state variable (ISV), the conjugate force function which equals to zero is a convex surface which divide the stress space into two parts, called equilibrium and nonequilibrium stress space respectively. The concept of effective stress of non-equilibrium has been proposed, in which the stress difference (d - d ) are the effective driving forces of non-equilibrium evolution process with hardening and damage effective. The formulations of inelastic strain rate and energy dissipation rate have been recast in the form of effective stresses. The inelastic complementary energy (ICE) is proposed to measure intensity of non-equilibrium evolution of material system. The specific equations of conjugated forces and kinetic equations of internal state variables are presented. The ICE and rate of energy dissipation are all calculated to be used to evaluate long-term stability of geo-engineering. The results show that rate of energy dissipation may be more convenience and flexible than ICE for internal state variable constitutive equation of rock mass.

Jan 1, 2015

By M. Walawalkar, C. K. Nichol, G. Azimi

Rare earth elements (REEs) are essential to the production of technologically advanced products, including wind turbines, electric and hybrid cars, automotive catalytic convertors, hand held electronic devices, and military defense systems. Their application in environmentally beneficial technologies puts them at a high priority for ameliorating environmental concerns. Because of their increasingly high demand and limiting resources, it is highly beneficial to find new sources for them. The present study is focused on leaching REEs from a phosphate based fertilizer by-product, known as phosphogypsum (PG). Global estimated world production of PG is about 240 million tonnes that is stored in the form of wet stacks adjacent to fertilizer production facilities, occupying vast amounts of land. PG also contains trace metals, which are detrimental to the environment, and hence it is crucial to monitor these stacks. We obtained our PG sample from Agrium Inc. (a fertilizer manufacturing company in Canada). We performed a comprehensive characterization of the PG sample and found that it predominantly contains Y, La, Ce, Nd, and Sm. For leaching purposes, we used hydrometallurgical acid leaching techniques. We investigated three acids over a range of operating conditions such as acid concentration, temperature, solid-to-liquid ratio, and residence time. We achieved an average leaching of 81% with hydrochloric acid as leachant using operating conditions that can be easily replicated in industrial settings. If our process is implemented, not only will it reduce the size of the existing PG stacks (mitigating their environmental risks), but also enable the development of new innovative routes for extracting REEs, making countries like Canada less vulnerable to supply disruptions.

Jan 1, 2016

By David Stone

North American Palladium is the owner/operator of the Lac des Iles mine located 115 km north of Thunder Bay, Ontario. The mine has produced an average of 65,000 ounces of palladium since 1994. Lac des Iles is currently undergoing the first phase of a major expansion that will increase ore production from 2,400 tons per day to 15,000 tons per day. A future expansion is being contemplated for 30,000 tons per day. The current expansion will increase palladium output to 250,000 ounces per year. The mine plan for the ultimate pit results in an excavation 1 km long by 800 m wide and over 450 m deep. This makes the ultimate pit at Lac des Iles one of the ten deepest open pits in the world, and the largest in Canada. This paper outlines the geotechnical data collection and design methodology for the ultimate open pit at Lac des Iles.

May 1, 2001

By M. A. Perras, H. A. Kasani, G. Walton

The development of an underground repository for the storage of spent nuclear fuel relies on a multi-barrier system to isolate the used fuel from the environment – ultimately, the biosphere. These barriers include engineered copper and steal canisters, dense bentonite, and the natural rock as well. In order to reduce the potential flow pathway along the excavation induced damage zone, various slots perpendicular to the excavation axis have been proposed. These slots are called cut-off seals and generally are filled with bentonite as it has high swelling potential and very low permeability. Many different geometries have been proposed, but clear guidelines on the selection of the optimum geometry have yet to be developed. In order to assess the optimum geometry, three shapes have been selected for this numerical study, including a trapezoid, a triangle, and a thin rectangular slot. The cut-offs were modelled for a conceptual shaft constructed in the Queenston Formation from Ontario over a stress to strength ratio that has been found at many underground research laboratories. It was found that the excavation of a triangular and thin slot cut-off had the smallest increase in the damage as quantified based on the extensile volumetric strain within the plastic yield zone. Both shapes showed an increase in the damaged zone by 13 cm at a stress to strength ratio of 2.75. The trapezoidal cut-off damage zone increased by 20 cm at the same stress to strength ratio. Further study is required to better define the optimum cut-off shape for wider variety of conditions and other geological formations.

Jan 1, 2020

By G. Hu, P. Chang

"As hydraulic fracturing is forbidden in some countries due to its accompanying environmental pollutions, microwave fracturing was proposed as a substitute technology to enhance coal seam permeability. Its mechanism includes: on the one hand, thermal stress caused by non-uniform microwave heating leads to fracture development in coal body, which increases coal seam permeability for gas; on the other hand, moisture within the coal body that hinders gas desorption and diffusion is removed during this process. To study the effect of moisture content on fracture development under microwave irradiation, high resolution photographs were taken before and after microwave treatment. These photographs were then binarized and skeletonized. Based on the processed photographs, the total area were calculated using ImageJ. The results showed that coal samples with 0% and 25% water saturation have more variation in fracture development as compared to that of samples with 50%, 75% and 100% water saturation. This study can be very useful in gas drainage and gas control in coal mines. INTRODUCTION Coalbed methane (CBM) is a clean energy source as well as a natural hazard to underground mining (Kurnia, Sasmito, & Mujumdar, 2014; Ranathunga, Perera, & Ranjith, 2014). One of the main goals of mine ventilation is diluting methane and removing underground contaminants (Xu, Huang, Nie, Chalmers, & Yang, 2017). Coalbed methane (CBM) has to be extracted prior and during underground mining to avoid methane outburst and explosion accidents. However, with the increase of mining depth, the coal permeability decreases significantly, which makes it hard to extract the methane (Huang, Xu, Hu, Kizil, & Chen, 2018). To enhance the coal seam permeability, many CBM recovery methods have been applied on-site, such as hydraulic fracturing and enhanced coalbed methane (ECBM). However, hydraulic fracturing is restricted due to environmental concerns and ECBM is only suitable for the coal seam with relative large permeability (Bustin, Bustin, Chikatamarla, Downey, & Mansoori, 2016; Huang, Xu, Chen, & Chen, 2018; Huang, Xu, Hu, et al., 2018)."

Jan 1, 2019