Suitable region of dynamic optimal interpolation for efficiently altimetry sea surface height mapping

The dynamic optimal interpolation(DOI)method is a technique based on quasi-geostrophic dynamics for merging multi-satellite altimeter along-track observations to generate gridded absolute dynamic topography(ADT).Compared with the linear optimal interpolation(LOI)method,the DOI method can improve the accuracy of gridded ADT locally but with low computational efficiency.Consequently,considering both computational efficiency and accuracy,the DOI method is more suitable to be used only for regional applications.In this study,we propose to evaluate the suitable region for applying the DOI method based on the correlation between the absolute value of the Jacobian operator of the geostrophic stream function and the improvement achieved by the DOI method.After verifying the LOI and DOI methods,the suitable region was investigated in three typical areas:the Gulf Stream(25°N-50°N,55°W-80°W),the Japanese Kuroshio(25°N-45°N,135°E-155°E),and the South China Sea(5°N-25°N,100°E-125°E).We propose to use the DOI method only in regions outside the equatorial region and where the absolute value of the Jacobian operator of the geostrophic stream function is higher than1×10^(-11).

Research on Interpolation Method for Missing Electricity Consumption Data

Missing value is one of the main factors that cause dirty data.Without high-quality data,there will be no reliable analysis results and precise decision-making.Therefore,the data warehouse needs to integrate high-quality data consistently.In the power system,the electricity consumption data of some large users cannot be normally collected resulting in missing data,which affects the calculation of power supply and eventually leads to a large error in the daily power line loss rate.For the problem of missing electricity consumption data,this study proposes a group method of data handling(GMDH)based data interpolation method in distribution power networks and applies it in the analysis of actually collected electricity data.First,the dependent and independent variables are defined from the original data,and the upper and lower limits of missing values are determined according to prior knowledge or existing data information.All missing data are randomly interpolated within the upper and lower limits.Then,the GMDH network is established to obtain the optimal complexity model,which is used to predict the missing data to replace the last imputed electricity consumption data.At last,this process is implemented iteratively until the missing values do not change.Under a relatively small noise level(α=0.25),the proposed approach achieves a maximum error of no more than 0.605%.Experimental findings demonstrate the efficacy and feasibility of the proposed approach,which realizes the transformation from incomplete data to complete data.Also,this proposed data interpolation approach provides a strong basis for the electricity theft diagnosis and metering fault analysis of electricity enterprises.

Integer multiple quantum image scaling based on NEQR and bicubic interpolation

As a branch of quantum image processing,quantum image scaling has been widely studied.However,most of the existing quantum image scaling algorithms are based on nearest-neighbor interpolation and bilinear interpolation,the quantum version of bicubic interpolation has not yet been studied.In this work,we present the first quantum image scaling scheme for bicubic interpolation based on the novel enhanced quantum representation(NEQR).Our scheme can realize synchronous enlargement and reduction of the image with the size of 2^(n)×2^(n) by integral multiple.Firstly,the image is represented by NEQR and the original image coordinates are obtained through multiple CNOT modules.Then,16 neighborhood pixels are obtained by quantum operation circuits,and the corresponding weights of these pixels are calculated by quantum arithmetic modules.Finally,a quantum matrix operation,instead of a classical convolution operation,is used to realize the sum of convolution of these pixels.Through simulation experiments and complexity analysis,we demonstrate that our scheme achieves exponential speedup over the classical bicubic interpolation algorithm,and has better effect than the quantum version of bilinear interpolation.

Wavelet Multi-Resolution Interpolation Galerkin Method for Linear Singularly Perturbed Boundary Value Problems

In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.

Comparison of Methods for Nitrate Interpolation in Wells in Aguascalientes, Mexico

The accuracy of interpolation models applied to groundwater depends, among other factors, on the interpolation method chosen. Therefore, it is necessary to compare different approaches. For this, different methods of interpolation of nitrate concentrations were contrasted in sixty-seven wells in an aquifer in Aguascalientes, Mexico. Four general interpolation methods were used in ArcGIS 10.5 to make the maps: IDW, Kriging, Natural Neighbor and Spline. In the modeling, only method type was varied. The input parameters (location, temporality, and nitrate concentration) were the same in the four interpolations;despite this, different maximum and minimum values were obtained for each interpolation method: for IDW, 0.2 to 22.0 mg/l, for Kriging, 3.5 to 16.5 mg/l, for Natural Neighbor, 0.3 to 21.7 mg/l and for Spline −30.8 to 37.2 mg/l. Finally, an assessment of the maps obtained was conducted by comparing them with the Official Mexican Standard (OMS), where 24 of the 67 wells were found outside the 10 mg/l that the OMS establishes as maximum permissible limit for human consumption. Taking as a starting point the measured values of nitrates (0.25 to 22.12 mg/l), as well as the spatial distribution of the interpolated values, it was determined that the Krigging method best fitted the data measured in the wells within the studied aquifer.

Comparison of Spatial Interpolation Methods of Precipitation Data in Central Macedonia, Greece

The purpose of this paper is to investigate the spatial interpolation of rainfall variability with deterministic and geostatic inspections in the Prefecture of Kilkis (Greece). The precipitation data where recorded from 12 meteorological stations in the Prefecture of Kilkis for 36 hydrological years (1973-2008). The cumulative monthly values of rainfall were studied on an annual and seasonal basis as well as during the arid-dry season. In the deterministic tests, the I.D.W. and R.B.F. checks were inspected, while in the geostatic tests, Ordinary Kriging and Universal Kriging respectively. The selection of the optimum method was made based on the least Root Mean Square Error (R.M.S.E.), as well as on the Mean Error (M.E.), as assessed by the cross validation analysis. The geostatical Kriging also considered the impact of isotropy and anisotropy across all time periods of data collection. Moreover, for Universal Kriging, the study explored spherical, exponential and Gaussian models in various combinations. Geostatistical techniques consistently demonstrated greater reliability than deterministic techniques across all time periods of data collection. Specifically, during the annual period, anisotropy was the prevailing characteristic in geostatistical techniques. Moreover, the results for the irrigation and seasonal periods were generally comparable, with few exceptions where isotropic methods yielded lower (R.M.S.E.) in some seasonal observations.

Quantum color image scaling based on bilinear interpolation

As a part of quantum image processing, quantum image scaling is a significant technology for the development of quantum computation. At present, most of the quantum image scaling schemes are based on grayscale images, with relatively little processing for color images. This paper proposes a quantum color image scaling scheme based on bilinear interpolation, which realizes the 2^(n_(1)) × 2^(n_(2)) quantum color image scaling. Firstly, the improved novel quantum representation of color digital images(INCQI) is employed to represent a 2^(n_(1)) × 2^(n_(2)) quantum color image, and the bilinear interpolation method for calculating pixel values of the interpolated image is presented. Then the quantum color image scaling-up and scaling-down circuits are designed by utilizing a series of quantum modules, and the complexity of the circuits is analyzed.Finally, the experimental simulation results of MATLAB based on the classical computer are given. The ultimate results demonstrate that the complexities of the scaling-up and scaling-down schemes are quadratic and linear, respectively, which are much lower than the cubic function and exponential function of other bilinear interpolation schemes.

OPTIMAL BIRKHOFF INTERPOLATION AND BIRKHOFF NUMBERS IN SOME FUNCTION SPACES

This paper investigates the optimal Birkhoff interpolation and Birkhoff numbers of some function spaces in space L∞[-1,1]and weighted spaces Lp,ω[-1,1],1≤p<∞,with w being a continuous integrable weight function in(-1,1).We proved that the Lagrange interpolation algorithms based on the zeros of some polynomials are optimal.We also show that the Lagrange interpolation algorithms based on the zeros of some polynomials are optimal when the function values of the two endpoints are included in the interpolation systems.

Missing interpolation model for wind power data based on the improved CEEMDAN method and generative adversarial interpolation network

Randomness and fluctuations in wind power output may cause changes in important parameters(e.g.,grid frequency and voltage),which in turn affect the stable operation of a power system.However,owing to external factors(such as weather),there are often various anomalies in wind power data,such as missing numerical values and unreasonable data.This significantly affects the accuracy of wind power generation predictions and operational decisions.Therefore,developing and applying reliable wind power interpolation methods is important for promoting the sustainable development of the wind power industry.In this study,the causes of abnormal data in wind power generation were first analyzed from a practical perspective.Second,an improved complete ensemble empirical mode decomposition with adaptive noise(ICEEMDAN)method with a generative adversarial interpolation network(GAIN)network was proposed to preprocess wind power generation and interpolate missing wind power generation sub-components.Finally,a complete wind power generation time series was reconstructed.Compared to traditional methods,the proposed ICEEMDAN-GAIN combination interpolation model has a higher interpolation accuracy and can effectively reduce the error impact caused by wind power generation sequence fluctuations.

Interpolation Technique for the Underwater DEM Generated by an Unmanned Surface Vessel

High-resolution underwater digital elevation models(DEMs)are important for water and soil conservation,hydrological analysis,and river channel dredging.In this work,the underwater topography of the Panjing River in Shanghai,China,was measured by an unmanned surface vessel.Five different interpolation methods were used to generate the underwater DEM and their precision and applicability for different underwater landforms were analyzed through cross-validation.The results showed that there was a positive correlation between the interpolation error and the terrain surface roughness.The five interpolation methods were all appropriate for the survey area,but their accuracy varied with different surface roughness.Based on the analysis results,an integrated approach was proposed to automatically select the appropriate interpolation method according to the different surface roughness in the surveying area.This approach improved the overall interpolation precision.The suggested technique provides a reference for the selection of interpolationmethods for underwater DEMdata.

Two Families of Multipoint Root-Solvers Using Inverse Interpolation with Memory

In this paper, a general family of derivative-free n + 1-point iterative methods using n + 1 evaluations of the function and a general family of n-point iterative methods using n evaluations of the function and only one evaluation of its derivative are constructed by the inverse interpolation with the memory on the previous step for solving the simple root of a nonlinear equation. The order and order of convergence of them are proved respectively. Finally, the proposed methods and the basins of attraction are demonstrated by the numerical examples.

Failure mechanisms and destruction characteristics of cemented coal gangue backfill under compression effect of non-uniform load

Backfill mining is one of the most important technical means for controlling strata movement and reducing surface subsidence and environmental damage during exploitation of underground coal resources. Ensuring the stability of the backfill bodies is the primary prerequisite for maintaining the safety of the backfilling working face, and the loading characteristics of backfill are closely related to the deformation and subsidence of the roof. Elastic thin plate model was used to explore the non-uniform subsidence law of the roof, and then the non-uniform distribution characteristics of backfill bodies’ load were revealed. Through a self-developed non-uniform loading device combined with acoustic emission (AE) and digital image correlation (DIC) monitoring technology, the synergistic dynamic evolution law of the bearing capacity, apparent crack, and internal fracture of cemented coal gangue backfills (CCGBs) under loads with different degrees of non-uniformity was deeply explored. The results showed that: 1) The uniaxial compressive strength (UCS) of CCGB increased and then decreased with an increase in the degree of non-uniformity of load (DNL). About 40% of DNL was the inflection point of DNL-UCS curve and when DNL exceeded 40%, the strength decreased in a cliff-like manner;2) A positive correlation was observed between the AE ringing count and UCS during the loading process of the specimen, which was manifested by a higher AE ringing count of the high-strength specimen. 3) Shear cracks gradually increased and failure mode of specimens gradually changed from “X” type dominated by tension cracks to inverted “Y” type dominated by shear cracks with an increase in DNL, and the crack opening displacement at the peak stress decreased and then increased. The crack opening displacement at 40% of the DNL was the smallest. This was consistent with the judgment of crack size based on the AE b-value, i. e., it showed the typical characteristics of “small b-value-large crack and large b-value-small crack”. The research results are of significance for preventing the instability and failure of backfill.

Study of the pressure transient behavior of directional wells considering the effect of non-uniform flux distribution

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.

Effect of non-uniform swelling on coal multiphysics during gas injection: The triangle approach

In current dual porosity/permeability models,there exists a fundamental assumption that the adsorption-induced swelling is distributed uniformly within the representative elementary volume (REV),irrespective of its internal structures and transient processes.However,both internal structures and transient processes can lead to the non-uniform swelling.In this study,we hypothesize that the non-uniform swelling is responsible for why coal permeability in experimental measurements is not only controlled by the effective stress but also is affected by the adsorption-induced swelling.We propose a concept of the swelling triangle composed of swelling paths to characterize the evolution of the non-uniform swelling and serve as a core link in coupled multiphysics.A swelling path is determined by a dimensionless volumetric ratio and a dimensionless swelling ratio.Different swelling paths have the same start and end point,and each swelling path represents a unique swelling case.The swelling path as the diagonal of the triangle represents the case of the uniform swelling while that as the two perpendicular boundaries represents the case of the localized swelling.The paths of all intermediate cases populate inside the triangle.The corresponding relations between the swelling path and the response of coal multiphysics are established by a non-uniform swelling coefficient.We define this method as the triangle approach and corresponding models as swelling path-based ones.The proposed concept and models are verified against a long-term experimental measurement of permeability and strains under constant effective stress.Our results demonstrate that during gas injection,coal multiphysics responses have a close dependence on the swelling path,and that in both future experiments and field predictions,this dependence must be considered.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

An improved non-uniform fast Fourier transform method for radio imaging of coronal mass ejections

Radioheliographs can obtain solar images at high temporal and spatial resolution,with a high dynamic range.These are among the most important instruments for studying solar radio bursts,understanding solar eruption events,and conducting space weather forecasting.This study aims to explore the effective use of radioheliographs for solar observations,specifically for imaging coronal mass ejections(CME),to track their evolution and provide space weather warnings.We have developed an imaging simulation program based on the principle of aperture synthesis imaging,covering the entire data processing flow from antenna configuration to dirty map generation.For grid processing,we propose an improved non-uniform fast Fourier transform(NUFFT)method to provide superior image quality.Using simulated imaging of radio coronal mass ejections,we provide practical recommendations for the performance of radioheliographs.This study provides important support for the validation and calibration of radioheliograph data processing,and is expected to profoundly enhance our understanding of solar activities.

Approximate solution of Volterra-Fredholm integral equations using generalized barycentric rational interpolant

It is well-known that interpolation by rational functions results in a more accurate approximation than the polynomials interpolation.However,classical rational interpolation has some deficiencies such as uncontrollable poles and low convergence order.In contrast with the classical rational interpolants,the generalized barycentric rational interpolants which depend linearly on the interpolated values,yield infinite smooth approximation with no poles in real numbers.In this paper,a numerical collocation approach,based on the generalized barycentric rational interpolation and Gaussian quadrature formula,was introduced to approximate the solution of Volterra-Fredholm integral equations.Three types of points in the solution domain are used as interpolation nodes.The obtained numerical results confirm that the barycentric rational interpolants are efficient tools for solving Volterra-Fredholm integral equations.Moreover,integral equations with Runge’s function as an exact solution,no oscillation occurrs in the obtained approximate solutions so that the Runge’s phenomenon is avoided.

A Remark on Pl Type Interpolation on Non-Uniformly Distributed Nodes on the Unit Circle

In this paper we study the problem of explicit representation and convergence of Pal type （0;1） interpolation and its converse, with some additional conditions, on the non-uniformly distributed nodes on the unit circle obtaIned by projecting the interlaced zeros of Pn （x） and Pn′ （x） on the unit circle. The motivation to this problem can be traced to the recent studies on the regularity of Birkhoff interpolation and Pal type interpolations on non-uniformly distributed zeros on the unit circle.

The Comparison of Spatial Interpolation Methods on Temperature and Precipitation of Sanjiangyuan Area

In order to get the spatial grid data of monthly precipitation and monthly average temperature of Sanjiangyuan area, the Co-Kriging (COK) and thin plate smoothing splines(TPS) interpolation methods were applied by using the climate data during 1971-2000 of 58 meteorological stations around Qinghai Province and the 3 arc-second digital elevation model (DEM) data. The performance was evaluated by the smallest statistical errors by general cross validation (GCV). Root-mean-squared predicted errors (RMSE) and mean absolute errors (MAE) were used to compare the performance of the two methods. The results showed that: 1) After combing covariates into the models, both methods performed better; 2) The performance of TPS was significantly better than COK: for monthly average temperature, the RMSE derived from TPS was 69.48% higher than COK, as MAE increased by 70.56%. And for monthly precipitation, the RMSE derived from TPS was 28.07% higher than COK, as MAE increased by 29.06%.

Analysis of radial basis function interpolation approach

The radial basis function （RBF） interpolation approach proposed by Freedman is used to solve inverse problems encountered in well-logging and other petrophysical issues. The approach is to predict petrophysical properties in the laboratory on the basis of physical rock datasets, which include the formation factor, viscosity, permeability, and molecular composition. However, this approach does not consider the effect of spatial distribution of the calibration data on the interpolation result. This study proposes a new RBF interpolation approach based on the Freedman＇s RBF interpolation approach, by which the unit basis functions are uniformly populated in the space domain. The inverse results of the two approaches are comparatively analyzed by using our datasets. We determine that although the interpolation effects of the two approaches are equivalent, the new approach is more flexible and beneficial for reducing the number of basis functions when the database is large, resulting in simplification of the interpolation function expression. However, the predicted results of the central data are not sufficiently satisfied when the data clusters are far apart.