沉积体系非均质性对致密油气储集层控制及其定量评价方法.pdf

非 常规油气 资源 矿物岩石地球化学通报 B u l l e t in o f Min e r a l o g y，P e t r o l o g y a n d Ge o ch e mis t r y Vo 1 3 4 No 1，J a n ，2 01 5 沉积体 系非均质性对致密油气储 集层 控 制及 其定量评价方法 黄秀 2 ， 刘可禹1 , 2 , 3 ， 吴松涛 4 ， 桂丽黎 2 1 中国石油天然气股份有 限公 司 勘探开发研究 院 ， 北京 1 0 0 0 8 3 ； 2 中国石油天然气股份有限公司 盆地构 造与油气成藏重点实验室 ， 北京 1 0 0 0 8 3 ； 3 C S I R O， P O B o x 1 1 3 0 ， B e n t l e y ， wA 6 1 0 2 ， A u s t r a l i a ； 4 中国石油天然气股份有 限公 司 油气储层重点实验室 ， 北京 1 0 0 0 8 3 摘 要 ： 以鄂 尔多斯盆地 延长组长 8 一 长 6致 密储 集层为例 ， 采用 地质 统计数 学( L e v y S i m) 和 沉积地层 正演模 拟对 沉积非均 质 性进行研 究。结果表 明， 地 质统计学 L e v y函数 能够有效地表 征非均质性参数 ， 不仅 可以表征油藏尺度( 岩心 、 露头和测井 曲线 尺度 ) ， 而且也 可以有效刻画岩石微观尺度 的非均质性。而三维正演地层模拟方法 S e d s im则能够有效地模拟浅水三角 洲砂 体 的空间展 布。沉积地层正演模拟 与非传统地 质统计学相结合 ， 能够 更真实地定量 刻画不 同尺度 的沉积 非均质性 ， 该方 法不仅 有效地解 决 了不 同尺度沉积储集层 、 烃源岩及盖层非均质 性定 量表 征的难题 ， 而且为非 常规 油气成 藏研究提供 了新思路 与新 的定 量评价 方法。 关 键 词 ： 沉积 非均质性； 浅水 三角洲； 鄂 尔多斯盆地 ； 地质统计 学， 三维定量正演模拟 中图分类号： P 6 1 8 1 3 0 2 文章编号 ： 1 0 0 7 2 8 0 2 ( 2 0 1 5 ) O 1 0 o 6 1 0 9 d o i： 1 0 3 9 6 9 j i s s n 1 0 0 7 2 8 0 2 2 0 1 5 0 1 0 0 7 Qu a n t it a t iv e An a l y s is o f S e d ime n t a r y He t e r 0 g e n e it ie s HUANG Xi u 。 一 ， LI U Ke y u 。 ， ， ” ，W U S o ng t a o ，GUI L i 1 i ， 1 R e s e a r ch I n s t i t u t e o f P e t r o l e u m E x p l o r a t i o n a n d D e v e l o p m e n t ， P e t r o C h i n a ，B e i j in g 1 0 0 0 8 3 ，C h in a ；2 C N P C K e y L a b o r a t o r y f o r Ba s in S t r u ct u r e a n d H y d r o ca r b o n Accu mu l a t io n ，Be lti n g 1 0 0 0 8 3 ，C h i n a；3 C S I R O E n e r g y F l a g s h i p，P 0 B o x 1 1 3 0，B e n t l e y，WA 6 1 0 2 ， A u s t r a l i a ； 4 C N P C K e y L a b o r a t o r y f o r O i l a n d G a s R e s e r v o i r ， B e n g 1 0 0 0 8 3 ，C h i n a Ab s t r a ct ：S e d ime n t a r y h e t e r o g e n e it ie s，wh ich e x is t in n a t u r e a t a l l s ca l e s ，c a n b e q u a n t i t a t iv e l y c h a r a ct e r i z e d b y u s in g g e o s t a t is t ic s a n d n u me r i c a l s imu l a t io n s I n t h is s t u d y，a t h r e e d i me n s io n a l s t r a t i g r a p h ic f o r wa r d mo d e l l in g p r o g r a m，S e d s - i m ，wa s u s e d t o s imu l a t e t h e d e v e l o pme n t a n d e v o l u t io n o f t h e s h a l l o w- wa t e r d e l t a s o f t h e C h a n g 8 - 6 me mb e r s o f t h e Ya h - ch a n g F o r ma t io n i n t h e Od o r s Ba s in Th e r e s u l t s s h o we d t h a t b a t h g e o s t a t is t ics a n d s t r a t i g r a p h i c f o r wa r d s imu l a t io n s ca n ch a r a ct e r iz e t h e s e d ime n t a r y h e t e r o g e n e i t ie s we l 1 F o r i n s t a n ce，t h e g e o s t a t is t ica l a p p r o a ch，s u ch a s L e v y S im ，ca n e ffe c t i v e l y c h a r a c t e r iz e t h e h e t e r o g e n e it y p a r a me t e r s o f s e d ime n t a ry f a cie s wit h o u t t e d io u s s u b d i v i s io n，wh il e n u me r ica l s imu l a t io n s ，s u ch a s s t r a t ig r a p h ic f o rwa r d mo d e l l in g t o o l s( S e d s im)， ca n b e u s e d t o r a p id l y t e s t mu l t ip l e w o r k i n g h y p o t h e s e s o f t h e r e l a t i o ns h ip s b e t we e n t h e h y d r o d y na mi c p r o ce s s e s a n d d e p o s it i o n a l f a cie s a t v a r io u s s ca l e s T h e co mb i n a t io n o f t h e n o n Ga u s s i a n g e o s t a t is t ica l me t h o d a n d s t r a t ig r a p h ic f o r wa r d mo d e l in g a p p r o a c h e s n o t o n l y ca n e ff e ct i v e l y ch a r a ct e r iz e s e d ime n t a ry h e t e r o g e n e it ie s a t d if f e r e n t s ca l e s b u t a l s o ca n p r o v id e a n e w wa y f o r t h e q u a n t i t a t iv e e v a l u a t io n o f u n co n v e n t i o na l r e s er v oir s Ke y wo r d s：s e d i me n t a ry h e t e r o g e n e it y；s h a l l o w wa t e r d e l t a s ；Or d o s B a s in；g e o s t a t i s t i c s ；s t r a t i g r a p h ic f o rw a r d mo d e l 收稿 日期 ： 2 0 1 4 - 0 6 3 0收到 ， 2 0 1 4 0 7 - 1 5改回 基金项 目： 国家重点基础研究计划 ( 2 0 1 4 C B 2 3 9 0 0 4 ) ； 中国石油天然气股份有 限公 司科技管理部项 目( 2 0 1 1 A 一 0 2 0 7 ) ；中国石油长庆油 田石 油分 公司示范工程项 目( 2 0 1 1 Z X0 5 0 4 4 ) 第一作 者简介 ： 黄秀 ( 1 9 8 1 一 ) ， 女 ， 博士研究生 ， 研究方 向： 盆地分析数值模拟 E m a il ：h u a n g x i u 1 9 8 3 p e t r o ch i n a C O B cn 通讯作者简介 ： 刘可禹 ( 1 9 6 3 一 ) ， 男 ， 教授 ， 研究方向 ： 含油气系统分析 E m a i l ：k e y u l i u p e t r o e h i n a e o m e n 6 2 储集 层非 均质 性 主要 与 不 同性 质 ( 岩性 ， 粒度 ， 分 选性 ) 沉 积 物 及 不 同 尺 度 界 面 的 空 间 叠 置 有 关 ( Ha l d o r s e n ，1 9 8 6 ) 。沉 积后 过 程 如构 造 应 力 、 压 实 和流体演化等， 均作用于初始的沉积格架和非均质 性 上 ， 最 终产 生 了不 同尺 度 的非 均 质 性 储集 层 。对 于基质与裂缝 ， 沉积相结合成岩和应力历史为控制 储集层流体演 化的基本要素 ( N a is h a n d L iu ， 2 0 0 0 ； Ky s e r a n d Hia t t ， 2 0 0 3； Eic h h u b l e t a 1 ， 2 00 4； Be i t l e r e t a 1 ， 2 0 0 5；Ab o u e s s a a n d Mo r a d，2 0 0 9；Mo r a d e t a 1 ， 2 0 1 0 ； C a r a c c io l o e t a 1 ， 2 0 1 3 ) 。不 同沉 积 环境形 成 的 砂 体 具有不 同的特 征 ， 早 期成 岩 作 用 也 受 到沉 积 环 境 的影响 ， 从 而影 响进 一 步 成 岩作 用 的类 型 、 强 度 ， 并对 砂岩 的孔 隙演 化 起 一 定 的 控 制 作 用 。尽 管 成 岩作用对储集层砂体进行强烈 的改造 ， 但是成岩作 用基于沉积作用 的基础上对原始孑 L 隙度进行改造 ， 对于大多数储集层而言 ， 原始沉积环境依然是 控制 储集 层 发 育 的 主 要 因 素 ( B l o c h a n d M c G o w e n ， 1 9 9 4 ) 。因此 ， 了解 沉 积 演 化 动 态 过 程 是 更 好 预 测 储 集层 位置 、 储集层 几何 结构 和渗 透率 隔层 的关 键 。 目前 ， 对 处理 非 均 质性 表 征 和 尺 度升 级 或 粗 化 的传 统 方 法 主 要 有 指 数 律 平 均 法 ( Wa r r e n a n d P r i c e ， 1 9 6 0 ； J o u r n e l a n d H u i j b r e g t s ， 1 9 7 8 ) 、 重 正 规 化 技 术 ( K in g ，1 9 8 9 ；1 9 9 6 ) 、 解 流 动 压 力 方 程 法 ( C h r is t ie a n d C l i f f o r d ， 1 9 9 7 ) 、 矢 量 法 ( P i c k u p e t a 1 ， 1 9 9 4 ) 和拟 函数 法 ( J o n e s e t a 1 ， 1 9 9 3 ) 等 5类 。这 些 传 统方 法 主 要 集 中在 对 网格 单 元 里 各 物性 参 数 的 粗化方法方 面， 属于纯物理 数学方法 ， 其共 同的特 点是倾 向于抹平一些在空 间分 布的渗透率极值条 带 ， 如 ， 泥 岩 夹 层 和 裂 缝 。然 而 对 非 常 规 油 气 藏 的 生产实践和模拟实验都表明， 这些渗透率极值条带 对非 常规 油气 资 源 评 价 影 响极 大 。随 着 非 常规 油 气勘 探 和开发 的进 展 ， 数 字 岩 心 技术 得 以广 泛 应 用 于刻 画致 密 储 集 层 的微 纳米 结 构 和 油 气 储 存 的 有 效储集空间 ， 但如何解决微纳米尺度储集层物性的 代表 性及 在 油气 藏 尺 度 的 贡 献 是 一 个 具 有 挑 战 性 的科 学 难题 。 因此 ， 为了 能有 效 地 解决 沉 积 非 均 质性 定 量 表 征 和 尺度升 级 ， 本 次 研究 主要 以鄂 尔 多 斯 盆 地长 8 一 长 6 为 例 ， 运 用 非 传 统 地 质 统 计 学 L e v y 函 数 ( P a in t e r a n d P a t e r s o n ， 1 9 9 4 ) 与 沉 积 地 层 正 演 模 拟 ( G r i ff it h s e t a 1 ， 2 0 0 1 ) 相结合 ， 真实地定量刻画沉积 非均 质性 ， 从 而 解决 不 同 尺度 ( 岩心 、 测 井 和地 震 尺 度 ) 沉 积储 集 层 、 烃 源 岩 及 盖 层 非 均 质 性 定 量 表 征 的难 题 。L e v y统 计 学 方 法 可 以用 于 表 征 储 集 层 微 纳米尺度到油藏尺度的非均质性 ， 而 S e d s im可以模 黄秀等 ： 沉积体系非均质性对致密油气储集层控制及其定量评价方法 拟厘 米 尺 度 到 数 百 公 里 ( 盆 地 ) 尺 度 的 非 均 质 性 。 本 文 旨在 为非 常 规 油 气 成藏 研 究 提 供 了 新 思 路 与 新 的定量 评价 方法 。 1 沉积体 系非均质性定量表征方法 沉 积非 均 质 性 存 在 于 自然 界 沉 积 体 系 的 任 何 尺度( Mia l l ， 1 9 9 0 ； L iu e t a 1 ， 2 0 0 2 ) ， 如 自小到大从纳 米 、 孔喉 、 岩 心 、 露 头 、 地 震 到 盆 地 尺度 ， 非 均 质 性 均 有不同的特征( 图 1 ) 。砂体的发育主要受到沉积环 境 的控 制 ， 不 同沉 积 环 境形 成 的砂 体 其特 征 各 不 相 同 ， 不 同沉 积过程 决 定 了沉积 体 系 的非 均质 性 和 自 相似性 程 度 及 可 粗化 区 间 ( L i u e t a 1 ， 2 0 0 2 ) 。这 些 不 同沉 积 体 系 的非 均 质 性 可 以通 过 地 质 统 计 学 和 沉积模 拟 来定量 表 征 。因此 ， 本 次 研 究 主要 是 利用 L e v y S im地 质统计 学 和 S e d s im三 维地 层 正演 模 拟方 法来对沉积体系非均质性进行定量表征。 1 1 L e v y函数 为 了更好 地 解 决 非 均 质 性 粗 化 的 问题 ， P a in t e r 提 出测井 数 据 和 储 集 层 非 均 质 性 具 有 列 维 ( L e v y ) 分 布 函数 的 统 计 特 征 ( P a in t e r a n d P a t e r s o n ， 1 9 9 4 ； P a in t e r ， 1 9 9 5 ； 1 9 9 6 ； P a in t e r e t a 1 ， 1 9 9 7 ) 。L e v y函数 是分数列维稳定运动 的形式， L e v y稳定概率密度 函 数 见S a mo r o d n i t s k y 和T a q q u( 1 9 9 4)和 L i u 等 ( 2 0 0 2 ) 。该方法与传统方法不 同， 其考虑了沉积环 境 对 非均质 性 的控 制 作用 ， 认 为储 集 层 非 均 质性 并 不是完全服从高斯分布 ， 高斯分布只是 L e v y分布的 一 种特殊分布特征。与高斯稳定分布相 比较 ， L e v y 分布函数更能够捕捉更多的地质信息 ， 从而能够更 真 实地反 映沉 积非 均 质性 。如 图 2所 示 ， 在 对 砂 岩 显微 x 光 C T照片 分 析结 果 中 ， L e v y分 布 函数 比高 斯稳 定分 布更 能反 映原 始 砂 岩显 微 照片 中的 信 息 。 而 且 ， L e v y函数使 用 L e v y指数 可 以有 效地 反 映沉 积非均质性 的程度 ， 值越 小， 其非 均质性程 度越 强 ， L e v y函数 中的间隔参数 日可以作为一个尺度参 数反映沉积环境的 自相似尺度区间， 粗化指数为最 大尺度值与最小尺度值 的比值 ， 与沉积 环境相关 ， 不 同 的 沉 积 环 境 具 有 不 同 的 粗 化 指 数 值 ( P a in t e r a n d Pa t e r s o n，1 99 4；Pa in t e r ， 1 9 9 5；1 99 6；Liu a nd Mo l z ， 1 9 9 7； Pa in t e r e t a 1 ， 1 9 9 7； Liu e t a 1 ， 2 0 0 2； Mo l z e t a 1 ， 2 0 0 4 ) ， 从而避免了人为地解释与描述 。 1 2 S e d s im 原 理 S e d s i m 能 够 模 拟 不 同尺 度 的 沉 积 非 均 质 性 。 S e d s i m主要以流体动力学方程( N a v ie r S t o k e s e q u a t i o n s ) 为核心。完整 的 N a v ie r - S t o k e s 方程 描述的是 流体在三维空间的运动 ， 然而 由于计算机速度的限 6 8 中的 自相似尺度区间约为 0 2 4 1 9 1 In， 粗化指数 为 8 ； 湖相沉积环境中的自相似尺度区间约为 0 2 4 6 0 i n ， 粗 化指 数为 2 5 。 ( 3 ) 地 质 统计学 L e v y函数 能够 有效地 表 征非均 质性参数 ， 不仅可以表征油藏尺度 ( 岩心 、 露头和测 井曲线尺度) 而且也可以有效刻 画岩石微观尺度的 非 均质性 。而 三维 正 演 地 层 模 拟 方 法 S e d s im 则 能 够 有效 地 模 拟 浅 水 三 角 洲 砂 体 的 空 间 展 布 。 由 于 在 时 间插 值 和 网 格 大 小 可 根 据 研 究 的 需 要 进 行 设 置， 模拟得到的垂向分辨率 比实际地震分辨率要高 得多 ， 横 向砂体 对 比比测 井 间对 比具有更 高 的准 确 度 。 致谢 ： 感谢 袁 选 俊 教 授 、 中 国 石 油 长 庆 油 田 石 油分公 司姚 泾利教授 以及 中国地质大学( 北京 ) 邓 宏文教授 对本 文的 悉 心指 导 ； 及 澳 大利 亚新 南威 尔 士 大学于 沐汉提供 图 2数据 。 参考文献 ( R e f e r e n ce s ) ： Ab o u e s s a A，Mo r a d S 2 0 0 9 An in t e g r a t e d s t u d y o f d ia g e n e s is a n d d e p o s it io n a l f a e ie s in t id a l s a n d s t o n e s ：Ha wa z f o r ma t io n( Mid d l e Or d o v i cl a n ) ，Mu r z u q B a s i n ，L i b y a J o u r n a l o f P e t r o l e u m G e o l o g y ，3 2 ( 1 ) ： 3 9 6 5 B e it l e r B，Pa r r yW T，Ch a n M A，Be it l e r B，P a r r y W T，Ch a n M A 2 0 0 5 F in g e r pr in t s o f flu id fl o w ：Ch e mica l d ia g e n e t ic h is t o r y o f t h e J u r a s s ic N a v a j o S a n d s t o n e ，s o u t h e r n U t a h ，U S A J o u r n a l o f S e d i m e n t a r y R e s e a r ch ，7 5 ( 4 ) ： 5 4 7 5 6 1 Bh a t t a ch a rya J P2 0 1 1 Pr a ct ic a l p r o b l e ms in the a p p l ica t io n o f t h e s e q u ci a e e s t r a t ig r a p h ie me t h o d a n d k e y s u r f a ce s ： I n t e g r a t ing o b s e r v a t i o n s f r o m a n ci e n t fl u v i a l d e h a i c w e d g e s w i t h Q u a t e rna r y a n d mo d e l l in g s t u d ie s S e d ime n t o l o g y，6 3 2 6 3 3 Bl o c h S，Mc Go we n J H1 9 9 4I n flu e n ce o f d e p o s it io n a l e n v ir o n me n t o n r e s e r v o ir q u Mit y p r e d ic t io nI n：W il s o n M D，e d R e s e rvo ir q u a l it y a s s e me n t a n d p r e d ict io n in e l a s t ic r o ck s S E PM s h o t Co u r s e，3 0： 41 57 Ca r a c cio l o L，Ar r ib a s J，I n g e r s o l l R V，T o r t o s a A 2 01 3 T h e d ia g e n e t ic d e s t r u c t io n o f p o r o s it y in p l u t o n ie l a s t ic p e t r o f a cie s ：T h e mio ce n e d il ig e n c ia a n d e o ce n e ma n io b r a f o r ma t io n s ， Or o e o p ia Mo u n t a in s ， s o u t h e r n Ca l if o r n ia ， US A Ge o l o g ic a l S o cie t y ， Lo n d o n， S p e cia l P u b l ic a t io n s ，3 8 6 Ch r is t ie M A，Cl if f o r d P J 1 9 9 7 A f a s t p r o ce d u r e f o r u p s ca l in g in C O B- p o s it io n a l s imu l a t io n S P E，3 7 9 8 6 Eich h u b l P，T a y l o r W L，Po l l a r d D D ，Ay d in A 2 0 0 4 P a l a e o fl u id flo w a n d d e f o rm a t io n in t h e Az t e c S a n d s t o n e a t t h e Va l l e y o f Fir e，Ne v a d a Ev id e n ce f o r t h e co u p l in g o f h y d r o g e o l o g ic ， d ia g e n e t ic ， a n d t e ct o n ic p r o ce s s e s G e o l o g ica l S o cie t y o f A m e rica B u l l e t in ，1 1 6 ( 9 1 O ) ：1 1 2 0 1 1 3 6 Gr iffit h s C M ，Dy t C，Pa r a s ch iv o iu E，L iu K 2 0 01 S e d s im in h y d r o ca r b o n e x p l o r a t io n I n ： Me r r ia m D， D a v is J C，( e d s ) Ge o l o g ic Mo d e l 一 黄秀等：沉积体系非均质性对致密油气储集层控制及其定量评 价方法 l in g a n d S imu l a t io n Ne w Yo r k：Kl u we r Aca d e mic Ha l d o r s e n H H1 9 8 6 S imu l a t o r p a r a me t e r a s s i g n me n t a n d t h e p r o b l e m o f s ca l e in r e s e rvo ir e ng in e e ri n g I n：L a k e L W ，C a r r o l l H B，( e d s ) Re s e r v o ir Ch a r a ct e r iz a t io n Ac a d e mic Pr e s s ，Or l a nd o2 9 3 3 4 0 Hu a n g X ，Dy t C，Gr if f it h s C，S a l l e s T 2 0 1 2 Nu me r ic a l f o r wa r d mo d e l l in g o f fl u x o t u r b id it e fl ume e x p e ri me n t s us in g S e d s imMa r in e a n d P e t r o l e u m G e o l o g y ， 3 5 ( 1 )：1 9 0 2 0 0 J o n e s A， D o y l e J ， J a co b s e n T K j f n s v i k D 1 9 9 3 Wh ich S u b s e i s mi c h e r e r o g e n e it ie s in flu e n ce wa t e r f l o o d p e rfo r ma n c e ?A c a s e s t u d y o f a l o w n e t t o - g r o s s fl u v ia l r e s e rvo ir P a p e r p r e s e n t e d a t t h e 7 t h Eu r o p e a n I OR S y mp o s iu m in MOS COW ，Ru s s ia ，2 72 9 J o u me l A G，H u i j b r e g t s C J 1 9 7 8 Mi n i n g G e o s t a t i s t i e s ：L o n d o n ，A ca d e mic Pr e s s，6 00 Kin g P R1 9 8 9 T h e u s e o f r e n o r ma l is a tio n for ca l cu l a t in g e f f e c t iv e p e r - me a b il it y T r a n s p o r t Po r o u s Me d ia，4：3 7 5 8 Kin g P R1 9 9 6Up s ca l in g p e r me a b il ity：Er r o r a n a l y s is f o r r e n o r ma l iz a t io n Tr a n s p o r t in Po r o u s Me d ia ，2 3： 3 3 7 3 5 4 Ky s e r K，Hia t t E E F l u id s in s e d ime n t a ry b a s in s ：An in t r o d u c t io n J o u r - h a l o f G e o ch e m ica l E x p l o r a t i o n， 2 0 0 3， 8 0 ( 2 ) ：1 3 9 1 4 9 L iu H H Mo l z F J 1 9 9 7 C o mme n t o n “ Ev id e n ce for n o n Ga u s s ia n s e a 1 in g b e h a v io r in h e t e r o g e n e o u s s e d ime n t a ry f o r ma t io n s ” b y S co t t W a t e r Re s o u r Re s ，3 3： 9 0 79 0 8 L i u K Y， P a t e r s o n L，Wo n g P，9 i D A A 2 0 0 2 s e d i me n t 0 l o g i c a l a p p r o a c b t o u p s e a l in$ T r a n s p o r t in Po r o u s，Me d ia 4 6：2 8 5 3 1 0 Mia l l A D1 9 9 0 P r in cip l e s o f s e d ime n t a r y b a s in a n a l y s is ，2 n d e d，Ne w Yo r k：S p ri n g e r Ve r l a g，6 6 8 Mo l z F J ， R a j a r a m H， L u S 2 0 0 4 S t o ch a s t ic f r a ct a l b a s e d m o d e l s o f h e t - e r o g e n e it y in s u b s u r f a ce h y d r o l o g y ： o ri g in s ， a p p l ic a t io n s ， l imit a t io n s，a n d f u t u r e r e s e a r c h q u e s t io n s Re v ie w s o f Ge o p h y s ic s ， 4 2RG1 0 0 2 d o i ：1 0 1 O 2 9 2 0 0 3 RGO 0 0l 2 6 Mo r a d S，A1 一 Ra ma d a n K ，Ke t z e r J M ，De Ro s L F 2 01 0 T h e imp a ct o f d ia g e n e s is o n t h e h e t e r o g e n e it y o f s a n d s t o n e r e s e r v o ir s ：A r e v ie w o f t h e r o l e o f d e p o s it io n a l f a cie s a n d s e q u e n ce s t r a t ig r a p h y AAPG b u l l e t in ， 9 4 ( 8 ) ：1 2 6 7 1 3 0 9 Na is h T， Liu K 2 0 0 0 Re v ie w o f “Co mp u t e ri z e d Mo d e l l in g o f S e d ime n t a ry S y s t e m”b y H a r f f e t a 1 ， ( 1 9 9 9 ) S e d ime n t a ry G e o l o g y ， 1 3 0：1 4 5 1 4 7 P a in t e r S，Pa t e r s o n L，Bo u h P1 9 9 7I mp r o v e d t e ch n iq u e f o r s t o cha s t ic in t e rpr e t a t io n o f r e s e rvo ir p r o p e rtie s ，S PE Pa p e r 3 0 5 9 9，S PE J 2 ( 1 ) ， 4 8 5 7 P a in t e r S， Pa t e r s o n L 1 9 9 4F r a c t io n a l L e v y mo t io n a s a mo d e l for s p a t ia l v a r ia b il it y in s e d ime n t a ry r o c k， Ge o p h y s Re s Le t t 2 1： 2 8 5 7 2 8 6 0 Pa in t e r S 1 9 9 5 Ra n d o m f r a ct a l mo d e l s o f h e t e r o g e n e it y： Th e L e v y s t a b l e a p p r o a c h，Ma t h e ma t ic a l Ge o l o g y，2 7：81 3 8 3 0 Pa in t e r S 1 9 9 6 S t o ch a s t ic in t e rp r e t a t i o n o f a q u if e r p r o p e rtie s u s in g f r a c t io n a l L e v y mo t io n，W a t e r R e s o u r Re s ，3 2：1 3 2 3 1 3 3 2 Pick u p G E， Rin g r o s e P S，Fo r e s t e r M M ，J e n s e n J L， S o r b ic K S 1 9 9 4Pe r me a b i l it y t e n s o r s for s e d ime n t a r y s t r u c t u r e s ，Ma t h e ma t ica l Ge o l o g y，2 6：2 2 7 2 5 0 S a l l e s T，Ma r e h 6 s E，Dy t C，Ha n q u ie z V，Mu l d e r T 2 01 0 S imu l a t io n o f t h e in t e r a ct io n s b e t we e n g r a v it y p r o ce s s e s a n d co n t o u r c u r r e n t s o n t h e Al g a r v e Ma r g in，S o u t h Po rtu g a l u s in g t h e s t r a t ig r a p h ic f o rw a r d 矿物岩石地球化学通报2 0 1 5 ， 3 4 ( 1 ) mo d e l S e d s imS e d ime n t a r y Ge o l o g y，2 2 9：9 51 0 9 S a mo r o d n i t s k y G，Ta q q u M S 1 9 9 4 S t a b l e n o n Ga u s s ia n r a n d o m p r o e e s s e s ： S t o c h a s t ic mo d e l s wit h in fi n it e v a ri a n ce Ne w Yo r k： Ch a p ma n a n d Ha l l ，6 3 2 Te t z l a f f D M Ha r b a ug h J W 1 9 8 9 S imu l a t in g e l a s t ic s e d ime n t a t io n c o mp u t e r me t h o d s i n t h e g e o s cie n ce s Ne w Yo r k：Va n No s t r a n d Re - in h o l d W a r r e n J E，P ri ce H S1 9 6 0 Fl o w in h e t e r o g e n e o u s p o r o u s me d ia Pr e s - e n t e d a t 3 5 t h An n u a l Fa l l Me e t in g o f S P E， De n v e r ， 1 0 Zo u C N，Z h a n gX Y，L u o P，W a n g L，L u o Z，Liu L H 2 0 1 0 S h a l l o w- l a cu s t ri n e s a n d- r ich d e l t a ic d e