APEX JASON
StatMod
只适用于三维工区。 主要算法：
1. Sequential Gaussian Simulation （序贯高斯模拟） SGS 2. Sequential Gaussian Collocated Co-simulation (序贯高斯配置协模拟) SGCCS 3. Sequential Gaussian Co-simulation (序贯高斯协模拟) SGCS 4. Sequential Indicator Simulation （序贯指示模拟,岩性模拟） SIS
5. Lithology masks (遮挡岩性指示模拟)
6. Sequential Indicator Simulation with a trend (带趋势的序贯指示模拟) SIS with trend 7. Stochastic inversion (随机反演)
一 、Sequential Gaussian Simulation (序贯高斯模拟） :
----只用井数据作一种曲线的随机模拟
例如对某一层（layer)，用多井的孔隙度，进行 孔隙度的随机模拟。 一般需要15-20口井才会有效果 分析：porosity first histogram first variogram (一) StatMod ? StatMod analysis …
1、Input :
(1) Variograms / transforms file … ? ( 键入文件名）simu.var (2) Trace gate …
(3) Time / Depth mode … ? Time
(4) Solid model … ? tdc1 ( 选EarthModel的结果)
(5) Layers … ? 选层(layer),例如：b.c层
(6) Data for histograms and transforms ? (作直方图分析） Well logs data … ? Model data … Horizon data …
a) 选井文件：（选多个时域井文件）
Log type : ( 选单一井曲线，如：porosity ) 出现 StatMod Transform 窗口 b) Histogram ?????? ? Nr of intervals : 10 (直方图条数）
Build the histogram
histogram
? 出现直方图?????Done
c) ????????? Use automatic fit Transform ???????????? Gaussion ? Type Log – Ganssion （选任一种来拟 By table 合分布函数）
d) Save as first … ? porosity_hist ( 存入simu.var )
? Assign current transform to selected layers ? Ok ? Done
e) Line width （置线条宽度）
StatMod 1
APEX JASON
Title （写图名） Done
(7)Data for variogram sampling and modeling --（作变异图分析）
? Primary data … ? ? Well log data … ? Secondary data … Model data … Horizon data …
a) 选多个时域井文件
选井曲线( porosity )，出现 StatMod variogram 窗口 b) Sample vario …
? Function Variogram ? Cross-vario ? Madogram Rodogram 注意： Threshold # of lag Lag intervel Indicator
Build the sample variogram ? (出现变异图）
c) 可以扫描岩性的横向与纵向变化： Select search parameters : Anisotropy : Isotropic （各向同性） Anisotropic （各向异性）
Scan ? 扫描储层的大致走向方向） Lateral Direction ( Azimuth )
# of angles Angle Tolerance Bandwidth (m) 5 10 15 5000
Vertical Direction ( Dip ):
Tolerance Bandwidth ( m )
15 1000
d) Anisotropy : ? Anisotropic Build the sample variogram e) ? Ok f)
Model vario …
Func Type: Nugget Gaussian 可选 Gaussian 或 Spherical Spherical ? 来拟合数据点 Expo Gen-Expo Rel-Quad StatMod 2 （
APEX JASON
填写 sill 值与 x.y.z ( lag 值),可以反复调。在 Gaussian simulation中，Z是主要的。 Sill Gaussion Z X Y
其中： Z: 纵向砂体的大致厚度
X: 在走向方向上砂体的大致分布范围 (根据地质分析) Y: 在倾向方向上砂体的大致分布范围 (根据地质分析) 操作步骤：
?填写 sill 值与 x.y.z ( lag 值) ? Add
?若不合适，点亮已有值（在最上边） ?修改 sill 值与 x.y.z ( lag 值) ? Change
? 反复调整
g) Direction : Anisotropic h) Save as first … ? porosity_var (存入 saimu.var )
? Done
? Assign current variogram to selected layers ? Done ? Ok
!! 此时： 直方图 porosity_hist
simu.var 变异图 porosity_var i) Line width … 线条宽度
Title 图头 Done
2、File ? Save & Exit
(二) StatMod ? StatMod modeling …
1、Input :
(1) Simulation mode …
? Run property simulation Mode : Gaussian simulation ? Gaussian cosimulation Gaussian collocated cosimulation Threshold indicator simulation Done
StatMod 3
APEX JASON
Time / Depth mode … ? Time
Solid model … ? 用EarthModel的结果
Layers … ? 选层，例如：b.c层,给采样间隔 Primary input data … ? Well log data … 选多个时域井文件
选单一井曲线 (如：porosity)
（6) Variogram / transforms file … ? 选 simu.var
(7) Trace gate … 2、Edit:
(1) Realization parameters … （实现参数）
Nr of realization : 10 ----输出多个等概率的随机模拟结果（例如：10个）
Output sample interval (s) : 0.001
? Perform kriging also （简单克里金）
? Use ordinary kriging （普通克里金）
(2) Search parameters …
Search radius (m) : 1000 ----扫描半径，要给大一些， 大于井间距，否则不计算。 While kriging
Number of neighbors per octant : 4
? Use irregular search on primary data
While simulation
Number of neighbors per octant : 2 Seed of random generator : 1 Use same path Always ?
For fine grid only Never ? Snap well to grid for fast performance
(3) Transforms … ?
Transform …
选层（例如：b.c)与对应的直方图参数 porosity_hist
(4) Variogram … ? Variogram …
? 层与对应的变异图参数 porosity_var
(6) Super codes …
Layer Super code b b c c
3、Output ? Generate …
? Run now
Output type : Model ? Horizon
Primary output
Primary file name : ssm_bc_porosity_simu.mod
? Generate primary kriging error : ssm_bc_porosity_simu_error.mod
并产生文件：（由于 Nr of realizations : 10 所以产生以下十个文件） ssm_bc_porosity_simu__01.mod ssm_bc_porosity_simu__02.mod StatMod 4 (2)
(3) (4) (5)
APEX JASON
ssm_bc_porosity_simu__03.mod 对孔隙度的10次随机 …… 模拟结果（等概率） ssm_bc_porosity_simu__10.mod
4、File ? Save & Exit
二、Sequential Gaussian Collocated Co-Simulation (序贯高斯配置协模拟): ------用井数据和阻抗体协同进行随机模拟，并与井匹配
1、First ( porosity) -------- Histogram porosity_hist 2、 First ( porosity) -------- Variogram porosity_var 3、 Second (AI) --------- Histogram aitm_hist 4、 两类数据的相关系数 (通过交会图分析)
不需作互变异图(cross variogram)，软件通过两类数据的相关系数实现协模拟。 (一) Statmod Analysis …
第一类数据的直方图与变异图已分析过了，可以应用，不必再分析。 下边做第二类数据的分析：
1、Input :
(1) Variograms / transforms file … ? ( 选文件名）simu.var (2) Trace gate …
(3) Data for histogram and transforms ? Model data …
(a) ? aitm_ai.mod ? None (b) Histograms … Build the histograms ?
? Done
(c) Transforms …
Type: ? Use automatic fit Gaussian
(d) Save as first … ? aitm_hist ? Ok (存入simu.var) ? Assign current transform to selected layers ? Done
(e) Done
2、File ? Save & Exit (二）StatMod modeling …
1、Input :
(1) Simulation mode …
? Run property simulation Mode : Gaussian simulation Gaussian cosimulation Gaussian collocated cosimulation ? Threshold indicator simulation (2) Time / Depth mode … ? Time (3) Solid model … ? tdc1
StatMod 5
APEX JASON
(4) Layers … ? 选层与采样间隔
(5) Primary input data … ? Well log data … ? 选井与曲线(porosity) (6) Secondary input data … ? Model data … ? aitm_ai.mod (7) Variogram / transforms file … ? simu.var (8) Trace gate …
2、Edit :
(1) Realization parameters …
Nr of realization : 10
Output sample interval (s) : 0.001
? Perform kriging also
? Use ordinary kriging
? Use cokriging second non-bias condition
(2) Search parameters … Global
Search radius (m) : 10000 While kriging
Number of neighbora per octant : 4
? Use irregular search on primary data
While simulation
Number of neighbors per octant : 2 Seed of random generator : 1 Use same path Always ?
For fine grid only Never ? Snap well to grid for fast performance
(3) Transforms … ? ? porosity_hist
Transform …
Second Transform … ? ? aitm_hist
(4) Variogram … ? ? porosity_var
Variogram …
(5) Super codes …
Layer Super code b b c c
(7) Correlation coefficients …
Layer Value
b 0.7 相关系数可通过交汇图大致给出 c 0.7
3、Output ? Generate … ? Run now
Output type : Model ? Horizon
Primary output
Primary file name : ssm_bc_SGCCS_first.mod
? Generate primary kriging error : ssm_bc_ SGCCS _first_error.mod 并产生文件：（由于 Nr of realizations : 10 ,所以产生以下十个文件）
StatMod 6
APEX JASON
ssm_bc_ SGCCS _first_01.mod ssm_bc_ SGCCS _first_02.mod 对孔隙度的配置协模拟结果
ssm_bc_ SGCCS _first_03.mod …… Secondary Output
? Generate secondary data : ssm_bc_ SGCCS _second.mod
? Generate primary kriging error : ssm_bc_ SGCCS _second_error.mod 并产生以下十个文件 :
ssm_bc_SGCCS_second_01.mod
ssm_bc_SGCCS_second _02.mod 对阻抗的配置协模拟结果
ssm_bc_SGCCS_second _03.mod ……
4、File ? Save & Exit
三、Sequential Gaussian Co-Simulation : (序贯高斯协模拟)
Transforms first : porosity (well) Variogram Transforms Cross Variogram second : aitm (cube) Variogram
(一) SdatMod analysis …
第一类数据的直方图与变异图和第二类数据的直方图均已分析完毕， 下边需做第二类数据的变异图与两类数据的交会变异图。 1、第二类数据的变异图（略） 2、Input :
(1) Variograms / transforms file … ? ( 选文件名）simu.var
(2) Data for variogram sampling and modeling ? Primary data … ? Well data …
(a) *.wll (b)
Sample
vario … … Function :
Cross variogram
First data … ? porosity_hist ? Ok Second data … ? aitm_hist ? Ok
(c)
Show the Cross Plot
ai
porosity StatMod 7
APEX JASON
Fit ?Edit polygons
画 polygon ? Save as … ? porosity_aitm.tab ? Done (d) Builder the sample variogram ? Ok
(e) Model Vario …
Func Type: Gaussian (可选择球形、指数） 填写sill （注意负值）, lag ( x, y, z) ， Save as first … ? porosity_aitm_cross_var (simu.var) ? Assign current variogram to selected laters ? Ok
1、File ? Save &Exit
(二) Statmod modeling …
1、Input :
(1) Simulation mode … ?
? Run Run property simulation Mode : Gaussian simulation Gaussian cosimulation ?
Gaussian collocated cosimulation Threshold indicator simulation (2)Time /Depth mode … ? Time (3) Solid model … ? tdc1 (4) Layrts … ? b,c (5) Primary input data … ? Well log data … ? porosity (6) Secondary input data … ? Model data … ? aitm_ai.mod (7) Variograms / transforms file … ? simu.var (8) Trace gate …
2. Edit :
(1) Realization parameters …
Nr of realization : 10
Output sample interval (s) : 0.001 ? Perform kriging also
? Use ordinary kriging
? Use cokriging second non-bias conditio
(2) Search parameters …
Global
Search radius (m) : 10000 While kriging
Number of neighbora per octant : 4 ? Use irregular search on primary dat
While simulation
Number of neighbors per octant : 2 Seed of random generator : 1
StatMod 8
APEX JASON
Use same path Always ?
For fine grid only Never ? Snap well to grid for fast performance
Transform … (3) Transforms … ? ? porosity_hist
Second Transform … ? ? aitm_hist
Variogram … (4) Variogram … ? ? porosity_var
Second Variograms … ? aitm_var ?
Cross variogram … ? ? porosity_aitm_cross_var
(5) Super codes …
Layer Super code
b b c c
3、Output ? Generate …
? Run now
Output type : Model ? Horizon
Primary output
Primary file name : ssm_bc_cosimu_first.mod
? Generate primary kriging error : ssm_bc_ cosimu _first_error.mod
并产生文件：（由于 Nr of realizations : 10 ,所以产生以下十个文件） ssm_bc_ cosimu _first_01.mod ssm_bc_ cosimu _first_02.mod
ssm_bc_ cosimu _first_03.mod …… Secondary Output
? Generate secondary data : ssm_bc_ cosimu _second.mod
? Generate primary kriging error : ssm_bc_ cosimu _second_error.mod 并产生以下十个文件 :
ssm_bc_ cosimu _second_01.mod ssm_bc_ cosimu _second _02.mod
ssm_bc_ cosimu _second _03.mod
……
4、File ? Save & Exit
四、 Sequential Indicator Simulation :(序贯指示模拟） 把井上的测井曲线数据，离散化为岩性数据， 作岩性指示变异图分析（ indication variogram） (一) 建立岩性类型：
Utilities ? Project management ? Project parameters …
StatMod 9
APEX JASON
Edit ? Type and units …
Add … lithology Data type :
? Lithologic data type : Edit litho-types …
Add … Lithotype : shale
Integerrepresentation : 注意：序号从0开始， 0 间隔为1。 Lthotype :
sand
Integerrepresentation : 1
……
(数据存入文件usertypelist.txt)
（二） 应用FunctionMod模块,通过测井曲线转换为岩性曲线（略）； （三）重新做 Earthmodel ,其中包括岩性曲线； （四）作岩性指示变异图分析：
1、 Iuput ： (作岩性指示变异图分析） (1) Time /Depth mode … ? Time (2) Trace gate … (3） Solid model … ? tdc1 (4) Layers … ? b,c,
（5） Data for variogram sampling and modeling ? Primary data ? well data
Data Well file : *.wll Model type :
lithology Done
( 出现StatMod variogram窗口)
Sample vario … （a) Function :
Indicator(Lithology) (b)Lithotype : shale (任选一种)
(c)
Build the sample Variogram … (d) Model vario …
Variogram type : Indicator(lithology) 1 st model type : lithology lithotype : shale Func Type : Gaussian 填写: sill , lag ( x, y, z) （Z:地层平均厚度, X*Y:地质体大小) Z是主要的!
StatMod 10
APEX JASON
(e) ? litho_var Save the first …
返回步骤(b)： Lithotype : sand 重复(c)-----(e) ? litho_var
……(注意：对每一种岩性分别作指示变异图分析,都存入同一个文件 litho_var）
2、File ? Save & Exit
(二) StatMod modeling …
1、Input ? Simulation mode … (打开两项） ? Property simulation Mode :
Gaussian Collocated cosimulation
? Run lithology simulation Mode :
Indicator simulation ? Indicator simulation with a trend
Done
其它类似。 2、Edit :
(1) Realization parameters … (2) Search parameyters … (3) Transforms …
Transform ? porosity_hist Second transform ? aitm_hist (4) Variograms… ?
Variogram ? porosity_var
Indicator variogram ? litho_var (5) Super codes …
(6) Lithology probabilities …(给各层岩性概率值）
3、Output ? Generate …
? Run now ……
Lithology Output Lithology file name :
Litho_sis_ind.mod
4、File ? Save & Exit
五、遮挡岩性指示模拟：(Lithology masks) 主要步骤与岩性模拟一样，不同的是：
遮挡某种岩性后，进行岩性指示模拟。
StatMod 11
APEX JASON
1、Statmod analysis :
Input ?Lithology masks …
? Enable liyhology masking
Lithology
Lithology
Select lithotypes for masking : (被照亮者为不被遮档的岩性，即所要的岩性）
Sand
Shale
Lithology data
Use : ? Model ? Horizon Model file :
Ind_aitm_ai.mod List …
其它同岩性指示模拟类似，部分岩性被遮档后，对未遮档的岩性部分 进行岩性的指示变异图分析。 2、Statmod modeling …
Input ? Lithology masks …
六、 Sequential Indicator Simulation with a trend (带趋势的序贯指示模拟)： （一）分别作砂岩与泥岩的指示变异图分析，
根据阻抗数据体，分别作砂岩与泥岩的分布趋势模型。
（二） StatMod modeling :
其它各项与指示模拟类似， 关键；
1、Input ? Simulation mode … （必须打开两项） ? Run property simulation Mode :
Gaussian Collocated cosimulation ? Run lithology simulation
Mode : Indicator simulation
Indicator simulation with a trend ?
Done
2、 Input ? Input Lithology probability trends ? 选 ai_to_sand_prob.mod ai_to_shale_prob.mod
3、Edit ：
（1） Realization parameters …
(2) Seach parameters … (3) Transform .. Transform Second transform (4) Variogram … ? Variogram Indicator variogram ? litho_test (5) Super codes …
StatMod 12 List …
APEX JASON
4、Output ? Generate … Lithology file name Simu_ssm_ind_trend.mod
5、File ? Save & Exit
七、 Stochastic inversion (随机反演):
（一） 需要如下分析数据
1、井上阻抗曲线的直方图分析， 2、 井上阻抗曲线的变异图分析。
3. 地震数据 4. 子波
(二) StatMod modeling:
1. Input :
(1) Simulation mode … (需打开两项） ? Run property simulation
Mode : Gaussian simulation ?
Gaussian cosimulation Gaussian collocated cosimulation Threshold indicator simulation
? Run lithology simulation
Mode : Indicator simulation
Indicator simulation with a trend
? Run Inversion ?
（2） Time /Depth mode :? Time (3) Solid mode : ? tdc1 (4) Layers … ? (5) Primary input data … ? Well log data … ? p_impedance (6) Variogram / transforms file … ? simu.var (7) Trace gate … (8) Data for inversion : (a) Seismic data … ? seis.mod (b) Wavelet : ? 选子波 2. Edit : (1) Realization parameters … (2) Search parameters … (3) Transforms … Transform (4) Variogram … Variogrem (5) Super codes (6) Inversion parameters … (a) Time gate … (b) Other inversion parameters … Inversion algorithm : Greedy ? Simulation annealing (模拟退火） Max. acceltance factor : 0.21 Min. acceltance factor : 0.01 StatMod 13 Max. seismic freq.(HZ) : 100 Nr. Iterations : 10 Convergence criterion : 0.1
APEX JASON
3. Output ? Generate … ?Run now Primary file name :
Inversion Output
? Generate synthetic data
? Generate residual data
? Generate synthetic/seismic corr.
simu_stochastic_inversion.mod stochastic_inversion_synthetic.mod stochastic_inversion_residuals.mod stochastic_inversion_sscor.mod StatMod 14