jason随机反演的中文手册.docx
《jason随机反演的中文手册.docx》由会员分享，可在线阅读，更多相关《jason随机反演的中文手册.docx（23页珍藏版）》请在冰豆网上搜索。
jason随机反演的中文手册
第五期用户通讯
随机模拟与随机反演
地震属性随机模拟、岩性指示随机模拟与随机反演：
1.需要输入的数据：
1、wells（logsandtops）
2、seismic.mod（=4ms…）
3、horizons
4、wavelets（=4ms…）
5、aitm.mod（=4ms…）
6、solidmodel（withTDC）
2.主要工作流程：
1、建立岩性类型与单位
2、利用density与gamma曲线建立岩性曲线（lithology）
3、建立岩性概率模型（shale_prob.mod,sand_prob.mod）
4、直方图与变异图的分析
并通过交会图确定两类数据的相关系数
5、随机模拟
6、计算多次孔隙度随机模拟结果的平均值
7、计算多次模拟的标准偏差
8、计算小层中的砂岩累计厚度
9、随机反演
（一）建立岩性类型与单位：
UtilitiesProjectmanagementProjectparameters
EditTypeandunits…
Add…
Lothologicdatatype
Editlithp-types…
sand/shale
Datatype:
Unit:
none
Add…
Lithologytype:
shale
0
Add…
Lithologytype:
sand
1
Ok（存入文件usertypelist.txt）
（二）利用density与gamma曲线建立岩性曲线（深度域、lithology）：
例如：
工区有六口井：
1.wll2.wll3.wll4.wll5.wll6.wll
经过分析：
当GR<85,density<2300kg/m^3时为sand
当GR>85,density.2300kg/m^3时为shale
1、建立一个函数（利用density与gamma曲线建立岩性曲线）：
FunctionModEditFunctionsNew
sand_shale_from_gamma_and_density
Functionname:
键入
List…
Sand/shale
Outputtype:
none
Outputunit:
Parameters:
Add…
Singlemodel,horizonorwell-log
Parameterskind选
gamma
Parametersname
List…
gamma
Parameterstype选
Edit…
gAPI
Parametersunit
Ok
Add…
Singlemodel,horizonorwell-log
Parameterskind选
density
Parametersname
List…
density
Parameterstype选
Edit…
Kg/m^3
Parametersunit
Ok
Usedomainrangedefinitions
Edit…
Add…
Range:
Gamma<85anddensity<2300
Case键入
Add…
Assignments:
sand_shale_from_gamma_and_density
Vriablename
1
Expressiontoassign键入
OkOk
Edit…
Default
Range:
Add
Assignments
sand_shalefrom_gamma_and_density
Vriablename
0
Expressiontoassign键入
OkOk
Oksaveas…sand_shale_from_gamma_and_density.fun
2、Time/DepthmodeDepth
InputVerticalgateofinputwheredataavaiable
Verticalgateofoutputwheredataavaiable
Horiontalgatewheredataavaiable
3、EditAssignments…
（1）FiletypeWell
Sampleinterval:
0.125m
Clear
Outputfilesettings:
sand/shale
Datatype
Add…
List…
Filename:
1.wllOk
Add…
2.wll3.wll4.wll5.wll6.wllOk
Add…
（六口井全进入窗口）
Assign
点亮1.wll
sand_shale_from_gamma_and_density
（2）Function:
选
Assignparameters…
点亮gamma（sfile）…
1.wll
WellOk
Assignparameters…
点亮density（sfile）…
1.wll
WellOk
Ok
Copy
点亮1.wll点亮2.wll3.wll…6.wll
Ok
Ok
4、EditExistingfiles…Appendanddon’toverwriteoverlap
5、OutputGenerater…Ok
6、ApplicationsDisplay&editsWelleditor（观察岩性曲线）
M
Color:
AttributesPlotattributes…（调颜色）
7、EarthModelbuilderwithTDC（建立时转换的岩性曲线）tdc
（三）建立岩性概率模型（shale_prob.mod,sand_prob.mod）
!
!
应针对某个小层（layer）
1、通过交会图了解砂泥岩的阻抗分布范围：
ApplicationsDisplayandEditors…Welledit...
选井，选曲线：
Porosity、impedance
CrossPlot
在交会图上了解砂泥岩在阻抗域的大致分布范围。
例如：
5e+67e+6
砂砂、泥泥
1100
2、通过交会图来建立岩性与阻抗转换的概率表：
Utilitiescreate/EditTabls…
（1）建立由阻抗转换成砂岩概率的表（ai_to_sand_pro.tab）：
X_nameY_name
AI
sand_prob
X_typeY_type
probability
impedance
5e+61
7e+60
Saveas…ai_to_sand_prob.tabOk
（2）建立由阻抗转换成泥岩概率的表（ai_to_shale_prob.tab）:
X_nameY_name
shale_prob
AI
X_typeY_type
probability
impedance
5e+60
7e+61
Saveas…ai_to_shale_prob.tabOk
3、生成岩性概率模型，作为随机指示模拟的趋势约束：
在FunctionMod模块中完成
（1）InputTime/Depthmode…Time
VerticalgateofinputSolidmodel
Selectsolidb,c,d
VerticalgateofoutputSolidmodeltdc
HorizontalgateTracegate…
Settracegate…
（2）EditAssignments…
Clear
Seismic/Property
Filetype:
Sampleinterval:
0.001s
Datatype:
probability
sand_prob.mod
Filename:
键入
Add…
抹亮sand_prob.mod
Assign
b,c,d
Assign
抹亮cai_to_sand_prob.tabOk
Assignparameters..
抹亮x（sfile…）
Seismic/Propertyaitm_ai.mod
Copy
Ok
抹亮cb,dOk
OkOk
（3）OutputGenerate…Ok
（4）同理产生：
shale_prob.mod
（5）剖面显示:
砂岩概率剖面（sand_pro.mod）、
泥岩概率剖面（shale_prob.mod），
作为后边随机指示模拟的趋势约束。
（四）直方图与变差图的分析（针对一个小层layer）:
例如：
对c层作分析：
（存入文件：
analysis.var）
=0.0005s
如果作：
collocatedco-simulationb
indicatesimulationwithtrend
需要作如下直方图与变异图的分析：
1、井数据c_sand_porosity_hist（15~35）
井曲线上c层段、砂岩的孔隙度直方图，数据范围限制在c层段
砂岩孔隙度的范围内（例如：
15~35）
2、井数据c_shale_porosity_hist（1~20）
井曲线上c层段、泥岩的孔隙度直方图，数据范围限制在c层段
泥岩孔隙度的范围内（例如：
1~20）
3、数据体c_sand_aitm_hist（1e+06~8e+06）
阻抗数据体上c层段、砂岩的孔隙度直方图，数据范围限制在c
层段砂岩孔隙度的范围内（例如：
1e+06~8e+06）
4、数据体c_shale_aitm_hist（4e+06~1e+07）
阻抗数据体上c层段、泥岩的孔隙度直方图，数据范围限制在c
层段泥岩孔隙度的范围内（例如：
4e+06~1e+07）
5、井数据c_sand_porosity_var
井曲线上c层段、砂岩的孔隙度变差图，数据范围限制在c层段
砂岩孔隙度的范围内（例如：
15~35）
6、井数据c_shale_porosity_var
井曲线上c层段、泥岩的孔隙度变差图，数据范围限制在c层段
泥岩孔隙度的范围内（例如：
1~20）
7、井数据litho_varsand
shale
井曲线上c层段的岩性指示变差图，砂岩与泥岩要分别做，但要
存在同一个litho_var中
8、两类数据的相关系数
通过交会图确定两类数据的相关系数。
利用井曲线做porosity与
impedance的交会图，用该图来确定两类数据的相关系数。
ApplicationsDisplay&EditorsWelledit..
InputWell…选井
EditDisplay&Edit选曲线（porosity、impedance）
CrossPlot
porosity
sand
shale
impedance
负相关正相关
相关性不好相关性好
（也可用井上的孔隙度与阻抗数据做交会图）
（五）随机模拟：
StatModStatModModeling
Input1.Simulationmode…
Runpropertysimulation
Mode:
Gaussiancollocatedco-simulation
Runlithologysimulation
Mode:
Indicatorsimulationwithtrend
2.Time/Depthmode…time
Solidmodel…tdc
3.Layerc0.0005s
4.PrimaryinoutdataWelllogdatatdc*.wll
impedance
5.SecondaryinputdataModeldataaitm_ai.mod
6.Lithologyinputdatawelllogdatatdc*.wll
sand/shale
7.Variograms/transformsfiles…analysis.var
8.Lithologyprobabilitytrends…
Trendtype:
Model
Selecttrendmodels:
shaleshale_prob.mod
sandsand_prob.mod
9.Tracegate
Edit
1、Realiationparameters:
10
Outputsampleinterval（s）:
0.0005
Performkrigingalso
Useordinarykriging
Usecokrigingsecondnon-biascondition
（克里金方法可选）
2、Searchparameters:
1000（视井距而定）
3、Transforms…
Transforms…
cshalec_shale_porosity_hist
csandc_sand_porosity_hist
SecondTransforms…
cshalec_shale_aitm_hist
csandc_sand_aitm_hist
4、Variograms…
Variograms…
cshalec_shale_porosity_var
csandc_sand_porosity_var
Indicatorvariograms…
clitho_var（包括shale、sand）
5、Supercodes
LayernameSupercode
cc
6、Correlationcoefficients…
cshale-0.7
csand-0.7
OutputGenerate…
Primaryfilename:
Primaryfilename:
ssm_first_c_poro.mod
Generateprimarykigingerror:
ssm_first_c_poro_error.mod
Secondaryfilename:
Generatesecondarydata:
ssm_second.mod
Generatesecondarykrigingerror:
ssm_second_error.mod
Lithologyfilename
Lithologyfilename:
ssm_sis.mod
Ok
（六）计算多次孔隙度随机模拟结果的平均值（mean_poro.mod）:
FunctionMod
InputTime/Depthmode…Time
VerticalgateofinputSolid…
VerticalgateofoutputSolid…
HoriontalgateTraceegate…
EditAssignments…
Filetype:
Seismic/Property
Sampleinterval:
0.0005s
clear
Outputfilesettings:
Datatype:
porosity
Filename:
mean_poro.modAdd
点亮mean_poro.modAssign
Mean_predefined_mfile
Function:
点亮v（mfile）…Assignmentsparameters…
Seismic/Property
File:
（选多次孔隙度模拟的结果）
ssm_c_first_poro_01.mod
ssm_c_first_poro_02.mod
……
ssm_c_first_poro_10.mod
AddOkOk
Ok
OutputGenerateOk
（七）计算多次模拟的标准偏差：
FunctionMod
InputTime/Depthmode…Time
VerticalgateofinputSolid…
VerticalgateofoutputSolid…
HoriontalgateTracegate…
EditAssignments…
Filetype:
Seismic/Property
Sampleinterval:
0.0005s
clear
Outputfilesettings:
none
Datatype:
Filename:
std_poro.mod
Add
点亮std_poro.modAssign
std_predefined_mfile
Function:
Ok
点亮v（mfile）…Assignmentsparameters…
porosity
Type:
Seismic/Property
File:
（选多次孔隙度模拟的结果）
ssm_c_first_poro_01.mod
ssm_c_first_poro_02.mod
……
ssm_c_first_poro_10.mod
AddOkOk
Ok
OutputGenerateOk
（八）计算小层中的砂岩累计厚度：
FunctionMod
InputTime/Depthmode…Time
VerticalgateofinputSolidmodel…
Slectsolidmodel
tdc
c（layer）
VerticalgateofoutputSolidmodel
HoriontalgateTracegate
EditAssignments…
Filetype:
Horion
clear
Outputfilesettings:
Datatype:
none
Filenamesand_thick_time_samp.horAdd
Horion:
sand_c
Add...点亮sand-thick_time_samp.hor/sand_c
Assign
Count_predefined_trace
Function:
点亮v（mfile）…Assignmentsparameters…
Ttype:
none
Seismic/Property
File:
（选岩性模拟的结果）
ssm_sis_c.mod
Ok
点亮I（value）Assignparameters…
Modeltype:
sand/shale
sand
Value
OkOkOk
OutputGenerateOk
（九）随机反演：
需要准备的分析数据：
1、井曲线impedance的直方图（c_impedance_hist）
2、井曲线impedance的变差图（c_impedance_var）
3、子波
4、地震数据
StatModStatModmodeling
InputSimulationmode…
InputSimulationmode…
Runpropertysimulation
Mode:
Gaussiansimulation
Runinversion
Time/Depthmode…time
Solidmodel…tdc
Layerc0.0005s
PrimaryinoutdataWelllogdatatdc*.wll
impedance
Variograms/transformsfiles…analysis.var
Tracegate
Dataforinversionseismicdataseismic.mod
Wavelets..选子波
Edit
1、Realiationparameters:
10
Outputsampleinterval（s）:
0.0005
Performkrigingalso
Useordinarykriging
Usecokrigingsecondnon-biascondition
（克里金方法可选）
2、Searchparameters:
1000（视井距而定）
3、Transforms…
Transforms…
cc_impedance_hist
4、Variograms…
Variograms…
cc_impedance_var
5、Supercodes
LayernameTimegate..
Supercode
cc
升级会员 