xuepso

函数基本用法

matlab

psokinv('配置文件.cfg') % 基本用法，只需要配置文件

完整的使用步骤

1. 首先生成配置文件

matlab

% 使用 sim_invconfig 函数生成配置文件模板 sim_invconfig('my_config.cfg');

2. 编辑配置文件

• 打开生成的.cfg文件（如my_config.cfg）

• 根据您的需求修改参数，特别是断层模型相关设置

• 可能需要设置 InSAR 数据路径、反演参数等

3. 运行反演函数

matlab

% 方法1：仅使用配置文件（UTMZONE 从数据中自动获取） psokinv('my_config.cfg'); % 方法2：指定 UTM 区域 psokinv('my_config.cfg', 50); % 假设 UTM 区域为 50 % 方法3：指定 UTM 区域和 CPU 核心数 psokinv('my_config.cfg', 50, 8); % UTM 区域 50，使用 8 个 CPU 核心

参数说明：

• <configure_file>.cfg: 必需参数，配置文件路径

• [UTMZONE]: 可选，UTM 投影区域号

  • 可省略（函数会尝试从 InSAR 数据文件名中提取）

  • 或直接指定数字（如 50 表示 UTM 50 区）

• [numcpus]: 可选，使用的 CPU 核心数

  • 用于并行计算加速

  • 省略时可能使用默认值

使用示例：

matlab

% 1. 生成配置文件 sim_invconfig('inv_setup.cfg'); % 2. 手动编辑 inv_setup.cfg 文件（用文本编辑器） % 设置数据路径、反演参数、断层约束等 % 3. 运行反演 psokinv('inv_setup.cfg', 48, 4); % UTM 48区，4核并行

注意：配置文件中的断层模型参数通常需要根据具体研究区域进行编辑才能得到合理的结果。

这个配置文件是为 PS-OKINV 断层反演程序准备的。让我详细解释如何根据您的需求修改这个配置文件：

配置文件结构解析

1.弹性参数部分

cfg

# elastic constant: lambda 3.23e+10 # elastic constant: mu 3.23e+10

• 作用：定义地壳的弹性参数

• 建议：根据研究区域的地壳结构调整

  • 常用值范围：λ = μ ≈ 30-40 GPa

  • 对于一般地壳：3.23e10 Pa（约32.3 GPa）

2.位移单位和投影

cfg

# unit of displacements: m,cm or mm m # utmzone of displacements, e.g. 19Q 30Q

• 位移单位：根据您的 InSAR 数据单位选择m、cm或mm

• UTM 区域：30Q表示 UTM 30区（西经6°到0°），Q表示北半球

3.反演算法参数

cfg

# number of particles 200 # PSO: maximum restart number 1 # PSO: total iteration number (default: 25) 25

• 粒子数：粒子群优化算法中的粒子数量（建议：100-500）

• 迭代次数：反演迭代次数（复杂问题可增加到 50-100）

4.输入文件设置

cfg

#No Inv(0/1) A&&B C Weight VCM MC(0/1) PathFileName 1 0 0 0 1.000 NULL 0 insar.inp

• 关键修改：insar.inp需要改为您的实际 InSAR 数据文件路径

• 格式：InSAR 数据文件通常包含：东坐标、北坐标、高程、视线向位移等

5.断层参数设置（最重要部分）

您目前看到的都是零值，需要根据地质资料设置合理的初始值和约束范围：

cfg

# Value MinV MaxV Inv(0/1) Symbol Parameters of fault 1 --------------------------------------------------------------------------- 0.00 0.00 0.00 1 f(1,1) x-start(km) 0.00 0.00 0.00 1 f(1,2) y-start(km) 0.00 0.00 0.00 1 f(1,3) strike(degree) 0.00 0.00 0.00 1 f(1,4) dip(degree) 0.00 0.00 0.00 1 f(1,5) depth(km) 0.00 0.00 0.00 1 f(1,6) width(km) 0.00 0.00 0.00 1 f(1,7) length(km) 90.00 -135.00 45.00 1 RakeCons Rake(degree) 0.00 0.00 0.00 0 MagnCons MW(mag)

如何修改断层参数示例

假设您要反演一个走滑断层，参数如下：

• 位置：UTM 东坐标 500000 m，北坐标 4000000 m

• 走向：30°，倾角：90°（垂直断层）

• 深度：5 km，宽度：10 km，长度：20 km

• 滑动角：-10°（右旋走滑）

修改后的配置示例：

cfg

# Value MinV MaxV Inv(0/1) Symbol Parameters of fault 1 --------------------------------------------------------------------------- 500.00 490.00 510.00 1 f(1,1) x-start(km) 4000.00 3990.00 4010.00 1 f(1,2) y-start(km) 30.00 0.00 60.00 1 f(1,3) strike(degree) 90.00 60.00 90.00 1 f(1,4) dip(degree) 5.00 2.00 15.00 1 f(1,5) depth(km) 10.00 5.00 20.00 1 f(1,6) width(km) 20.00 10.00 30.00 1 f(1,7) length(km) -10.00 -30.00 10.00 1 RakeCons Rake(degree) 0.00 0.00 0.00 0 MagnCons MW(mag)

修改步骤

1. 确定 InSAR 数据文件

   • 准备 InSAR 数据文件（如my_insar_data.inp）

   • 格式参考：x(E) y(N) height disp_LoS uncertainty

2. 更新输入文件路径

   cfg

#No Inv(0/1) A&&B C Weight VCM MC(0/1) PathFileName 1 0 0 0 1.000 NULL 0 /path/to/my_insar_data.inp

1. 设置断层初始值

   • 根据地震震中、余震分布、地质资料等

   • 设置合理的Value（初始值）

   • 设置合理的MinV和MaxV（搜索范围）

2. 设置反演参数

   • Inv(0/1)：1 表示反演该参数，0 表示固定该参数

   • 通常除矩震级（MW）外，其他参数都设为 1

3. 调整算法参数

   • 复杂模型：增加粒子数（300-500）和迭代次数（50-100）

   • 简单模型：使用默认值即可

运行示例

修改完成后：

matlab

% 运行反演 psokinv('my_config.cfg'); % 或者指定CPU核心数 psokinv('my_config.cfg', [], 4); % 使用4个CPU核心

注意：配置文件中的COORTYPE: 0表示使用 UTM 坐标系统，坐标单位为公里。

10