LogoRefATT

Input Parameters

Input parameters are specified in a YAML file. The following is full list of input parameters:

#################################################
#                generate mesh3D                #
#################################################
mesh:
  min_max_x: [30, 32] # x in m
  min_max_y: [30, 32] # y in m
  min_z: -50000 # depth in m
  n_xy: [61, 61] # number of nodes in x,y,dep direction
  interface_model_path: ./models/interface.h5 
  ngrid_layers: [61] # number of grid at each layer in z direction [n_bottom, n_top]
  suppress_utm: False
  utm_projection_zone: 36

#################################################
#              generate database                #
#################################################
database:
  model_type: 2 # model type: 0 for 3d model with uniform grid in hdf5 or 1 for 1d model in txt, 2 for 3d model with fitted grid
  input_model_path: ./models/initial_model.h5
  # input_model_path: ./input/vp1d.vel
  z_interface: [1000] # interface depth in meter for 1d model

#################################################
#            traveltime data file path          #
#################################################
source:
  src_rec_file: ./OUTPUT_FILES/src_rec_file_forward.dat # source receiver file path
  swap_src_rec: true # swap source and receiver
  source_depth_in_km: false

#################################################
#            parallel computation settings      #
#################################################
parallel: # parameters for parallel computation
  n_sims: 8 # number of simultanoues runs (parallel the sources)
  ndiv_xy: [1, 1] # number of subdivision on each direction (parallel the computional domain)

############################################
#            output file setting           #
############################################
output_setting:
  output_dir: OUTPUT_FILES/ # path to output director (default is ./OUTPUT_FILES/)
  output_initial_model:    true # output initial model (initial_model.h5) or not.                                   
  output_final_model:      true # output merged final model (final_model.h5) or not.                                
  output_in_process_model:       false # output model at each inv iteration or not.                       
  output_in_process_data:  true # output src_rec_file at each inv iteration or not. 
  logger_level: 1       # 0 for debug, 1 for info, 2 for warning, 3 for error, 4 for critical             

#################################################
#          inversion or forward modeling        #
#################################################
# run mode
# 0 for forward simulation only,
# 1 for inversion
run_mode: 1

###################################################
#          model update parameters setting        #
###################################################
model_update:
  n_inversion_grid: 5 # number of inversion grid sets
  n_inv_xyz: [11, 11, 12]
  max_iterations: 10 # maximum number of inversion iterations
  conv_tol: 0.0001
  optim_method: 1 # 0 for steepest descent, 1 for L-BFGS
  step_length: 0.02 # initial step length of model perturbation. 0.01 means maximum 1% perturbation for each iteration.
  max_shrink: 0.618 # shrink factor for step length (0 - 1)
  max_sub_iterations: 20
  c1: 0.001 # parameter for Armijo condition
  c2: 0.9 # parameter for curvature condition