An inverse problem methodology to identify flow channels in fractured media using synthetic steady-state head and geometrical data

被引:24
作者
Le Goc, R. [1 ,2 ]
de Dreuzy, J.-R. [2 ]
Davy, P. [2 ]
机构
[1] Itasca Consultants SAS, F-69134 Ecully, France
[2] Univ Rennes 1, CNRS, UMR 6118, F-35042 Rennes, France
关键词
Inverse problem; Simulated annealing; Fractured media; Flow channels; NUMERICAL INVERSION; FLUID-FLOW; IDENTIFICATION; TRANSPORT; ZONE; OPTIMIZATION; CALIBRATION; VALIDATION; PARAMETERS; NETWORK;
D O I
10.1016/j.advwatres.2010.04.011
中图分类号
TV21 [水资源调查与水利规划];
学科分类号
081501 ;
摘要
We present a methodology for identifying highly-localized flow channels embedded in a significantly less permeable medium using steady-state head and geometrical data. This situation is typical of fractured media where flows are often strongly channeled at the scales of interest (10 m-1 km). The objective is to identify both geometrical and hydraulic characteristics of the conducting structures. Channels are identified in decreasing order of importance by successive optimizations of an objective function. The identification strategy takes advantage of the hierarchical flow organization to restrict the dimension of the solution space of each individual optimization step. The characteristics of the secondary channels are strongly determined by the main flow channels. The latter are slightly modified by the secondary channels through the addition of a regularization term to the main channel characteristics in the objective function. As the objective function is strongly non-convex with numerous local minima, inversion is performed using a stochastic algorithm (simulated annealing). We assess the possibilities of the hierarchical identification strategy on simple synthetic steady-state flow configurations where hydraulic data are made up of 25 regularly spaced heads and of the boundary conditions. Those flow structures that are dominated by at most two simple channels can be identified with these head data only. Configurations comprising up to three complex and interconnected channels can still be identified with additional geometrical information including the distances of piezometers to their closest channel. The capabilities of the hierarchical identification strategy are limited to flow structures dominated by at most three equivalent flow channels. We finally discuss the perspectives of application of the method to transient-state data obtained on a more restricted number of piezometers. (C) 2010 Elsevier Ltd. All rights reserved.
引用
收藏
页码:782 / 800
页数:19
相关论文
共 56 条
  • [1] [Anonymous], 1987, METHODS DATA FITTING
  • [2] A travel time based hydraulic tomographic approach
    Brauchler, R
    Liedl, R
    Dietrich, P
    [J]. WATER RESOURCES RESEARCH, 2003, 39 (12) : SBH201 - SBH2012
  • [3] An inverse problem from 2D ground-water modelling
    Bruckner, G
    Handrock-Meyer, S
    Langmach, H
    [J]. INVERSE PROBLEMS, 1998, 14 (04) : 835 - 851
  • [4] Using the migration of the induced seismicity as a constraint for fractured Hot Dry Rock reservoir modelling
    Bruel, Dominique
    [J]. INTERNATIONAL JOURNAL OF ROCK MECHANICS AND MINING SCIENCES, 2007, 44 (08) : 1106 - 1117
  • [5] BRUYELLE J, 2009, P 11 ANN C GEN EV CO
  • [6] BRUYELLE J, 2009, 71 EAGE C EXH 8 11 J
  • [7] MODELING FRACTURE FLOW WITH A STOCHASTIC DISCRETE FRACTURE NETWORK - CALIBRATION AND VALIDATION .1. THE FLOW MODEL
    CACAS, MC
    LEDOUX, E
    DEMARSILY, G
    TILLIE, B
    BARBREAU, A
    DURAND, E
    FEUGA, B
    PEAUDECERF, P
    [J]. WATER RESOURCES RESEARCH, 1990, 26 (03) : 479 - 489
  • [8] Inverse problem in hydrogeology
    Carrera, J
    Alcolea, A
    Medina, A
    Hidalgo, J
    Slooten, LJ
    [J]. HYDROGEOLOGY JOURNAL, 2005, 13 (01) : 206 - 222
  • [9] ESTIMATION OF AQUIFER PARAMETERS UNDER TRANSIENT AND STEADY-STATE CONDITIONS .3. APPLICATION TO SYNTHETIC AND FIELD DATA
    CARRERA, J
    NEUMAN, SP
    [J]. WATER RESOURCES RESEARCH, 1986, 22 (02) : 228 - 242
  • [10] CHANG LC, 2005, WATER RESOUR RES, P41