The loss of hydrocarbon production caused by the dynamic behavior of the inner boundary and propped fractures under long-term production conditions has been widely reported in recent studies.However,the quantitative r...The loss of hydrocarbon production caused by the dynamic behavior of the inner boundary and propped fractures under long-term production conditions has been widely reported in recent studies.However,the quantitative relationships for the variations of the inner boundary and propped fractures have not been determined and incorporated in the semi-analytical models for the pressure and rate transient analysis.This work focuses on describing the variations of the inner boundary and propped fractures and capturing the typical characteristics from the pressure transient curves.A generalized semi-analytical model was developed to characterize the dynamic behavior of the inner boundary and propped fractures under long-term production conditions.The pressure-dependent length shrinkage coefficients,which quantify the length changes of the inner zone and propped fractures,are modified and incorporated into this multi-zone semi-analytical model.With simultaneous numerical iterations and numerical inversions in Laplace and real-time space,the transient solutions to pressure and rate behavior are determined in just a few seconds.The dynamic behavior of the inner boundary and propped fractures on transient pressure curves is divided into five periods:fracture bilinear flow(FR1),dynamic PFs flow(FR2),inner-area linear flow(FR3),dynamic inner boundary flow(FR4),and outer-area dominated linear flow(FR5).The early hump during FR2 period and a positive upward shift during FR4period are captured on the log-log pressure transient curves,reflecting the dynamic behavior of the inner boundary and propped fractures during the long-term production period.The transient pressure behavior will exhibit greater positive upward trend and the flow rate will be lower with the shrinkage of the inner boundary.The pressure derivative curve will be upward earlier as the inner boundary shrinks more rapidly.The lower permeability caused by the closure of un-propped fractures in the inner zone results in greater upward in pressure derivative curves.If the permeability loss for the dynamic behavior of the inner boundary caused by the closure of un-propped fractures is neglected,the flow rate will be overestimated in the later production period.展开更多
A 30 mm electrothermal-chemical(ETC) gun experimental system is employed to research the burning rate characteristics of 4/7 high-nitrogen solid propellant. Enhanced gas generation rates(EGGR) of propellants during an...A 30 mm electrothermal-chemical(ETC) gun experimental system is employed to research the burning rate characteristics of 4/7 high-nitrogen solid propellant. Enhanced gas generation rates(EGGR) of propellants during and after electrical discharges are verified in the experiments. A modified 0D internal ballistic model is established to simulate the ETC launch. According to the measured pressure and electrical parameters, a transient burning rate law including the influence of EGGR coefficient by electric power and pressure gradient(dp/dt) is added into the model. The EGGR coefficient of 4/7 high-nitrogen solid propellant is equal to 0.005 MW-1. Both simulated breech pressure and projectile muzzle velocity accord with the experimental results well. Compared with Woodley's modified burning rate law, the breech pressure curves acquired by the transient burning rate law are more consistent with test results. Based on the parameters calculated in the model, the relationship among propellant burning rate, pressure gradient(dp/dt) and electric power is analyzed. Depending on the transient burning rate law and experimental data, the burning of solid propellant under the condition of plasma is described more accurately.展开更多
This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of com...This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.展开更多
Although the single-particle model enhanced with electrolyte dynamics(SPMe)is simplified from the pseudo-twodimensional(P2D)electrochemical model for lithium-ion batteries,it is difficult to solve the partial differen...Although the single-particle model enhanced with electrolyte dynamics(SPMe)is simplified from the pseudo-twodimensional(P2D)electrochemical model for lithium-ion batteries,it is difficult to solve the partial differential equations of solid–liquid phases in real-time applications.Moreover,working temperatures have a heavy impact on the battery behavior.Hence,a thermal-coupling SPMe is constructed.Herein,a lumped thermal model is established to estimate battery temperatures.The order of the SPMe model is reduced by using both transfer functions and truncation techniques and merged with Arrhenius equations for thermal effects.The polarization voltage drop is then modified through the use of test data because its original model is unreliable theoretically.Finally,the coupling-model parameters are extracted using genetic algorithms.Experimental results demonstrate that the proposed model produces average errors of about 42 mV under 15 constant current conditions and 15 mV under nine dynamic conditions,respectively.This new electrochemicalthermal coupling model is reliable and expected to be used for onboard applications.展开更多
During the production,the fluid in the vicinity of the directional well enters the wellbore with different rates,leading to non-uniform flux distribution along the directional well.However,in all existing studies,it i...During the production,the fluid in the vicinity of the directional well enters the wellbore with different rates,leading to non-uniform flux distribution along the directional well.However,in all existing studies,it is oversimplified to a uniform flux distribution,which can result in inaccurate results for field applications.Therefore,this paper proposes a semi-analytical model of a directional well based on the assumption of non-uniform flux distribution.Specifically,the direction well is discretized into a carefully chosen series of linear sources,such that the complex well trajectory can be captured and the nonuniform flux distribution along the wellbore can be considered to model the three-dimensional flow behavior.By using the finite difference method,we can obtain the numerical solutions of the transient flow within the wellbore.With the aid of Green's function method,we can obtain the analytical solutions of the transient flow from the matrix to the wellbore.The complete flow behavior of a directional well is perfectly represented by coupling the above two types of transient flow.Subsequently,on the basis of the proposed model,we conduct a comprehensive analysis of the pressure transient behavior of a directional well.The computation results show that the flux variation along the direction well has a significant effect on pressure responses.In addition,the directional well in an infinite reservoir may exhibit the following flow regimes:wellbore afterflow,transition flow,inclined radial flow,elliptical flow,horizontal linear flow,and horizontal radial flow.The horizontal linear flow can be observed only if the formation thickness is much smaller than the well length.Furthermore,a dip region that appears on the pressure derivative curve indicates the three-dimensional flow behavior near the wellbore.展开更多
The paper is focused on the coupling effect between film boiling heat transfer and evaporation drag around a hot-particle in cold liquid. Based on the continuity, momentum and energy equations of the vapor film, a tra...The paper is focused on the coupling effect between film boiling heat transfer and evaporation drag around a hot-particle in cold liquid. Based on the continuity, momentum and energy equations of the vapor film, a transient two-dimensional single particle model has been established. This paper contains a detailed description of HPMC (High-temperature Particle Moving in Coolant) model for studying some aspects of the premixing stage of fuel-coolant interactions (FCIs). The transient process of high-temperature particles moving in coolant can be simu-lated. Comparisons between the experiment results and the calculations using HPMC model demonstrate that HPMC model achieves a good agreement in predicting the time-varying characteristic of high-temperature spheres moving in coolant.展开更多
Accuracy of time-depth conversion in data processing of transient electromagnetic prospecting always affects the accu- rate positioning of water bodies in coal mines. In order to improve the accuracy of time-depth con...Accuracy of time-depth conversion in data processing of transient electromagnetic prospecting always affects the accu- rate positioning of water bodies in coal mines. In order to improve the accuracy of time-depth conversion, we established a mathe- matical model of time-depth conversion for a transient electromagnetic method based on the theory of "double smoke ring effect" of full space transient electromagnetic field transmission. Using a 3-layer as well as a 4-layer geo-electric model for roadway floors, we performed the time-depth conversion of theoretical curves of apparent resistance varying over time. In these curves, the depth corresponding to extreme value points is nearly the same as the depth of a geo-electric model. The position of water body deter- mined by our time-depth conversion method agrees well with the result of borehole drilling, indicating that the established time-depth conversion model can clearly improve the accuracy of spatial positioning of water bodies in coal mines.展开更多
Severe well interference through complex fracture networks(CFNs)can be observed among multi-well pads in low permeability reservoirs.The well interference analysis between multi-fractured horizontal wells(MFHWs)is vit...Severe well interference through complex fracture networks(CFNs)can be observed among multi-well pads in low permeability reservoirs.The well interference analysis between multi-fractured horizontal wells(MFHWs)is vitally important for reservoir effective development.Well interference has been historically investigated by pressure transient analysis,while it has shown that rate transient analysis has great potential in well interference diagnosis.However,the impact of complex fracture networks(CFNs)on rate transient behavior of parent well and child well in unconventional reservoirs is still not clear.To further investigate,this paper develops an integrated approach combining pressure and rate transient analysis for well interference diagnosis considering CFNs.To perform multi-well simulation considering CFNs,non-intrusive embedded discrete fracture model approach was applied for coupling fracture with reservoir models.The impact of CFN including natural fractures and frac-hits on pressure and rate transient behavior in multi-well system was investigated.On a logelog plot,interference flow and compound linear flow are two new flow regimes caused by nearby producers.When both NFs and frac-hits are present in the reservoir,frac-hits have a greater impact on well#1 which contains frac-hits,and NFs have greater impact on well#3 which does not have frac-hits.For all well producing circumstances,it might be challenging to see divergence during pseudosteady state flow brought on by frac-hits on the logelog plot.Besides,when NFs occur,reservoir depletion becomes noticeable in comparison to frac-hits in pressure distribution.Application of this integrated approach demonstrates that it works well to characterize the well interference among different multi-fractured horizontal wells in a well pad.Better reservoir evaluation can be acquired based on the new features observed in the novel model,demonstrating the practicability of the proposed approach.The findings of this study can help for better evaluating well interference degree in multi-well systems combing PTA and RTA,which can reduce the uncertainty and improve the accuracy of the well interference analysis based on both field pressure and rate data.展开更多
基金financial funding of National Natural Science Foundation of China (No.52004307)China National Petroleum Corporation (No.ZLZX2020-02-04)the Science Foundation of China University of Petroleum,Beijing (No.2462018YJRC015)。
文摘The loss of hydrocarbon production caused by the dynamic behavior of the inner boundary and propped fractures under long-term production conditions has been widely reported in recent studies.However,the quantitative relationships for the variations of the inner boundary and propped fractures have not been determined and incorporated in the semi-analytical models for the pressure and rate transient analysis.This work focuses on describing the variations of the inner boundary and propped fractures and capturing the typical characteristics from the pressure transient curves.A generalized semi-analytical model was developed to characterize the dynamic behavior of the inner boundary and propped fractures under long-term production conditions.The pressure-dependent length shrinkage coefficients,which quantify the length changes of the inner zone and propped fractures,are modified and incorporated into this multi-zone semi-analytical model.With simultaneous numerical iterations and numerical inversions in Laplace and real-time space,the transient solutions to pressure and rate behavior are determined in just a few seconds.The dynamic behavior of the inner boundary and propped fractures on transient pressure curves is divided into five periods:fracture bilinear flow(FR1),dynamic PFs flow(FR2),inner-area linear flow(FR3),dynamic inner boundary flow(FR4),and outer-area dominated linear flow(FR5).The early hump during FR2 period and a positive upward shift during FR4period are captured on the log-log pressure transient curves,reflecting the dynamic behavior of the inner boundary and propped fractures during the long-term production period.The transient pressure behavior will exhibit greater positive upward trend and the flow rate will be lower with the shrinkage of the inner boundary.The pressure derivative curve will be upward earlier as the inner boundary shrinks more rapidly.The lower permeability caused by the closure of un-propped fractures in the inner zone results in greater upward in pressure derivative curves.If the permeability loss for the dynamic behavior of the inner boundary caused by the closure of un-propped fractures is neglected,the flow rate will be overestimated in the later production period.
文摘A 30 mm electrothermal-chemical(ETC) gun experimental system is employed to research the burning rate characteristics of 4/7 high-nitrogen solid propellant. Enhanced gas generation rates(EGGR) of propellants during and after electrical discharges are verified in the experiments. A modified 0D internal ballistic model is established to simulate the ETC launch. According to the measured pressure and electrical parameters, a transient burning rate law including the influence of EGGR coefficient by electric power and pressure gradient(dp/dt) is added into the model. The EGGR coefficient of 4/7 high-nitrogen solid propellant is equal to 0.005 MW-1. Both simulated breech pressure and projectile muzzle velocity accord with the experimental results well. Compared with Woodley's modified burning rate law, the breech pressure curves acquired by the transient burning rate law are more consistent with test results. Based on the parameters calculated in the model, the relationship among propellant burning rate, pressure gradient(dp/dt) and electric power is analyzed. Depending on the transient burning rate law and experimental data, the burning of solid propellant under the condition of plasma is described more accurately.
基金This work was supported by China Postdoctoral Science Foundation(No.2022M723391)the Science and Technology Innovation Project of Higher Education in Shanxi Province(No.2019L0754)+1 种基金the Central Guiding Local Science and Technology Development Fund Project(No.YDZJSX2021B021)Shanxi Province Basic Research Plan General Project(No.202203021221294).
文摘This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.
基金the financial support from the National Key Research and Development Program of China(Grant No.2021YFF0601101)。
文摘Although the single-particle model enhanced with electrolyte dynamics(SPMe)is simplified from the pseudo-twodimensional(P2D)electrochemical model for lithium-ion batteries,it is difficult to solve the partial differential equations of solid–liquid phases in real-time applications.Moreover,working temperatures have a heavy impact on the battery behavior.Hence,a thermal-coupling SPMe is constructed.Herein,a lumped thermal model is established to estimate battery temperatures.The order of the SPMe model is reduced by using both transfer functions and truncation techniques and merged with Arrhenius equations for thermal effects.The polarization voltage drop is then modified through the use of test data because its original model is unreliable theoretically.Finally,the coupling-model parameters are extracted using genetic algorithms.Experimental results demonstrate that the proposed model produces average errors of about 42 mV under 15 constant current conditions and 15 mV under nine dynamic conditions,respectively.This new electrochemicalthermal coupling model is reliable and expected to be used for onboard applications.
基金the financial support provided by the National Natural Science Foundation of China(No.52104043)。
文摘During the production,the fluid in the vicinity of the directional well enters the wellbore with different rates,leading to non-uniform flux distribution along the directional well.However,in all existing studies,it is oversimplified to a uniform flux distribution,which can result in inaccurate results for field applications.Therefore,this paper proposes a semi-analytical model of a directional well based on the assumption of non-uniform flux distribution.Specifically,the direction well is discretized into a carefully chosen series of linear sources,such that the complex well trajectory can be captured and the nonuniform flux distribution along the wellbore can be considered to model the three-dimensional flow behavior.By using the finite difference method,we can obtain the numerical solutions of the transient flow within the wellbore.With the aid of Green's function method,we can obtain the analytical solutions of the transient flow from the matrix to the wellbore.The complete flow behavior of a directional well is perfectly represented by coupling the above two types of transient flow.Subsequently,on the basis of the proposed model,we conduct a comprehensive analysis of the pressure transient behavior of a directional well.The computation results show that the flux variation along the direction well has a significant effect on pressure responses.In addition,the directional well in an infinite reservoir may exhibit the following flow regimes:wellbore afterflow,transition flow,inclined radial flow,elliptical flow,horizontal linear flow,and horizontal radial flow.The horizontal linear flow can be observed only if the formation thickness is much smaller than the well length.Furthermore,a dip region that appears on the pressure derivative curve indicates the three-dimensional flow behavior near the wellbore.
基金Supported partially by National Natural Science Fundation of China (No.50046026 and No.50376036) and 985 Engineering Fund of Shanghai Jiaotong University.
文摘The paper is focused on the coupling effect between film boiling heat transfer and evaporation drag around a hot-particle in cold liquid. Based on the continuity, momentum and energy equations of the vapor film, a transient two-dimensional single particle model has been established. This paper contains a detailed description of HPMC (High-temperature Particle Moving in Coolant) model for studying some aspects of the premixing stage of fuel-coolant interactions (FCIs). The transient process of high-temperature particles moving in coolant can be simu-lated. Comparisons between the experiment results and the calculations using HPMC model demonstrate that HPMC model achieves a good agreement in predicting the time-varying characteristic of high-temperature spheres moving in coolant.
基金Project 2007CB209400 supported by the National Basic Research Program of China
文摘Accuracy of time-depth conversion in data processing of transient electromagnetic prospecting always affects the accu- rate positioning of water bodies in coal mines. In order to improve the accuracy of time-depth conversion, we established a mathe- matical model of time-depth conversion for a transient electromagnetic method based on the theory of "double smoke ring effect" of full space transient electromagnetic field transmission. Using a 3-layer as well as a 4-layer geo-electric model for roadway floors, we performed the time-depth conversion of theoretical curves of apparent resistance varying over time. In these curves, the depth corresponding to extreme value points is nearly the same as the depth of a geo-electric model. The position of water body deter- mined by our time-depth conversion method agrees well with the result of borehole drilling, indicating that the established time-depth conversion model can clearly improve the accuracy of spatial positioning of water bodies in coal mines.
基金The authors are grateful to the financial support from China Postdoctoral Science Foundation(2022M712645)Opening Fund of Key Laboratory of Enhanced Oil Recovery(Northeast Petroleum University),Ministry of Education(NEPU-EOR-2021-03).
文摘Severe well interference through complex fracture networks(CFNs)can be observed among multi-well pads in low permeability reservoirs.The well interference analysis between multi-fractured horizontal wells(MFHWs)is vitally important for reservoir effective development.Well interference has been historically investigated by pressure transient analysis,while it has shown that rate transient analysis has great potential in well interference diagnosis.However,the impact of complex fracture networks(CFNs)on rate transient behavior of parent well and child well in unconventional reservoirs is still not clear.To further investigate,this paper develops an integrated approach combining pressure and rate transient analysis for well interference diagnosis considering CFNs.To perform multi-well simulation considering CFNs,non-intrusive embedded discrete fracture model approach was applied for coupling fracture with reservoir models.The impact of CFN including natural fractures and frac-hits on pressure and rate transient behavior in multi-well system was investigated.On a logelog plot,interference flow and compound linear flow are two new flow regimes caused by nearby producers.When both NFs and frac-hits are present in the reservoir,frac-hits have a greater impact on well#1 which contains frac-hits,and NFs have greater impact on well#3 which does not have frac-hits.For all well producing circumstances,it might be challenging to see divergence during pseudosteady state flow brought on by frac-hits on the logelog plot.Besides,when NFs occur,reservoir depletion becomes noticeable in comparison to frac-hits in pressure distribution.Application of this integrated approach demonstrates that it works well to characterize the well interference among different multi-fractured horizontal wells in a well pad.Better reservoir evaluation can be acquired based on the new features observed in the novel model,demonstrating the practicability of the proposed approach.The findings of this study can help for better evaluating well interference degree in multi-well systems combing PTA and RTA,which can reduce the uncertainty and improve the accuracy of the well interference analysis based on both field pressure and rate data.
