The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LW...The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.展开更多
The vertical two-dimensional non-hydrostatic pressure models with multiple layers can make prediction more accurate than those obtained by the hydrostatic pres- sure assumption. However, they are time-consuming and un...The vertical two-dimensional non-hydrostatic pressure models with multiple layers can make prediction more accurate than those obtained by the hydrostatic pres- sure assumption. However, they are time-consuming and unstable, which makes them unsuitable for wider application. In this study, an efficient model with a single layer is developed. Decomposing the pressure into the hydrostatic and dynamic components and integrating the x-momentum equation from the bottom to the free surface can yield a horizontal momentum equation, in which the terms relevant to the dynamic pressure are discretized semi-implicitly. The convective terms in the vertical momentum equation are ignored, and the rest of the equation is approximated with the Keller-box scheme. The velocities expressed as the unknown dynamic pressure are substituted into the continuity equation, resulting in a tri-diagonal linear system solved by the Thomas algorithm. The validation of solitary and sinusoidal waves indicates that the present model can provide comparable results to the models with multiple layers but at much lower computation cost.展开更多
The efficiency of solving computationally partial differential equations can be profoundly highlighted by the creation of precise,higher-order compact numerical scheme that results in truly outstanding accuracy at a g...The efficiency of solving computationally partial differential equations can be profoundly highlighted by the creation of precise,higher-order compact numerical scheme that results in truly outstanding accuracy at a given cost.The objective of this article is to develop a highly accurate novel algorithm for two dimensional non-linear Burgers Huxley(BH)equations.The proposed compact numerical scheme is found to be free of superiors approximate oscillations across discontinuities,and in a smooth ow region,it efciently obtained a high-order accuracy.In particular,two classes of higherorder compact nite difference schemes are taken into account and compared based on their computational economy.The stability and accuracy show that the schemes are unconditionally stable and accurate up to a two-order in time and to six-order in space.Moreover,algorithms and data tables illustrate the scheme efciency and decisiveness for solving such non-linear coupled system.Efciency is scaled in terms of L_(2) and L_(∞) norms,which validate the approximated results with the corresponding analytical solution.The investigation of the stability requirements of the implicit method applied in the algorithm was carried out.Reasonable agreement was constructed under indistinguishable computational conditions.The proposed methods can be implemented for real-world problems,originating in engineering and science.展开更多
This paper present an implementation of＂modified cubic B-spline differential quadrature method （MCB-DQM）＂ proposed by Arora ＆ Singh （Applied Mathematics and Computation Vol. 224（1） （2013） 161-177） for numer...This paper present an implementation of＂modified cubic B-spline differential quadrature method （MCB-DQM）＂ proposed by Arora ＆ Singh （Applied Mathematics and Computation Vol. 224（1） （2013） 161-177） for numerical computation of Fokker-Planck equations. The modified cubic B-splines are used as set of basis functions in the differential quadrature to compute the weighting coefficients for the spatial derivatives, which reduces Fokker-Planck equation into system of first-order ordinary differential equations （ODEs）, in time. The well known SSP-RK43 scheme is then applied to solve the resulting system of ODEs. The efficiency of proposed method has been confirmed by three examples having their exact solutions. This shows that MCB-DQM results are capable of achieving high accuracy. Advantage of the scheme is that it can be applied very smoothly to solve the linear or nonlinear physical problems, and a very less storage space is required which causes less accumulation of numerical errors.展开更多
基金supported by the Major Research Project on Scientific Instrument Development of the National Natural Science Foundation of China(42327901)National Natural Science Foundation of China(42030806,42074120,41904104,423B2405).
文摘The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.
基金Project supported by the Specialized Research Fund for the Doctoral Program of Higher Education(No. 20110142110064)the Ministry of Water Resources’ Science and Technology Promotion Plan Program (No. TG1316)
文摘The vertical two-dimensional non-hydrostatic pressure models with multiple layers can make prediction more accurate than those obtained by the hydrostatic pres- sure assumption. However, they are time-consuming and unstable, which makes them unsuitable for wider application. In this study, an efficient model with a single layer is developed. Decomposing the pressure into the hydrostatic and dynamic components and integrating the x-momentum equation from the bottom to the free surface can yield a horizontal momentum equation, in which the terms relevant to the dynamic pressure are discretized semi-implicitly. The convective terms in the vertical momentum equation are ignored, and the rest of the equation is approximated with the Keller-box scheme. The velocities expressed as the unknown dynamic pressure are substituted into the continuity equation, resulting in a tri-diagonal linear system solved by the Thomas algorithm. The validation of solitary and sinusoidal waves indicates that the present model can provide comparable results to the models with multiple layers but at much lower computation cost.
文摘The efficiency of solving computationally partial differential equations can be profoundly highlighted by the creation of precise,higher-order compact numerical scheme that results in truly outstanding accuracy at a given cost.The objective of this article is to develop a highly accurate novel algorithm for two dimensional non-linear Burgers Huxley(BH)equations.The proposed compact numerical scheme is found to be free of superiors approximate oscillations across discontinuities,and in a smooth ow region,it efciently obtained a high-order accuracy.In particular,two classes of higherorder compact nite difference schemes are taken into account and compared based on their computational economy.The stability and accuracy show that the schemes are unconditionally stable and accurate up to a two-order in time and to six-order in space.Moreover,algorithms and data tables illustrate the scheme efciency and decisiveness for solving such non-linear coupled system.Efciency is scaled in terms of L_(2) and L_(∞) norms,which validate the approximated results with the corresponding analytical solution.The investigation of the stability requirements of the implicit method applied in the algorithm was carried out.Reasonable agreement was constructed under indistinguishable computational conditions.The proposed methods can be implemented for real-world problems,originating in engineering and science.
文摘This paper present an implementation of＂modified cubic B-spline differential quadrature method （MCB-DQM）＂ proposed by Arora ＆ Singh （Applied Mathematics and Computation Vol. 224（1） （2013） 161-177） for numerical computation of Fokker-Planck equations. The modified cubic B-splines are used as set of basis functions in the differential quadrature to compute the weighting coefficients for the spatial derivatives, which reduces Fokker-Planck equation into system of first-order ordinary differential equations （ODEs）, in time. The well known SSP-RK43 scheme is then applied to solve the resulting system of ODEs. The efficiency of proposed method has been confirmed by three examples having their exact solutions. This shows that MCB-DQM results are capable of achieving high accuracy. Advantage of the scheme is that it can be applied very smoothly to solve the linear or nonlinear physical problems, and a very less storage space is required which causes less accumulation of numerical errors.
