Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media withi...Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media within these rocks.Faced with the challenge of calculating the three-dimensional fractal dimensions of rock porosity,this study proposes an innovative computational process that directly calculates the three-dimensional fractal dimensions from a geometric perspective.By employing a composite denoising approach that integrates Fourier transform(FT)and wavelet transform(WT),coupled with multimodal pore extraction techniques such as threshold segmentation,top-hat transformation,and membrane enhancement,we successfully crafted accurate digital rock models.The improved box-counting method was then applied to analyze the voxel data of these digital rocks,accurately calculating the fractal dimensions of the rock pore distribution.Further numerical simulations of permeability experiments were conducted to explore the physical correlations between the rock pore fractal dimensions,porosity,and absolute permeability.The results reveal that rocks with higher fractal dimensions exhibit more complex pore connectivity pathways and a wider,more uneven pore distribution,suggesting that the ideal rock samples should possess lower fractal dimensions and higher effective porosity rates to achieve optimal fluid transmission properties.The methodology and conclusions of this study provide new tools and insights for the quantitative analysis of complex pores in rocks and contribute to the exploration of the fractal transport properties of media within rocks.展开更多
Fractal dimensions of a terrain quantitatively describe the self-organizedstructure of the terrain geometry. However, the local topographic variation cannot be illustrated bythe conventional box-counting method. This ...Fractal dimensions of a terrain quantitatively describe the self-organizedstructure of the terrain geometry. However, the local topographic variation cannot be illustrated bythe conventional box-counting method. This paper proposes a successive shift box-counting method,in which the studied object is divided into small sub-objects that are composed of a series of gridsaccording to its characteristic scaling. The terrain fractal dimensions in the grids are calculatedwith the successive shift box-counting method and the scattered points with values of fractaldimensions are obtained. The present research shows that the planar variation of fractal dimensionsis well consistent with fault traces and geological boundaries.展开更多
The Laguerre spectral and pseudospectral methods are investigated for multidimensional nonlinear partial differential equations. Some results on the modified Laguerre orthogonal approximation and interpolation are est...The Laguerre spectral and pseudospectral methods are investigated for multidimensional nonlinear partial differential equations. Some results on the modified Laguerre orthogonal approximation and interpolation are established, which play important roles in the related numerical methods for unbounded domains. As an example, the modified Laguerre spectral and pseudospectral methods are proposed for two-dimensional Logistic equation. The stability and convergence of the suggested schemes are proved. Numerical results demonstrate the high accuracy of these approaches.展开更多
Stochastic fractional differential systems are important and useful in the mathematics,physics,and engineering fields.However,the determination of their probabilistic responses is difficult due to their non-Markovian ...Stochastic fractional differential systems are important and useful in the mathematics,physics,and engineering fields.However,the determination of their probabilistic responses is difficult due to their non-Markovian property.The recently developed globally-evolving-based generalized density evolution equation(GE-GDEE),which is a unified partial differential equation(PDE)governing the transient probability density function(PDF)of a generic path-continuous process,including non-Markovian ones,provides a feasible tool to solve this problem.In the paper,the GE-GDEE for multi-dimensional linear fractional differential systems subject to Gaussian white noise is established.In particular,it is proved that in the GE-GDEE corresponding to the state-quantities of interest,the intrinsic drift coefficient is a time-varying linear function,and can be analytically determined.In this sense,an alternative low-dimensional equivalent linear integer-order differential system with exact closed-form coefficients for the original highdimensional linear fractional differential system can be constructed such that their transient PDFs are identical.Specifically,for a multi-dimensional linear fractional differential system,if only one or two quantities are of interest,GE-GDEE is only in one or two dimensions,and the surrogate system would be a one-or two-dimensional linear integer-order system.Several examples are studied to assess the merit of the proposed method.Though presently the closed-form intrinsic drift coefficient is only available for linear stochastic fractional differential systems,the findings in the present paper provide a remarkable demonstration on the existence and eligibility of GE-GDEE for the case that the original high-dimensional system itself is non-Markovian,and provide insights for the physical-mechanism-informed determination of intrinsic drift and diffusion coefficients of GE-GDEE of more generic complex nonlinear systems.展开更多
In this paper, the box-counting dimension is used to derive an explicit formula for the dimension of a fractal constructed using several contractions or by combining fractals. This dimension agrees with the Hausdorff ...In this paper, the box-counting dimension is used to derive an explicit formula for the dimension of a fractal constructed using several contractions or by combining fractals. This dimension agrees with the Hausdorff dimension in the particular case when the scales factors considered are all the same. A more general sufficient condition for the box-counting dimension and the Hausdorff dimension to be the same is given. It is also shown that the dimension of the fractal obtained by combining two fractals is the weighted average of the dimensions of the two fractals.展开更多
In this paper, a new iterated function system consisting of non-linear affinemaps is constructed. We investigate the fractal interpolation functions generated bysuch a system and get its differentiability, its box dim...In this paper, a new iterated function system consisting of non-linear affinemaps is constructed. We investigate the fractal interpolation functions generated bysuch a system and get its differentiability, its box dimension, its packing dimension,and a lower bound of its Hansdorff dimension.展开更多
Two of the main challenges in optimal control are solving problems with state-dependent running costs and developing efficient numerical solvers that are computationally tractable in high dimensions.In this paper,we p...Two of the main challenges in optimal control are solving problems with state-dependent running costs and developing efficient numerical solvers that are computationally tractable in high dimensions.In this paper,we provide analytical solutions to certain optimal control problems whose running cost depends on the state variable and with constraints on the control.We also provide Lax-Oleinik-type representation formulas for the corresponding Hamilton-Jacobi partial differential equations with state-dependent Hamiltonians.Additionally,we present an efficient,grid-free numerical solver based on our representation formulas,which is shown to scale linearly with the state dimension,and thus,to overcome the curse of dimensionality.Using existing optimization methods and the min-plus technique,we extend our numerical solvers to address more general classes of convex and nonconvex initial costs.We demonstrate the capabilities of our numerical solvers using implementations on a central processing unit(CPU)and a field-programmable gate array(FPGA).In several cases,our FPGA implementation obtains over a 10 times speedup compared to the CPU,which demonstrates the promising performance boosts FPGAs can achieve.Our numerical results show that our solvers have the potential to serve as a building block for solving broader classes of high-dimensional optimal control problems in real-time.展开更多
Phase space reconstruction is the first step to recognizing the chaos from observed time series. On the basis of differential entropy, this paper introduces an efficient method to estimate the embedding dimension and ...Phase space reconstruction is the first step to recognizing the chaos from observed time series. On the basis of differential entropy, this paper introduces an efficient method to estimate the embedding dimension and the time delay simultaneously. The differential entropy is used to characterize the disorder degree of the reconstructed attractor. The minimum value of the differential entropy corresponds to the optimum set of the reconstructed parameters. Simulated experiments show that the original phase space can be effectively reconstructed from time series, and the accuracy of the invariants in phase space reconstruction is greatly improved. It provides a new method for the identification of chaotic signals from time series.展开更多
The finite dimension of the global attractors for the systems of the perturbed and unperturbed dissipative Hamiltonian amplitude equations governing modulated wave are investigated.An interesting result is also obtain...The finite dimension of the global attractors for the systems of the perturbed and unperturbed dissipative Hamiltonian amplitude equations governing modulated wave are investigated.An interesting result is also obtained that the upper bound of the dimension of the global attractor for the perturbed equation is independent of ε.展开更多
Phase space reconstruction is the first step of recognizing the chaotic time series.On the basis of differential entropy ratio method,the embedding dimension opt m and time delay t are optimal for the state space reco...Phase space reconstruction is the first step of recognizing the chaotic time series.On the basis of differential entropy ratio method,the embedding dimension opt m and time delay t are optimal for the state space reconstruction could be determined.But they are not the optimal parameters accepted for prediction.This study proposes an improved method based on the differential entropy ratio and Radial Basis Function(RBF)neural network to estimate the embedding dimension m and the time delay t,which have both optimal characteristics of the state space reconstruction and the prediction.Simulating experiments of Lorenz system and Doffing system show that the original phase space could be reconstructed from the time series effectively,and both the prediction accuracy and prediction length are improved greatly.展开更多
A set is called regular if its Hausdorff dimension and upper box-counting dimension coincide. In this paper,we prove that the random self-conformal set is regular almost surely. Also we determine the dimensions for a ...A set is called regular if its Hausdorff dimension and upper box-counting dimension coincide. In this paper,we prove that the random self-conformal set is regular almost surely. Also we determine the dimensions for a class of random self-conformal sets.展开更多
A class of multi-point boundary value problems are studied.Easily verified suffcient conditions to guarantee the existence of at least one solutions of above mentioned BVPs are established.The examples are presented t...A class of multi-point boundary value problems are studied.Easily verified suffcient conditions to guarantee the existence of at least one solutions of above mentioned BVPs are established.The examples are presented to illustrate the main results.展开更多
Increase in application fields of video has boosted the demand to analyze and organize video libraries for efficient scene analysis and information retrieval. This paper addresses the detection of shot transitions, wh...Increase in application fields of video has boosted the demand to analyze and organize video libraries for efficient scene analysis and information retrieval. This paper addresses the detection of shot transitions, which plays a crucial role in scene analysis, using a novel method based on fractal dimension (FD) that carries information on roughness of image intensity surface and textural structure. The proposed method is tested on sport videos including soccer and tennis matches that contain considerable amount of abrupt and gradual shot transitions. Experimental results indicate that the FD based shot transition detection method offers promising performance with respect to pixel and histogram based methods available in the literature.展开更多
We study a new algorithm for solvingparabolic partial differential equations(PDEs)and backward stochastic differential equations(BSDEs)in high dimension,which is based on an analogy between the BSDE and reinforcement ...We study a new algorithm for solvingparabolic partial differential equations(PDEs)and backward stochastic differential equations(BSDEs)in high dimension,which is based on an analogy between the BSDE and reinforcement learning with the gradient of the solution playing the role of the policy function,and the loss function given by the error between the prescribed terminal condition and the solution of the BSDE.The policy function is then approximated by a neural network,as is done in deep reinforcement learning.Numerical results using TensorFlow illustrate the efficiency and accuracy of the studied algorithm for several 100-dimensional nonlinear PDEs from physics and finance such as the Allen–Cahn equation,the Hamilton–Jacobi–Bellman equation,and a nonlinear pricing model for financial derivatives.展开更多
In the research of fractal cities, the fractal dimension is very important. It is used to describe the fractal character of the city. The authors have designed two approaches to calculate the fractal dimension by the ...In the research of fractal cities, the fractal dimension is very important. It is used to describe the fractal character of the city. The authors have designed two approaches to calculate the fractal dimension by the box-counting method through an example of Beijing, which are called the vector method and the grid method, respectively. The former calculates the fractal dimension through an intersecting analysis in ArcView; and the latter is carried out by programming in Matlab. They are compared from three aspects: the calculating process, the limits in use, and the results. As a result, the conclusion is made that there are merits and faults on both methods, and they should be chosen to use properly in practical situation.展开更多
In this paper we extend the theory of Grbner bases to difference-differential modules and present a new algorithmic approach for computing the Hilbert function of a finitely generated difference-differential module eq...In this paper we extend the theory of Grbner bases to difference-differential modules and present a new algorithmic approach for computing the Hilbert function of a finitely generated difference-differential module equipped with the natural filtration. We present and verify algorithms for construct-ing these Grbner bases counterparts. To this aim we introduce the concept of "generalized term order" on Nm ×Zn and on difference-differential modules. Using Grbner bases on difference-differential mod-ules we present a direct and algorithmic approach to computing the difference-differential dimension polynomials of a difference-differential module and of a system of linear partial difference-differential equations.展开更多
基金supported by the National Natural Science Foundation of China (Nos.52374078 and 52074043)the Fundamental Research Funds for the Central Universities (No.2023CDJKYJH021)。
文摘Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media within these rocks.Faced with the challenge of calculating the three-dimensional fractal dimensions of rock porosity,this study proposes an innovative computational process that directly calculates the three-dimensional fractal dimensions from a geometric perspective.By employing a composite denoising approach that integrates Fourier transform(FT)and wavelet transform(WT),coupled with multimodal pore extraction techniques such as threshold segmentation,top-hat transformation,and membrane enhancement,we successfully crafted accurate digital rock models.The improved box-counting method was then applied to analyze the voxel data of these digital rocks,accurately calculating the fractal dimensions of the rock pore distribution.Further numerical simulations of permeability experiments were conducted to explore the physical correlations between the rock pore fractal dimensions,porosity,and absolute permeability.The results reveal that rocks with higher fractal dimensions exhibit more complex pore connectivity pathways and a wider,more uneven pore distribution,suggesting that the ideal rock samples should possess lower fractal dimensions and higher effective porosity rates to achieve optimal fluid transmission properties.The methodology and conclusions of this study provide new tools and insights for the quantitative analysis of complex pores in rocks and contribute to the exploration of the fractal transport properties of media within rocks.
文摘Fractal dimensions of a terrain quantitatively describe the self-organizedstructure of the terrain geometry. However, the local topographic variation cannot be illustrated bythe conventional box-counting method. This paper proposes a successive shift box-counting method,in which the studied object is divided into small sub-objects that are composed of a series of gridsaccording to its characteristic scaling. The terrain fractal dimensions in the grids are calculatedwith the successive shift box-counting method and the scattered points with values of fractaldimensions are obtained. The present research shows that the planar variation of fractal dimensionsis well consistent with fault traces and geological boundaries.
基金the Science Foundation of the Science and Technology Commission of Shanghai Municipality(No.075105118)the Shanghai Leading Academic Discipline Project(No.T0401)the Fund for E-institute of Shanghai Universities(No.E03004)
文摘The Laguerre spectral and pseudospectral methods are investigated for multidimensional nonlinear partial differential equations. Some results on the modified Laguerre orthogonal approximation and interpolation are established, which play important roles in the related numerical methods for unbounded domains. As an example, the modified Laguerre spectral and pseudospectral methods are proposed for two-dimensional Logistic equation. The stability and convergence of the suggested schemes are proved. Numerical results demonstrate the high accuracy of these approaches.
基金The supports of the National Natural Science Foundation of China(Grant Nos.51725804 and U1711264)the Research Fund for State Key Laboratories of Ministry of Science and Technology of China(SLDRCE19-B-23)the Shanghai Post-Doctoral Excellence Program(2022558)。
文摘Stochastic fractional differential systems are important and useful in the mathematics,physics,and engineering fields.However,the determination of their probabilistic responses is difficult due to their non-Markovian property.The recently developed globally-evolving-based generalized density evolution equation(GE-GDEE),which is a unified partial differential equation(PDE)governing the transient probability density function(PDF)of a generic path-continuous process,including non-Markovian ones,provides a feasible tool to solve this problem.In the paper,the GE-GDEE for multi-dimensional linear fractional differential systems subject to Gaussian white noise is established.In particular,it is proved that in the GE-GDEE corresponding to the state-quantities of interest,the intrinsic drift coefficient is a time-varying linear function,and can be analytically determined.In this sense,an alternative low-dimensional equivalent linear integer-order differential system with exact closed-form coefficients for the original highdimensional linear fractional differential system can be constructed such that their transient PDFs are identical.Specifically,for a multi-dimensional linear fractional differential system,if only one or two quantities are of interest,GE-GDEE is only in one or two dimensions,and the surrogate system would be a one-or two-dimensional linear integer-order system.Several examples are studied to assess the merit of the proposed method.Though presently the closed-form intrinsic drift coefficient is only available for linear stochastic fractional differential systems,the findings in the present paper provide a remarkable demonstration on the existence and eligibility of GE-GDEE for the case that the original high-dimensional system itself is non-Markovian,and provide insights for the physical-mechanism-informed determination of intrinsic drift and diffusion coefficients of GE-GDEE of more generic complex nonlinear systems.
文摘In this paper, the box-counting dimension is used to derive an explicit formula for the dimension of a fractal constructed using several contractions or by combining fractals. This dimension agrees with the Hausdorff dimension in the particular case when the scales factors considered are all the same. A more general sufficient condition for the box-counting dimension and the Hausdorff dimension to be the same is given. It is also shown that the dimension of the fractal obtained by combining two fractals is the weighted average of the dimensions of the two fractals.
基金Supported by the National Natural Science Foundation of China and the Natural Science Foundtion of Zhejiang Province.
文摘In this paper, a new iterated function system consisting of non-linear affinemaps is constructed. We investigate the fractal interpolation functions generated bysuch a system and get its differentiability, its box dimension, its packing dimension,and a lower bound of its Hansdorff dimension.
基金supported by the DOE-MMICS SEA-CROGS DE-SC0023191 and the AFOSR MURI FA9550-20-1-0358supported by the SMART Scholarship,which is funded by the USD/R&E(The Under Secretary of Defense-Research and Engineering),National Defense Education Program(NDEP)/BA-1,Basic Research.
文摘Two of the main challenges in optimal control are solving problems with state-dependent running costs and developing efficient numerical solvers that are computationally tractable in high dimensions.In this paper,we provide analytical solutions to certain optimal control problems whose running cost depends on the state variable and with constraints on the control.We also provide Lax-Oleinik-type representation formulas for the corresponding Hamilton-Jacobi partial differential equations with state-dependent Hamiltonians.Additionally,we present an efficient,grid-free numerical solver based on our representation formulas,which is shown to scale linearly with the state dimension,and thus,to overcome the curse of dimensionality.Using existing optimization methods and the min-plus technique,we extend our numerical solvers to address more general classes of convex and nonconvex initial costs.We demonstrate the capabilities of our numerical solvers using implementations on a central processing unit(CPU)and a field-programmable gate array(FPGA).In several cases,our FPGA implementation obtains over a 10 times speedup compared to the CPU,which demonstrates the promising performance boosts FPGAs can achieve.Our numerical results show that our solvers have the potential to serve as a building block for solving broader classes of high-dimensional optimal control problems in real-time.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.50775198 and 60102002)
文摘Phase space reconstruction is the first step to recognizing the chaos from observed time series. On the basis of differential entropy, this paper introduces an efficient method to estimate the embedding dimension and the time delay simultaneously. The differential entropy is used to characterize the disorder degree of the reconstructed attractor. The minimum value of the differential entropy corresponds to the optimum set of the reconstructed parameters. Simulated experiments show that the original phase space can be effectively reconstructed from time series, and the accuracy of the invariants in phase space reconstruction is greatly improved. It provides a new method for the identification of chaotic signals from time series.
基金Supported Partially by the National Natural Science Foundation of China ( 1 0 1 31 0 5 0 ) ,theEducation Ministry of China and Shanghai Science and Technology Committee( 0 3QMH1 40 7)Supported by the National Natural Science Foundation of China( 1 986
文摘The finite dimension of the global attractors for the systems of the perturbed and unperturbed dissipative Hamiltonian amplitude equations governing modulated wave are investigated.An interesting result is also obtained that the upper bound of the dimension of the global attractor for the perturbed equation is independent of ε.
基金Supported by the Key Program of National Natural Science Foundation of China(Nos.61077071,51075349)Program of National Natural Science Foundation of Hebei Province(Nos.F2011203207,F2010001312)
文摘Phase space reconstruction is the first step of recognizing the chaotic time series.On the basis of differential entropy ratio method,the embedding dimension opt m and time delay t are optimal for the state space reconstruction could be determined.But they are not the optimal parameters accepted for prediction.This study proposes an improved method based on the differential entropy ratio and Radial Basis Function(RBF)neural network to estimate the embedding dimension m and the time delay t,which have both optimal characteristics of the state space reconstruction and the prediction.Simulating experiments of Lorenz system and Doffing system show that the original phase space could be reconstructed from the time series effectively,and both the prediction accuracy and prediction length are improved greatly.
文摘A set is called regular if its Hausdorff dimension and upper box-counting dimension coincide. In this paper,we prove that the random self-conformal set is regular almost surely. Also we determine the dimensions for a class of random self-conformal sets.
基金Supported by the Science Foundation of Educational Committee of Hunan Province(08C794)
文摘A class of multi-point boundary value problems are studied.Easily verified suffcient conditions to guarantee the existence of at least one solutions of above mentioned BVPs are established.The examples are presented to illustrate the main results.
文摘Increase in application fields of video has boosted the demand to analyze and organize video libraries for efficient scene analysis and information retrieval. This paper addresses the detection of shot transitions, which plays a crucial role in scene analysis, using a novel method based on fractal dimension (FD) that carries information on roughness of image intensity surface and textural structure. The proposed method is tested on sport videos including soccer and tennis matches that contain considerable amount of abrupt and gradual shot transitions. Experimental results indicate that the FD based shot transition detection method offers promising performance with respect to pixel and histogram based methods available in the literature.
文摘We study a new algorithm for solvingparabolic partial differential equations(PDEs)and backward stochastic differential equations(BSDEs)in high dimension,which is based on an analogy between the BSDE and reinforcement learning with the gradient of the solution playing the role of the policy function,and the loss function given by the error between the prescribed terminal condition and the solution of the BSDE.The policy function is then approximated by a neural network,as is done in deep reinforcement learning.Numerical results using TensorFlow illustrate the efficiency and accuracy of the studied algorithm for several 100-dimensional nonlinear PDEs from physics and finance such as the Allen–Cahn equation,the Hamilton–Jacobi–Bellman equation,and a nonlinear pricing model for financial derivatives.
文摘In the research of fractal cities, the fractal dimension is very important. It is used to describe the fractal character of the city. The authors have designed two approaches to calculate the fractal dimension by the box-counting method through an example of Beijing, which are called the vector method and the grid method, respectively. The former calculates the fractal dimension through an intersecting analysis in ArcView; and the latter is carried out by programming in Matlab. They are compared from three aspects: the calculating process, the limits in use, and the results. As a result, the conclusion is made that there are merits and faults on both methods, and they should be chosen to use properly in practical situation.
基金supported by the National Natural Science Foundation of China (Grant No. 60473019)the KLMM (Grant No. 0705)
文摘In this paper we extend the theory of Grbner bases to difference-differential modules and present a new algorithmic approach for computing the Hilbert function of a finitely generated difference-differential module equipped with the natural filtration. We present and verify algorithms for construct-ing these Grbner bases counterparts. To this aim we introduce the concept of "generalized term order" on Nm ×Zn and on difference-differential modules. Using Grbner bases on difference-differential mod-ules we present a direct and algorithmic approach to computing the difference-differential dimension polynomials of a difference-differential module and of a system of linear partial difference-differential equations.