Numerical Simulation of Oil-Water Two-Phase Flow in Low Permeability Tight Reservoirs Based on Weighted Least Squares Meshless Method

In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.

Legendre-Weighted Residual Methods for System of Fractional Order Differential Equations

The numerical approach for finding the solution of fractional order systems of boundary value problems (BPVs) is derived in this paper. The implementation of the weighted residuals such as Galerkin, Least Square, and Collocation methods are included for solving fractional order differential equations, which is broadened to acquire the approximate solutions of fractional order systems with differentiable polynomials, namely Legendre polynomials, as basis functions. The algorithm of the residual formulations of matrix form can be coded efficiently. The interpretation of Caputo fractional derivatives is employed here. We have demonstrated these methods numerically through a few examples of linear and nonlinear BVPs. The results in absolute errors show that the present method efficiently finds the numerical solutions of fractional order systems of differential equations.

Evaluation of soil erosion vulnerability in Hubei Province of China using RUSLE model and combination weighting method

Soil erosion has been recognized as a critical environmental issue worldwide.While previous studies have primarily focused on watershed-scale soil erosion vulnerability from a natural factor perspective,there is a notable gap in understanding the intricate interplay between natural and socio-economic factors,especially in the context of spatial heterogeneity and nonlinear impacts of human-land interactions.To address this,our study evaluates the soil erosion vulnerability at a provincial scale,taking Hubei Province as a case study to explore the combined effects of natural and socio-economic factors.We developed an evaluation index system based on 15 indicators of soil erosion vulnerability:exposure,sensitivity,and adaptability.In addition,the combination weighting method was applied to determine index weights,and the spatial interaction was analyzed using spatial autocorrelation,geographical temporally weighted regression and geographical detector.The results showed an overall decreasing soil erosion intensity in Hubei Province during 2000 and 2020.The soil erosion vulnerability increased before 2000 and then.The areas with high soil erosion vulnerability were mainly confined in the central and southern regions of Hubei Province(Xiantao,Tianmen,Qianjiang and Ezhou)with obvious spatial aggregation that intensified over time.Natural factors(habitat quality index)had negative impacts on soil erosion vulnerability,whereas socio-economic factors(population density)showed substantial spatial variability in their influences.There was a positive correlation between soil erosion vulnerability and erosion intensity,with the correlation coefficients ranging from-0.41 and 0.93.The increase of slope was found to enhance the positive correlation between soil erosion vulnerability and intensity.

WEIGHTED LEAST SQUARE METHOD FOR S-N CURVE FITTING

An S-N curve fitting approach is proposed based on the weighted least square method, and the weights are inversely proportional to the length of mean confidence intervals of experimental data sets. The assumption coincides with the physical characteristics of the fatigue life scatter. Two examples demonstrate the method. It is shown that the method has better accuracy and reasonableness compared with the usual least square method.

A Simple Weighted Integration Method for Calculating Surface Tension Force to Suppress Parasitic Flow in the Level Set Approach

Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded interface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.

Method of Semi Analytic Perturbation Weighted Residuals for Nonlinear Bending of Shallow Shells

The semi? analytic perturbation weighted residuals method was used to solve the nonlinear bending problem of shallow shells, and the fifth order B spline was taken as trial function to seek an efficient method for nonlinear bending problem of shallow shells. The results from the present method are in good agreement with those derived from other methods. The present method is of higher accuracy, lower computing time and wider adaptability. In addition, the design of computer program is simple and it is easy to be programmed.

ADAPTIVE FUSION ALGORITHMS BASED ON WEIGHTED LEAST SQUARE METHOD

Weighted fusion algorithms, which can be applied in the area of multi-sensor data fusion, are advanced based on weighted least square method. A weighted fusion algorithm, in which the relationship between weight coefficients and measurement noise is established, is proposed by giving attention to the correlation of measurement noise. Then a simplified weighted fusion algorithm is deduced on the assumption that measurement noise is uncorrelated. In addition, an algorithm, which can adjust the weight coefficients in the simplified algorithm by making estimations of measurement noise from measurements, is presented. It is proved by emulation and experiment that the precision performance of the multi-sensor system based on these algorithms is better than that of the multi-sensor system based on other algorithms.

A Simple Weighted Integration Method for Calculating Surface Tension Force to Suppress Parasitic Flow in the Level Set Approach

Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded in- terface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.

Preconditioned iterative methods for solving weighted linear least squares problems

A class of preconditioned iterative methods, i.e., preconditioned generalized accelerated overrelaxation （GAOR） methods, is proposed to solve linear systems based on a class of weighted linear least squares problems. The convergence and comparison results are obtained. The comparison results show that the convergence rate of the preconditioned iterative methods is better than that of the original methods. Furthermore, the effectiveness of the proposed methods is shown in the numerical experiment.

Gas kinetic flux solver based finite volume weighted essentially non-oscillatory scheme for inviscid compressible flows

A high-order gas kinetic flux solver(GKFS)is presented for simulating inviscid compressible flows.The weighted essentially non-oscillatory(WENO)scheme on a uniform mesh in the finite volume formulation is combined with the circular function-based GKFS(C-GKFS)to capture more details of the flow fields with fewer grids.Different from most of the current GKFSs,which are constructed based on the Maxwellian distribution function or its equivalent form,the C-GKFS simplifies the Maxwellian distribution function into the circular function,which ensures that the Euler or Navier-Stokes equations can be recovered correctly.This improves the efficiency of the GKFS and reduces its complexity to facilitate the practical application of engineering.Several benchmark cases are simulated,and good agreement can be obtained in comparison with the references,which demonstrates that the high-order C-GKFS can achieve the desired accuracy.

A redundant subspace weighting procedure for clock ensemble

A redundant-subspace-weighting(RSW)-based approach is proposed to enhance the frequency stability on a time scale of a clock ensemble.In this method,multiple overlapping subspaces are constructed in the clock ensemble,and the weight of each clock in this ensemble is defined by using the spatial covariance matrix.The superimposition average of covariances in different subspaces reduces the correlations between clocks in the same laboratory to some extent.After optimizing the parameters of this weighting procedure,the frequency stabilities of virtual clock ensembles are significantly improved in most cases.

A robust triaxial localization method of AE source using refraction path

Acoustic emission(AE)localization algorithms based on homogeneous media or single-velocity are less accurate when applied to the triaxial localization experiments.To the end,a robust triaxial localization method of AE source using refraction path is proposed.Firstly,the control equation of the refraction path is established according to the sensor coordinates and arrival times.Secondly,considering the influence of time-difference-of-arrival(TDOA)errors,the residual of the governing equation is calculated to estimate the equation weight.Thirdly,the refraction points in different directions are solved using Snell’s law and orthogonal constraints.Finally,the source coordinates are iteratively solved by weighted correction terms.The feasibility and accuracy of the proposed method are verified by pencil-lead breaking experiments.The simulation results show that the new method is almost unaffected by the refraction ratio,and always holds more stable and accurate positioning performance than the traditional method under different ratios and scales of TDOA outliers.

A weight's agglomerative method for detecting communities in weighted networks based on weight's similarity

This paper proposes the new definition of the community structure of the weighted networks that groups of nodes in which the edge＇s weights distribute uniformly but at random between them. It can describe the steady connections between nodes or some similarity between nodes＇ functions effectively. In order to detect the community structure efficiently, a threshold coefficient t~ to evaluate the equivalence of edges＇ weights and a new weighted modularity based on the weight＇s similarity are proposed. Then, constructing the weighted matrix and using the agglomerative mechanism, it presents a weight＇s agglomerative method based on optimizing the modularity to detect communities. For a network with n nodes, the algorithm can detect the community structure in time O（n2 log~）. Simulations on networks show that the algorithm has higher accuracy and precision than the existing techniques. Furthermore, with the change of t~ the algorithm discovers a special hierarchical organization which can describe the various steady connections between nodes in groups.

A WEIGHTED RESIDUAL METHOD FOR ELASTIC-PLASTIC ANALYSIS NEAR A CRACK TIP AND THE CALCULATION OF THE PLASTIC STRESS INTENSITY FACTORS

In this paper, a weighted residual method for the elastic-plastic analysis near a crack tip is systematically given by taking the model of power-law hardening under plane strain condition as a sample. The elastic-plastic solutions of the crack lip field and an approach based on the superposition of the nonlinear finite element method on the complete solution in the whole crack body field, to calculate the plastic stress intensity factors, are also developed. Therefore, a complete analvsis based on the calculation both for the crack tip field and for the whole crack body field is provided.

Indoor air quality risk assessment on the Canadian campus by a multilevel integrated weighted average method

Increasing incidents of indoor air quality(IAQ) related complaints lead us to the fact that IAQ has become a significant occupational health and environmental issue. However, how to effectively evaluate IAQ under different scale of multiple indicators is still a challenge. The traditional single-indicator method is subjected to uncertainties in assessing IAQ due to different subjectivity on good or bad quality and scalar differences of data set. In this study, a multilevel integrated weighted average IAQ method including initial walking through assessment(IWA) and two-layers weighted average method are developed and applied to evaluate IAQ of the laboratory building at the University of Regina in Canada. Some important chemical parameters related to IAQ in terms of volatile organic compounds(VOCs), methanol(HCHO), carbon dioxide(CO2), and carbon monoxide(CO) are evaluated based on 5 months continuous monitoring data. The new integrated assessment result can not only indicates the risk of an individual parameter, but also able to quantify the overall IAQ risk on the sampling site. Finally, some recommendations based on the result are proposed to address sustainable IAQ practices in the sampling area.

Criteria for Weighted Moving-Mean Method

The moving-mean method is one of the conventional approaches for trend-extraction from a data set. It is usually applied in an empirical way. The smoothing degree of the trend depends on the selections of window length and weighted coefficients, which are associated with the change pattern of the data. Are there any uniform criteria for determining them? The present article is a reaction to this fundamental problem. By investigating many kinds of data, the results show that: 1) Within a certain range, the more points which participate in moving-mean, the better the trend function. However, in case the window length is too long, the trend function may tend to the ordinary global mean. 2) For a given window length, what matters is the choice of weighted coefficients. As the five-point case concerned, the local-midpoint, local-mean and global-mean criteria hold. Among these three criteria, the local-mean one has the strongest adaptability, which is suggested for your usage.

Dynamic Equivalent Method of Motor Loads for Power Systems Based on the Weighted

Dynamic equivalence can not only largely reduce the system size and the computation time but also stress the dominant features of the system [1]-[3]. This paper firstly recommends the basic concept of dynamic equivalent and the status of both domestic and abroad development in this area. The most existing equivalent methods usually only deal with static load models and neglect the dynamic characteristics of loads such as induction motors. In addition, the existing polymerization method which is based on the frequency domain algorithm of induction electric machines parameters takes a long time to equivalent for the large system, then the new method based on the weighted is proposed. Then, the basic steps for dynamic equivalence with the weighted method are introduced as follows. At first, the clustering criterion of motor loads based on time domain simulation is given. The motors with similar dynamic characteristics are classified into one group. Then, the simplication of the buses of motors in same group and network is carried out. Finally, parameters of the equivalent motor are calculated and the equivalent system is thus obtained based on the weighted. This aggregation method is applied to the simple distribution system of 4 generators. Simulation results show that the method can quickly obtain polymerization parameters of generator groups and the aggregation model retains the dynamic performance of the original model with good accuracy, the active and reactive power fitting error is smaller as well.

Free Energy Level Correction by Monte Carlo Resampling with Weighted Histogram Analysis Method

Free energy calculations may provide vital information for studying various chemical and biological processes.Quantum mechanical methods are required to accurately describe interaction energies,but their computations are often too demanding for conformational sampling.As a remedy,level correction schemes that allow calculating high level free energies based on conformations from lower level simulations have been developed.Here,we present a variation of a Monte Carlo(MC)resampling approach in relation to the weighted histogram analysis method(WHAM).We show that our scheme can generate free energy surfaces that can practically converge to the exact one with sufficient sampling,and that it treats cases with insufficient sampling in a more stable manner than the conventional WHAM-based level correction scheme.It can also provide a guide for checking the uncertainty of the levelcorrected surface and a well-defined criterion for deciding the extent of smoothing on the free energy surface for its visual improvement.We demonstrate these aspects by obtaining the free energy maps associated with the alanine dipeptide and proton transfer network of the KillerRed protein in explicit water,and exemplify that the MC resampled WHAM scheme can be a practical tool for producing free energy surfaces of realistic systems.

Evaluating gravity gradient components based on a reweighted inversion method

In gravity gradient inversion,to choose an appropriate component combination is very important,that needs to understand the function of each component of gravity gradient in the inversion.In this paper,based on the previous research on the characteristics of gravity gradient components,we propose a reweighted inversion method to evaluate the influence of single gravity gradient component on the inversion resolution The proposed method only adopts the misfit function of the regularized equation and introduce a depth weighting function to overcome skin effect produced in gravity gradient inversion.A comparison between different inversion results was undertaken to verify the influence of the depth weighting function on the inversion result resolution.To avoid the premise of introducing prior information,we select the depth weighting function based on the sensitivity matrix.The inversion results using the single-prism model and the complex model show that the influence of different components on the resolution of inversion results is different in different directions,however,the inversion results based on two kind of models with adding different levels of random noise are basically consistent with the results of inversion without noises.Finally,the method was applied to real data from the Vinton salt dome,Louisiana,USA.

Gaussian Weighted Trajectory Method Ⅳ： No Rainbow Effect in Practice

The Gaussian weighted trajectory method （GWTM） is a practical implementation of classical S matrix theory （CSMT） in the random phase approximation, CSMT being the first and simplest semi-classical approach of molecular collisions, developped in the early seventies. Though very close in spirit to the purely classical description, GWTM accounts to some extent for the quantization of the different degrees-of-freedom involved in the processes. While CSMT may give diverging final state distributions, in relation to the rainbow effect of elastic scattering theory, GWTM has never led to such a mathematical catastrophe. The goal of the present note is to explain this finding.