Self-oscillating systems abound in the natural world and offer substantial potential for applications in controllers,micro-motors,medical equipments,and so on.Currently,numerical methods have been widely utilized for ...Self-oscillating systems abound in the natural world and offer substantial potential for applications in controllers,micro-motors,medical equipments,and so on.Currently,numerical methods have been widely utilized for obtaining the characteristics of self-oscillation including amplitude and frequency.However,numerical methods are burdened by intricate computations and limited precision,hindering comprehensive investigations into self-oscillating systems.In this paper,the stability of a liquid crystal elastomer fiber self-oscillating system under a linear temperature field is studied,and analytical solutions for the amplitude and frequency are determined.Initially,we establish the governing equations of self-oscillation,elucidate two motion regimes,and reveal the underlying mechanism.Subsequently,we conduct a stability analysis and employ a multi-scale method to obtain the analytical solutions for the amplitude and frequency.The results show agreement between the multi-scale and numerical methods.This research contributes to the examination of diverse self-oscillating systems and advances the theoretical analysis of self-oscillating systems rooted in active materials.展开更多
The face stability problem is a major concern for tunnels excavated in rock masses governed by the Hoek-Brown strength criterion.To provide an accurate prediction for the theoretical solution of the critical face pres...The face stability problem is a major concern for tunnels excavated in rock masses governed by the Hoek-Brown strength criterion.To provide an accurate prediction for the theoretical solution of the critical face pressure,this study adopts the piecewise linear method(PLM)to account for the nonlinearity of the strength envelope and proposes a new multi-horn rotational mechanism based on the Hoek-Brown strength criterion and the associative flow rule.The analytical solution of critical support pressure is derived from the energy-work balance equation in the framework of the plastic limit theorem;it is formulated as a multivariable nonlinear optimization problem relying on 2m dependent variables(m is the number of segments).Meanwhile,two classic linearized measures,the generalized tangential technique(GTT)and equivalent Mohr-Coulomb parameters method(EMM),are incorporated into the analysis for comparison.Surprisingly,the parametric study indicates a significant improvement in support pressure by up to 13%compared with the GTT,and as expected,the stability of the tunnel face is greatly influenced by the rock strength parameters.The stress distribution on the rupture surface is calculated to gain an intuitive understanding of the failure at the limit state.Although the limit analysis is incapable of calculating the true stress distribution in rock masses,a rough approximation of the stress vector on the rupture surface is permitted.In the end,sets of normalized face pressure are provided in the form of charts for a quick assessment of face stability in rock masses.展开更多
The problem of linear time-varying(LTV) system modal analysis is considered based on time-dependent state space representations, as classical modal analysis of linear time-invariant systems and current LTV system mo...The problem of linear time-varying(LTV) system modal analysis is considered based on time-dependent state space representations, as classical modal analysis of linear time-invariant systems and current LTV system modal analysis under the "frozen-time" assumption are not able to determine the dynamic stability of LTV systems. Time-dependent state space representations of LTV systems are first introduced, and the corresponding modal analysis theories are subsequently presented via a stabilitypreserving state transformation. The time-varying modes of LTV systems are extended in terms of uniqueness, and are further interpreted to determine the system's stability. An extended modal identification is proposed to estimate the time-varying modes, consisting of the estimation of the state transition matrix via a subspace-based method and the extraction of the time-varying modes by the QR decomposition. The proposed approach is numerically validated by three numerical cases, and is experimentally validated by a coupled moving-mass simply supported beam exper- imental case. The proposed approach is capable of accurately estimating the time-varying modes, and provides anew way to determine the dynamic stability of LTV systems by using the estimated time-varying modes.展开更多
This paper focuses on the methodology analysis for the stability and the corresponding tracking performance of a closed-loop digital jump linear control system with a stochastic switching signal. The method is applied...This paper focuses on the methodology analysis for the stability and the corresponding tracking performance of a closed-loop digital jump linear control system with a stochastic switching signal. The method is applied to a flight control system. A distributed recoverable platform is implemented on the flight control system and subject to independent digital upsets. The upset processes are used to stimulate electromagnetic environments. Specifically, the paper presents the scenarios that the upset process is directly injected into the distributed flight control system, which is modeled by independent Markov upset processes and independent and identically distributed (IID) processes. A theoretical performance analysis and simulation modelling are both presented in detail for a more complete independent digital upset injection. The specific examples are proposed to verify the methodology of tracking performance analysis. The general analyses for different configurations are also proposed. Comparisons among different configurations are conducted to demonstrate the availability and the characteristics of the design.展开更多
Linear Stability Analysis(LSA)of parallel shear flows,via local and global approaches,is presented.The local analysis is carried out by solving the Orr-Sommerfeld(OS)equation using a spectral-collocation method based ...Linear Stability Analysis(LSA)of parallel shear flows,via local and global approaches,is presented.The local analysis is carried out by solving the Orr-Sommerfeld(OS)equation using a spectral-collocation method based on Chebyshev polynomials.A stabilized finite element formulation is employed to carry out the global analysis using the linearized disturbance equations in primitive variables.The local and global analysis are compared.As per the Squires theorem,the two-dimensional disturbance has the largest growth rate.Therefore,only two-dimensional disturbances are considered.By its very nature,the local analysis assumes the disturbance field to be spatially periodic in the streamwise direction.The global analysis permits a more general disturbance.However,to enable a comparison with the local analysis,periodic boundary conditions,at the inlet and exit of the domain,are imposed on the disturbance.Computations are carried out for the LSA of the Plane Poiseuille Flow(PPF).The relationship between the wavenumber,a,of the disturbance and the streamwise extent of the domain,L,in the global analysis is explored for Re=7000.It is found that a and L are related by L=2pn/a,where n is the number of cells of the instability along the streamwise direction within the domain length,L.The procedure to interpret the results from the global analysis,for comparison with local analysis,is described.展开更多
Based on the hydrodynamic stability theory of distorted laminar flow and the kind of distortion profiles on the mean velocity in parallel shear flow given in paper [1], this paper investigates the linear stability beh...Based on the hydrodynamic stability theory of distorted laminar flow and the kind of distortion profiles on the mean velocity in parallel shear flow given in paper [1], this paper investigates the linear stability behaviour of parallel shear flow, presents unstable results of plane Couette flow and pipe Poiseuille flow to two-dimensional or axisymmetric disturbances for the first time, and obtains neutral curves of these two motions under certain definition.展开更多
In this paper,a linear/nonlinear switching active disturbance rejection control(SADRC)based decoupling control approach is proposed to deal with some difficult control problems in a class of multi-input multi-output(M...In this paper,a linear/nonlinear switching active disturbance rejection control(SADRC)based decoupling control approach is proposed to deal with some difficult control problems in a class of multi-input multi-output(MIMO)systems such as multi-variables,disturbances,and coupling,etc.Firstly,the structure and parameter tuning method of SADRC is introduced into this paper.Followed on this,virtual control variables are adopted into the MIMO systems,making the systems decoupled.Then the SADRC controller is designed for every subsystem.After this,a stability analyzed method via the Lyapunov function is proposed for the whole system.Finally,some simulations are presented to demonstrate the anti-disturbance and robustness of SADRC,and results show SADRC has a potential applications in engineering practice.展开更多
An analysis method based on the fuzzy Lyapunov functions is presented to analyze the stability of the continuous affine fuzzy systems. First, a method is introduced to deal with the consequent part of the fuzzy local ...An analysis method based on the fuzzy Lyapunov functions is presented to analyze the stability of the continuous affine fuzzy systems. First, a method is introduced to deal with the consequent part of the fuzzy local model. Thus, the stability analysis method of the homogeneous fuzzy system can be used for reference. Stability conditions are derived in terms of linear matrix inequalities based on the fuzzy Lyapunov functions and the modified common Lyapunov functions, respectively. The results demonstrate that the stability result based on the fuzzy Lyapunov functions is less conservative than that based on the modified common Lyapunov functions via numerical examples. Compared with the method which does not expand the consequent part, the proposed method is simpler but its feasible region is reduced. Finally, in order to expand the application of the fuzzy Lyapunov functions, the piecewise fuzzy Lyapunov function is proposed, which can be used to analyze the stability for triangular or trapezoidal membership functions and obtain the stability conditions. A numerical example validates the effectiveness of the proposed approach.展开更多
The aim of this paper is to study the determination of the stability regions for continuous-time systems subject to actuator saturation. Using an affine saturation-dependent Lyapunov function, a new method is proposed...The aim of this paper is to study the determination of the stability regions for continuous-time systems subject to actuator saturation. Using an affine saturation-dependent Lyapunov function, a new method is proposed to obtain the estimation of the domain of attraction of the closed-loop system. A family of linear matrix inequalities (LMIs) that provides sufficient conditions for the existence of this type of Lyapunov function are presented. The results obtained in this paper can reduce the conservativeness compared with the existing ones. Numerical examples are given to illustrate the effectiveness of the proposed results.展开更多
This paper uses the Oseen transformation to solve the differential equations governing motion of the vertical linear gradient flow distribution close to a wall surface. The Navier-Stokes equations are used to consider...This paper uses the Oseen transformation to solve the differential equations governing motion of the vertical linear gradient flow distribution close to a wall surface. The Navier-Stokes equations are used to consider the inertia term along the flow direction. A novel contour integral method is used to solve the complex Airy function. The boundary conditions of linear gradient flow distribution for finite problems are determined. The vorticity function, the pressure function, and the turbulent velocity profiles are provided, and the stability of particle trajectories is studied. An Lx-function form of the third derivative circulation is used to to simplify the solution. Theoretical results are compared with the experimental measurements with satisfactory agreement.展开更多
The issue of the stability and controller design of Takagi-Sugeno(T-S) fuzzy control systems with time-delay is investigated under imperfect premise matching when the T-S fuzzy time-delay model and fuzzy controller ...The issue of the stability and controller design of Takagi-Sugeno(T-S) fuzzy control systems with time-delay is investigated under imperfect premise matching when the T-S fuzzy time-delay model and fuzzy controller do not share the same membership functions.A new stability criterion which contains the information of membership functions is derived.The new stability criterion is less conservative,and enhances the design flexibility.Two numerical examples are presented to illustrate the conservativeness and effectiveness of the proposed method.展开更多
为了提升永磁直线同步电机(permanent magnet linear synchronous motor,PMLSM)在负载变化、参数摄动和其他不确定因素下的抗扰性能和速度跟踪性能,提出一种基于级联线性-非线性自抗扰控制器的PMLSM速度控制策略。首先,建立考虑负载扰...为了提升永磁直线同步电机(permanent magnet linear synchronous motor,PMLSM)在负载变化、参数摄动和其他不确定因素下的抗扰性能和速度跟踪性能,提出一种基于级联线性-非线性自抗扰控制器的PMLSM速度控制策略。首先,建立考虑负载扰动和参数失配的PMLSM数学模型;其次,设计级联线性-非线性扩张状态观测器来实时估计和补偿系统所受的不确定扰动,前级线性扩张状态观测器保证系统在大扰动下保持稳定,后级非线性扩张状态观测器利用非线性机制进一步提高系统对扰动的估计精度,从而将线性自抗扰控制和非线性自抗扰控制的优势相结合,以此提升系统的速度跟踪性能和抗扰动能力;并且,对所提控制器提出基于劳斯判据的稳定性分析方法,并对系统的抗扰性能和噪声抑制性能进行了频域分析;最后,对基于PI控制、级联线性自抗扰控制、非线性自抗扰控制和级联线性-非线性自抗扰控制的永磁直线同步电机系统进行仿真和实验对比,验证所提方法的优越性。展开更多
Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operation...Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.展开更多
基金Project supported by the National Natural Science Foundation of China (No.12172001)the Anhui Provincial Natural Science Foundation of China (No.2208085Y01)+1 种基金the University Natural Science Research Project of Anhui Province of China (No.2022AH020029)the Housing and Urban-Rural Development Science and Technology Project of Anhui Province of China (No.2023-YF129)。
文摘Self-oscillating systems abound in the natural world and offer substantial potential for applications in controllers,micro-motors,medical equipments,and so on.Currently,numerical methods have been widely utilized for obtaining the characteristics of self-oscillation including amplitude and frequency.However,numerical methods are burdened by intricate computations and limited precision,hindering comprehensive investigations into self-oscillating systems.In this paper,the stability of a liquid crystal elastomer fiber self-oscillating system under a linear temperature field is studied,and analytical solutions for the amplitude and frequency are determined.Initially,we establish the governing equations of self-oscillation,elucidate two motion regimes,and reveal the underlying mechanism.Subsequently,we conduct a stability analysis and employ a multi-scale method to obtain the analytical solutions for the amplitude and frequency.The results show agreement between the multi-scale and numerical methods.This research contributes to the examination of diverse self-oscillating systems and advances the theoretical analysis of self-oscillating systems rooted in active materials.
基金supported by Fundamental Research Funds for the central universities of Central South University(No.2022ZZTS0153).
文摘The face stability problem is a major concern for tunnels excavated in rock masses governed by the Hoek-Brown strength criterion.To provide an accurate prediction for the theoretical solution of the critical face pressure,this study adopts the piecewise linear method(PLM)to account for the nonlinearity of the strength envelope and proposes a new multi-horn rotational mechanism based on the Hoek-Brown strength criterion and the associative flow rule.The analytical solution of critical support pressure is derived from the energy-work balance equation in the framework of the plastic limit theorem;it is formulated as a multivariable nonlinear optimization problem relying on 2m dependent variables(m is the number of segments).Meanwhile,two classic linearized measures,the generalized tangential technique(GTT)and equivalent Mohr-Coulomb parameters method(EMM),are incorporated into the analysis for comparison.Surprisingly,the parametric study indicates a significant improvement in support pressure by up to 13%compared with the GTT,and as expected,the stability of the tunnel face is greatly influenced by the rock strength parameters.The stress distribution on the rupture surface is calculated to gain an intuitive understanding of the failure at the limit state.Although the limit analysis is incapable of calculating the true stress distribution in rock masses,a rough approximation of the stress vector on the rupture surface is permitted.In the end,sets of normalized face pressure are provided in the form of charts for a quick assessment of face stability in rock masses.
基金Supported by the China Scholarship Council,National Natural Science Foundation of China(Grant No.11402022)the Interuniversity Attraction Poles Programme of the Belgian Science Policy Office(DYSCO)+1 种基金the Fund for Scientific Research–Flanders(FWO)the Research Fund KU Leuven
文摘The problem of linear time-varying(LTV) system modal analysis is considered based on time-dependent state space representations, as classical modal analysis of linear time-invariant systems and current LTV system modal analysis under the "frozen-time" assumption are not able to determine the dynamic stability of LTV systems. Time-dependent state space representations of LTV systems are first introduced, and the corresponding modal analysis theories are subsequently presented via a stabilitypreserving state transformation. The time-varying modes of LTV systems are extended in terms of uniqueness, and are further interpreted to determine the system's stability. An extended modal identification is proposed to estimate the time-varying modes, consisting of the estimation of the state transition matrix via a subspace-based method and the extraction of the time-varying modes by the QR decomposition. The proposed approach is numerically validated by three numerical cases, and is experimentally validated by a coupled moving-mass simply supported beam exper- imental case. The proposed approach is capable of accurately estimating the time-varying modes, and provides anew way to determine the dynamic stability of LTV systems by using the estimated time-varying modes.
基金Project supported by the Young Scientists Fund of the National Natural Science Foundation of China(Grant No.61403395)the Natural Science Foundation of Tianjin,China(Grant No.13JCYBJC39000)+2 种基金the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education Ministry,Chinathe Tianjin Key Laboratory of Civil Aircraft Airworthiness and Maintenance in Civil Aviation of China(Grant No.104003020106)the Fund for Scholars of Civil Aviation University of China(Grant No.2012QD21x)
文摘This paper focuses on the methodology analysis for the stability and the corresponding tracking performance of a closed-loop digital jump linear control system with a stochastic switching signal. The method is applied to a flight control system. A distributed recoverable platform is implemented on the flight control system and subject to independent digital upsets. The upset processes are used to stimulate electromagnetic environments. Specifically, the paper presents the scenarios that the upset process is directly injected into the distributed flight control system, which is modeled by independent Markov upset processes and independent and identically distributed (IID) processes. A theoretical performance analysis and simulation modelling are both presented in detail for a more complete independent digital upset injection. The specific examples are proposed to verify the methodology of tracking performance analysis. The general analyses for different configurations are also proposed. Comparisons among different configurations are conducted to demonstrate the availability and the characteristics of the design.
文摘Linear Stability Analysis(LSA)of parallel shear flows,via local and global approaches,is presented.The local analysis is carried out by solving the Orr-Sommerfeld(OS)equation using a spectral-collocation method based on Chebyshev polynomials.A stabilized finite element formulation is employed to carry out the global analysis using the linearized disturbance equations in primitive variables.The local and global analysis are compared.As per the Squires theorem,the two-dimensional disturbance has the largest growth rate.Therefore,only two-dimensional disturbances are considered.By its very nature,the local analysis assumes the disturbance field to be spatially periodic in the streamwise direction.The global analysis permits a more general disturbance.However,to enable a comparison with the local analysis,periodic boundary conditions,at the inlet and exit of the domain,are imposed on the disturbance.Computations are carried out for the LSA of the Plane Poiseuille Flow(PPF).The relationship between the wavenumber,a,of the disturbance and the streamwise extent of the domain,L,in the global analysis is explored for Re=7000.It is found that a and L are related by L=2pn/a,where n is the number of cells of the instability along the streamwise direction within the domain length,L.The procedure to interpret the results from the global analysis,for comparison with local analysis,is described.
文摘Based on the hydrodynamic stability theory of distorted laminar flow and the kind of distortion profiles on the mean velocity in parallel shear flow given in paper [1], this paper investigates the linear stability behaviour of parallel shear flow, presents unstable results of plane Couette flow and pipe Poiseuille flow to two-dimensional or axisymmetric disturbances for the first time, and obtains neutral curves of these two motions under certain definition.
基金supported by the Scientific Research Innovation Development Foundation of Army Engineering University((2019)71).
文摘In this paper,a linear/nonlinear switching active disturbance rejection control(SADRC)based decoupling control approach is proposed to deal with some difficult control problems in a class of multi-input multi-output(MIMO)systems such as multi-variables,disturbances,and coupling,etc.Firstly,the structure and parameter tuning method of SADRC is introduced into this paper.Followed on this,virtual control variables are adopted into the MIMO systems,making the systems decoupled.Then the SADRC controller is designed for every subsystem.After this,a stability analyzed method via the Lyapunov function is proposed for the whole system.Finally,some simulations are presented to demonstrate the anti-disturbance and robustness of SADRC,and results show SADRC has a potential applications in engineering practice.
基金Specialized Research Fund for the Doctoral Program of Higher Education ( No. 20090092110051)the Key Project of Chinese Ministry of Education ( No. 108060)the National Natural Science Foundation of China ( No. 51076027, 51036002, 51106024)
文摘An analysis method based on the fuzzy Lyapunov functions is presented to analyze the stability of the continuous affine fuzzy systems. First, a method is introduced to deal with the consequent part of the fuzzy local model. Thus, the stability analysis method of the homogeneous fuzzy system can be used for reference. Stability conditions are derived in terms of linear matrix inequalities based on the fuzzy Lyapunov functions and the modified common Lyapunov functions, respectively. The results demonstrate that the stability result based on the fuzzy Lyapunov functions is less conservative than that based on the modified common Lyapunov functions via numerical examples. Compared with the method which does not expand the consequent part, the proposed method is simpler but its feasible region is reduced. Finally, in order to expand the application of the fuzzy Lyapunov functions, the piecewise fuzzy Lyapunov function is proposed, which can be used to analyze the stability for triangular or trapezoidal membership functions and obtain the stability conditions. A numerical example validates the effectiveness of the proposed approach.
基金supported by the National Creative Research Groups Science Foundation of China (No.60721062)the National High Technology Research and Development Program of China (863 Program) (No.2006AA04 Z182)National Natural Science Foundation of China (No.60736021)
文摘The aim of this paper is to study the determination of the stability regions for continuous-time systems subject to actuator saturation. Using an affine saturation-dependent Lyapunov function, a new method is proposed to obtain the estimation of the domain of attraction of the closed-loop system. A family of linear matrix inequalities (LMIs) that provides sufficient conditions for the existence of this type of Lyapunov function are presented. The results obtained in this paper can reduce the conservativeness compared with the existing ones. Numerical examples are given to illustrate the effectiveness of the proposed results.
基金Project supported by the National Natural Science foundation of China(No.51079095)the Science Fund for Creative Research Groups of the National Natural Science Foundation of China(No.51021004)
文摘This paper uses the Oseen transformation to solve the differential equations governing motion of the vertical linear gradient flow distribution close to a wall surface. The Navier-Stokes equations are used to consider the inertia term along the flow direction. A novel contour integral method is used to solve the complex Airy function. The boundary conditions of linear gradient flow distribution for finite problems are determined. The vorticity function, the pressure function, and the turbulent velocity profiles are provided, and the stability of particle trajectories is studied. An Lx-function form of the third derivative circulation is used to to simplify the solution. Theoretical results are compared with the experimental measurements with satisfactory agreement.
基金Supported by the National Natural Science Foundation of China(60874084)the Academy of Finland(135225,127299)
文摘The issue of the stability and controller design of Takagi-Sugeno(T-S) fuzzy control systems with time-delay is investigated under imperfect premise matching when the T-S fuzzy time-delay model and fuzzy controller do not share the same membership functions.A new stability criterion which contains the information of membership functions is derived.The new stability criterion is less conservative,and enhances the design flexibility.Two numerical examples are presented to illustrate the conservativeness and effectiveness of the proposed method.
文摘为了提升永磁直线同步电机(permanent magnet linear synchronous motor,PMLSM)在负载变化、参数摄动和其他不确定因素下的抗扰性能和速度跟踪性能,提出一种基于级联线性-非线性自抗扰控制器的PMLSM速度控制策略。首先,建立考虑负载扰动和参数失配的PMLSM数学模型;其次,设计级联线性-非线性扩张状态观测器来实时估计和补偿系统所受的不确定扰动,前级线性扩张状态观测器保证系统在大扰动下保持稳定,后级非线性扩张状态观测器利用非线性机制进一步提高系统对扰动的估计精度,从而将线性自抗扰控制和非线性自抗扰控制的优势相结合,以此提升系统的速度跟踪性能和抗扰动能力;并且,对所提控制器提出基于劳斯判据的稳定性分析方法,并对系统的抗扰性能和噪声抑制性能进行了频域分析;最后,对基于PI控制、级联线性自抗扰控制、非线性自抗扰控制和级联线性-非线性自抗扰控制的永磁直线同步电机系统进行仿真和实验对比,验证所提方法的优越性。
基金funded by the joint fund of the National Key Research and Development Program of China(No.2021YFC2902101)National Natural Science Foundation of China(Grant No.52374084)+1 种基金Open Foundation of National Energy shale gas R&D(experiment) center(2022-KFKT-12)the 111 Project(B17009)。
文摘Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.