It is known that there is a discrepancy between field data and the results predicted from the previous equations derived by simplifying three-dimensional（3-D） flow into two-dimensions（2-D）.This paper presents a ne...It is known that there is a discrepancy between field data and the results predicted from the previous equations derived by simplifying three-dimensional（3-D） flow into two-dimensions（2-D）.This paper presents a new steady-state productivity equation for horizontal wells in bottom water drive gas reservoirs.Firstly,the fundamental solution to the 3-D steady-state Laplace equation is derived with the philosophy of source and the Green function for a horizontal well located at the center of the laterally infinite gas reservoir.Then,using the fundamental solution and the Simpson integral formula,the average pseudo-pressure equation and the steady-state productivity equation are achieved for the horizontal section.Two case-studies are given in the paper,the results calculated from the newly-derived formula are very close to the numerical simulation performed with the Canadian software CMG and the real production data,indicating that the new formula can be used to predict the steady-state productivity of such horizontal gas wells.展开更多
This article establishes the precise asymptotics Eu^m(t, x)(t → ∞ or m → ∞) for the stochastic heat equation ?u/?t(t, x) =1/2?u(t, x) + u(t, x)(t, x)?W/?t(t, x) with the time-derivative Gaussian noise W?/?t(t, x) ...This article establishes the precise asymptotics Eu^m(t, x)(t → ∞ or m → ∞) for the stochastic heat equation ?u/?t(t, x) =1/2?u(t, x) + u(t, x)(t, x)?W/?t(t, x) with the time-derivative Gaussian noise W?/?t(t, x) that is fractional in time and homogeneous in space.展开更多
Based on the measurement mechanism of mobility in pressure measurement while drilling, through analyzing a large number of mobility data, it is found that under the condition of water-based mud drilling, the product o...Based on the measurement mechanism of mobility in pressure measurement while drilling, through analyzing a large number of mobility data, it is found that under the condition of water-based mud drilling, the product of mobility from pressure measurement while drilling and the viscosity of mud filtrate is infinitely close to the water phase permeability under the residual oil in relative permeability experiment. Based on this, a method converting the mobility from pressure measurement while drilling to core permeability is proposed, and the permeability based on Timur formula has been established. Application of this method in Penglai 19-9 oilfield of Bohai Sea shows:(1) Compared with the permeability calculated by the model of adjacent oilfields, the permeability calculated by this model is more consistent with the permeability calculated by core analysis.(2) Based on the new model, the correlation between the calculated mobility of well logging and the actual drilling specific productivity index bas been established. Compared with the relationship established by using the permeability model of an adjacent oilfield, the correlation of the new model is better.(3) Productivity of four directional wells was predicted, and the prediction results are in good agreement with the actual production after drilling.展开更多
To transform the exponential traveling wave solutions to bilinear differential equations, a sufficient and necessary condition is proposed. Motivated by the condition, we extend the results to the(2+1)-dimensional Kad...To transform the exponential traveling wave solutions to bilinear differential equations, a sufficient and necessary condition is proposed. Motivated by the condition, we extend the results to the(2+1)-dimensional Kadomtsev–Petviashvili(KP) equation, the(3+1)-dimensional generalized Kadomtsev–Petviashvili(g-KP) equation, and the B-type Kadomtsev–Petviashvili(BKP) equation. Aa a result, we obtain some new resonant multiple wave solutions through the parameterization for wave numbers and frequencies via some linear combinations of exponential traveling waves. Finally, these new resonant type solutions can be displayed in graphs to illustrate the resonant behaviors of multiple wave solutions.展开更多
The fractional Feynman-Kac equations describe the distributions of functionals of non-Brownian motion, or anomalous diffusion, including two types called the forward and backward fractional Feynman-Kac equations, wher...The fractional Feynman-Kac equations describe the distributions of functionals of non-Brownian motion, or anomalous diffusion, including two types called the forward and backward fractional Feynman-Kac equations, where the nonlocal time-space coupled fractional substantial derivative is involved. This paper focuses on the more widely used backward version. Based on the newly proposed approximation operators for fractional substantial derivative, we establish compact finite difference schemes for the backward fractional Feynman-Kac equation. The proposed difference schemes have the q-th(q = 1, 2, 3, 4) order accuracy in temporal direction and fourth order accuracy in spatial direction, respectively. The numerical stability and convergence in the maximum norm are proved for the first order time discretization scheme by the discrete energy method, where an inner product in complex space is introduced. Finally, extensive numerical experiments are carried out to verify the availability and superiority of the algorithms. Also, simulations of the backward fractional Feynman-Kac equation with Dirac delta function as the initial condition are performed to further confirm the effectiveness of the proposed methods.展开更多
In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,t...In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,the original Boussinesq system is transformed into an equivalent one.Then we discretize it using the second-order backward di erentiation formula(BDF2)and Crank-Nicolson(CN)to obtain two second-order time-advanced schemes.In both numerical schemes,a pressure-correction method is employed to decouple the velocity and pressure.These two schemes possess the desired property that they can be fully decoupled with satisfying unconditional stability.We rigorously prove both the unconditional stability and unique solvability of the discrete schemes.Furthermore,we provide detailed implementations of the decoupling procedures.Finally,various 2D numerical simulations are performed to verify the accuracy and energy stability of the proposed schemes.展开更多
We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditi...We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition （without the incompressibility condition）, which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in R^n（n ≥ 3） with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.展开更多
In this article, we study the nonlinear stochastic heat equation in the spatial domain R^d subject to a Gaussian noise which is white in time and colored in space. The spatial correlation can be any symmetric, nonnega...In this article, we study the nonlinear stochastic heat equation in the spatial domain R^d subject to a Gaussian noise which is white in time and colored in space. The spatial correlation can be any symmetric, nonnegative and nonnegative-definite function that satisfies Dalang's condition. We establish the existence and uniqueness of a random field solution starting from measure-valued initial data. We find the upper and lower bounds for the second moment. With these moment bounds, we first establish some necessary and sufficient conditions for the phase transition of the moment Lyapunov exponents, which extends the classical results from the stochastic heat equation on Z^d to that on R^d.Then, we prove a localization result for the intermittency fronts, which extends results by Conus and Khoshnevisan [9] from one space dimension to higher space dimension. The linear case has been recently proved by Huang et al [17] using different techniques.展开更多
基金financial support from the Open Fund(PLN1003) of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation(Southwest Petroleum University)the National Science and Technology Major Project in the l lth Five-Year Plan(Grant No.2008ZX05054)
文摘It is known that there is a discrepancy between field data and the results predicted from the previous equations derived by simplifying three-dimensional（3-D） flow into two-dimensions（2-D）.This paper presents a new steady-state productivity equation for horizontal wells in bottom water drive gas reservoirs.Firstly,the fundamental solution to the 3-D steady-state Laplace equation is derived with the philosophy of source and the Green function for a horizontal well located at the center of the laterally infinite gas reservoir.Then,using the fundamental solution and the Simpson integral formula,the average pseudo-pressure equation and the steady-state productivity equation are achieved for the horizontal section.Two case-studies are given in the paper,the results calculated from the newly-derived formula are very close to the numerical simulation performed with the Canadian software CMG and the real production data,indicating that the new formula can be used to predict the steady-state productivity of such horizontal gas wells.
基金Research partially supported by the “1000 Talents Plan” from Jilin University,Jilin Province and Chinese Governmentby the Simons Foundation(244767)
文摘This article establishes the precise asymptotics Eu^m(t, x)(t → ∞ or m → ∞) for the stochastic heat equation ?u/?t(t, x) =1/2?u(t, x) + u(t, x)(t, x)?W/?t(t, x) with the time-derivative Gaussian noise W?/?t(t, x) that is fractional in time and homogeneous in space.
基金Supported by the China National Science and Technology Major Project(2016ZX058-001)the CNOOC Scientific and Technological Project(CNOOC-KJ135-ZDXM36-TJ).
文摘Based on the measurement mechanism of mobility in pressure measurement while drilling, through analyzing a large number of mobility data, it is found that under the condition of water-based mud drilling, the product of mobility from pressure measurement while drilling and the viscosity of mud filtrate is infinitely close to the water phase permeability under the residual oil in relative permeability experiment. Based on this, a method converting the mobility from pressure measurement while drilling to core permeability is proposed, and the permeability based on Timur formula has been established. Application of this method in Penglai 19-9 oilfield of Bohai Sea shows:(1) Compared with the permeability calculated by the model of adjacent oilfields, the permeability calculated by this model is more consistent with the permeability calculated by core analysis.(2) Based on the new model, the correlation between the calculated mobility of well logging and the actual drilling specific productivity index bas been established. Compared with the relationship established by using the permeability model of an adjacent oilfield, the correlation of the new model is better.(3) Productivity of four directional wells was predicted, and the prediction results are in good agreement with the actual production after drilling.
基金Project supported by the Yue-Qi Scholar of the China University of Mining and Technology(Grant No.102504180004)the 333 Project of Jiangsu Province,China(Grant No.BRA2018320)
文摘To transform the exponential traveling wave solutions to bilinear differential equations, a sufficient and necessary condition is proposed. Motivated by the condition, we extend the results to the(2+1)-dimensional Kadomtsev–Petviashvili(KP) equation, the(3+1)-dimensional generalized Kadomtsev–Petviashvili(g-KP) equation, and the B-type Kadomtsev–Petviashvili(BKP) equation. Aa a result, we obtain some new resonant multiple wave solutions through the parameterization for wave numbers and frequencies via some linear combinations of exponential traveling waves. Finally, these new resonant type solutions can be displayed in graphs to illustrate the resonant behaviors of multiple wave solutions.
基金Project supported by the National Natural Science Foundation of China(Grant No.11471262)Henan University of Technology High-level Talents Fund,China(Grant No.2018BS039)
文摘The fractional Feynman-Kac equations describe the distributions of functionals of non-Brownian motion, or anomalous diffusion, including two types called the forward and backward fractional Feynman-Kac equations, where the nonlocal time-space coupled fractional substantial derivative is involved. This paper focuses on the more widely used backward version. Based on the newly proposed approximation operators for fractional substantial derivative, we establish compact finite difference schemes for the backward fractional Feynman-Kac equation. The proposed difference schemes have the q-th(q = 1, 2, 3, 4) order accuracy in temporal direction and fourth order accuracy in spatial direction, respectively. The numerical stability and convergence in the maximum norm are proved for the first order time discretization scheme by the discrete energy method, where an inner product in complex space is introduced. Finally, extensive numerical experiments are carried out to verify the availability and superiority of the algorithms. Also, simulations of the backward fractional Feynman-Kac equation with Dirac delta function as the initial condition are performed to further confirm the effectiveness of the proposed methods.
基金Supported by Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)+2 种基金Basic Research Plan of Shanxi Province(202203021211129)Shanxi Province Natural Science Research(202203021212249)Special/Youth Foundation of Taiyuan University of Technology(2022QN101)。
文摘In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,the original Boussinesq system is transformed into an equivalent one.Then we discretize it using the second-order backward di erentiation formula(BDF2)and Crank-Nicolson(CN)to obtain two second-order time-advanced schemes.In both numerical schemes,a pressure-correction method is employed to decouple the velocity and pressure.These two schemes possess the desired property that they can be fully decoupled with satisfying unconditional stability.We rigorously prove both the unconditional stability and unique solvability of the discrete schemes.Furthermore,we provide detailed implementations of the decoupling procedures.Finally,various 2D numerical simulations are performed to verify the accuracy and energy stability of the proposed schemes.
基金supported in part by the National Science Foundation under Grants DMS-0807551, DMS-0720925, and DMS-0505473the Natural Science Foundationof China (10728101)supported in part by EPSRC grant EP/F029578/1
文摘We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition （without the incompressibility condition）, which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in R^n（n ≥ 3） with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.
基金supported by the National Research Foundation of Korea (NRF-2017R1C1B1005436)the TJ Park Science Fellowship of POSCO TJ Park Foundation
文摘In this article, we study the nonlinear stochastic heat equation in the spatial domain R^d subject to a Gaussian noise which is white in time and colored in space. The spatial correlation can be any symmetric, nonnegative and nonnegative-definite function that satisfies Dalang's condition. We establish the existence and uniqueness of a random field solution starting from measure-valued initial data. We find the upper and lower bounds for the second moment. With these moment bounds, we first establish some necessary and sufficient conditions for the phase transition of the moment Lyapunov exponents, which extends the classical results from the stochastic heat equation on Z^d to that on R^d.Then, we prove a localization result for the intermittency fronts, which extends results by Conus and Khoshnevisan [9] from one space dimension to higher space dimension. The linear case has been recently proved by Huang et al [17] using different techniques.
