Marion Parquer

Na Yan

Luca Colombera

Nigel P Mountney

Pauline Collon

Guillaume Caumon

Modelling fluvial meander-belt deposits

Keywords: fluvial, channel, point bar, facies, abandoned meander, simulation

The sedimentary record of meandering rivers contains a diverse and complex set of lithological heterogeneities, which impact natural resource management. Different methods exist to model such accumulated successions present in the subsurface by integrating knowledge of system evolutionary behaviour and geometries visible on seismic time or stratal slices. With reference to case-study examples, we review, discuss and employ two of these methods: (i) ChaRMigS generates possible scenarios for channel evolution and meander cut-offs by a reverse migration process; (ii) PB-SAND is a forward stratigraphic model which simulates fluvial point-bar geometry and facies distributions from known palaeo-channel geometries. We introduce a workflow to demonstrate how these two methods can be applied in combination to predict fluvial meander-belt facies distributions, using a subsurface dataset on a Pleistocene succession from the Gulf of Thailand where abandoned channels are visible on seismic time slices, but for which bar-accretion geometries and the exact timing of channel abandonment are unclear. Results show the value of a combined modelling approach to automate the stochastic generation of facies distributions constrained by seismic interpretations.

Introduction

The deposits of meandering river systems are important hosts for natural resources including water, hydrocarbons and ore minerals [START_REF] Fitterman | Electromagnetic mapping of buried paleochannels in eastern Abu Dhabi Emirate, UAE[END_REF][START_REF] García-Gil | GIS-supported mapping of low-temperature geothermal potential taking groundwater flow into account[END_REF][START_REF] Mayall | Turbidite channel reservoirs -Key elements in facies prediction and effective development[END_REF]. In subsurface successions, reflection seismic imaging is a primary tool used to document the architecture of reservoirs composed of the deposits of meandering channel belts [START_REF] Maynard | One Million Years from the Upper Arang Formation, West Natuna Basin, Implications for Reservoir Distribution and Facies Variation in Fluvial Deltaic Deposits[END_REF][START_REF] Reijenstein | Seismic geomorphology and highresolution seismic stratigraphy of inner-shelf fluvial, esturarine, deltaic, marine sequences, Gulf of Thaïland[END_REF]. Depending on their resolution, seismic images commonly reveal (i) the last (pre-avulsion) channel path, (ii) remnant meander-loop fragments delineated by abandoned channel elements [START_REF] Lonsdale | Cut-offs at an abyssal meander south of Iceland[END_REF][START_REF] Erskine | Alluvial cutoffs as indicators of former channel conditions[END_REF], and (iii) remnant point-bar fragments, some with associated accretion patterns evident in planview (e.g., [START_REF] Maynard | One Million Years from the Upper Arang Formation, West Natuna Basin, Implications for Reservoir Distribution and Facies Variation in Fluvial Deltaic Deposits[END_REF][START_REF] Hubbard | Seismic geomorphology and sedimentology of a tidally influenced river deposit, Lower Cretaceous Athabasca oil sands, Alberta, Canada[END_REF]Durkin et al., 2018a). Abandoned channel paths are commonly imaged clearly on seismic time or stratal slices, since they tend to be infilled with fine-grained sediment (mudstone) that exhibits a marked acoustic-impedance contrast relative to adjacent sand-rich point-bar deposits [START_REF] Posamentier | Seismic geomorphology and stratigraphy of depositional elements in deep-water settings[END_REF][START_REF] Veeken | Seismic stratigraphy, basin analysis and reservoir characterisation[END_REF][START_REF] Issautier | Impacts of fluvial reservoir heterogeneity on connectivity : Implications in estimating geological storage capacity for CO 2[END_REF][START_REF] Miall | Fluvial depositional systems[END_REF]. Assessing the juxtaposition of preserved architectural elements is key to understand the history of channel-belt evolution. The nature of channel cut-offs and the partial overprinting of abandoned channel and bar elements can commonly be used to reconstruct the relative chronology of the evolution of a channel belt (Durkin et al., 2018a;[START_REF] Parquer | Integration of channel meander abandonment age uncertainty into a stochastic channelized system reconstruction method[END_REF]. Where architectural features of meandering channel belts -for example architectural elements composed of palaeo-channel loops -can be observed in planform, it is possible to use such observations to condition models for prediction of facies architectures. This can be done in the frame of stochastic modelling.

The deposits of meandering fluvial systems are characterised by strong lithological heterogeneities at different spatial scales, which makes prediction of subsurface sedimentary architecture and natural resource characterization challenging [START_REF] Fustic | Differential entrapment of charged oil -New insights on McMurray Formation oil trapping mechanisms[END_REF]. At the scale of architectural elements, sandprone point-bar elements are commonly partitioned by mud-prone abandoned channel deposits [START_REF] Musial | Subsurface and outcrop characterization of large tidally influenced point bars of the Cretaceous McMurray Formation (Alberta, Canada)[END_REF][START_REF] Colombera | Geometry and compartmentalization of fluvial meander-belt reservoirs at the bar-form scale: Quantitative insight from outcrop, modern and subsurface analogues[END_REF]. These types of architectural elements are themselves composed internally of various lithofacies, which commonly occur in a predictable order. Importantly, for these types of succession, mud-prone packages of strata, or individual mud-prone beds that may only comprise a volumetrically minor part of an overall point-bar element, can significantly partition sand-dominated bar deposits [START_REF] Nanson | Point bar and floodplain formation of the meandering Beatton River, northeastern British Columbia, Canada[END_REF][START_REF] Thomas | Inclined heterolithic stratification-Terminology, description, interpretation and significance[END_REF][START_REF] Pyrcz | ALLUVSIM : A program for event-based stochastic modeling of fluvial depositional systems[END_REF][START_REF] Cojan | Process-based reservoir modelling in the example of meandering channel[END_REF][START_REF] Jackson | Three-dimensional modeling of a shoreface-shelf parasequence reservoir analog: Part 2. Geologic controls on fluid flow and hydrocarbon recovery[END_REF][START_REF] Labrecque | Sedimentology and stratigraphic architecture of a point bar deposit, Lower Cretaceous McMurray Formation, Alberta, Canada[END_REF][START_REF] Colombera | Seismic-driven geocellular modeling of fluvial meander-belt reservoirs using a rule-based method[END_REF]. This is notably the case where thin but laterally continuous and inclined mudstone beds are draped over sandy lateral-accretion units within point-bar elements [START_REF] Yan | Three-dimensional modelling of fluvial point-bar architecture and facies heterogeneity using analogue data and associated analysis of intra-bar static connectivity: application to humid coastal-plain and dryland fluvial systems[END_REF]. Conversely, mud-dominated abandoned channel elements commonly have thin, sand-rich components in their basal-most part, which permit the flow of fluids through overall low-permeability channel-plug deposits [START_REF] Donselaar | Connectivity of fluvial point-bar deposits: An example from the Miocene Huesca fluvial fan, Ebro Basin, Spain[END_REF]. Consequently, the internal distribution of lithofacies inside these bodies exerts a major control on porous fluid-flow pathways [START_REF] Bridge | A revised mathematical model and FORTRAN IV program to predict flow, bed topography, and grain size in open-channel bends[END_REF][START_REF] Brownlie | Flow depth in sand-bed channels[END_REF]Sun et al., 2001a[START_REF] Sun | Meander migration and the lateral tilting of floodplains[END_REF]Willis & Sech, 2019a, 2019b). Knowledge of the expected internal facies arrangement within architectural elements is therefore an essential requirement for effective groundwater or hydrocarbon reservoir evaluation (e.g., [START_REF] Jackson | Effect of discontinuous shales on reservoir performance during horizontal waterflooding[END_REF][START_REF] Jackson | Three-dimensional modeling of a shoreface-shelf parasequence reservoir analog: Part 2. Geologic controls on fluid flow and hydrocarbon recovery[END_REF]. However, in many subsurface settings, internal facies arrangements and proportions are at best only poorly constrained by sparse wells that provide wireline or core data [START_REF] Miall | Reconstructing fluvial macroform architecture from two-dimensional outcrops; examples from the Castlegate Sandstone, Book Cliffs, Utah[END_REF][START_REF] Kostic | Sedimentary architecture and 3D ground-penetrating radar analysis of gravelly meandering river deposits (Neckar Valley, SW Germany)[END_REF][START_REF] Musial | Subsurface and outcrop characterization of large tidally influenced point bars of the Cretaceous McMurray Formation (Alberta, Canada)[END_REF]. In such cases, databases are commonly employed to provide information about the distributions of facies known from analogous modern rivers and outcropping successions (e.g., Colombera et al., 2012aColombera et al., , 2012b[START_REF] Colombera | A quantitative approach to fluvial facies models: methods and example results[END_REF]Colombera & Mountney, 2019). Such databases are used to inform reservoir models that seek to depict the expected distribution of facies and elements within subsurface successions [START_REF] Arnold | Uncertainty Quantification in Reservoir Prediction: Part 1 Model Realism in History Matching Using Geological Prior Definitions[END_REF].

In general, these static reservoir models represent the internal anatomy, lithological configuration and petrophysical properties of hydrocarbon reservoirs as cell-based corner-point reservoir grids [START_REF] Farmer | Geological modelling and reservoir simulation[END_REF][START_REF] Ringrose | Reservoir model design[END_REF]. For each reservoir, or sector thereof, the uncertainty relating to geological and petrophysical features is usually addressed by generating multiple equiprobable models built using stochastic modelling methods [START_REF] Caumon | Geological Objects and Physical Parameter Fields in the Subsurface: A Review[END_REF]. The approaches that are typically used to model the distributions of meander-belt facies associations in subsurface fluvial successions consist of general-purpose geostatistical modelling tools, such as those based on two-point statistics (e.g., sequential indicator simulations; [START_REF] Journel | Focusing on spatial connectivity of extreme-valued attributes: stochastic indicator models of reservoir heterogeneities[END_REF][START_REF] Deutsch | GSLIB : Geostatistical Software Library and User's Guide[END_REF] or on multi-point statistics (e.g., [START_REF] Srivastava | An overview of stochastic methods for reservoir characterization[END_REF][START_REF] Strebelle | Conditional simulation of complex geological structures using multiple-point statistics[END_REF][START_REF] Straubhaar | Addressing Conditioning Data in Multiple-Point Statistics Simulation Algorithms Based on a Multiple Grid Approach[END_REF][START_REF] Arnold | Uncertainty Quantification in Reservoir Prediction: Part 1 Model Realism in History Matching Using Geological Prior Definitions[END_REF][START_REF] Calderon | Geological Facies Recovery Based on Weighted ℓ 1 -Regularization[END_REF]. However, pixel-based geostatistical methods commonly struggle to reproduce the geometry and connectivity of complex facies configurations in three dimensions [START_REF] Renard | Conditioning Facies Simulations with Connectivity Data[END_REF][START_REF] Rongier | Comparing Connected Structures in Ensemble of Random Fields[END_REF].

Alternatively, object-based methods can be used to reproduce complex sedimentary architectures through the distribution of objects that represent sedimentary bodies within a modelled volume, for example channel bodies of given geometry (e.g. width, length, sinuosity) [START_REF] Deutsch | Hierarchical Object-Based Stochastic Modeling of Fluvial Reservoirs[END_REF][START_REF] Keogh | The development of fluvial stochastic modelling in the Norwegian oil industry: A historical review, subsurface implementation and future directions[END_REF][START_REF] Hassanpour | Improved geostatistical models of inclined heterolithic strata for McMurray Formation, Alberta, Canada[END_REF][START_REF] Ruiu | Modeling Channel Forms and Related Sedimentary Objects Using a Boundary Representation Based on Non-Uniform Rational B-Splines[END_REF]Rongier et al., 2017a;Colombera et al., 2019). In these methods, each object is randomly sampled and modelled independently from the others combining statistical point processes and Monte Carlo sampling of shape parameters. In these approaches, conditioning to well or seismic data is notably challenging when the observation density is large relative to the sizes of the simulated objects [START_REF] Wang | Conditioning 3D object-based models to dense well data[END_REF]. In classical object-based simulation, architectural elements are also modelled independently one from another, meaning that natural depositional processes are not captured by the simulation algorithm. This shortcoming motivates another family of methods that seeks to mimic the geological processes producing specific types of sedimentary deposits; these are so-called dynamic, event-based or rule-based approaches [START_REF] Pyrcz | Stratigraphic rule-based reservoir modeling[END_REF]. Several variants of this approach have been proposed for modelling the architecture of meandering channel systems, such as those of fluvial and submarine channels. Among these, process-based methods seek to model physical deposition by a forward geometric evolution depending on curvature of meanders (e.g., [START_REF] Ikeda | Bend theory of river meanders. Part 1. Linear development[END_REF][START_REF] Cojan | Process-based reservoir modelling in the example of meandering channel[END_REF][START_REF] Lemay | Geomorphic variability of submarine channelized systems along continental margins: Comparison with fluvial meandering channels[END_REF][START_REF] Sylvester | High curvatures drive river meandering[END_REF]. Descriptive-based methods, by contrast, are concerned with the reproduction of the geometry generated by the sedimentary processes [START_REF] Labourdette | Three-dimensional modelling of stacked turbidite channels in West Africa : impact on dynamic reservoir simulations[END_REF][START_REF] Pyrcz | ALLUVSIM : A program for event-based stochastic modeling of fluvial depositional systems[END_REF][START_REF] Ruiu | Modeling Channel Forms and Related Sedimentary Objects Using a Boundary Representation Based on Non-Uniform Rational B-Splines[END_REF][START_REF] Rongier | A geostatistical approach to the simulation of stacked channels[END_REF]. For these dynamic methods, stratal configurations are modelled with consideration of how the system evolves through time. However, a common challenge for meandering channelized systems is the difficulty in making the models fit closely to fragmentary observations of the channel remnants, such as abandoned meanders that are evident on seismic images [START_REF] Labourdette | Three-dimensional modelling of stacked turbidite channels in West Africa : impact on dynamic reservoir simulations[END_REF][START_REF] Ruiu | Semiautomatic interpretation of 3D sedimentological structures on geologic images: An object-based approach[END_REF], Rongier et al., 2017a; the ChaRMigS method of [START_REF] Parquer | Reconstruction of channelized systems through a conditioned reverse migration method[END_REF]; the PB-SAND method of [START_REF] Yan | A 3D forward stratigraphic model of fluvial meander-bend evolution for prediction of point-bar lithofacies architecture[END_REF]. A solution to this is offered by inverse (i.e., reverse) stratigraphic modelling methods (e.g., ChaRMigS method of [START_REF] Parquer | Reconstruction of channelized systems through a conditioned reverse migration method[END_REF], which honour observations of abandoned meanders by integrating them into a vector-based modelling structure that can be used to determine the temporal evolution of a river channel. In all cases, the question of uncertainty assessment remains central and calls for considering multiple realizations of reservoir geometry [START_REF] Pyrcz | Stratigraphic rule-based reservoir modeling[END_REF][START_REF] Caumon | Geological Objects and Physical Parameter Fields in the Subsurface: A Review[END_REF][START_REF] Arnold | Uncertainty Quantification in Reservoir Prediction: Part 1 Model Realism in History Matching Using Geological Prior Definitions[END_REF]. Overall, these rule-based approaches have been shown to be a promising way to address the shortcomings of pixel-based and object-based methods, in particular towards honouring all quantitative spatial observations while maintaining consistent relationships between the simulated architectural elements. However, only a limited number of these methods (e.g., PB-SAND method of [START_REF] Yan | A 3D forward stratigraphic model of fluvial meander-bend evolution for prediction of point-bar lithofacies architecture[END_REF] permit simulation of the distribution of facies inside the reconstructed architectures.

The aim of this study is to highlight how two numerical modelling approaches -ChaRMigS [START_REF] Parquer | Reconstruction of channelized systems through a conditioned reverse migration method[END_REF] and PB-SAND [START_REF] Yan | A 3D forward stratigraphic model of fluvial meander-bend evolution for prediction of point-bar lithofacies architecture[END_REF], can be jointly employed to model the internal architecture of the accumulated deposits of meandering fluvial successions, with consideration of how this architecture varies in relation to the likely evolutionary scenarios of the formative river system. The modelling approach seeks to account for deposition and erosion over time and space as a fundamental control on preserved sedimentary architecture. Specific research objectives are as follows: (i) to compare the contexts of application, the inputs and outputs, and the relative merits and limitations of two modelling approaches used to predict the sedimentary architecture of meandering fluvial channel belts, based on inverse and forward modelling, respectively; (ii) to illustrate how the two modelling approaches can be linked and used sequentially to obtain several fully stochastic realizations of the facies architecture of large fluvial meander-belt successions whose channel fills are imaged in three-dimensional seismic datasets.

Methods

The ChaRMigS and PB-SAND models both use geometric, vector-based approaches [START_REF] Parquer | Reconstruction of channelized systems through a conditioned reverse migration method[END_REF][START_REF] Yan | A 3D forward stratigraphic model of fluvial meander-bend evolution for prediction of point-bar lithofacies architecture[END_REF]. Characteristics of the two modelling methods are described below and summarized in Table 1.

ChaRMigS approach

The Channel Reverse Migration Simulation (ChaRMigS) model [START_REF] Parquer | Reconstruction of channelized systems through a conditioned reverse migration method[END_REF] uses an inverse stratigraphic modelling approach to propose a wide range of possible reconstructed temporal evolutions of meandering channel planforms, which could have originated via several possible migration patterns. These migration patterns can be the products of meander-bend expansion, translation and rotation [START_REF] Daniel | Channel movement of meandering Indiana streams[END_REF][START_REF] Hickin | The development of meanders in natural river-channels[END_REF][START_REF] Holbrook | Expansion to contraction, translation to retreat, osculation or toggle, and degree of wobble, and the many means by which a meander may move[END_REF], and these processes can be punctuated by channel abandonment, or by neck or chute cut-off of individual meander loops.

ChaRMigS simulates channel migration patterns in the form of time-lapse channel trajectories. Multiple stochastic outputs of temporal channel evolution are generated for each set of input parameters, consisting of the centreline of the youngest active channel and, optionally, of centrelines associated with multiple oxbow-shaped abandoned meander loops, as well as probability distributions of migration offset and cut-off rates. Both distributions can be inferred based on analogues or knowledge of the system under study. The spatial organization of the abandoned meanders is used to determine likely histories of planform evolution of the channel belt, with consideration of gradual channel migration and meanderbend cut-off processes. Abandoned meander loops are integrated into the output scenarios of channel evolution in different orders, which are consistent with potentially different bend cut-off chronologies.

The global chronology of bend cut-off emerges from the simulation process: at every simulation time step, the probability of integrating a particular abandoned meander loop in the channel trajectory is computed, based on the current floodplain geometry; the older channel trajectory corresponding to the next simulation time step then reflects gradual migration and integrates loops which have been randomly selected by Monte Carlo sampling [START_REF] Parquer | Integration of channel meander abandonment age uncertainty into a stochastic channelized system reconstruction method[END_REF]. A history of events can also be imposed when the chronology of abandonment is known [START_REF] Parquer | Reconstruction of channelized systems through a conditioned reverse migration method[END_REF].

The channel trajectories simulated by ChaRMigS appear naturally sinuous, therefore tend to result in irregular streamwise variations in the magnitude of accretion and erosion. In ChaRMigS, temporal and spatial variations in the thickness of channel deposits resulting from along-channel bathymetric variations are included as model parameters, but are not reproduced in the outputs presented herein for the sake of simplicity. Channel-belt aggradation can be modelled as a global or vertical offset of channel centrelines between accretion steps. Establishing a close fit to accretion geometries (e.g., as evidenced by scroll bars) may be achieved in principle by interactive NURBS editing [START_REF] Ruiu | Semiautomatic interpretation of 3D sedimentological structures on geologic images: An object-based approach[END_REF], but has not been considered in the ChaRMigS workflow. Facies simulation is not incorporated in this modelling approach.

PB-SAND approach

The Point-Bar Sedimentary Architecture Numerical Deduction (PB-SAND) model [START_REF] Yan | A 3D forward stratigraphic model of fluvial meander-bend evolution for prediction of point-bar lithofacies architecture[END_REF] uses a forward geometric-based modelling approach. This model has been devised as a tool for the reconstruction and prediction of the complex spatio-temporal evolution of fluvial meanders, their generated 3D lithofacies distributions and resulting heterogeneity. The PB-SAND model permits the reconstruction of point-bar geometries and internal sedimentary architectures using a deterministic approach to simulate accretion patterns as they evolve over a series of time steps [START_REF] Yan | Three-dimensional modelling of fluvial point-bar architecture and facies heterogeneity using analogue data and associated analysis of intra-bar static connectivity: application to humid coastal-plain and dryland fluvial systems[END_REF]. The input trajectories that control the planform morphology of point bars can be digitized from seismic images, from remotely sensed images of modern systems, or devised based on field observations of ancient outcropping successions. Output channel trajectories are linearly interpolated and extrapolated between two input trajectories, and are therefore deterministic (i.e., a single solution is obtained for each set of input trajectories). The approach can be used to model up to three consecutive meander bends in a single run. Multiple model outputs can then be chained together to model larger river reaches [START_REF] Yan | Three-dimensional forward stratigraphic modelling of the sedimentary architecture of meandering-river successions in evolving half-graben rift basins[END_REF]. The accretion patterns between two successive time steps are extrapolated based on the type of meander transformation and a specified number of accretion increments. As a result, different accretion patterns can be generated, arising from any combination of meander expansion, rotation, and downstream translation (cf. [START_REF] Daniel | Channel movement of meandering Indiana streams[END_REF][START_REF] Ghinassi | Plan-form evolution of ancient meandering rivers reconstructed from longitudinal outcrop sections[END_REF][START_REF] Holbrook | Expansion to contraction, translation to retreat, osculation or toggle, and degree of wobble, and the many means by which a meander may move[END_REF].

The PB-SAND model additionally permits the stochastic simulation of the internal facies arrangement of point bars developed through different modes of meander transformation (expansion, translation, rotation). This is done using a rule-based method that is derived from and constrained by known examples of fluvial meandering systems [START_REF] Yan | Three-dimensional modelling of fluvial point-bar architecture and facies heterogeneity using analogue data and associated analysis of intra-bar static connectivity: application to humid coastal-plain and dryland fluvial systems[END_REF]. Trends of fining upwards, fining outwards through meander expansion, and fining downstream beyond a meander apex and in concave-bank areas, can all be reproduced, for example (see details in [START_REF] Yan | A 3D forward stratigraphic model of fluvial meander-bend evolution for prediction of point-bar lithofacies architecture[END_REF]. Parameters that control the geometry of facies transitions within accretion packages (e.g., inclination of boundaries) are stochastically drawn from predefined probability curves (e.g., Gaussian distributions; see Fig. A.1 in Yan et al., 2017). Variations in the thickness of channel deposits might result from streamwise bathymetric variations (from riffle to pool) and temporal changes (in sinuosity or hydraulic geometry). Aggradation of the channel floor can also be modelled [START_REF] Yan | Three-dimensional forward stratigraphic modelling of the sedimentary architecture of meandering-river successions in evolving half-graben rift basins[END_REF]. In PB-SAND, incorporation of abandoned cut-off meander loops is achieved through integration of the outputs of separate modelling steps. It is difficult to incorporate abandoned meander loops specifically in a temporal sequence, because the evolutionary trajectories of the point-bar elements (for example, as sometimes recorded in accretion patterns observed in seismic slices) commonly have been heavily overprinted by more recently developed point-bar elements. The modelling approach of PB-SAND can therefore be used to generate stochastic models of the internal facies architecture of meander belts, but only in cases where the planview accretion geometries of the channel belt can be inferred with confidence, as is sometimes the case with 3D seismic datasets [START_REF] Colombera | Seismic-driven geocellular modeling of fluvial meander-belt reservoirs using a rule-based method[END_REF][START_REF] Yan | Three-dimensional modelling of fluvial point-bar architecture and facies heterogeneity using analogue data and associated analysis of intra-bar static connectivity: application to humid coastal-plain and dryland fluvial systems[END_REF].

FAKTS database

The Fluvial Architecture Knowledge Transfer System (FAKTS; Colombera et al., 2012b) is a relational database storing descriptions of the sedimentary architecture of many modern river systems and ancient fluvial successions, themselves classified on several parameters (e.g., climate, age, drainage area). FAKTS currently contains data from over 300 sedimentological case studies, and incorporates data derived from both original field studies and the published literature. FAKTS provides a quantitative characterization of sedimentary architectures at multiple scales (e.g., lithofacies, point bars, channel belts). This database serves to provide input parameters to PB-SAND, relating to the geometry of architectural elements (e.g., bar and channel-fill dimensions) and to the proportion and geometry of lithofacies (e.g., proportion, thickness and length of mud drapes).

Applied comparisons of ChaRMigS and PB-SAND

Channel belt of the Oulanka River (Finland)

To demonstrate an application of the two modelling approaches considered in this work, ChaRMigS and PB-SAND were applied separately to model the planform evolution of a stretch of the meandering channel belt of the present-day Oulanka River (Finland), which occupies a pre-glacial valley [START_REF] Koutaniemi | Some aspects of the palaeohydrology connected with the development of the relief in the Oulanka river valley, NE Finland[END_REF].

The planform evolution can be inferred from a cursory analysis of the valley topography, as seen in a LiDAR acquisition from National Land Survey of Finland and of satellite imagery (Figure 1A). In addition to the current river path and an abandoned meander loop, the high-resolution topographic dataset allows the identification and digitization of scroll topographies (see the red arrows 1 and 2 and 3 in Figure 1A).

Figure 1 A-LiDAR data of the Oulanka River, Finland. One abandoned meander shape can be clearly identified and used to constrain the reverse migration. Three areas with marked scroll topographies (arrows 1, 2 and 3) are clearly observable, but not used as input constraints in the ChaRMigS method (https://tiedostopalvelu.maanmittauslaitos.fi/tp/kartta?lang=en). B-C-D-Three stochastic equiprobable realizations obtained using only the ChaRMigS method in 30 time steps with a probability distribution of horizontal offset equal to a Gaussian probability distribution G (36,2.62 [m]), and a probability distribution of cut-off equal to triangular probability distribution (min = 0, mode = 1, max = 2 at a given time step). The relative age is expressed backwards, from 0 (present day; shown in darker blue) to 30 (oldest; shown in white).

Three equiprobable ChaRMigS reconstructions are presented (Figure 1B-D), each showing different possible river-migration histories. The counter-point bars present on the western part of the abandoned meander (Figure 1A, red arrow number 1) are reconstructed in part by the realizations in Figure 1C,D (green arrow 1). However, none of the reconstructed geometries fits to the scroll-bar patterns of point bars observed on the eastern side of the abandoned meander (Figure 1A, red arrow number 3). Furthermore, on the northern part of the study area, the observed scroll patterns (Figure 1A, red arrow number 2) are also not well reproduced by the model. For this area, only one of the three realizations (Figure 1C, green arrow 2) reproduces channel migration by translation. These two limitations are due to the reverse-migration vectors, which attract the main channel path upstream and towards the centreline in order to reproduce the most likely migration patterns [START_REF] Parquer | Reconstruction of channelized systems through a conditioned reverse migration method[END_REF]. In the ChaRMigS approach, the point bars observed on the eastern side (red arrow number 3) could only be reproduced by significantly shifting or enlarging the range of the migration offset probability distribution, or by generating an additional abandoned meander bordering the point bars.

In parallel to the ChaRMigS application on this dataset, and in order to enable a comparison with ChaRMigS outputs, PB-SAND was used to reproduce the accretion history of the north-western part of the channel belt, on the basis of three manually digitized channel trajectories (corresponding to the current river course and to two earlier river paths inferred from observed scroll geometries). Downstreamtranslating point bars are modelled. This point-bar type is typical of meandering rivers developed in confined valleys where bend expansion is inhibited by confinement by the valley walls [START_REF] Smith | Counter point bar deposits: lithofacies and reservoir significance in the meandering modern Peace River and ancient McMurray Formation, Alberta, Canada[END_REF][START_REF] Nicoll | Planform geometry and channel migration of confined meandering rivers on the Canadian prairies[END_REF]. The modelled scroll-bar patterns match closely those that can be identified in the LiDAR image (Figure 2A). This case study highlights how the two modelling approaches provide different insights in the channelized system in terms of both geometry and dynamics. ChaRMigS proposes multiple scenarios of how a meander is abandoned and its relation to the chronology of cut-offs during channel evolution, the style of meander-bend transformation, and the influence of these events on the planview geometries of meander-belt deposits. By contrast, PB-SAND produces a deterministic model of planform channel evolution that aligns with observations of local point-bar accretion geometries arising from specified histories of channel migration, so that it becomes possible to predict the expected lithological heterogeneity (e.g., sandier point bars vs muddier counter-point bars; cf. [START_REF] Smith | Counter point bar deposits: lithofacies and reservoir significance in the meandering modern Peace River and ancient McMurray Formation, Alberta, Canada[END_REF][START_REF] Durkin | Predicting heterogeneity in meandering fluvial and tidal-fluvial deposits: The point-bar to counter point-bar transition[END_REF] from styles of accretion (e.g., meander expansion vs downstream translation).

Pleistocene succession of the Gulf of Thailand

In a second case-study example, ChaRMigS and PB-SAND were again applied separately, in this case to model the planform evolution of a meandering channel belt seen on a time slice from a 3D reflectionseismic survey of a Pleistocene succession from the Gulf of Thailand [START_REF] Reijenstein | Seismic geomorphology and highresolution seismic stratigraphy of inner-shelf fluvial, esturarine, deltaic, marine sequences, Gulf of Thaïland[END_REF]. FAKTS was used to provide input parameters on channel-bar scaling relationships and on the proportions of lithologies, for modelling the geometry and facies distributions of point-bar and channel-fill deposits. The seismic data clearly reveals the planform bar geometry of the channel belt, within which the path of the last abandoned (avulsed) channel can be recognized, together with abandoned meanders that have undergone cut-offs. The internal accretion geometry of the point bars can be mapped to infer the paths taken through time by the river channel. Some of the inputs of ChaRMigS and PB-SAND are the same; for example, the input channel trajectory immediately prior to avulsion, digitized directly from the seismic image, as well as the channel width. randomly sampled from these distributions, as described in [START_REF] Parquer | Reconstruction of channelized systems through a conditioned reverse migration method[END_REF]. These realizations vary with regards to the direction and magnitude of channel migration (offset between trajectories) between any two time steps, which thus implies variable downstream channel migration relative to the degree of bend expansion. In some realizations, counter-point bars develop, as observed on the seismic image (see red circle 1 in B and model output in C, D and E), but also in areas where they are not visible (see red circle 2 in B and model output in C, D and E). This emerges from the stochastic nature of the ChaRMigS method, which samples the reverse migration rate randomly for each half-meander. The two abandoned meanders (Figure 3) are integrated for each realization in the same relative order but not necessarily at the same modelling time steps. The chronology according to which abandoned loops are integrated in the model is due to the integration probability inferred by the geometry and organization of the channel belt and updated at each time step [START_REF] Parquer | Integration of channel meander abandonment age uncertainty into a stochastic channelized system reconstruction method[END_REF]: on the seismic slice, meander loop A (Figure 3B) displays a channel trajectory that is more perpendicular to that of meander loop B, and is located closer to the main channel path than meander loop B (Figure 3B). This makes loop A a better candidate for integration in the reverse migration process than meander loop B, in the first stages of reconstruction. This is because the probability of integration is directly related to the angular separation between the main channel path and the meander loop, as well as the distance between the two [START_REF] Parquer | Integration of channel meander abandonment age uncertainty into a stochastic channelized system reconstruction method[END_REF]. Therefore, consistency in the relative chronology of cut-off occurrence reflects limited uncertainty in the actual relative timing. The absolute timing of cut-off integration is determined stochastically in the algorithm used to simulate the reverse migration of the main channel path for each time step. The probability of integration of an abandoned meander is recomputed at each time step and controlled by the random sampling migration offset. Therefore, the relative location of the reverse migrated channel path will be different for each realization due to difference in the integration probability. The PB-SAND model has also been employed separately to reconstruct the same channelized system by using the point-bar accretion geometries evident in the seismic imagery and which reflect the evolutionary history of the formative river channel, together with the trajectory of the last (pre-avulsion) channel path (Figure 4A). An overall fining-upward facies arrangement -as is common for the deposits of these types of meandering fluvial systems [START_REF] Nanson | Point bar and floodplain formation of the meandering Beatton River, northeastern British Columbia, Canada[END_REF] -is used to depict the internal facies distribution within the modelled lateral-accretion packages, as shown in cross sections (Figure 4C). Contrasting facies proportions have been selected for deposits of point-bars (gravel: 10%, coarse sand: 20%, medium sand: 40%, fine sand: 10%, silt: 10%, clay-rich mud: 10%) and counter-point-bars (gravel: 10%, coarse sand: 10%, medium sand: 10%, fine sand: 20%, silt: 40%, clay-rich mud: 10%). These values of facies proportions are taken from the FAKTS database of geological analogues and are representative of meander belts of this scale [START_REF] Colombera | A quantitative approach to fluvial facies models: methods and example results[END_REF]. The proportion for each facies varies systematically between -5 % to +5%, to mimic inherent natural variability [START_REF] Yan | A 3D forward stratigraphic model of fluvial meander-bend evolution for prediction of point-bar lithofacies architecture[END_REF]. The modelled facies architecture can be rendered in cross sections (Figure 4C) that permit visualisation of the sedimentary heterogeneity created by the different migration patterns. For example, the finer counter-point-bar deposits in the left-hand part of cross section A-A', and the transition from sandier point-bar deposits to finer counter-point-bar deposits in cross section F-F' are reconstructed in a similar way to that observed on the seismic image and reproduced by ChaRMigS in Figures 3D andE. The 3D geometries of lateral accretion packages and the facies distribution of point-bar deposits is similar to those of comparable outcrop successions (e.g., [START_REF] Ghinassi | Plan-form evolution of ancient meandering rivers reconstructed from longitudinal outcrop sections[END_REF][START_REF] Ghinassi | Downstream-migrating fluvial point bars in the rock record[END_REF][START_REF] Durkin | Predicting heterogeneity in meandering fluvial and tidal-fluvial deposits: The point-bar to counter point-bar transition[END_REF].

Considerations on the modelling approaches

Based on the separate application of the two modelling techniques to the two datasets presented above, the weaknesses and strengths of the two approaches can be discussed. The case-study example presented above highlights how PB-SAND is able to model meander-belt sedimentary architectures in 3D using input channel trajectories that can be inferred from clearly identifiable point-bar accretion geometries, in situations where these geometries can be seen in 3D seismic surveys. However, integration of abandoned meanders cannot be achieved directly in PB-SAND, and the trajectories of these meander loops cannot be used to condition the model outputs. This can, however, be achieved using ChaRMigS, which produces a range of possible planform geometries associated with equiprobable scenarios of channel evolution while being conditioned to abandoned meander loops. As such, the two modelling methods can be used in a complementary manner as demonstrated in the combined workflow that is presented next.

4. Application of the combined workflow: a case study (Pleistocene, Gulf of Thailand)

Workflow and dataset

The ChaRMigS and PB-SAND modelling approaches each have advantages and limitations. This naturally leads us to propose a combined use of both methods following the workflow illustrated in Figure 5. The main idea of this workflow is to generate scenarios of palaeo-channel evolution from incomplete observations on a seismic slice using ChaRMigS, then to use the channel trajectories in PB-SAND to generate facies models, using parameters gathered from the FAKTS database. We employed this approach to generate 3D models of the internal facies architecture of meander-belt deposits imaged in seismic data from the Pleistocene of the Gulf of Thailand. In this work, sets of channel trajectories representing stochastic outputs of ChaRMigS have been smoothed, trimmed and loaded as input into PB-SAND, together with analogue data relating to the geometry of element and facies proportions from the FAKTS database (see below). The proposed workflow is applied to show how the ChaRMigS and PB-SAND modelling approaches, augmented and supported by information from the FAKTS database, can be used to reveal complementary insights into the evolution of fluvial systems and their preserved meander-belt architectures. This case study was chosen to demonstrate how the two modelling approaches can be used in combination to provide a novel workflow for subsurface modelling (Figure 5). Such an approach is especially useful for situations in which channel belts of meandering fluvial systems (including associated channel fills) are recognizable in 3D seismic surveys, but where internal point-bar accretion geometries cannot be discerned. In this case, ChaRMigS is applied to generate multiple stochastic realizations of the internal accretion history of a meander belt in a way that is compatible with the observable geometries of abandoned channel fills seen on a seismic time or stratal slice. The ChaRMigS model outputs channel trajectories of different stochastic realizations, which can then be transferred as input to PB-SAND for simulation of the facies heterogeneities. Such a situation is exemplified in a 3D reflection seismic time slice from a meander belt sitting inside a valley fill in the Pleistocene succession of the Gulf of Thailand [START_REF] Reijenstein | Seismic geomorphology and highresolution seismic stratigraphy of inner-shelf fluvial, esturarine, deltaic, marine sequences, Gulf of Thaïland[END_REF]. In the considered time slice (Figure 6A), the internal anatomy of individual pointbar elements cannot be resolved from the seismic data alone. Therefore, channel trajectories with which to attempt the direct application of PB-SAND cannot be inferred with confidence. In this case, five abandoned meanders that have undergone cut-off prior to channel-belt abandonment (Figure 6A) have been used to condition ChaRMigS to generate possible planform evolution histories. This particular time slice has been selected for illustrating this combined workflow for two specific reasons: (i) it reveals substantial architectural complexity of the channelized system (i.e., high number of abandoned meander loops with highly sinuous channel fills); and (ii) it exemplifies a situation in which this approach would be most useful (i.e., a case where internal point-bar accretion geometries cannot be seen directly in the seismic data). ). For both presented realizations, all five observed abandoned meanders have been integrated. However, the chronology of integration is markedly different. This is due to the difference in the migration pattern of the bar deposits in each reconstruction: the migration pattern is mainly expansional in Figure 6B and mainly translational in Figure 6C. This highlights the ability of ChaRMigS to develop models of channel-belt architecture with markedly contrasting accretion geometries for a given set of starting conditions.

ChaRMigS

PB-SAND

Figure 6 A-Seismic image of a channelized system from the Pleistocene of the Gulf of Thailand [START_REF] Reijenstein | Seismic geomorphology and highresolution seismic stratigraphy of inner-shelf fluvial, esturarine, deltaic, marine sequences, Gulf of Thaïland[END_REF]. B-C-Two stochastic and equiprobable realizations obtained with ChaRMigS within 150 time steps. In parts B and C, the numbers to which integration chronologies are associated (1 to 5)

refer to meander loops labelled in part

A. The white boxes denote the area that has been modelled in 3D using PB-SAND, as shown in figures 7 and 8.

Both the relative chronology of integration of the abandoned meanders and the channel migration patterns vary significantly across the two realizations (see the arrow that indicates the temporal order of integration on the top of Figure 6B,C). The abandoned meander loop 2 (Figure 6A) is the first one to be integrated in the realization in Figure 6C (i.e., it corresponds to the most recent cut off in this scenario), but the last one to be integrated in realization Figure 6B (i.e., it corresponds to the oldest cut off in this scenario).

The integration of ChaRMigS and PB-SAND for the creation of static reservoir models effectively requires the use of the stochastic inverse modelling outputs from ChaRMigS, (i.e., channel trajectories) as inputs to the forward modelling runs of PB-SAND, which is used to generate realizations of the expected internal architecture of selected parts of the modelled channel belt.

For the channel belt depicted within the white boxes in Figure 6, PB-SAND has been applied to model the southern sector (Figure 7 and Figure 8). The deterministic planform realizations created by PB-SAND were conditioned using the ChaRMigS output channel trajectories resulting from both equiprobable realizations. In each realization, a reduced number of ChaRMigS output channel trajectories have been manually chosen, that are sufficient to control the transformation styles of point bars at their different evolutional stages (Figure 6): 15 and 12 bends, respectively. PB-SAND runs for multiple meander bends and associated point-bar and channel-fill elements have been incorporated in two resulting planform architectures, one for each ChaRMigS realization (Figures 6B,C), as presented in Figure 7E and Figure 8E. These outputs from PB-SAND embody smoother accretion geometries compared to the ChaRMigS realizations (Figure 6). Indeed, input trajectories from ChaRMigS were smoothed to meet the requirements of PB-SAND, which requires that each input trajectory defining a point bar has the same number of points (definition nodes). The accretion histories and styles predicted using ChaRMigS are maintained overall, however. This is recognizable when comparing, for example, the reconstruction of the lateral migration of the abandoned meander 2 in Figure 7E and the downstream migration of abandoned meander 1 in Figure 8E. The two planform realizations produced by PB-SAND (Figures 7 and8) have been used as frameworks for modelling cross-sectional geometries and facies distributions. The internal facies architecture of the modelled point-bar elements has been simulated using data on the proportions of facies in a range of successions that act as geological analogues contained in the FAKTS database (Colombera et al., 2012b[START_REF] Colombera | A quantitative approach to fluvial facies models: methods and example results[END_REF], as presented in Figure 9. The proportions of each facies (fines, sand, and gravel) are randomly sampled between the first and third quartile in such way that they add up to 100 %, using the approach adopted in [START_REF] Yan | Three-dimensional modelling of fluvial point-bar architecture and facies heterogeneity using analogue data and associated analysis of intra-bar static connectivity: application to humid coastal-plain and dryland fluvial systems[END_REF]. This enables the sampled facies proportions to be representative whilst ensuring reasonable variability between samples. As the modelled point-bar elements are relatively sand-rich overall [START_REF] Feng | An investigation of fluvial geomorphology in the Quaternary of the Gulf of Thailand, with implications for river classification[END_REF][START_REF] Reijenstein | Seismic geomorphology and highresolution seismic stratigraphy of inner-shelf fluvial, esturarine, deltaic, marine sequences, Gulf of Thaïland[END_REF], differences between samples are limited, and therefore only one facies association is used in the model presented here.

Input facies proportions are as follows: 3% gravel, 81% sand, and 16% fines. To represent grain-size trends more clearly, sands are differentiated into three classes (coarse, medium and fine) of equal proportion (each 27% of the total), whereas fines are divided into two classes (silt and clay-prone) of equal proportion (each 8% of the total). To simulate and represent lithological contrasts in planview, the proportion of fines in counter-point-bar deposits is increased at the expense of sand classes (clay: 24 %; silt: 24 %; fine sand: 16 %; medium sand: 16 %; coarse sand: 17 %; gravel: 3 %). Abandoned-channel fills are modelled as homogeneous mud plugs. The thickness of the channel fills and of the associated point-bar deposits reflect scaling relationships between channel-fill width -which is observed in the seismic data -and thickness, as based on geological analogues from the FAKTS database. For the sake of simplicity, temporal variations in river depth that match observed changes in channel-fill width have not been modelled in this case, although such functionality is incorporated in PB-SAND.

The facies architecture of the modelled meander belt is presented on cross-sections, which have been traced through both realizations in Figure 7E and Figure 8E at corresponding positions (Figure 10 and 11). Green arrows point to contacts between accretion packages in point-bar units with contrasting accretion style and lithological make-up.

Figure 11 Cross-sections made in the PB-SAND realization Figure 8E. The differences with the realization in Figure 7E are highlighted by these cross-sections (compare the blue dashed boxes here with those shown in Figure 10). Green arrows point to contacts between accretion packages in point-bar units with contrasting accretion style and lithological make-up.

Discussion

Reservoir-modelling applications

Seismic data are routinely used for constraining subsurface reservoir models. For example, seismic attributes are employed as soft-data (e.g., facies probability) for conditioning geostatistical models (e.g., [START_REF] Strebelle | Conditional simulation of complex geological structures using multiple-point statistics[END_REF][START_REF] Yao | Integrating seismic data for lithofacies modeling: A comparison of sequential indicator simulation algorithms[END_REF]Rongier et al., 2017a), or can be used to determine outputs of process-based forward models that best honour subsurface observations [START_REF] Koneshloo | A workflow for static reservoir modelling guided by seismic data in a fluvial system[END_REF]. However, geological insight that can be drawn from geometries of sedimentary units seen in 3D seismic datasets is typically underused in reservoir-modelling practice. The presented workflow leverages on planview observations of meander-belt architectural elements, i.e. point bars and channel fills, to build stochastic models of facies heterogeneity that honour the distribution of these elements, through the sequential application of two algorithms: ChaRMigS and PB-SAND.

ChaRMigS permits the creation of different equiprobable 3D geometries resulting for meandering-channel evolution, through an inverse modelling approach, based on observations of the youngest channel trajectory and of abandoned meander loops. PB-SAND allows the creation of stochastic facies distributions that consider the type of planform evolution of a meander belt, using a forward modelling approach. PB-SAND complements ChaRMigS in that it can be used to model lithologies, whereas ChaRMigS complements PB-SAND in its ability to generate stochastic planforms that incorporate abandoned meanders. These two complementary methods can therefore be integrated effectively to build 3D geocellular models of subsurface meander-belt successions that are seen in 3D reflection seismic surveys. This approach is particularly applicable to cases in which the internal accretion geometry of pointbar elements cannot be resolved from the available data, but the trajectory of channel fills can be identified in horizontal slices, for example in pre-drill situations. The integration of these two methods, which is exemplified in the applications described herein and summarized in Figure 5, represents a novel approach to reservoir modelling. The application to the subsurface of the Gulf of Thailand (above) shows how this integrated workflow can provide additional insight to subsurface characterization in data-limited settings, and how the combined modelling approach significantly enhances our ability to predict lithological heterogeneity and facies distributions in a way that is consistent with possible depositional history. The proposed workflow can be readily applied in pre-drill situations where well conditioning is not necessary. It can also be applied in conjunction with pixel-based geostatistical modelling algorithms to reservoir units that have been drilled for which well conditioning is required [START_REF] Colombera | Seismic-driven geocellular modeling of fluvial meander-belt reservoirs using a rule-based method[END_REF].

Future research

The effective combination of the two modelling methods can be further improved and refined by (i) increasing the internal consistency of channel trajectories, facies distributions and the level of integration of both methods, and (ii) improving the automation in matching all available observations. Both aspects would lead to a reduction of uncertainties and would further improve the predictive ability of the produced reservoir models.

In ChaRMigS, improvements can be made in the morphodynamic realism of trajectories that are produced in some circumstances. For example, better interpolation between channel centrelines is needed to fill gaps determined in certain cases by the integration of abandoned meanders, particularly in association with closely spaced meander bends. Also, the way in which ChaRMigS incorporates abandoned meander loops in the model outputs may, at times, result in river planforms that are not realistic, because they would represent situations whereby the inception of a channel path following cut-off results in a very highly sinuous trajectory (see, for example, the evolution of meander 2 of Figure 6C as modelled in Figure 8E). The same abandoned meander highlights limitations of PB-SAND, with regards to its ability to handle sinuous input trajectories and to its requirement for centreline smoothing. Such smoothing may have an impact on the reproduction of architectures that are delivered by ChaRMigS, such as abandoned loops integrated into the main channel path. In terms of method integration, modelling the architecture of bars and channel fills associated with meander belts containing several abandoned meanders is currently a cumbersome process. Indeed, PB-SAND models short stretches of a channel reach, with up to three meander bends, and therefore requires that the results of multiple modelling steps be merged when a complex amalgamated meander belt is modelled. The application of the combined workflow to the Pleistocene Gulf of Thailand clearly demonstrates the value of implementing a tighter integration between ChaRMigS and PB-SAND.

In terms of facies prediction, PB-SAND, uses geological rules that link facies to accretion styles (e.g., downstream fining to counter-point-bar regions, Figure 7E and Figure 8E; fining upward, Figure 10 and Figure 11). ChaRMigS fulfils the need to generate planform geometries of meandering reaches, which are required to model these 3D facies architectures by PB-SAND. Future research goals to complete the 3D consideration of reservoirs include the introduction of additional components of lithological heterogeneity within the facies using geostatistics, incorporation of the effects of vertical sand-body stacking, plus incorporation of additional architectural elements, such as levees and crevasse splays. Such model features could be used to account for external forcing parameters (see, for example, [START_REF] Yan | Three-dimensional modelling of fluvial point-bar architecture and facies heterogeneity using analogue data and associated analysis of intra-bar static connectivity: application to humid coastal-plain and dryland fluvial systems[END_REF].

A second general avenue for further research is to better and more effectively honour features seen on seismic sections and stratal or time slices or other types of subsurface data. One way of achieving this more effectively is by finding ways to infer model parameters to generate more models that automatically honour all observations. In the present state of the methods, some model results need to be eliminated by experts. For example, when applying ChaRMigS on the Pleistocene Gulf of Thailand data set, we observed that several realizations could not incorporate all abandoned loops. This may occur in particular when the channel path crosses an abandoned meander, as it may then need to migrate for a large number of time steps before the integration can be considered. In a sense, this represents a "missed opportunity" for meander loop integration, which can hardly be recovered as the interpreted older loop should then no longer be visible. Therefore, when this happens, the considered realization should be rejected as incompatible with observations. In this case, a possible strategy to explore the uncertainties and obtain a wide range of realizations requires changing the migration rate distributions. More generally, parameter inference for the proposed method is clearly an aspect which calls for further research.

A second way towards honouring all observations is to augment the input of the methods to automatically allow for new types of observations. For example, some of the stochastic channel trajectories simulated by ChaRMigS sit outside the confines of the channel belt imaged in the seismic data (i.e., they are placed on the valley interfluves where in reality the formative river had never wandered) (Figures 6B-C). To prevent problems of this type, it would be ideal to be able to define spatial domains within which the model would operate. Along the same lines, more types of hard and soft data should be incorporated for model conditioning (e.g., channel-belt limits, boreholes, vertical facies proportions). The difficulty in honouring heterogeneous data sets is a limitation that affects this class of modelling methods more widely [START_REF] Pyrcz | Stratigraphic rule-based reservoir modeling[END_REF]. The evaluation of model outputs is currently based on visual screening: a quantitative and automatic approach for screening of model results would be preferable.

Conclusions

This work has presented applications of inverse (ChaRMigS) and forward (PB-SAND) models for the architecture of fluvial meander belts; these models have been applied to modern and ancient fluvial systems for scopes of comparison and to demonstrate how they can be linked sequentially in a new workflow for geocellular reservoir modelling.

Based on input trajectories of the planform of channel fills, including those of oxbow-shaped abandoned meander loops, ChaRMigS generates several equiprobable migration patterns that can be used to assess uncertainty about internal reservoir architectures. PB-SAND instead uses input channel trajectories, such as those produced by ChaRMigS, to generate a geometric framework of meander-belt accretion and to populate this framework by facies of different types, stochastically, based on suitable analogue data, for example as provided by FAKTS, a database of real-world case-study examples.

When working with a 3D seismic dataset in which abandoned meander scars and avulsed channel fills are observed, but in which internal point-bar growth geometries are below seismic resolution, the proposed workflow, which combines ChaRMigS with PB-SAND, makes it possible to explore a range of meander-belt facies architectures that are realistic with respect to expected styles of accretion and associated lithological heterogeneity. This novel approach is valuable for quantification of subsurface uncertainty.

members listed on http://ring.georessources.univ-lorraine.fr/consortium/members) for funding that has enabled the development of ChaRMigS, and Emerson for donating SKUA-GOCAD License. We thank an anonymous reviewer who provided constructive comments that improved the paper.

Figure 2

 2 Figure 2 Application of the PB-SAND method to the LiDAR data of the Oulanka River on the basis of three observed channel stages detailed in A. Realization detailed in B for point bars stage 1 (t1 -t2) and C for point bars stage 2 (t1-t2-t3).

Figure 3 A

 3 Figure3A-Seismic image of a palaeo-channelized system from the Gulf of Thailand[START_REF] Reijenstein | Seismic geomorphology and highresolution seismic stratigraphy of inner-shelf fluvial, esturarine, deltaic, marine sequences, Gulf of Thaïland[END_REF]. B-Input trajectories of the last avulsed channel path and of two abandoned meander loops. Two counter point-bar domains are clearly observable in the seismic image (labelled 1 and 2) and are used to discuss realization quality (see text). C-, D-, E-Three stochastic equiprobable realizations of reconstruction of the palaeo-channelized system made using only the ChaRMigS inverse model. Results from ChaRMigS are shown in Figure 3. The three reconstructions of the channelized system (Figure 3C-E) represent different stochastic realizations generated with ChaRMigS, each in 30 steps with a Gaussian probability distribution of offset of horizontal migration (mean = 333 m, standard deviation = 166 m) and a triangular probability distribution of channel cut-off (minimum = 0, mode = 1, maximum = 2 [number of cut-offs per time step], at a given time step). For each realization, the actual values are

Figure 4

 4 Figure4PB-SAND modelling results of the seismic image from the Gulf of Thailand[START_REF] Reijenstein | Seismic geomorphology and highresolution seismic stratigraphy of inner-shelf fluvial, esturarine, deltaic, marine sequences, Gulf of Thaïland[END_REF] presented in Figure3A. A-Input lines for PB-SAND reconstruction. B-Modelled accretion surfaces between two consecutive input trajectories are denoted from blue to red in a temporal sequence. C-Cross-sections of the reconstruction shown in B. The reconstruction here is based on four channel-path centrelines. Then, a simulation of the internal facies architecture of the point-bar deposits has been undertaken to reveal the expected facies distribution within the reconstructed system. The cross sections depict the expected spatial distribution of six different facies within the point-bar elements. Facies distributions are determined from distributions from known natural point-bar examples stored in the FAKTS database.

Figure 5

 5 Figure5Workflow combining ChaRMigS and PB-SAND, together with inputs from 3D seismic data and FAKTS database analogues, to create stochastic models of meander-belt architectures.

  Figures6B-Cshow two examples of proposed stochastic ChaRMigS reconstructions of the planform geometry of channelized systems observed on the seismic image. These results have been obtained with a Gaussian probability distribution of horizontal migration offset (mean = 60, standard deviation = 15 [m]) and a triangular probability distribution of cut-off (minimum = 0, mode = 1, maximum = 2 [number of cutoffs per time step]). For both presented realizations, all five observed abandoned meanders have been integrated. However, the chronology of integration is markedly different. This is due to the difference in the migration pattern of the bar deposits in each reconstruction: the migration pattern is mainly expansional in Figure6Band mainly translational in Figure6C. This highlights the ability of ChaRMigS to develop models of channel-belt architecture with markedly contrasting accretion geometries for a given set of starting conditions.

Figure 7 .

 7 Figure 7. Parts A-D show input channel paths for PB-SAND realization 1, corresponding to fifteen selected time steps of the ChaRMigS realization of Figure 6B (inside the white box). Part E shows PB-SAND results for the ChaRMigS simulation presented in Figure 6B. The yellow, dark-grey, and light-grey areas denote relatively sandy point-bar deposits, finer counter-point-bar deposits, and a transitional zone between point bar and counter-point bars where downstream fining takes place, respectively. Grey and blue channel fills have been abandoned by cut off and avulsion, respectively. Thin lines of different colours denote accretion surfaces at different time steps. Pink lines indicate cross-sections shown in Figure 10.

Figure 8

 8 Figure 8 Parts A-D show input channel paths for PB-SAND realization 2, corresponding to twelve selected time steps of the ChaRMigS realization of Figure 6C (inside the white box). Part E shows PB-SAND results for the ChaRMigS simulation presented in Figure 6C. The yellow, dark-grey, and light-grey areas denote relatively sandy point-bar deposits, finer counter-point-bar deposits, and a transitional zone between point bar and counter-point bars where downstream fining takes place, respectively. Grey and blue channel fills have been abandoned by cut off and avulsion, respectively. Thin lines of different colours denote accretion surfaces at different time steps. Pink lines indicate cross-sections shown in Figure 11.

Figure 9

 9 Figure 9 Facies proportions extracted from the FAKTS database for point-bar elements in meandering fluvial systems. The part of each distribution between the first and second quartiles is denoted in red; the part between the second and third quartiles is denoted in blue; the white bar inside each box denotes the median; the whiskers on both sides denote the minimum and maximum values excluding outliers shown as black crosses (outside of 1.5 times the interquartile range). Data relate to 99 point bars containing 593 facies units. "Fine" refers to clay and silt sediment fractions.

Figure 11 ,

 11 Figure11, respectively). Together with the planview model outputs (Figure7Eand Figure8E) these crosssections convey the variability in facies architecture that results from the contrasting accretion styles associated with the stochastic planform evolutions provided by ChaRMigS. It is evident that local abrupt changes of migration direction induce point-bar reactivation and juxtaposition of sandier and muddier sediment volumes (e.g., cross-section B-B' in Figures7E and 10). It is also evident, for example, how facies distributions vary across realizations in relation to downstream-fining trends observed within point-bar deposits and through their transition to counter-point-bar deposits (e.g., cross section C-C' in Figure8Eand Figure11). These architectural styles can cause significant petrophysical heterogeneity and can locally compartmentalize part of the meander-belt reservoir (see the green arrows in Figures10 and 11).

Figure 10

 10 Figure 10 Cross-sections made in the PB-SAND realization of Figure 7E. The differences with the realization in Figure 8E are highlighted by these cross-sections (compare the blue dashed boxes here with those shown in Figure11). Green arrows point to contacts between accretion packages in point-bar units with contrasting accretion style and lithological make-up.

Table 1 :

 1 Summary of main characteristics and differences between the two modelling approaches used in this work.

ChaRMigS (Parquer et al. 2017) PB-SAND (Yan et al. 2017) Approach

		Inverse temporal model, geometric (not	Forward temporal model, geometric (not
		process-oriented)	process-oriented)
	Framework	Vector-based, grid-free	Vector-based, grid-free
	Coding Language	C++		Matlab
	Channel evolution and accretion	•	Stochastic output	•	Deterministic output
	geometries	•	Channel trajectories inferred from	•	Channel trajectories inferred from
			channel fills		observed point-bar accretion
		•	Ability to model outcomes of bend		geometries
			expansion, translation, and rotation	•	Ability to model outcomes of bend
		•	Magnitude of inner-bank accretion		expansion, translation, and rotation
			may vary irregularly along stream	•	Magnitude of inner-bank accretion
			and between time steps		varies linearly along stream and is
		•	Loci of intra-channel belt erosion		constant between time steps of two
			may vary along stream between		input trajectories
			pairs of time steps	•	Intra-channel belt erosion limited by
					constraints imposed by centreline
					intersection

Number of meander bends modelled in a single run Unlimited 3 (but able to link model outputs together to model an unlimited number) Incorporation of oxbow-shaped abandonments

		•	Cut-off rate probability distribution		trajectories of channel-paths at key
		•	Abandoned meander loop paths		times
				•	Scroll-bar number and migration rate
					between two consecutive
					trajectories
				•	Facies associations
				•	Facies trends
				See Fig. A.1 in Yan et al. (2017) for a fuller
				account of input parameters.
	Outputs	Multiple scenarios of reconstructed 3D	3D modelled geometry of fluvial meander-
		channelized system geometry	belt and internal facies composition of
				point-bar and related architectural
				elements in 3D
		Directly by the model; trajectories used	Achieved by overlaying results of separate
		for conditioning model	modelling runs
	Variations in deposit thickness	Modelled by a local vertical migration rate	Resulting from streamwise bathymetric
		parameter, but assumed constant in	variations and temporal changes
		applications to-date.
	Channel-belt aggradation	As vertical offset between accretion steps	As vertical offset between accretion steps
	Lithologies	•	Not currently modelled	•	Stochastic output
			automatically	•	Modelled using a rule-based,
		•	Can be modelled manually by		stochastic approach
			creating a grid in each channel	•	Not conditioned on well data
	Ideal modelling applications	Sector models of meander-belt reservoirs	Sector models of meander-belt reservoirs
		of any scale, with channel fills and oxbow-	at any scale, but especially suitable for
		shaped abandonments imaged in 3D	cases where internal accretion geometries
		seismic datasets	are clearly imaged in 3D seismic datasets
	Main Inputs	•	Youngest preserved channel path	•	Youngest preserved channel path
		•	Migration rate probability	•	Numerous remnant point-bar
			distributions		deposits used to determine

Acknowledgements

The authors thank Hernan Reijenstein for authorizing the use of seismic data from Reijenstein et al. (2011). MP was supported by an award from the International Association of Sedimentologists (IAS). NY, LC and NPM acknowledge sponsors of the Fluvial & Eolian Research Group (Aker BP, Anadarko, Areva [now Orano], BHP Billiton, Cairn India [Vedanta], Chevron, ConocoPhillips, Equinor, Murphy Oil, Nexen-CNOOC, Occidental, Saudi Aramco, Shell, Tullow Oil, Woodside, and YPF) and our partner Petrotechnical Data Systems (PDS) for financial support. LC has been supported by the Natural Environment Research Council (Catalyst Fund award NE/M007324/1, Follow-on Fund NE/N017218/1). PC, GC and MP acknowledge the support of the RING Consortium managed by ASGA (ANDRA, CCR, Chevron, Equinor, ExxonMobil, IFPEN, Emerson, Petrobras, PDS, Repsol, Saudi Aramco, Schlumberger, Total, WintershallDEA and the associate


Sylvain Tourancheau

Patrick Le Callet
email: patrick.lecallet@univ-nantes.fr

Kjell Brunnström
email: kjell.brunnstrom@acreo.se

Dominique Barba

DISPLAY AWARENESS IN SUBJECTIVE AND OBJECTIVE VIDEO QUALITY EVALUATION

Most of studies on video quality assessment are focused on impact of coding distortion or transmission error. In this paper, display is considered. Regarding technology, subjective experiments suggest that there are differences in term of quality between LCD and CRT. CRT provides better quality when viewing HD video content while it is LCD when viewing still colour images. One explanation of this behaviour is explained through motion blur. From a motion blur perception model, an efficient metric of the quality loss due to this effect on LCD is proposed. Finally, subjective experiments results using SD video content are consistent with a motion blur based model and point out that the effect of display technology is linked with video resolution.

INTRODUCTION

The incoming of the high-definition new visual experience at home have boosted new display technologies, since they enable the increase of screen size necessary to sense immersion, impact and immediacy as in a movie theatre [START_REF] Fujio | Future broadcasting and high definition television[END_REF]. For these reasons, these new displays, and particularly liquid crystal displays (LCD), will soon replace the old mature CRT technology.

Liquid crystal displays have many differences with CRT displays. Subjective preference tests between this two types of displays have highlighted a high preference for CRT displays concerning moving pictures [START_REF] Itu | Report on results of comparative subjective picture quality assessment test between CRT and LCD[END_REF]. Many defects have been counted by viewers, such as colour differences, degradations in dark areas and de-interlacing artifacts for interlaced sequences. But among all these defects, motion blur seems to be the most annoying one, particularly in sequences with significant movements. On the other hand, CRT displays suffer from several shortcomings too. Flickering can be annoying in certain conditions, and the small luminance range can led to flat pictures with dirty colours.

In the recent years, subjective and objective quality assessment becomes a research topic of interest. Activities of the Video Quality Experts Group (VQEG) are a good example of this interest. Previous works [START_REF] Vqeg | Final report from the video quality experts group on the validation of objective models of video quality assessment[END_REF] and work in progress [START_REF] Vqeg | Multimedia group test plan[END_REF] are mainly related to coding or transmitting purpose at a given resolution, e.g. coding artifacts and transmission errors. Considering the whole chain, quality assessment should be able to manage dependency to other technology issues.

In this paper, display is studied. As a consequence, only the high part of the quality range is considered, using video and pictures with no (uncompressed) or very few coding distortions. In order to know the impact of distortions due to flat panel displays technology with respect to the old mature CRT technology, subjective quality assessment of HDTV sequences is performed both on LCD and CRT displays.

Assuming that LCD motion blur is the most annoying artifacts when displaying moving pictures on LCD, its perception is described and its magnitude measured. This leads to the design of an objective metric which enables the prediction of the perceived quality on LCD with respect to the CRT one. Such a metric could be used in order to evaluate LCD improvements introduced by manufacturers to reduce technology artifacts.

In the last section, the impact of LCD technology on perceived quality regarding the video resolution is explored with subjective experiments on SD sequences.

SUBJECTIVE QUALITY ASSESSMENT TESTS

Tests conditions and equipment

Subjective quality assessment tests have been performed in a specific showroom, with lighting conditions and display parameters precisely measured and adjusted according to BT.500-11 and BT.710-4 ITU recommendations. Tests have been conducted both on a HDTV CRT display JVC DT-V 1910CG and then on a HDTV LCD Philips T370 HW01. Surrounding conditions and display parameters for each screen are presented in Tables 1 and2 quality of processed sequences (resp. pictures) both between them and with the quality of the explicit reference sequence (resp. picture). This leads to precise and reliable measures of the quality [START_REF] Blin | New quality evaluation method suited to multimedia context: Samviq[END_REF]. Notation scale is continuous, each score can take a value between 0 and 100.

Observers

Observers were mostly male students between 20 and 25. All were familiar with television and cinema but not with HDTV. Acuity and colour blindness have been checked for each observer, respectively with Monoyer's plates and Ishihara colour vision test. Observers with at least on error in Ishihara test or with an acuity less than 9/10 are rejected. After the tests have been completed by all the observers, a rejection technique from the EBU [START_REF] Ebu | SAMVIQ -Subjective assessment methodology for video quality[END_REF] is applied. This process verifies the level of consistency of the scores of one observer according to mean score of all observers. Following the application of this rejection process, 15 valid subjects should be retained at minimum.

VIDEO QUALITY ASSESSMENT

Material

In order to measure the difference of quality between the two types of displays for moving pictures, nine 1080i sequences with significant movements have been chosen. These videos have been supplied by European broadcasters SVT and Euro1080.

Each of them contains 250 frames which corresponds to a 10-second duration. Each reference (uncompressed) has been distorted with the H.264 reference coder at seven different bit-rates in a range of quality from . Sequences are received in 1080i format by the two displays. They're displayed in interlaced format on the CRT but not on the LCD which de-interlaces them since the flat panel matrix are only able to display progressive format.

Results

Mean opinion scores (MOS) of the observers for the ten reference sequences and for the two types of displays are presented in Table 3. ∆MOS is the difference of MOS from CRT and LCD:

∆MOS = MOS CRT -MOS LCD (1)
The perceived quality of moving pictures displayed on LCD is globally lower than the perceived quality of moving pictures displayed on CRT display. It's interesting to notice that this loss of quality is quite important for sequences with quick movements such as Concert, Parkrun, Foot and Voile. This loss of quality on LCD seems to be related to the quantity and/or the fastness of the movements in the sequence. To validate this hypothesis, it has been decided to conduct the same experiment with still pictures.

Séquence

STILL PICTURES QUALITY ASSESSMENT

Material

Five images have been chosen in order to measure the difference of quality between the two types of screens for still images. They contain specific contents such as natural textures, flesh colours, oriented contours, water reflection, characters, etc. Each of them has been distorted with two types of process: jpeg compression and down-scaling/up-scaling filtering. Pictures are displayed in interlaced format on the CRT and in progressive format on the LCD in order to repeat exactly the same conditions as those of the video quality assessment.

For this test, the same group of observers has been used for the two displays. The group has been split in two parts: the observers of the first part have passed the test on CRT first, the observers of the second part have passed the test on LCD first.

Results

MOS of the observers for the five reference pictures and for the two displays are presented in Table 4 It can be observed that for still pictures the quality on LCD is globally preferred. For the pictures House and Map the difference between the two types of displays is largely in favour of LCD (with a ∆MOS of about a third of the quality scale). This can be explained by the presence of horizontally fine oriented contours which make the flickering of the CRT display more noticeable.

Overall, shortcomings of CRT displays such as flickering and limited range of luminance seems to lead to a lower feeling of natural and sense of immersion. LCD is brighter, vivid and colourful and the perceived quality of still pictures is clearly higher on it (∆MOS mean = -17.6). However, with the exact same viewing conditions and displays parameters, the perceived quality of moving pictures is higher on CRT (∆MOS mean = 6.4). It's assumed that this difference must be due to moving artifacts such as LCD motion blur which are not present with still pictures.

LCD MOTION BLUR

Results described in the previous part lead to the statement that the excellent perceived quality on LCD with still pictures is strongly reduced with moving pictures. Moving artifacts due to LCD technology, and particularly LCD motion blur, seem to be responsible for this loss of quality in video.

In this part, LCD motion blur is described. The perception of this motion blur is then measured and a model of perception is proposed. Finally, this perception model is used to design an objective metric which enables the prediction of the loss of quality on LCD with respect the perceived quality on CRT.

Description

LCD motion blur has been widely studied in recent works [START_REF] Michiel | The temporal MTF of displays and related video signal processing[END_REF][START_REF] Hao Pan | LCD motion blur modeling and analysis[END_REF][START_REF] Feng | LCD motion-blur analysis, perception, and reduction using synchronized backlight flashing[END_REF]]. It's mainly caused by the hold-type LCD's displaying method: the light intensity is maintained on the screen for the duration of the frame, whereas on CRT light intensity is a pulse which fades over the frame duration.

The main difference happens when the eyes of the observer are tracking a moving object on the LCD screen: for a given frame, the picture is sustained on the screen while the eyes are still moving slightly anticipating the movement of the object. Edges of this object are displaced on the retina resulting in a blur [START_REF] Kurita | Moving picture quality improvement for hold-type AM-LCDs[END_REF].

Model of motion blur perception

In order to measure the relation between the motion velocity and the magnitude of perceived blur in the simple case of scrolling bars, psycho-physics measurements have been designed [START_REF] Tourancheau | Prediction of perceived quality differences between CRT and LCD displays based on motion blur[END_REF]. Results of these experiments are presented in Figure 1, they lead to the following linear relation :

W = aV ; (2)
The width W (in pixels) of motion blur that appears on the edges of a moving object is proportional to its velocity V (in pixels per frame). Pan et al. have developed a theoretical model of LCD motion blur perception [START_REF] Hao Pan | LCD motion blur modeling and analysis[END_REF] and they obtain the same relation. Their model permits to identify the parameter a, which depends on the temporal function of the display.

Prediction of ∆MOS based on motion blur

Since subjective assessment on HD video content provides a value of ∆MOS between CRT and LCD per sequence, an objective metric is developed in order to predict this quality loss ∆MOS p from the motion blur perception model. This metric is made in several steps. First, a spatio-temporal classification is done in two passes. First pass is corresponding to a block based motion estimation that leads to the construction of tubes which are the sets of blocks positions along the direction of motion. Second pass is the classification of each tube according to its spatial content. Since motion blur is only visible at sufficient contrast [START_REF] Laird | Spatio-velocity CSF as a function of retinal velocity using unstabilized stimuli[END_REF], only tubes categorised as textures and edges are selected. An average motion vector is computed from all the vectors of the remaining tubes. Norm of this global vector is used to compute the width of perceived motion blur according to Equation 2. This value W is an indicator of the average magnitude of perceived blur along the sequence. Finally, ∆MOS p is computed from a function of W . This function is non linear since there is no influence on perceived quality below a threshold of W , and the quality difference saturates for high values (cf. Figure 2).

An estimation of the subjective quality scores on LCD from the subjective quality scores on CRT can be made using the following relation: The quality of the model can be measured by the linear correlation coefficient (CC) and the root mean square error (RMSE) between estimated LCD scores and actual LCD scores. Values of 0.958 for CC and 1.30 for RMSE are obtained. These performances can be compared with those ob-tain comparing results between CRT and LCD: 0.241 for CC and 6.75 for RMSE.

MOS LCD est = MOS CRT -∆MOS p . (3)
Therefore, the proposed prediction model is able to evaluate solutions to LCD motion blur reduction from obtained quality on CRT.

IMPACT OF LCD MOTION BLUR WITH RESPECT TO THE DISPLAY RESOLUTION

As it has been shown in the previous section, a high correlation exists between the magnitude of perceived blur in a sequence and the loss of quality observed on LCD when displaying this sequence. Furthermore, the magnitude of perceived blur depends on display characteristics and resolution. More precisely in Equation 2, the velocity V of motion would be less important if the resolution is lower. As a consequence, the perceived blur would be smaller and should be less annoying as the resolution decreases.

Standard definition video quality assessment

In order to assess the impact of LCD motion blur relatively to the resolution, the same video quality assessment tests have been realised with sequences at a standard definition (SD). Four sequences of the previous experiment have been chosen and reduced to SD resolution by computing the HD versions through a half-band filtering followed by a down-sampling by a factor of 2 (both along horizontal and vertical directions). This processing is performed on each field of the interlaced 1080i sequences, with the advantage that it doesn't necessitate any interpolation.

Tests have been led both on the CRT display and then on the LCD. Viewing conditions and display parameters are the same as those described in Section 2.1, except for the viewing distances which have been set to six times the pictures'height (6H). SD sequences have been displayed inserted in a HD resolution grey level sequence in order to suit the displays native HD resolution.

Each reference (uncompressed) has been distorted with the H.264 reference coder at seven different bit-rates. The impact of technology on the perceived quality is assessed. As a result, only the sequences with few or no coding distortions are taken in account: the mean of MOS for sequences whith good to excellent quality is computed.

Results

Table 5 shows the mean of MOS for good-to-excellent quality HD sequences, on the two displays, and ∆MOS is computed. Same results for SD sequences is shown in Table 6. As expected, the loss of quality on LCD is strongly smaller in SD resolution with respect to HD resolution. Moreover, the loss of quality on SD sequences is not significant with regards to intervals of confidence.

As explained before, the perception of LCD motion blur closely depends on the display resolution. The quantity of perceived blur linearly depends on the velocity of motion which is twice smaller in SD than in HD (since the resolution is divided by two). As a result, the perceived motion blur should be less annoying in SD and the advantages of LCD such as colourfulness and large luminance range seem to tower over this artifact leading to a better global quality on LCD display.

CONCLUSION

Subjective quality evaluation of HD moving pictures have shown that the perceived quality is better on CRT display than on LCD. This loss of quality ∆MOS seems to be due to the flat panel technology. Actually, new artifacts such as motion blur are very annoying with quick movements. Benefits of LCD, for instance colourfulness and larger luminance range, have been highlighted with subjective assessment of still pictures: they lead to a better perceived quality on LCD than on CRT. However, in video they don't achieve to compensate the loss of quality due to moving artifacts. LCD motion blur have been studied and a mathematical model is used to measure its magnitude as a function of the quantity of movements. This magnitude depends on the display resolution in the same way as the velocity of moving objects. A high correlation has been highlighted between this motion blur magnitude and the loss of quality on LCD which enables the prediction of the loss of quality ∆MOS between CRT and LCD.

Consequently, the loss of quality ∆MOS due to LCD motion blur would depend on the display resolution and should be weaker with lower resolutions than HD. Same video quality assessment tests with SD sequences have confirmed this. At a lower resolution, the moving artifacts due to LCD technology are less annoying. Benefits of LCD such as colourfulness and larger luminance range seem to tower over this defects since the perceived quality is better on LCD in SD resolution.

In this paper, it's shown that the new LCD technology leads to new shortcomings when displaying moving pictures. These new artifacts are not significant at low resolutions, which validate the use of LCD for subjective video quality assessment in the Multimedia Testplan of the VQEG [START_REF] Vqeg | Multimedia group test plan[END_REF]. However, when increasing the resolution, artifacts due to flat panel technology become more annoying and have an important impact on the perceived quality. It comes that the subjective video quality assessment at high resolution (HDTV for example) should be led very carefully on LCD since a significant part of perceived distortions could be due to the display.

Figure 1 :

 1 Figure 1: Perceived blur width W as a function of motion velocity V .

Figure 2 :

 2 Figure 2: Prediction of the loss of quality ∆MOS p from the average blur magnitude.

  . Viewing distances have been set to 3H, where H is the height of the displayed pictures.

	Background luminance	35 cd/m 2
	Chromaticity of background	D 65
	Pictures height (H)	46 cm
	Viewing distance	138 cm (3H)
	Display black luminance	0.641 cd/m 2
	Display peak luminance	471 cd/m 2
		Background luminance	7 cd/m 2
		Chromaticity of background	D 65
		Picture height (H)	20.5 cm
		Viewing distance	61.5 cm (3H)
		Display black luminance	0.531 cd/m 2
		Display peak luminance	70.9 cd/m 2
		Table 1: Viewing conditions and display parameters for the
		CRT display JVC DT-V 1910CG.
		2.2 Protocol
		These tests have been performed with sequences and pictures
		of rather good to excellent quality. As a consequence, the
		used protocol should enable the quality discrimination. A
		well known stable method for this purpose is the SAMVIQ
		protocol [5], developed by France Telecom RD and standard-
		ised by the EBU and the ITU.
		SAMVIQ is multi stimuli continuous quality scale pro-
		tocol. With this procedure, the observers can compare the

Table 2 :

 2 Viewing conditions and display parameters for the LC display Philips T370 HW01.

Table 3 :

 3 Mean opinion scores by sequences and displays.

		MOS CRT MOS LCD ∆MOS
	VOILE	83.9	77.7	6.2
	FOOT	82.8	76.3	6.5
	CONCERT	84.5	73.8	10.7
	SHOW	82.9	75.3	7.6
	CREDITS	83.1	79.1	4.0
	MOBCAL	81.4	81.0	0.4
	PARKRUN	87.6	80.2	7.4
	SHIELDS	86.7	78.2	8.5
	STOCKHOLM	86.1	82.3	3.8

Table 4 :

 4 . Mean opinion scores by pictures and displays.

	Séquence	MOS CRT MOS LCD ∆MOS
	FOOTBALL	66.3	79.4	-13.1
	HAND	73.6	80.3	-6.7
	HOUSE	51.8	81.8	-30.0
	LANDSCAPE	73.5	78.7	-5.2
	MAP	51.4	84.4	-33.0

Table 5 :

 5 Mean of the MOS for good-to-excellent quality sequences in HD resolution, on CRT and LCD.

	Séquence	MOS CRT MOS LCD ∆MOS
	MOBCAL	79.5	71.9	7.6
	PARKRUN	83.0	70.7	12.3
	SHIELDS	81.4	67.9	13.5
	STOCKHOLM	81.6	73.0	8.7
	Séquence	MOS CRT MOS LCD ∆MOS
	MOBCAL	71.6	68.0	3.6
	PARKRUN	77.6	72.9	4.7
	SHIELDS	75.5	73.5	2.0
	STOCKHOLM	75.7	73.2	2.5

Table 6 :

 6 Mean of the MOS for good-to-excellent quality sequences in SD resolution, on CRT and LCD.