The capillary pressure curves from upscaling interfacial and unsaturated flows in porous layers with vertical heterogeneity

We provide the capillary pressure curves p_(c)(s)as a function of the effective saturation s based on the theoretical framework of upscaling unsaturated flows in vertically heterogeneous porous layers proposed recently(Z.Zheng,Journal of Fluid Mechanics,950,A17,2022).Based on the assumption of vertical gravitational-capillary equilibrium,the saturation distribution and profile shape of the invading fluid can be obtained by solving a nonlinear integral-differential equation.The capillary pressure curves p_(c)(s)can then be constructed by systematically varying the injection rate.Together with the relative permeability curves k_(rn)(s)that are already obtained.One can now provide quick estimates on the overall behaviours of interfacial and unsaturated flows in vertically-heterogeneous porous layers.

结合UPSCALING技术与对象多特征的土地利用/覆盖信息提取研究

应用SPOT融合数据,以北京密云地区为例,提出了整合Upscaling技术与对象多特征方法的新思路,通过基于半变异函数的空间变异特征分析,建立了面向对象多特征与多分辨率数据集的多尺度分类决策树,并对自动分类效率进行了初步探讨。

A Method for Upscaling Genetic Parameters of CERES-Rice in Regional Applications

To upscale the genetic parameters of CERES-Rice in regional applications, Jiangsu Province, the second largest rice producing province in China, was taken as an example. The province was divided into four rice regions with different rice variety types, and five to six sites in each region were selected. Then the eight genetic parameters of CERES-Rice, particularly the four parameters related to the yield, were modified and validated using the Trial and Error Method and the local statistical data of rice yield at a county level from 2001 to 2004, combined with the regional experiments of rice varieties in the province as well as the local meteorological and soil data （Method 1）. The simulated results of Method 1 were compared with those of other three traditional methods upscaling the genetic parameters, i.e., using one-site experimental data from a local representative rice variety （Method 2）, using local long-term rice yield data at a county level after deducting the trend yield due to progress of science and technology （Method 3）, and using rice yield data at a super scale, such as provincial, ecological zone, country or continent levels （Method 4）. The results showed that the best fitness was obtained by using the Method 1. The coefficients of correlation between the simulated yield and the statistical yield in the Method 1 were significant at 0.05 or 0.01 levels and the root mean squared error （RMSE） values were less than 9% for all the four rice regions. The method for upscaling the genetic parameters of CERES-Rice presented is not only valuable for the impact studies of climate change, but also favorable to provide a methodology for reference in crop model applications to the other regional studies.

New dynamic permeability upscaling method for flow simulation under depletion drive and no-crossflow conditions

The main purpose of upscaling in reservoir simulation is to capture the dynamic behavior of fine scale models at the coarse scale. Traditional static or dynamic methods use assumptions about the boundary conditions to determine the upscaled properties, in this paper, we show that the upscaled properties are strongly dependent on the flow process observed at the fine scale. We use a simple no- crossflow depletion drive process and demonstrate that an upscaled property is not a constant value. Instead, if the goal is to match the performance of the fine scale model, the upscaled permeability changes with time. We provide an analytical solution to determine the upscaled permeability and present the value of upscaled permeability under limiting conditions. Our equation suggests that it is possible that upscaled value can fall outside the range of fine scale values under certain conditions. We show that for pseudo steady state flow, using common averaging methods like arithmetic or even geometric averaging methods can lead to optimistic results. We also show that the no-crossflow solution is significantly different than crossflow solution at late times. We validate our method by comparing the results of the method with flow simulation results in two and multi-layered models.

Data-Driven Upscaling of Orientation Kinematics in Suspensions of Rigid Fibres

Describing the orientation state of the particles is often critical in fibre suspension applications.Macroscopic descriptors,the so-called second-order orientation tensor(or moment)leading the way,are often preferred due to their low computational cost.Closure problems however arise when evolution equations for the moments are derived from the orientation distribution functions and the impact of the chosen closure is often unpredictable.In this work,our aim is to provide macroscopic simulations of orientation that are cheap,accurate and closure-free.To this end,we propose an innovative data-based approach to the upscaling of orientation kinematics in the context of fibre suspensions.Since the physics at the microscopic scale can be modelled reasonably enough,the idea is to conduct accurate offline direct numerical simulations at that scale and to extract the corresponding macroscopic descriptors in order to build a database of scenarios.During the online stage,the macroscopic descriptors can then be updated quickly by combining adequately the items from the database instead of relying on an imprecise macroscopic model.This methodology is presented in the well-known case of dilute fibre suspensions(where it can be compared against closure-based macroscopic models)and in the case of suspensions of confined or electrically-charged fibres,for which state-of-the-art closures proved to be inadequate or simply do not exist.

Simultaneous upscaling of two properties of reservoirs in one dimension using adaptive bandwidth in kernel function method

Upscaling of primary geological models with huge cells, especially in porous media, is the first step in fluid flow simulation. Numerical methods are often used to solve the models. The upscaling method must preserve the important properties of the spatial distribution of the reservoir properties. An grid upscaling method based on adaptive bandwidth in kernel function is proposed according to the spatial distribution of property. This type of upscaling reduces the number of cells, while preserves the main heterogeneity features of the original fine model. The key point of the paper is upscaling two reservoir properties simultaneously. For each reservoir feature, the amount of bandwidth or optimal threshold is calculated and the results of the upscaling are obtained. Then two approaches are used to upscaling two properties simultaneously based on maximum bandwidth and minimum bandwidth. In fact, we now have a finalized upscaled model for both reservoir properties for each approach in which not only the number of their cells, but also the locations of the cells are equal. The upscaling error of the minimum bandwidth approach is less than that of the maximum bandwidth approach.

Muddy interlayer forecasting and an equivalent upscaling method based on tortuous paths: A case study of Mackay River oil sand reservoirs in Canada

Based on the abundant core data of oil sands in the Mackay river in Canada,the termination frequency of muddy interlayers was counted to predict the extension range of interlayers using a queuing theory model,and then the quantitative relationship between the thickness and extension length of muddy interlayer was established.An equivalent upscaling method of geologic model based on tortuous paths under the effects of muddy interlayer has been proposed.Single muddy interlayers in each coarse grid are tracked and identified,and the average length,width and proportion of muddy interlayer in each coarse grid are determined by using the geological connectivity tracing algorithm.The average fluid flow length of tortuous path under the influence of muddy interlayer is calculated.Based on the Darcy formula,the formula calculating average permeability in the coarsened grid is deduced to work out the permeability of equivalent coarsened grid.The comparison of coarsening results of the oil sand reservoir of Mackay River with actual development indexes shows that the equivalent upscaling method of muddy interlayer by tortuous path calculation can reflect the blocking effect of muddy interlayer very well,and better reflect the effects of geological condition on production.

Sensitivity-based upscaling for history matching of reservoir models

Simulation of reservoir flow processes at the finest scale is computationally expensive and in some cases impractical.Consequently,upscaling of several fine-scale grid blocks into fewer coarse-scale grids has become an integral part of reservoir simulation for most reservoirs.This is because as the number of grid blocks increases,the number of flow equations increases and this increases,in large proportion,the time required for solving flow problems.Although we can adopt parallel computation to share the load,a large number of grid blocks still pose significant computational challenges.Thus,upscaling acts as a bridge between the reservoir scale and the simulation scale.However as the upscaling ratio is increased,the accuracy of the numerical simulation is reduced;hence,there is a need to keep a balance between the two.In this work,we present a sensitivity-based upscaling technique that is applicable during history matching.This method involves partial homogenization of the reservoir model based on the model reduction pattern obtained from analysis of the sensitivity matrix.The technique is based on wavelet transformation and reduction of the data and model spaces as presented in the 2Dwp-wk approach.In the 2Dwp-wk approach,a set of wavelets of measured data is first selected and then a reduced model space composed of important wavelets is gradually built during the first few iterations of nonlinear regression.The building of the reduced model space is done by thresholding the full wavelet sensitivity matrix.The pattern of permeability distribution in the reservoir resulting from the thresholding of the full wavelet sensitivity matrix is used to determine the neighboring grids that are upscaled.In essence,neighboring grid blocks having the same permeability values due to model space reduction are combined into a single grid block in the simulation model,thus integrating upscaling with wavelet multiscale inverse modeling.We apply the method to estimate the parameters of two synthetic reservoirs.The history matching results obtained using this sensitivity-based upscaling are in very close agreement with the match provided by fine-scale inverse analysis.The reliability of the technique is evaluated using various scenarios and almost all the cases considered have shown very good results.The technique speeds up the history matching process without seriously compromising the accuracy of the estimates.

Development of a convolutional neural network based geomechanical upscaling technique for heterogeneous geological reservoir

Geomechanical assessment using coupled reservoir-geomechanical simulation is becoming increasingly important for analyzing the potential geomechanical risks in subsurface geological developments.However,a robust and efficient geomechanical upscaling technique for heterogeneous geological reservoirs is lacking to advance the applications of three-dimensional(3D)reservoir-scale geomechanical simulation considering detailed geological heterogeneities.Here,we develop convolutional neural network(CNN)proxies that reproduce the anisotropic nonlinear geomechanical response caused by lithological heterogeneity,and compute upscaled geomechanical properties from CNN proxies.The CNN proxies are trained using a large dataset of randomly generated spatially correlated sand-shale realizations as inputs and simulation results of their macroscopic geomechanical response as outputs.The trained CNN models can provide the upscaled shear strength(R^(2)>0.949),stress-strain behavior(R^(2)>0.925),and volumetric strain changes(R^(2)>0.958)that highly agree with the numerical simulation results while saving over two orders of magnitude of computational time.This is a major advantage in computing the upscaled geomechanical properties directly from geological realizations without the need to perform local numerical simulations to obtain the geomechanical response.The proposed CNN proxybased upscaling technique has the ability to(1)bridge the gap between the fine-scale geocellular models considering geological uncertainties and computationally efficient geomechanical models used to assess the geomechanical risks of large-scale subsurface development,and(2)improve the efficiency of numerical upscaling techniques that rely on local numerical simulations,leading to significantly increased computational time for uncertainty quantification using numerous geological realizations.

Super-resolution GANs for upscaling unplanned urban settlements from remote sensing satellite imagery - the case of Chinese urban village detection

The semantic segmentation of informal urban settlements represents an essential contribution towards renovation strategies and reconstruction plans.In this context,however,a big challenge remains unsolved when dealing with incomplete data acquisitions from multiple sensing devices,especially when study areas are depicted by images of different resolutions.In practice,traditional methodologies are directed to downgrade the higher-resolution data to the lowest-resolution measure,to define an overall homogeneous dataset,which is however ineffective in downstream segmentation activities of such crowded unplanned urban environments.To this purpose,we hereby tackle the problem in the opposite direction,namely upscaling the lower-resolution data to the highest-resolution measure,contributing to assess the use of cutting-edge super-resolution generative adversarial network(SR-GAN)architectures.The experimental novelty targets the particular case involving the automatic detection of‘urban villages’,sign of the quick transformation of Chinese urban environments.By aligning image resolutions from two different data sources(Gaofen-2 and Sentinel-2 data),we evaluated the degree of improvement with regard to pixel-based landcover segmentation,achieving,on a 1 m resolution target,classification accuracies up to 83%,67%and 56%for 4x,8x,and 10x resolution upgrades respectively,disclosing the advantages of artificially-upscaled images for segmenting detailed characteristics of informal settlements.

RECENT ADVANCES OF UPSCALING METHODS FOR THE SIMULATION OF FLOW TRANSPORT THROUGH HETEROGENEOUS POROUS MEDIA

We review some of our recent efforts in developing upscaling methods for simulating the flow transport through heterogeneous porous media. In particular, the steady flow transport through highly heterogeneous porous media driven by extraction wells and the flow transport through unsaturated porous media will be considered.

Upscaling in situ albedo for validation of coarse scale albedo product over heterogeneous surfaces

One of the essential steps for satellite albedo validation is upscaling in situ measurements to corresponding pixel scale over relatively heterogeneous land surfaces.Although the multi-scale validation strategy is applicable for heterogeneous surfaces,the calibration of the high-resolution imagery during upscaling process is never perfect,and thus the upscaling results suffer from errors.The regression-kriging(RK)technique can compensate the calibration part by applying kriging to upscale residuals and produce more accurate upscaling results.In this paper,in situ measurements and high spatial resolution albedo imagery combined with RK technique was proposed.This method is illustrated by upscaling surface albedo from in situ measurements scale to the coarse pixel scale in the core experimental area of HiWATER,where 17 WSN nodes were deployed at heterogeneous area.The upscaling results of this method were compared with the upscaling results from multi-scale strategy.The results show that the upscaling method based on in situ measurements and high-resolution imagery combined with RK technique can capture the spatial characteristics of surface albedo better.Further,an attempt was made to expand this method in time series.Finally,a preliminary validation of the Moderate Resolution Imaging Spectroradiometer albedo product was performed as the tentative application.

Rock Plasticity from Microtomography and Upscaling

We present a workflow for upscaling of rock properties using microtomography and percolation theory. In this paper we focus on a pilot study for assessing the plastic strength of rocks from a digital rock image. Firstly, we determine the size of mechanical representative volume ele- ment （RVE） by using upper/lower bound dissipation computations in accordance with thermody- namics. Then the mechanical RVE is used to simulate the rock failure at micro-scale using FEM. Two cases of different pressures of linear Drucker-Prager plasticity of rocks are computed to compute the macroscopic cohesion and the angle of internal friction of the rock. We also detect the critical exponents of yield stress for sealing laws from a series of derivative models that are created by a shrinking/expanding algorithm. We use microtomographic data sets of two carbonate samples and compare the results with previous results. The results show that natural rock samples with irregular structures may have the critical exponent of yield stress different from random models. This unexpected result could have significant ramifications for assessing the stability of solid materials with internal structure. Therefore our pilot study needs to be extended to investigate the scaling laws of strength of many more natural rocks with irregular microstructure.

Forward a spatio-temporal trend surface for long-term groundmeasured albedo upscaling over heterogeneous land surface

Upscaling ground albedo is challenged by the serious discrepancy between the heterogeneity of land surfaces and the small number of ground-based measurements.Conventional ground-based measurements cannot provide sufficient information on the characteristics of surface albedo at the scale of coarse pixels over heterogeneous land surfaces.One method of overcoming this problem is to introduce high-resolution albedo imagery as ancillary information for upscaling.However,due to the low frequency of updating of high-resolution albedo maps,upscaling time series of ground-based albedo measurements is difficult.This paper proposes a method that is based on the idea of conceptual universal scaling methodology for establishing a spatiotemporal trend surface using very few high-resolution images and time series of ground-based measurements for spatial-temporal upscaling of albedo.The construction of the spatiotemporal trend surface incorporates the spatial information provided by auxiliary remote sensing images and the temporal information provided by long time series of ground observations.This approach was illustrated by upscaling ground-based fine-scale albedo measurements to a coarse scale over the core study area in HiWATER.The results indicate that this method can characterize the spatiotemporal variations in surface albedo well,and the overall correlation coefficient was 0.702 during the study period.

Upscaling the use of fallout radionuclides in soil erosion and sediment budget investigations:Addressing the challenge

The application of fallout radionuclides in soil erosion investigations and related sediment budget studies has provided a widely used tool for improving understanding of soil erosion and sediment transfer processes.However,most studies using fallout radionuclides undertaken to date have focussed on small areas.This focus on small areas reflects both the issues addressed and practical constraints associated with sample collection and analysis.Increasing acceptance of the important role of fine sediment in degrading aquatic habitats and in the transfer and fate of nutrients and contaminants within terrestrial and fluvial systems has emphasised the need to consider larger areas and the catchment or regional scale.The need to upscale existing approaches to the use of fallout radionuclides to larger areas represents an important challenge.This contribution provides a brief review of existing and potential approaches to upscaling the use of fallout radionuclides and presents two examples where such approaches have been successfully applied.These involve a national scale assessment of soil erosion rates in England and Wales based on 137Cs measurements and an investigation of the sediment budgets of three small/intermediate-size catchments in southern Italy.

Least Squares One-Class Support Tensor Machine

One-class classification problem has become a popular problem in many fields, with a wide range of applications in anomaly detection, fault diagnosis, and face recognition. We investigate the one-class classification problem for second-order tensor data. Traditional vector-based one-class classification methods such as one-class support vector machine (OCSVM) and least squares one-class support vector machine (LSOCSVM) have limitations when tensor is used as input data, so we propose a new tensor one-class classification method, LSOCSTM, which directly uses tensor as input data. On one hand, using tensor as input data not only enables to classify tensor data, but also for vector data, classifying it after high dimensionalizing it into tensor still improves the classification accuracy and overcomes the over-fitting problem. On the other hand, different from one-class support tensor machine (OCSTM), we use squared loss instead of the original loss function so that we solve a series of linear equations instead of quadratic programming problems. Therefore, we use the distance to the hyperplane as a metric for classification, and the proposed method is more accurate and faster compared to existing methods. The experimental results show the high efficiency of the proposed method compared with several state-of-the-art methods.

Numerical investigation of dual-porosity model with transient transfer function based on discrete-fracture model

Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model （DPM）. The comparisons results show that the DPM overestimates the inter-porosity coefficient ）, and the storage factor w. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density. The pressure propagation is slower in the direction of larger fracture density.

Recent progress in laser texturing of battery materials: a review of tuning electrochemical performances, related material development, and prospects for large-scale manufacturing

Traditional electrode manufacturing for lithium-ion batteries is well established,reliable,and has already reached high processing speeds and improvements in production costs.For modern electric vehicles,however,the need for batteries with high gravimetric and volumetric energy densities at cell level is increasing;and new production concepts are required for this purpose.During the last decade,laser processing of battery materials emerged as a promising processing tool for either improving manufacturing flexibility and product reliability or enhancing battery performances.Laser cutting and welding already reached a high level of maturity and it is obvious that in the near future they will become frequently implemented in battery production lines.This review focuses on laser texturing of electrode materials due to its high potential for significantly enhancing battery performances beyond state-of-the-art.Technical approaches and processing strategies for new electrode architectures and concepts will be presented and discussed with regard to energy and power density requirements.The boost of electrochemical performances due to laser texturing of energy storage materials is currently proven at the laboratory scale.However,promising developments in high-power,ultrafast laser technology may push laser structuring of batteries to the next technical readiness level soon.For demonstration in pilot lines adapted to future cell production,process upscaling regarding footprint area and processing speed are the main issues as well as the economic aspects with regards to CapEx amortization and the benefits resulting from the next generation battery.This review begins with an introduction of the three-dimensional battery and thick film concept,made possible by laser texturing.Laser processing of electrode components,namely current collectors,anodes,and cathodes will be presented.Different types of electrode architectures,such as holes,grids,and lines,were generated;their impact on battery performances are illustrated.The usage of high-energy materials,which are on the threshold of commercialization,is highlighted.Battery performance increase is triggered by controlling lithium-ion diffusion kinetics in liquid electrolyte filled porous electrodes.This review concludes with a discussion of various laser parameter tasks for process upscaling in a new type of extreme manufacturing.

Study on a Dual Embedded Discrete Fracture Model for Fluid Flow in Fractured Porous Media

Simulation of fluid flow in the fractured porous media is very important and challenging.Researchers have developed some models for fractured porous media.With the development of related research in recent years,the prospect of embedded discrete fracture model(EDFM)is more and more bright.However,since the size of the fractures in the actual reservoir varies greatly,a very fine grid should be used which leads to a huge burden to the computing resources.To address this challenge,in the present paper,an upscaling based model is proposed.In this model,the flow in large-scale fractures is directly described by the EDFM while that in the small-scale fractures is upscaled through local simulation by EDFM.The EDFM is used to simulate the large-and small-scale fractures independently two times,so the new model is called dual embedded discrete fracture model(D-EDFM).In this paper,the detailed implementation process of D-EDFM is introduced and,through test cases,it is found the proposed model is a feasible method to simulate the flow in fractured porous media.

Dendritic spine density changes and homeostatic synaptic scaling:a meta-analysis of animal studies

Mechanisms of homeostatic plasticity promote compensatory changes of cellular excitability in response to chronic changes in the network activity.This type of plasticity is essential for the maintenance of brain circuits and is involved in the regulation of neural regeneration and the progress of neurodegenerative disorders.One of the most studied homeostatic processes is synaptic scaling,where global synaptic adjustments take place to restore the neuronal firing rate to a physiological range by the modulation of synaptic receptors,neurotransmitters,and morphology.However,despite the comprehensive literature on the electrophysiological properties of homeostatic scaling,less is known about the structural adjustments that occur in the synapses and dendritic tree.In this study,we performed a meta-analysis of articles investigating the effects of chronic network excitation(synaptic downscaling)or inhibition(synaptic upscaling)on the dendritic spine density of neurons.Our results indicate that spine density is consistently reduced after protocols that induce synaptic scaling,independent of the intervention type.Then,we discuss the implication of our findings to the current knowledge on the morphological changes induced by homeostatic plasticity.