; Oh+'0HP
$TU Delft Repository search results0TU Delft Repository search results (max. 1000)TU Delft LibraryTU Delft Library@T̞@T̞՜.+,0HPX`hp
x
WorksheetFeuilles de calcul
B=%r8X"1Calibri1Calibri1Calibri1
Calibri 83ffff̙̙3f3fff3f3f33333f33333.bTU Delft Repositoryg euuidrepository linktitleauthorcontributorpublication yearabstract
subject topiclanguagepublication type publisherisbnissnpatent
patent statusbibliographic noteaccess restrictionembargo datefaculty
departmentresearch group programmeprojectcoordinates)uuid:60193e69d6b44ea8bef862f71ddcf7f0Dhttp://resolver.tudelft.nl/uuid:60193e69d6b44ea8bef862f71ddcf7f0[Adjointbased adaptive convergence control of the iterative finite volume multiscale method
de Zeeuw, W. (Student TU Delft); Jesus de Moraes, R. (TU Delft Petroleum Engineering; Petrobras Research & Development Center); Heemink, A.W. (TU Delft Mathematical Physics); Jansen, J.D. (TU Delft Civil Engineering and Geosciences; TU Delft Geoscience and Engineering)Klie, Hector (editor)We propose a novel adaptive, adjointbased, iterative multiscale finite volume (iMSFV) method. The method aims to reduce the computational cost of the smoothing stage of the original iMSFV method by selectively choosing finescale subdomains (or subset of primary variables) to solve for. The selection of finescale primary variables is obtained from a goaloriented adjoint model. An adjointbased indicator is utilized as a criterion to select the primary variables having the largest errors. The Lagrange multipliers from the adjoint model can be interpreted as sensitivities of the objective function value with respect to deviations from the constraints. In case of adjoining the porous media flow equations with Lagrange multipliers, this implies that the multipliers are the sensitivities of the objective function with respect to the residuals of the flow equations, i.e., to the residual error that remains after approximately solving linear equations with the aid of an iterative solver. This allow us to recognize at which locations the solution contains more errors. More specifically, we propose a modification to the iMSFV method to adaptively reduce the size of the finescale system that must be smoothed. The aim is to make the finescale smoothing stage less computationally demanding. To that end, we introduce a goaloriented, adjointbased finescale system reduction criterion. We demonstrate the performance of our method via singlephase, incompressible flow simulation models with challenging geological settings and using a historymatching like misfit objective function as the goal. The performance of the newly introduced method is compared to the original iMSFV method. We investigate the adaptivity versus accuracy of the method and demonstrate how the solution accuracy varies by varying the number of unknowns selected to be smoothed. It is shown that the method can provide accurate solutions at reduced computational cost. The proofofconcept applications indicate that the method deserves further investigations.enconference paperSociety of Petroleum Engineers9781613996348EGreen Open Access added to TU Delft Institutional Repository You share, we take care! Taverne project https://www.openaccess.nl/en/yousharewetakecare Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public.
20191011)uuid:604d8ddd862046659acadb4e9c8451c4Dhttp://resolver.tudelft.nl/uuid:604d8ddd862046659acadb4e9c8451c4^An Efficient Robust Optimization Workflow using Multiscale Simulation and Stochastic GradientsJesus de Moraes, R. (TU Delft Petroleum Engineering; Petrobras); Fonseca, Rahul Mark (TNO); Helici, Mircea A. (TNO); Heemink, A.W. (TU Delft Mathematical Physics); Jansen, J.D. (TU Delft Geoscience and Engineering)We present an efficient workflow that combines multiscale (MS) forward simulation and stochastic gradient computation  MSStoSAG  for the optimization of well controls applied to waterflooding under geological uncertainty. A twostage iterative Multiscale Finite Volume (iMSFV), a mass conservative reservoir simulation strategy, is employed as the forward simulation strategy. MS methods provide the ability to accurately capture fine sc< ale heterogeneities, and thus the finescale physics of the problem, while solving for the primary variables in a more computationally efficient coarsescale simulation grid. In the workflow, the construction of the basis fuctions is performed at an offline stage and they are not reconstructed/updated throughout the optimization process. Instead, inaccuracies due to outdated basis functions are addressed by the iMSFV smoothing stage. The Stochastic Simplex Approximate Gradient (StoSAG) method, a stochastic gradient technique is employed to compute the gradient of the objective function using forward simulation responses. Our experiments illustrate that iMSFV simulations provide accurate forward simulation responses for the gradient computation, with the advantage of speeding up the workflow due to faster simulations. Speedups up to a factor of five on the forward simulation, the most computationally expensive step of the optimization workflow, were achieved for the examples considered in the paper. Additionally, we investigate the impact of MS parameters such as coarsening ratio and heterogeneity contrast on the optimization process. The combination of speed and accuracy of MS forward simulation with the flexibility of the StoSAG technique allows for a flexible and efficient optimization workflow suitable for largescale problems.YGradientbased optimization; multiscale methods; Robust optimization; Stochastic gradientjournal article)uuid:cadf667d49284c8cb0f940ffed1b0f49Dhttp://resolver.tudelft.nl/uuid:cadf667d49284c8cb0f940ffed1b0f49KIterative multiscale gradient computation for heterogeneous subsurface flow.Jesus de Moraes, R. (TU Delft Petroleum Engineering; Petrobras Research & Development Center); R. P. Rodrigues, Jos (Petrobras Research & Development Center); Hajibeygi, H. (TU Delft Petroleum Engineering); Jansen, J.D. (TU Delft Civil Engineering and Geosciences; TU Delft Geoscience and Engineering)yWe introduce a semianalytical iterative multiscale derivative computation methodology that allows for error control and reduction to any desired accuracy, up to finescale precision. The model responses are computed by the multiscale forward simulation of flow in heterogeneous porous media. The derivative computation method is based on the augmentation of the model equation and state vectors with the smoothing stage defined by the iterative multiscale method. In the formulation, we avoid additional complexity involved in computing partial derivatives associated to the smoothing step. We account for it as an approximate derivative computation stage. The numerical experiments illustrate how the newly introduced derivative method computes misfit objective function gradients that converge to finescale one as the iterative multiscale residual converges. The robustness of the methodology is investigated for test cases with high contrast permeability fields. The iterative multiscale gradient method casts a promising approach, with minimal accuracyefficiency tradeoff, for largescale heterogeneous porous media optimization problems.Adjoint method; Direct method; Gradient computation; Iterative multiscale finite volume; Multiscale methods; Subsurface flowAccepted Author Manuscript)uuid:31b1847ee32c482e9e8f286de866e751Dhttp://resolver.tudelft.nl/uuid:31b1847ee32c482e9e8f286de866e751OMultiscale Analytical Derivative Formulations for Improved Reservoir Management4Jesus de Moraes, R. (TU Delft Petroleum Engineering)qJansen, J.D. (promotor); Hajibeygi, H. (copromotor); Delft University of Technology (degree granting institution)The exploitation of subsurface resources is, inevitably, surrounded by uncertainty. Limited knowledge on the economical, operational, and geological setting are just a few instances of sources of uncertainty. From the geological point of view, the currently available technology is not able to provide the description of the fluids and rock properties at the necessary level of detail required by the mathematical models utilized in the exploitation decisionmaking process. However, even if a fu< ll, accurate description of the subsurface was available, the outcome of such hypothetical mathematical model would likely be computationally too expensive to be evaluated considering the currently available computational power, hindering the decision making process.<br/><br/>Under this reality, geoscientists are consistently making effort to improve the mathematical models, while being inherently constrained by uncertainty, and to find more efficient ways to computationally solve these models.<br/><br/>Closedloop Reservoir Management (CLRM) is a workflow that allows the continuous update of the subsurface models based on production data from different sources. It relies on computationally demanding optimization algorithms (for the assimilation of production data and control optimization) which require multiple simulations of the subsurface model. One important aspect for the successful application of the CLRM workflow is the definition of a model that can both be run multiple times in a reasonable timespan and still reasonably represent the underlying physics. <br/><br/>Multiscale (MS) methods, a reservoir simulation technique that solves a coarser simulation model, thus increasing the computational speed up, while still utilizing the finescale representation of the reservoir, figures as an accurate and efficient simulation strategy.<br/><br/>This thesis focuses on the development of efficient algorithms for subsurface models optimization by taking advantage of multiscale simulation strategies. It presents (1) multiscale analytical derivative computation strategies to efficiently and accurately address the optimization algorithms employed in the CLRM workflow and (2) novel strategies to handle the mathematical modeling of subsurface management studies from a multiscale perspective. On the latter, we specifically address a more fundamental multiscale aspect of data assimilation studies: the assimilation of observations from a distinct spatial representation compared to the simulation model scale.<br/><br/>As a result, this thesis discusses in detail the development of mathematical models and algorithms for the derivative computation of subsurface model responses and their application into gradientbased optimization algorithms employed in the data assimilation and lifecycle optimization steps of CLRM. The advantages are improved computational efficiency with accuracy maintenance and the ability to address the subsurface management from a multiscale view point not only from the forward simulation perspective, but also from the inverse modeling side.tmultiscale simulation; analytical derivative computation; adjoint method; lifecycle optimization; data assimilationdoctoral thesis9789461869906Petroleum Engineering)uuid:dc72b5467dc14c0f9352e2455c7811eeDhttp://resolver.tudelft.nl/uuid:dc72b5467dc14c0f9352e2455c7811eeJComputing derivative information of sequentially coupled subsurface modelsJesus de Moraes, R. (TU Delft Petroleum Engineering; Petrobras); Rodrigues, Jos R.P. (Petrobras); Hajibeygi, H. (TU Delft Petroleum Engineering); Jansen, J.D. (TU Delft Civil Engineering and Geosciences; TU Delft Geoscience and Engineering)A generic framework for the computation of derivative information required for gradientbased optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic waterflooding model, where the forward model is constructed based on a sequentially coupled flowtransport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and lifecycle optimization. Results are compared with the classical fully c< oupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradientbased optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.{Adjoint method; Data assimilation; Direct method; Gradientbased optimization; Lifecycle optimization; Sequential coupling!Civil Engineering and GeosciencesGeoscience and Engineering)uuid:0dd36cc6797c4d8d928e7d755ca17220Dhttp://resolver.tudelft.nl/uuid:0dd36cc6797c4d8d928e7d755ca17220)A Multiscale Method For Data AssimilationJesus de Moraes, R. (TU Delft Petroleum Engineering); Hajibeygi, H. (TU Delft Petroleum Engineering); Jansen, J.D. (TU Delft Civil Engineering and Geosciences; TU Delft Geoscience and Engineering)Gunasekera, D. (editor)In data assimilation problems, various types of data are naturally linked to different spatial resolutions (e.g. seismic and electromagnetic data), and these scales are usually not coincident to the subsurface simulation model scale. Alternatives like down/upscaling of the data and/or the simulation model can be used, but with potential loss of important information. To address this issue, a novel Multiscale (MS) data assimilation method is introduced. The overall idea of the method is to keep uncertain parameters and observed data at their original representation scale, avoiding down/upscaling of any quantity. The method relies on a recently developed mathematical framework to compute adjoint gradients via a MS strategy. The finescale uncertain parameters are directly updated and the MS grid is constructed in a resolution that meets the observed data resolution. The advantages of the technique are demonstrated in the assimilation of data represented at a coarser scale than the simulation model. The misfit objective function is constructed to keep the MS nature of the problem. The regularization term is represented at the simulation model (fine) scale, whereas the data misfit term is represented at the observed data (coarse) scale. The performance of the method is demonstrated in synthetic models and compared to down/upscaling strategies. The experiments show that the MS strategy provides advantages 1) on the computational side expensive operations are only performed at the coarse scale; 2) with respect to accuracy the matched uncertain parameter distribution is closer to the truth ; and 3) in the optimization performance faster convergence behaviour due to faster gradient computation. In conclusion, the newly developed method is capable of providing superior results when compared to strategies that rely on the up/downscaling of the response/observed data, addressing the scale dissimilarity via a robust, consistent MS strategy.EAGE
20190303)uuid:ace246d6efa747e8994214d4d75ef04cDhttp://resolver.tudelft.nl/uuid:ace246d6efa747e8994214d4d75ef04cFMultiscale gradient computation for flow in heterogeneous porous mediaJesus de Moraes, R. (TU Delft Petroleum Engineering); Rodrigues, Jos R P (Petrobras); Hajibeygi, H. (TU Delft Petroleum Engineering); Jansen, J.D. (TU Delft Civil Engineering and Geosciences; TU Delft Geoscience and Engineering)An efficient multiscale (MS) gradient computation method for subsurface flow management and optimization is introduced. The general, algebraic framework allows for the calculation of gradients using both the Direct and Adjoint derivative methods. The framework also allows for the utilization of any MS formulation that can be algebraically expressed in terms of a restriction and a prolongation operator. This is achieved via an implicit differentiation formulation. The approach favors algorithms for multiplying the sensitivity matrix and its transpose with arbitrary vectors. This provides a flexible way of computing gradients in a form suitable for any given gradientbased optimization algorithm. No assumption w.r.t. the nature of the problem< or specific optimization parameters is made. Therefore, the framework can be applied to any gradientbased study. In the implementation, extra partial derivative information required by the gradient computation is computed via automatic differentiation. A detailed utilization of the framework using the MS Finite Volume (MSFV) simulation technique is presented. Numerical experiments are performed to demonstrate the accuracy of the method compared to a finescale simulator. In addition, an asymptotic analysis is presented to provide an estimate of its computational complexity. The investigations show that the presented method casts an accurate and efficient MS gradient computation strategy that can be successfully utilized in nextgeneration reservoir management studies.iAdjoint method; Automatic differentiation; Direct method; Gradientbased optimization; Multiscale methods
20170814)uuid:926df424834344d080d1255c1553f168Dhttp://resolver.tudelft.nl/uuid:926df424834344d080d1255c1553f168qImproving the computational efficiency of approximate gradients using a multiscale reservoir simulation frameworkJesus de Moraes, R. (TU Delft Petroleum Engineering; Petrobras); Fonseca, R.M. (TNO); Helici, M.; Heemink, A.W. (TU Delft Mathematical Physics); Jansen, J.D. (TU Delft Civil Engineering and Geosciences; TU Delft Geoscience and Engineering)0In this work, the application of tensor methodologies for computerassisted history matching of channelized reservoirs is explored. A tensorbased approach is used for the parameterization of petrophysical parameters to reduce the dimensionality of the parameter estimation problem. Building on the work of Afra and Gildin (2013); Afra et.al. (2014); Afra and Gildin (2016), permeability fields of multiple model realizations are collected in a tensor form which is subsequently decomposed to derive a lowdimensional representation of the dominant spatial structures in the models. This representation then is used to estimate an identifiable reduced set of parameters using an ensemble Kalman filter (EnKF) strategy. This approach is attractive for the parameter estimation of permeabilities because it increases the ability to represent channelized structures in the updates resulting in an improved predictive capacity of the historymatched models. In particular, channel continuity is better preserved than with a Principal Component Analysis (PCA) parameterization.SPE
20170822)uuid:c0120e5f68464d7faf22733347ebdee2Dhttp://resolver.tudelft.nl/uuid:c0120e5f68464d7faf22733347ebdee2CMultiscale gradient computation for multiphase flow in porous mediaJesus de Moraes, R. (TU Delft Petroleum Engineering); Rodrigues, J. R.P. (Petrobras); Hajibeygi, H. (TU Delft Petroleum Engineering); Jansen, J.D. (TU Delft Civil Engineering and Geosciences; TU Delft Geoscience and Engineering)wA multiscale gradient computation method for multiphase flow in heterogeneous porous media is developed. The method constructs multiscale primal and dual coarse grids, imposed on the given finescale computational grid. Local multiscale basis functions are computed on (dual) coarse blocks, constructing an accurate map (prolongation operator) between coarse and finescale systems. While the expensive operations involved in computing the gradients are performed at the coarse scale, sensitivities with respect to uncertain parameters (e.g., grid block permeabilities) are expressed in the fine scale via the partial derivatives of the prolongation operator. Hence, the method allows for updating of the geological model, rather than the dynamic model only, avoiding upscaling and the inevitable loss of information. The formulation and implementation are based on automatic differentiation (AD), allowing for convenient extensions to complex physics. An IMPES coupling strategy for flow and transport is followed, in the forward simulation. The flow equation is computed using a multiscale finite volume (MSFV) formulation and the transport equation is computed at the fine scale, after reconstruction of mass conservative velocity field. To< assess the performance of the method, a synthetic multiphase flow test case is considered. The multiscale gradients are compared against those obtained from a finescale reference strategy. Apart from its computational efficiency, the benefits of the method include flexibility to accommodate variables expressed at different scales, specially in multiscale data assimilation and reservoir management studies.)uuid:facd9e6db58a4aa09e0dff1edde1bfb1Dhttp://resolver.tudelft.nl/uuid:facd9e6db58a4aa09e0dff1edde1bfb1:Multiscale Gradient Computation for Subsurface Flow ModelsJesus de Moraes, R. (TU Delft Petroleum Engineering); Rodrigues, J.R.P. (Petrobras); Hajibeygi, H. (TU Delft Petroleum Engineering); Jansen, J.D. (TU Delft Civil Engineering and Geosciences; TU Delft Geoscience and Engineering)7We present an efficient multiscale (MS) gradient computation that is suitable for reservoir management studies involving optimization techniques for, e.g., computerassisted history matching or lifecycle production optimization. The general, algebraic framework allows for the calculation of gradients using both the Direct and Adjoint derivative methods. The framework also allows for the utilization of any MS formulation in the forward reservoir simulation that can be algebraically expressed in terms of a restriction and a prolongation operator. In the implementation, extra partial derivative information required by the gradient methods is computed via automatic differentiation. Numerical experiments demonstrate the accuracy of the method compared against those based on finescale simulation (industry standard).
20170301
*+&ffffff?'ffffff?(?)?"dXX333333?333333?U}}}}}}}}}} }
}}}
}}}}}}}}}}}}
@
!
"
#
$
%
&
'@
(
)
*
+
,

.@
/
0
*
1
2
3
4
5
6@
7
8
9
:
;
<
=
>
?@
@
A
*
B
C
;
D
E
F
G
H@
I
J
"
K
B
C
;
L
M
N
O@
P
Q
*
"
R
B
C
;
S
T
U
V@
W
X
"
Y
B
C
;
Z
[
\
] @
^
"
Y
B
C
;
_
`
a
b
@
c
J
"
d
B
C
;
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~>@dd
yKyKhttp://resolver.tudelft.nl/uuid:60193e69d6b44ea8bef862f71ddcf7f0yKyKhttp://resolver.tudelft.nl/uuid:604d8ddd862046659acadb4e9c8451c4yKyKhttp://resolver.tudelft.nl/uuid:cadf667d49284c8cb0f940ffed1b0f49yKyKhttp://resolver.tudelft.nl/uuid:31b1847ee32c482e9e8f286de866e751yKyKhttp://resolver.tudelft.nl/uuid:dc72b5467dc14c0f9352e2455c7811eeyKyKhttp://resolver.tudelft.nl/uuid:0dd36cc6797c4d8d928e7d755ca17220yKyKhttp://resolver.tudelft.nl/uuid:ace246d6efa747e8994214d4d75ef04cyKyKhttp://resolver.tudelft.nl/uuid:926df424834344d080d1255c1553f168 yKyKhttp://resolver.tudelft.nl/uuid:c0120e5f68464d7faf22733347ebdee2
yKyKhttp://resolver.tudelft.nl/uuid:facd9e6db58a4aa09e0dff1edde1bfb1gg
Root Entry FT̞T̞@SummaryInformation( F<Workbook FDocumentSummaryInformation8 F
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~