# # Simulation "A" from Barbour, Norbeck, and Rubinstein (2017, SRL) # # For use with poel version 2012 by R. Wang (GFZ) # e-mail: wang@gfz-potsdam.de # ############################################################## ## Conventions ## * Cylindrical coordinates (Z positive downwards!) are used ## * If not others specified, SI Unit System is used overall ## * Positive strain is extensional ## * Tilt is positive when the upper end of a borehole tiltmeter body moves away from the pumping well ## * Input file encoding: US-ASCII ## (use in unix, for example, to convert) ############################################################## # # SOURCE PARAMETERS A: SOURCE GEOMETRY # ==================================== #------------------------------------------------------------------------------- 1335.0 1920.0 |dble: s_top_depth, s_bottom_depth; 0.1 |dble: s_radius; #------------------------------------------------------------------------------- # # SOURCE PARAMETERS B: SOURCE TYPE # ================================ # 1 = injection (bvp) #------------------------------------------------------------------------------- 1 |int: sw_source_type; #------------------------------------------------------------------------------- # # SOURCE PARAMETERS C: SOURCE TIME HISTORY # ======================================== #------------------------------------------------------------------------------- 47 |int: no_data_lines; #------------------------------------------------------------------------------- # no time source_function (+ = injection, - = pumping) # [-] [s] [m^3/s] #------------------------------------------------------------------------------- 1 0 0.00e+00 2 0 2.56e-03 3 2678400 2.16e-02 4 5097600 2.45e-02 5 7776000 2.20e-02 6 10368000 1.42e-02 7 13046400 2.45e-02 8 15638400 3.40e-02 9 18316800 2.96e-02 10 20995200 3.35e-02 11 23587200 3.84e-02 12 26265600 2.62e-02 13 28857600 2.48e-02 14 31536000 3.08e-02 15 34214400 2.19e-02 16 36633600 2.77e-02 17 39312000 2.64e-02 18 41904000 2.38e-02 19 44582400 1.90e-02 20 47174400 1.88e-02 21 49852800 1.91e-02 22 52531200 9.89e-03 23 55123200 9.30e-03 24 57801600 9.78e-03 25 60393600 1.23e-02 26 63072000 3.99e-03 27 65750400 5.90e-03 28 68169600 5.59e-03 29 70848000 4.93e-03 30 73440000 6.78e-03 31 76118400 6.09e-03 32 78710400 3.83e-03 33 81388800 3.62e-03 34 84067200 3.77e-03 35 86659200 3.69e-03 36 89337600 3.87e-03 37 91929600 3.78e-03 38 94608000 7.92e-03 39 97286400 5.16e-03 40 99792000 7.78e-03 41 102470400 7.58e-03 42 105062400 5.90e-03 43 107740800 5.80e-03 44 110332800 3.83e-03 45 113011200 3.75e-03 46 115689600 3.24e-04 47 115689600 0.00e+00 # ################################################################################ # # RECEIVER PARAMETERS A: RECEIVER DEPTH SAMPLING # ============================================== #------------------------------------------------------------------------------- 1 |int: sw_receiver_depth_sampling; 36 |int: no_depths; 0.0 7000.0 |dble: zr_1,zr_n; or zr_1,zr_2,...,zr_n; #------------------------------------------------------------------------------- # # RECEIVER PARAMETERS B: RECEIVER DISTANCE SAMPLING # ================================================= #------------------------------------------------------------------------------- 1 |int: sw_receiver_distance_sampling; 41 |int: no_distances; 0.0 10000.0 |dble: r_1,r_n; or r_1,r_2,...,r_n; #------------------------------------------------------------------------------- # # RECEIVER PARAMETERS C: Time SAMPLING # ==================================== #------------------------------------------------------------------------------- 1.2096e+08 |dble: time_window; 100 |int: no_time_samples; ################################################################################ # # WAVENUMBER INTEGRATION PARAMETERS # ================================= # relative accuracy (0.01 for 1% error) for numerical wavenumber integration; #------------------------------------------------------------------------------- 0.015 |dble: accuracy; ################################################################################ # # OUTPUTS A: DISPLACEMENT TIME SERIES # =================================== # (change 0 to 1 to turn on output, and 1 to 0 to turn off output) #------------------------------------------------------------------------------- 1 1 |int: sw_t_files(1-3); 'uz.t' 'ur.t' |char: t_files(1-3); #------------------------------------------------------------------------------- # # OUTPUTS B: STRAIN TENSOR & TILT TIME SERIES # =========================================== #------------------------------------------------------------------------------- 1 1 1 1 0 |int: sw_t_files(3-7); 'ezz.t' 'err.t' 'ett.t' 'ezr.t' 'tlt.t' |char: t_files(3-7); #------------------------------------------------------------------------------- # # OUTPUTS C: PORE PRESSURE & DARCY VELOCITY TIME SERIES # ===================================================== #------------------------------------------------------------------------------- 1 0 0 |int: sw_t_files(8-10); 'pp.t' 'dvz.t' 'dvr.t' |char: t_files(8-10); #------------------------------------------------------------------------------- # # OUTPUTS D: SNAPSHOTS OF ALL OBSERVABLES # ======================================= #------------------------------------------------------------------------------- 1 |int: no_sn; 10.0 'snapshot0010.dat' |dble: sn_time(i),sn_file(i), i=1,2,... ################################################################################ # # GLOBAL MODEL PARAMETERS # ======================= # surfcon: 0 == whole space, 1 == unconfined free surface (p=0), 2 == confined free surface (dp/dz=0) #------------------------------------------------------------------------------- 1 |int: isurfcon 15 |int: no_model_lines; #------------------------------------------------------------------------------- # # MULTILAYERED MODEL PARAMETERS # ============================= # no depth[m] mu[Pa] nu nu_u B D[m^2/s] #------------------------------------------------------------------------------- 1 0.00 1.0E+10 0.30 0.40 0.75 0.1000 2 1100.00 1.0E+10 0.30 0.40 0.75 0.1000 3 1100.00 1.5E+10 0.30 0.40 0.75 0.0005 4 1300.00 1.5E+10 0.30 0.40 0.75 0.0005 5 1300.00 1.5E+10 0.25 0.40 0.75 1.0000 6 2000.00 1.5E+10 0.25 0.40 0.75 1.0000 7 2000.00 3.0E+10 0.25 0.38 0.75 1.0000 8 3300.00 3.0E+10 0.25 0.38 0.75 0.236 9 3300.00 3.0E+10 0.25 0.38 0.75 0.236 10 4400.00 3.0E+10 0.25 0.38 0.75 0.063 11 4400.00 3.0E+10 0.25 0.38 0.75 0.063 12 5600.00 3.0E+10 0.25 0.38 0.75 0.016 13 5600.00 3.0E+10 0.25 0.38 0.75 0.016 14 8000.00 3.0E+10 0.25 0.38 0.75 0.002 15 8000.00 3.0E+10 0.25 0.38 0.75 0.002 #------------------------------------------------------------------------------- ###########################end of all inputs####################################