DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 4. 有机质类型 不同类型的有机质，在热演化生烃的过程中，其产物的数量和性质都有较大的差异。研究有机质类型的方法较多 [ 9 ] ，该次利用干酪根的显微组分和碳同位素资料进行研究。干酪根显微组分是在显微镜下直接观测烃源岩中有机质显微组分的组成情况，腐殖型和腐泥型有机质显微组分差异明显，常用类型指数(I t )表示。研究区侏罗系各组12件烃源岩干酪根显微组分组成分析结果(表2)表明，主要以镜质组和腐泥组为主，整体镜质组略微占优，两者占到显微组分体积分数的95%以上。其次含极少量的惰质组，体积分数不高于5%；没有发现壳质组组分。I t 值分布范围：20.3~38.7，均值为3.5。结果表明，侏罗系各组烃源岩有机质类型为II 2 -III型，为偏腐殖型有机质。
system: [SYSTEM]

user: 摘 要 黄河口凹陷古近系储集砂体类型和成岩过程复杂，为储层优选提供依据，需要全面总结黄河口凹陷古近系砂岩储层的岩石学特征和成岩作用特征。通过大量的基础分析测试工作以及钻井薄片观察，将杂砂岩和净砂岩的特征进行区分，明确了净砂岩是古近系的主要储层，并将其划分至3类8种成分-成因储层类型。黄河口凹陷古近系砂岩储层以长石砂岩类为主，地层单元由下至上石英砂岩类比例呈增加趋势，岩屑砂岩类呈减少趋势。砂岩成岩过程中经历了压实压溶、胶结、交代、蚀变、溶解、黏土矿物转化等6种成岩作用，储层岩石学特征的差异影响了不同层位储层的成岩作用及其储层物性。 关键词 :砂岩储层，岩石学特征，成分-成因类型，成岩作用，古近系，黄河口凹陷 Copyright © 2017 by authors， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 黄河口凹陷位于渤海湾盆地东南部，覆盖面积约2370 km 2 ，整个凹陷坐落在古老中朝地台之上，受控于渤中34构造东界的走滑断层和渤南低凸起的南界断层。整体而言，黄河口凹陷内古近系发育较为齐全，由老至新发育有孔店组(E 1-2 k)、沙河街组(E 2-3 s)、东营组(E 3 d) (其中沙河街组三段上亚段( )由于遭受强烈剥蚀而不发育)。黄河口凹陷是渤海湾盆地内的一富油气凹陷，其古近系和新近系油气藏作为两大勘探层系，资源量十分可观 [ 1 ] - [ 7 ] 。鉴于古近系具有良好的生储组合和勘探前景，其储集层一直也是众多专家关注的重点 [ 8 ] - [ 13 ] 。已有学者对E 2-3 s等部分层段的砂岩储层特征进行了描述 [ 8 ] [ 9 ] ，针对储层的主控因素 [ 10 ] [ 11 ] 及平面预测 [ 11 ] [ 12 ] 等内容，也开展了一定的研究工作。但整体而言，黄河口凹陷储层的研究力度还需深入，尤其是从宏观方面针对整个古近系砂岩储层岩石学和成岩作用的研究较为欠缺，影响了未来优质储层的优选以及有利目标的确定，限制勘探的步伐。 该次研究基于基础资料，通过大量的基础分析测试以及岩屑、壁心和岩心的染色铸体薄片或普通薄片观察等方法 [ 14 ] [ 15 ] [ 16 ] [ 17 ] ，系统观察总结、详细统计分析，最终确定了砂岩储层的分类方案，并对不同类型储层进行划分、统计，明确了黄河口凹陷整个古近系砂岩储层的优势类别、变化趋势，分析了其微观成岩作用以及岩石学特征对成岩作用的影响，突破原有研究层系单一的局限，以期为日后的黄河口凹陷的储层评价、油气藏综合研究、勘探目标的优选等提供依据。
system: [SYSTEM]

user: 2. 砂岩储层岩石学特征 2.1. 砂岩储层成分-成因分类 为更清晰地掌握古近系砂岩储层类型，通过砂岩中杂基含量的差异，将砂岩样品划分为杂砂岩和净砂岩(图1)，并将净砂岩储层设定为研究目标，制定了成分-成因分类方案。根据石英、长石和岩屑的体积分数，对净砂岩进行砂岩类型的划分，将其分为石英砂岩、长石砂岩和岩屑砂岩3大类，并依据内部矿物体积分数的差异进一步划分为8种成分-成因类型(图2)。 图1. 杂砂岩与净砂岩微观特征 图2. E 1-2 k-E 2 s 4 (沙河街组四段)砂岩储层成分-成因分类 2.2. 砂岩储层成分-成因类型及分布规律 研究区出现的砂岩成分类型较多，8种成分-成因类型的砂岩均有出现，总体上以长石砂岩类为主，但不同层位砂岩的优势成分-成因类型及发育规律具有明显的差别。 2.2.1 . E 1-2 k- E 2 s 4 储层类型 E 1-2 k-E 2 s 4 为古近系下部地层单元，该次样品数据共有7个。就目前分析结果而言，E 1-2 k-E 2 s 4 中岩屑砂岩类占优势，并以长石质岩屑砂岩居多，所占比例约为57% (由于E 1-2 k和E 2 s 4 岩性地层难以区分，因此将二者做合并分析) (表1，图2)。
system: [SYSTEM]

user: 摘 要 白豹油田157区内主要的含油层位是长4 + 5油层组，通过分析岩性、沉积构造、测井及垂向剖面结构等沉积相标志表明，区内长4 + 5油层组为三角洲前缘亚相沉积。研究认为，长4 + 5油层组油藏的形成和分布主要受沉积相控制，为典型的岩性油藏，油气富集在水下分流河道、河口坝等有利相带中的砂体发育部位，以岩性上倾尖灭油藏为主。 关键词 :长4 + 5油层组，沉积微相，油藏分布 Copyright © 2017 by authors， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 白豹油田157区位于鄂尔多斯盆地东部斜坡带，总体为一平缓的西倾单斜，倾角不足1˚，断层和褶皱不发育，在局部形成起伏较小轴向近东西或北东的鼻状隆起，缺乏油气富集的二级构造带和构造圈闭，主要发育岩性圈闭 [ 1 ] [ 2 ] [ 3 ] 。区内生产已进入开发后期，产量下降快，水淹程度大，产能不断降低，剩余油分布不均。因此对各小层沉积相展布及油藏分布规律进行细致研究，对油田后期的开发极为重要 [ 4 ] [ 5 ] [ 6 ] 。
system: [SYSTEM]

user: 2. 沉积相类型及特征 2.1. 测井相模版 根据取心井观察，并结合相应的测井曲线分析，建立了测井相模版，认为研究区长4 + 5油层组主要沉积类型为三角洲前缘亚相，进一步可细分为水下分流河道、河口坝、席状砂、支流间湾和湖泥共5种微相类型。 1) 水下分流河道 是研究区强水动力沉积区域，沉积物以中-细砂岩、粉砂岩为主，多具正韵律特征。自然电位、自然伽马曲线呈较明显的箱形-钟形。 2) 河口坝 水下分流河道中所携带的泥砂物质，因流速下降从而在河口处堆积下来形成河口坝，主要是砂和粉砂质，沉积物的颗粒自下而上，由细变粗，呈现反韵律特征。测井曲线表现为明显的漏斗形，常与水下分流河道共生，储集物性良好。 3) 席状砂 河口坝或远砂坝在湖浪作用下经淘洗、颠簸，沉积物逐渐向河口两侧的湖岸带迁移从而形成前缘席状砂，常在河口坝两侧呈小规模片状分布，广泛发育交错层理。测井曲线呈低幅指形和粗齿形，在平面上位于河道前方及侧缘。 4) 支流间湾 位于水下分流河道之间，河流带来的悬浮质淤泥、粉砂质泥岩、泥质粉砂岩等沉积物在此汇集，岩石粒度较细、分选性中等，在剖面上，岩性呈复合正韵律或复合正、反韵律。测井曲线主要表现为低幅齿形。 5) 湖泥 位于三角洲前缘沉积的末端，是三角洲体系中分布较为广泛、沉积厚度最大的部位。测井曲线上表现为曲线低幅，且变化趋于平直，偶尔会出现细砂岩的小齿峰。 2.2. 单井相解释 依据白豹油田157区测井相模版，综合各种沉积微相岩性特征、沉积特征、电性特征及物性特征等建立了该区沉积微相综合柱状图(图1)。 图1. 白豹油田157区长4 + 5油层组沉积微相柱状图 图1中深度范围为1925~2015 m井段，自下而上表现为一个海退现象，总体为一个由弱-强的反韵律沉积旋回，其中水下分流河道、河口坝、席状砂及支流间湾等微相均有发育，砂体发育较好。
system: [SYSTEM]

user: 3. 油藏类型及油层分布特征 3.1. 油藏类型 鄂尔多斯盆地是典型的克拉通盆地，其演化较为稳定，具有持续沉降、整体抬升、坡降宽缓、背斜微弱、近于水平接触等特点，缺乏形成构造型油气聚集或富集的二级构造带 [ 7 ] [ 8 ] 。白豹油田157区长4 + 5油层组的油层、干层沿上倾方向岩性尖灭，具典型的岩性油藏特征(图2)。 图2. 白豹油田157区长4 + 5油层组油藏剖面图 3.2. 油层分布特征 区内长4 + 5油层组的油层厚度一般在4~16 m，生产层位中油层厚度在4~8 m约占49.43%，厚度为12~16 m的厚油层约占7.57%。从油层展布情况来看，油层主要分布在研究区水下分流河道及河口坝砂体较为发育处，厚度均较厚，含油连片性较好。而席状砂砂体内油层厚度较薄，且断续分布。
system: [SYSTEM]

user: 4. 结论 1) 白豹油田157区长4 + 5油层组主要发育三角洲前缘亚相，进一步可细分为水下分流河道、河口坝、席状砂、支流间湾、湖泥等微相类型，以水下分流河道沉积为主。 2) 区内长4 + 5油层组油藏未经明显的构造运动破坏，保存条件较好，长4 + 5油层组各小层在油层、干层沿上倾方向存在岩性尖灭，为典型的岩性油藏。
system: [SYSTEM]

user: 摘 要 砂体厚度(或含量)是油气勘探中的重要参数。以地震资料与测井解释成果为基础，研究储层砂体厚度预测方法，提出利用地震属性优化技术实现降维，建立敏感地震属性集；并将井旁地震道对应的敏感属性集与测井解释砂体厚度输入神经网络，通过网络训练使误差最小化；在该基础上，逐道输入敏感属性集，由网络输出对应的砂体厚度。在胜利油田H4区块应用上述方法预测砂岩厚度，相对误差基本小于20%，满足了勘探生产的需要。 关键词 :砂体厚度，地震属性优化，网络函数逼近，储层预测 Copyright © 2017 by author， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 地震资料解释与反演的目标就是从三维地震数据中提取更多的信息参数，该信息被称为地震波特征参数或地震属性，是将地震数据作不同的数学运算得到的一种新参数 [ 1 ] - [ 6 ] ，利用该参数可以解决地质构造、预测地层岩性和油气等油气勘探的关键问题。地震属性是指叠前或叠后地震数据经过数学变换而导出的有关地震波几何学、运动学、动力学和统计学的特征，是地震资料中可描述的、可定量化的特征。地震属性的种类较多，早期主要以振幅属性、时间属性、频率属性和地层吸收属性等为主。近年来，随着数学、图像与信号处理等相关领域技术的发展，计算地震属性的方法越来越多，提取的地震属性也越来越多。地震属性从不同角度反映了地层砂体厚度、孔隙及孔隙流体、渗透率、地层异常压力等地层信息。关于砂体厚度预测的方法较多，可归结为两大类：一是以波阻抗反演为基础，二是地震属性(叠前和叠后)综合应用。该次研究属于地震属性的综合应用，笔者在前人研究的基础上，将属性优化降维与网络函数逼近相结合，优选出对砂体敏感的地震属性集，进而研究储层砂体厚度分布。
system: [SYSTEM]

user: 2. 属性优化与神经网络函数逼近方法 2.1. 地震属性优化 储层参数预测是要将地震属性通过某种数学算法(如模式识别、神经网络)赋予地质含义。地震属性种类繁多，目前国内外用于储层预测的地震属性有300余种，这些地震属性有独立的，也有相互关联的；有些具有直接的地球物理意义，有些则没有。地震属性优化就是要在储层预测中消去不敏感地震属性的影响，减少有关联的地震属性影响，提升敏感的、独立的地震属性影响。 在利用地震属性进行储层预测(包括储层砂体厚度、孔隙度、饱和度等)中，通常提取一组属性集，并结合模式识别、多元统计、灰色系统或神经网络法等进行储层预测。但是，不同地区、不同层位、不同的储层参数(孔隙度、饱和度、砂体厚度等)，相应的地震属性敏感度不同，其中K-L变换(Karhunen-Loeve transform)能有效地压制相关属性，提高独立属性的作用。该次研究采用K-L变换法实现降维运算，提高独立属性在砂体厚度预测中的影响。 2.2. 神经网络函数逼近方法 对于神经网络函数逼近法作储层参数预测，关键是网络的选择，神经网络有40种左右，结构、性能各不相同，但它们都是由大量的基本处理单元广泛连接而成，基本单元称为神经元，也称为节点，节点是所有输入的加权和，并通过一个非线性函数输出结果。网络性质的决定要素有3个，即神经元特性、网络结构、训练方法。该次研究使用的神经元的激活函数是S型函数。神经网络结构是多层感知器，如图1所示，多层感知器是一种层状结构的前馈神经网络，它由输入层、输出层和一个或多个中间层构成；网络中的每个节点只与邻层(上、下层)的各个节点相连接，与同一层间的各节点不相连，各节点之间连接的权系数大小不同。 图1. 储层砂体厚度预测的网络结构 在网络训练过程中，首先由计算机自动生成一组小的随机数，并将其作为网络内部各个节点之间连接的权系数和各个节点上的阈值的初始值；然后，输入所有训练样本集，即K-L变换优化后得到的敏感地震属性集，计算网络输出结果，进而计算实际输出与期望输出的差值(实际输出是网络计算得到的输出结果，期望输出为测井解释的砂体厚度)；最后比较差值是否满足误差要求，若不满足误差要求，按照误差反传播算法(BP算法)的规则，修改节点间连接的权系数大小和节点上阈值的大小，不断反复上述过程，直至满足误差要求，此时估价函数降至可接收值，将节点间连接的权系数大小和节点上阈值的大小固定下来，作为下一步预测使用。 在利用网络函数逼近法作储层砂体厚度预测时，使用如图1所示的3层感知器作函数逼近，一个输入层、一个输出层和一个中间层，导出的关系式为： (1) 式中：Y k 为网络的输出层各节点；k为输出层节点数，该次研究只做砂体厚度预测，k取1；M为中间层的节点数，该次研究通过计算比较，取M = 9时收敛快、精度高；N为输入层的节点数，取敏感属性集的分量个数；W 2 为中间层节点与输出层节点间的连接权；W 1 为输入层节点与中间层节点间的连接权；q 1 为输入层到中间层节点上的阈值；q 2 为中间层到输出节点上的阈值。
system: [SYSTEM]

user: 3. 砂体厚度预测流程 1) 地震资料的构造精细解释。在层位标定的基础上，对目的层顶、底界面作地震构造精细解释，并提取目的层顶、底界面之间的地震数据。 2) 地震属性的提取，建立敏感属性集。对目的层段地震数据，作傅里叶变换、小波变换、希尔伯特变换、相关运算、谱比法吸收系数提取等不同的数学运算，得到目的层段的地震属性。 3) 地震属性的优化。对上面提取的地震属性集通过K-L变换优化，确立对砂体厚度敏感的属性集。 4) 神经网络的训练。设计神经网络的结构，根据井旁敏感的地震属性和测井解释砂体厚度或录井测得砂体厚度，建立神经网络训练集，输入设计的网络，按照BP算法原理进行网络训练，确定神经网络内部各个节点之间连接的权系数和节点上的阈值。 5) 储层砂体厚度分布预测。将敏感地震属性集逐道输入训练后的网络，输出即为该道对应的砂体厚度。
system: [SYSTEM]

user: 4. 应用效果与认识 在构造精细解释的基础上，以地震属性提取为基础，以K-L变换地震属性优化和神经网络函数逼近为手段，对胜利油田H4区块的沙河街组二段10亚段上( )、沙河街组二段10亚段下( )等层段的砂体厚度分布进行预测。在应用过程中，重点考虑K-L变换优化后的敏感地震属性集，选取合理的网络结构，并充分利用测井资料综合解释的砂岩厚度作为约束条件。图2和图3分别是 和 砂体厚度平面分布图，将2个层位近60口井的砂岩厚度预测结果与井资料进行对比分析可知，二者符合程度较高，相对误差基本小于20%，砂体平面展布规律与沉积相分析规律基本一致，满足了勘探生产的需要。
system: [SYSTEM]

user: 5. 结论 1) 不同的储层参数，通过地震属性优化分析能提取反映敏感的一组地震属性，优化后的属性集有利于提高储层参数预测精度。 2) 该次研究提出的基于地震属性优化与网络函数逼近储层砂体预测方法是可行的，在胜利油田H4区块的应用中，预测结果与井资料的符合率较高，相对误差均在20%以下。 图2. 研究区 预测砂岩厚度分布图 图3. 研究区 预测砂岩厚度分布图
system: [SYSTEM]

user: 摘 要 地震成果由时间域转换到深度域，涉及速度建场和变速成图技术。随着水平井的广泛应用，对转深精度要求更高。特别在复杂山地探区，地震资料品质和地层速度在横向和纵向上变化比较剧烈，加之基准面计算方法不合理，导致转深的区域构造形态与深度预测精度往往很低，使得水平井的钻探设计和成本控制变得十分困难。通过在建场过程中对各个环节减小误差，进而提高整体速度场精度。对引起速度场建场误差的关键环节进行分析后认为，从基于表层模型的CMP (共中心点)基准面转换、模型迭代层速度计算、地震反射层约束速度场建立、两步法层速度误差校正方法入手，降低各个环节的误差，能够提高适用于盆边山地过渡带的速度场精度。 关键词 :速度场建模，模型迭代层速度计算，两步法误差趋势面校正 Copyright © 2017 by author， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 涪陵地区地处四川盆地和盆边山地过渡地带，构造位置位于四川盆地川东南构造区川东高陡褶皱带，万县复向斜包鸾–焦石坝背斜带中的焦石坝背斜，构造呈北东向展布。境内地势以低山丘陵为主，横跨长江南北、纵贯乌江东西两岸。涪陵地区地震地质条件复杂，速度变化剧烈，传统速度场建场方法误差大，不能满足实际生产需求。在研究该地区表层结构、深层地质和速度资料特点的基础上，通过分析影响速度建场和变速成图精度的因素，从基于表层模型的速度谱CMP (共中心点)时间面到深度面的转换、模型迭代层速度计算、基于层位约束的两步法层速度误差趋势面构建等方法入手，建立了适用于复杂山地的高精度速度场建场技术。该技术在焦石坝地区的水平井钻井实践中，减少了钻探过程中井轨迹的调整次数，提高了水平井设计和钻探精度，降低了成本，取得了很好的效果，具有较高的生产应用价值。
system: [SYSTEM]

user: 2. 高精度速度场建场思路 高精度速度场建立，实际上就是提高速度精度、消除速度误差。通过对焦石坝三维地震工区速度误 差产生的原因分析，主要是在基准面转换、层速度计算、地震速度与测井速度差异、速度场校正方法等方面产生的误差。针对误差产生的特点，采取在建场过程中对各个环节进行误差校正的思路，通过降低每一环节的误差，使整体误差降低，进而提高速度场精度。具体方法为：采用基于表层模型的基准面转换技术，减小CMP基准面转换误差；采用模型迭代法计算层速度，用以减小直接使用Dix公式产生的误差；利用层约束速度建场技术减小层速度在纵横向空间上的误差；采用两步法误差趋势面校正技术减小井震之间的误差；利用钻探开发的时移性，用新的钻井信息对最终速度场进行局部误差校正，使速度场不断求精。高精度速度场建场流程如图1所示。 图1. 高精度速度场建场流程图
system: [SYSTEM]

user: 3. 关键技术方法研究 3.1. 基于表层模型的CMP基准面转换 在地震资料解释中，关键的一个环节是时深转换，目前时深转换选取的基准面均为速度分析基准面，即CMP面 [ 1 ] 。由于CMP面是速度谱的零线，是时间域的浮动面，是静校正量中的低频分量 [ 2 ] ，不能简单地作为时间域转换的基准面，其大小受表层厚度、速度、基准面高程以及充填速度影响。传统时深转换方法大多忽略了CMP基准面转换对深度误差的影响，特别是在复杂山地地区，仅用替换速度和CMP时间量计算转换量，会造成与实际高程值的不符。图2为建立的水平模型，其中(a)为表层深度模型，(b)为(a)的时间模型，当中包含了CMP时间面。从图2(a)与图2(b)的对比中可以看出，若不考虑低降速带的速度和厚度，计算出来的CMP基准面高程与实际CMP基准面(地表面)高程之间误差为35 m，该误差势必会造成基准面转换时的速度误差，从而影响时深转换的精度。 采用基于表层模型的浮动基准面转换技术，考虑表层模型对速度的影响，可以提高基准面转换的精度，使速度更趋向真实。改进的CMP基准面高程的计算方法为：考虑静校正的替换计算，采用反替换(或称为反静校)来求取CMP基准面高程。 常规方法： (1) 图2. 计算CMP基准面高程误差示意图 改进方法： (2) 式中：H cmp1 为常规方法计算得到的CMP基准面高程，m；H d 为固定基准面高程，m；T cmp 为CMP校正量，s；v r 为替换速度，m/s；H cmp2 为改进方法计算得到的CMP基准面高程，m；H g 为高速顶(低速带底)高程，m；v 0 为低速度带速度，m/s。 从改进公式(式(2))中可以看出，考虑了替换速度、低降速层速度和厚度的影响，计算得到更为接近实际情况的CMP面高程。基于表层模型的速度谱基准面转换，是使时间域CMP基准面向深度域转换误差更小的有效方法。换算后的CMP直接转深面不仅受地表高程影响，也受高速层顶界面起伏影响，使得换算后的基准面与地表面并非完全一致，而采用基于表层模型的CMP转深面更贴近地表面(图3)。 图3. 地表面、基于表层模型的CMP转深圆滑面、CMP直接转深面圆滑面对比 3.2. 模型迭代层速度计算 常规速度建场中，通常采用Dix公式将均方根速度转换成层速度或平均速度。但Dix公式是在地下水平均匀层状介质、射线垂直入射的假设前提下建立起来的数学关系式。对于速度各向异性较强和速度出现倒转现象的地区是不适合的。采用模型迭代方法进行反演是解决层速度精度问题的有效途径之一 [ 3 ] 。模型迭代方法为：首先由Dix公式给出初始层速度模型；然后利用T 0 计算的谱点对应的倾角，通过射线追踪法求出反射波时距曲线，得到模拟速度谱；再用理论曲线对其进行拟合，求出理论叠加速度；最后将该速度与速度谱上的叠加速度进行比较，当误差小到合理精度时，即认为所给模型是正确的，输出所求的层速度值。该方法不受模型是水平层状介质和速度不能翻转的限制，计算精度要高于Dix公式。模型迭代流程如图4所示。 图4. 模型迭代流程图 3.3. 反射层位约束层速度场 由地质人员解释得到的T 0 图，转换为能够定井位的构造图，需要平均速度作为桥梁。平均速度的精度直接影响构造图的精度。对平均速度的直接编辑、平滑，会影响到各个目的层，导致当前目的层精度提高，其他目的层精度下降。为了解决该问题，引入地质反射层约束技术，原理是利用叠加速度和若干地质反射层，建立按T 0 层位变化的层速度场(同时也可以按产状对地质分层进行细分小层)，误差的校正均在约束的层速度场上进行，对分层内层速度进行编辑、平滑，既可以保证当前层速度的纵向和横向速 度变化趋势，使层间速度更能代表实际速度的变化规律，又不影响其他层的速度变化，用该方法来提高速度场精度。图5~图7为间隔抽取的约束前、后的层速度对比图，可以看到，反射层约束前的层速度(图5)不能反映当前层速度的纵向和横向速度变化趋势，通过地震剖面解释的层位(图6)约束后的层速度(图7)更能反映层速度沿层的变化规律。 3.4. 两步法误差趋势面校正 速度误差产生的原因除了基准面转换、层速度计算等因素外，还有地震分辨率对速度谱影响产生的 图5. 反射层约束前层速度剖面图 图6. 地震解释层位剖面图 图7. 反射层约束后层速度剖面图 误差、岩性变化引起的误差、地震速度与井筒测量的速度间误差、追踪层位精度引起的误差、资料处理中速度陷阱等 [ 4 ] 。上述误差参杂在一起，反映到井点上也就是井点误差。以往对误差的校正主要采用数学插值网格化运算，以求出误差分布的趋势面，这与井的数量、分布均匀程度都有密切关系，其结果是在井点处精度较高，但其他部位容易引起变形扭曲，使得预测精度不高 [ 5 ] 。 通过研究认为，综合误差可以分解成规律误差和残余误差两部分。规律误差是与整体速度场相关的误差，为随时间、速度、深度等因素变化的函数，相对较平缓(图8)。残余误差是与具体位置相关的、局部的误差。两步法误差趋势面校正就是针对上述两种类型误差进行分别校正。该方法的优势是规律速度校正时不使用井数据，校正后基本保留了地震速度的变化趋势，能够较好地保持构造形态(图9)。 图8. 目的层规律速度校正前、后井点误差对比图 图9. 目的层两步法校正后的剩余误差分布图
system: [SYSTEM]

user: 4. 效果分析 利用上述方法在104口正钻井中进行应用，统计完钻的100口水平井，水平段共计152935.2 m，优质储层段钻遇率为85%，A靶点调整后的深度与实钻靶点深度误差小于20 m的占91.3%，水平井的入靶 精度和穿行效率均较高，取得了较好的应用效果。
system: [SYSTEM]

user: 5. 结语 笔者提出的引起速度场建场误差的几个关键环节，以实际勘探开发需求为驱动，改进了复杂山地地区传统建场方法和研究思路，通过实际应用，使速度场建场的精度得到了很大的提高。
system: [SYSTEM]

user: 摘 要 根据标准模型下的基本粒子结构，一种力可以确定为超对称互反费米子之间的引力。 关键词 :费米子，互易矢量，超对称，中微子，引力波 Copyright © 2017 by author and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 1) 有限加法群 设 为集合G的三个元素。即 具有下列二元运算： 对于这种具有二元运算的集合，它是适合下列条件的 a) 结合律成立 b) G中存在一个元0，使得 c) 对于G中任意元 存在逆元， d) 适合交换律： 因此，G为交换群(AbeI)。 2) 几百年前，伟大的科学家牛顿给出了万有引力定律，几十年前著名的物理学家爱因斯坦创建了引力场理论，为宇宙空间的研究做出了历史性的贡献。随着现代物理学的研究深入，万有引力量子化是目前理论物理学最前沿，也是最需要必须突破的课题。四种力的大统一，将会带来整个物理学革命性的变化。粒子物理学的研究，特别是规范场的建立，为我们解决这个难题提供了理论基础。规范场 [ 1 ] 告诉我们，电子，轻子等基本粒子都是不同的费米子通过强力，电磁力，弱力形成的。费米子之间是否存在超对称场，是否可以量子化处理，费米子之间的超对称场是引力场吗？
system: [SYSTEM]

user: 2. 场论分析 1) 对于两个费米子群 设定量子基态电荷相等，其中 ( 为度量系数 [ 2 ] )则 与 与 与 为互易矢量。如果存在下列两个矢量基态矩阵线性相关 (根据有限加法群) (1) (2) 则费米子 构成超对称互反费米子场。 (两个矩阵的线性相关很显然) 2) 如果将一个质点看出是费米子群，则这个质点的质量是费米子群的质量总和。设两组费米子S力矩，由万有引力定律可以得出 (K为引力系数) (3) 对于 当S为0时，G为强力且形成暗物质。 当S > 0时，G为质点之间的引力。 3) 如果两组费米子中有一对互易矢量形成强力，由于能量不守恒 [ 3 ] ，则构成四维绕轴旋转动态粒子模型(图1)，且分水平滚动式动能波–引力波，与垂直螺旋式势能波–中微子。 物理状态分析 垂直螺旋式势能波 (4) ( 中任一活动标架， 为角速度) (5) 水平滚动式动能波 (6) (7) 垂直螺旋式势能波粒子由于存在三个不同方向的活动标架，相应存在三个不同相位的中微子；水平滚动式动能波粒子与两个无关矢量(标架)相关。
system: [SYSTEM]

user: 3. 结束语 1) 很显然，这种超对称引力 与对称费米子数量成正比，与距离的平方成反比。结论不但适用于质点之间的引力结构量子化分析，还适用于微观结构，如原子核与电子之间的量子力学分析。更为我们了解宇宙结构提供了理论支持。 2) 动态粒子模型的构造，从理论上给出了中微子和引力波的量子化解析，特别是对引力系数的实质给出了量子化解释。 图1. 四维绕轴旋转动态粒子模型 引力量子化，是目前理论物理学研究的最大课题，愿本文给出的数学分析给物理学特别是粒子物理学与宇宙空间研究带来新的突破，也希望学术界特别是粒子物理学做进一步的验证。
system: [SYSTEM]

user: 摘 要 随钻伽马测井地质导向是提高水平井储层有效钻遇率的重要方法，但目前研究随钻伽马测井的方法计算速度有限。为实现实时地质导向，基于水平井地层中随钻伽马测井探测范围的空间特征，确立了随钻自然伽马测井的正演方法；利用该正演方法，建立了在不同地层中离边界不同距离情况下的随钻自然伽马测井响应图版；再利用该图版，建立普遍适用的随钻伽马测井快速反演储层界面的方法。在双层介质模型和随钻实例地层中，可以快速、准确地反演出地层界面的位置，表明该方法可用于随钻伽马测井的实时地质导向工作。 关键词 :随钻伽马测井，水平井，快速反演，储层界面，正演算法 Copyright © 2017 by author， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 由于水平井在单井产量、经济效益以及在非常规油气勘探中的优势明显，所以应用十分广泛。当前水平井钻井技术的挑战之一是实时地质导向 [ 1 ] [ 2 ] [ 3 ] [ 4 ] 。水平井钻前地质模型具有一定的不确定性，故随钻实时地质导向对于提高储层有效钻遇率显得极为必要，而实时地质导向的关键在于实时预测地层界面 [ 3 ] [ 4 ] 。国内常用随钻自然伽马测井和电阻率测井方法进行地质导向，随钻电阻率因其探测深度较深，常用于确定地层界面。但在低阻、低对比度储层或薄层中，因受到围岩影响，随钻电阻率用于指示地层界面的可信度降低，而自然伽马具有更高的分辨率，在上述储层中可以较为准确地指示岩性和界面 [ 5 ] [ 6 ] [ 7 ] [ 8 ] 。目前，常用蒙特卡洛等方法对随钻伽马测井进行正反演研究。考虑到该类方法的计算速度和时效性问题 [ 5 ] [ 9 ] [ 10 ] ，该次研究首先基于地层中自然伽马强度的分布原理，分析了常规伽马测井的响应模型，然后根据水平井中随钻伽马测井探测范围的空间特征，确立了水平井中随钻自然伽马测井的正演方法；利用该正演方法，确定了在不同伽马特征的地层中，离地层边界不同距离情况下的随钻自然伽马测井响应图版；然后利用该图版，建立了随钻伽马测井快速反演储层界面的方法。在不同伽马比例的双层介质模型中进行测试，该方法可快速反演出地层界面；将该方法应用在理论地层模型和随钻实例地层中，可以快速、准确地反演出地层界面的位置。因此，该方法可用于指导随钻伽马测井的实时地质导向工作。
system: [SYSTEM]

user: 2. 随钻伽马测井和正演模型 2.1. 随钻伽马测井 自然伽马测井是用于测量地层天然放射性伽马射线强度的测井方法 [ 10 ] 。在原理上，随钻自然伽马与常规伽马测井方法类似，常用于指示水平井井周的地层岩性和界面，在薄层、低阻储层的随钻地质导向中应用广泛 [ 5 ] [ 7 ] [ 8 ] [ 10 ] 。通常，随钻仪器中伽马测量短节会采用近钻头测量方式，近钻头测量能准确无滞后地测量钻头处的储层伽马特性，及时确定地层和边界位置，根据地质信息调整井眼轨迹，实现实时的地质导向 [ 8 ] [ 11 ] 。 2.2. 无限均匀介质自然伽马响应模型 考虑地层模型为无限均匀介质，设地层密度为ρ，每克岩石中含q克放射性物质，每克放射性物质每秒发射的伽马光子数为a，地层对伽马光子的吸收系数为μ，取球坐标系(见图1)，那么体积元dv在距离为r的观测点处产生的射线强度为 [ 5 ] [ 10 ] ： (1) 式中：J r 为自然伽马通量，s −1 ∙cm −2 ；a为伽马光子数，g −1 ∙s −1 ；q为每克中含岩石放射性物质的质量分数，g/g；ρ为地层密度，g/cm 3 ；r为观测点的探测距离，cm；μ为地层对伽马光子的吸收系数，cm −1 ；v为球坐标微元体积，cm 3 ；θ为球坐标天顶角，(˚)；φ为球坐标方位角，(º)。 对式(1)积分后可得均匀球体空间在球心观测点产生的自然伽马通量： (2) 根据式(2)，假设探测点处99%的自然伽马通量由半径为R的球体造成，则有1−e −μr = 0.99，由此可得μr = 4.065，从而根据不同地层的μ计算出相应的R。一般自然伽马测井仪器设计的探测深度约为0.3 m，如果考虑套管、钻杆、钻井液等因素影响，仪器实际探测深度一般为0.2 m左右 [ 5 ] [ 10 ] [ 12 ] ，该次研究设定探测深度R = 0.3 m。 图1. 自然伽马探测范围球坐标系积分示意图 2.3. 水平井层状地层自然伽马测井正演模型 邵才瑞等 [ 5 ] [ 6 ] [ 10 ] 分析了水平井中随钻伽马测井的探测范围与地层空间系(图2)，基于伽马射线在地层中的衰减规律，建立了随钻伽马仪器位于地层不同位置和不同地层模型中的快速正演算法。考虑到储层厚度通常是随钻伽马测井探测直径的数倍 [ 2 ] [ 6 ] [ 13 ] ，该次研究中薄层情况不予考虑，故图2中受到地层1和地层2的界面影响时，仪器探测范围内不同区域(A、B、C)的伽马射线强度为： (3) (4) (5) 式中：J A 、J B 和J C 分别代表图2中A、B和C三个体积域在观测点O处的伽马射线强度，s − 1 ∙cm −2 ；参数a、q、ρ、μ的下标1、2与图2中的地层序号1、2一致；r 0 为探测深度，cm；φ 0 为图2中地层界面与探测范围的夹角，(˚)；d为探测器到地层界面的距离，cm；f是金格函数。 上述算式中使用的金格函数 [ 12 ] ： (6) 式中：x为金格函数变量，1。 图2. 横向各向同性砂岩、泥岩层状地层伽马探测示意图
system: [SYSTEM]

user: 3. 随钻伽马测井快速反演储层界面方法 随钻伽马测井正演模型无法直接用于储层界面距离的反演，假设随钻伽马仪器由地层1 (砂岩)穿过界面进入地层2(泥岩)时(图2)，设界面处仪器离界面距离为0，在界面之上和之下时d分别为负值和正值，利用上述正演模型模拟出离界面距离不同时自然伽马强度的分布情况(图3)。由图3可以看到，当d大于0.3 m后，仪器响应几乎不受影响，因此在界面附近0.3 m处自然伽马响应才能反映出界面的影响；当界面上、下地层自然伽马比值一定时，随钻伽马仪器响应是该比值的确定函数；当地层1与地层2的自然伽马比值C为1:12时，离界面的距离是伽马不同比值的函数，当C为其他值时，可以看到随钻伽马测井具有类似的响应特征曲线。因此，图3可以作为确定随钻伽马离界面距离的图版使用。 为快速获取储层界面信息，利用该图版建立快速反演模型，当C = 1:12时，在 区域，可获得较好的拟合曲线(相关系数近似为1)，如图3中红色虚线所示，公式如下： (7) 类似地，可以获得不同C情况下d的函数，建立起离界面距离的快速反演方法。 图3. 在不同地层水平井中随钻伽马离界面不同位置时的响应特征
system: [SYSTEM]

user: 4. 界面快速反演测试 为测试该快速反演储层界面方法的正确性，分别利用理论地层界面模型和现场实例来验证该方法的有效性。 4.1. 双层介质模型界面反演 建立理论层状地层模型(图4)，水平井轨迹从上往下穿过泥岩层(150API)和砂岩层(30API)，依据水平井轨迹穿过该地层的情况，利用式(3)~(6)正演计算得到图中的自然伽马测井响应特征；再利用上述快速反演储层界面模型，反演得到图4中离界面的距离情况，可以看到，在随钻自然伽马测井探测范围内(0.3 m)，离界面的距离随水平位移变化曲线与井轨迹几乎重合，对比实际界面的距离绝对误差仅0.01 m，证明了快速反演方法的准确性。 图4. 理论层状地层模型中的响应特征 4.2. 实际地层界面反演 新疆油田某水平井S井，选取水平井轨迹靠近储层边界的深度段，为清晰显示界面反演深度段的对比情况，选择井轨迹起点为1000 m的垂直深度，位移起点为0 m，可以观察到在水平位移约20 m时，井轨迹钻遇界面，且上方泥岩层自然伽马为107API，下方砂岩层自然伽马为53API；考虑到随钻自然伽马的探测范围约0.3 m，利用快速反演储层界面方法反演了水平位移20~165 m的储层界面，共选取了5个界面关键点(如图5中A、B、C、D、E，e a 为绝对误差)；对比利用电法、声波和密度测井等曲线综合反演出的最终真实界面位置 [ 8 ] [ 11 ] [ 13 ] ，可以看到利用快速反演储层界面方法得到的储层界面准确性高，最大绝对误差仅0.06 m (A点)，平均绝对误差为0.03 m。考虑到地层横向上自然伽马特征的非均质性特点，误差分析结果表明快速反演储层界面方法可较为准确地反演储层界面，为实时地质导向工作提供参考。 图5. S井随钻伽马测井反演储层界面示意图
system: [SYSTEM]

user: 5. 结论 1) 利用随钻伽马快速正演算法，建立了不同伽马特征地层的界面与随钻伽马测井响应之间的关系图版，在该基础之上，建立了适应不同自然伽马特征地层的随钻伽马测井快速反演储层界面的方法。 2) 在建立的理论地层界面模型中，快速反演储层界面方法计算的界面绝对误差仅0.01 m；在S井随钻伽马测井现场实例中，该方法计算界面的最大绝对误差仅0.06 m，平均绝对误差为0.03 m。理论和实际地层界面模型反演结果表明，该方法可获得快速、准确的地层界面信息，可为随钻伽马测井的实时地质导向提供界面参考信息。
system: [SYSTEM]

user: 摘 要 胜北洼陷低幅度构造油藏受孔隙结构和油藏高度影响，油气选择性成藏特征明显，储层定量评价较困难。研究表明，储层特征与束缚孔隙度有一定相关性，利用束缚孔隙度建立了储层品质和饱和度评价新方法，较好地解决了低幅度构造油藏储层评价的难题，在胜北洼陷应用效果较好。 关键词 :低幅度构造，束缚孔隙，孔隙结构，饱和度 Copyright © 2017 by authors， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 研究区位于吐哈盆地台北凹陷胜北洼陷次凹中部，受平行断层切割，且位于油气运移充注的主力方向，形成了自北向东的胜北洼陷2号、3号及4号局部断背斜构造。该区自上而下在侏罗系 七克 台组、三间房组、西山窑组以及喀拉扎组(J 3 k)见到丰富油气显示。J 3 k以曲流河沉积为主，储层为曲流河河道砂体，厚度小。J 3 k油气藏闭合高度30 m，背斜两翼倾角小于3°。受低幅度构造和储层孔隙结构影响，J 3 k油藏油气选择性成藏特征较为明显，油水系统较复杂。
system: [SYSTEM]

user: 2. 储层特征 J 3 k储层岩性主要为棕色泥质粉砂岩、粉砂岩互层。依据胜北洼陷7口井512块岩样粒度分析数据统计(图1)表明，J 3 k岩性主要为细砂岩，粉砂岩次之，细砂岩和粉砂岩占总样品的95%，储层岩性粒径分布较广，分选较差。
system: [SYSTEM]

user: 3. 孔隙结构 经取心铸体薄片(图2)观察，胜北洼陷J 3 k储层孔隙以粒间孔隙为主，受粒间填系物充填影响，孔隙结构较为复杂，局部发育溶蚀粒间孔、粒内孔，总体孔隙结构连通性较好 [ 1 ] 。 图1. 胜北洼陷平均粒度(a)及岩性(b)分布图 图2. 胜北洼陷J 3 k铸体薄片 经扫描电镜(图3)观察，J 3 k见方沸石，局部见伊-蒙混层；伊-蒙混层以黏土桥形式充填孔隙，卷曲丝发状分布于颗粒表面，石英次生加大I级，颗粒溶蚀孔隙普遍发育，粒间孔发育为中等~好，微孔发育差~中等，胶结物内微孔、颗粒内溶蚀微孔发育。 图3. 胜北洼陷J 3 k扫描电镜 J 3 k储层岩石结构成熟度较低，碎屑粒径分布较广，导致孔隙结构差异较大，压汞测试中启动压力和最大进汞饱和度差异明显(图4)，储层孔隙结构相对较复杂 [ 2 ] 。孔隙度主要分布在6%~19%，平均13.95%；渗透率主要分布在1~100 mD。孔隙度与渗透率呈指数关系(图5)，表明储层渗透性主要受孔隙度影响，但相同孔隙度下渗透率差异能达到3个数量级，说明相同孔隙度下储层孔隙结构有较大差异。 图4. 胜北洼陷毛细管压力曲线 图5. 胜北洼陷J 3 k孔隙度-渗透率交会图
system: [SYSTEM]

user: 4. 储层特性评价 胜北洼陷油藏为构造岩性油藏，储层特征基本符合阿尔奇公式砂岩模型，阿尔奇公式能间接反映储层孔隙结构特征 [ 3 ] 。 (1) 式中：S w 为含水饱和度，1；a、b为岩性系数，1；m为胶结指数，1；n为饱和度指数，1；f为孔隙度，1；ρ w 、ρ t 分别为地层水电阻率、地层真电阻率，Ω×m。 对公式(1)双边取对数获得公式(2)，经过整理得到公式(3)： (2) (3) 为了保障钻井施工安全，泥浆柱压力一般大于地层压力，使得泥浆在钻井过程中向地层侵入，地层侵入液体主要为泥浆滤液。由于测井过程中浅电阻仪器探测范围小，探测范围内地层被泥浆滤液冲洗完全，测井值受储层所含流体属性影响小，主要受储层孔隙特征影响。岩石受泥浆冲洗部分孔隙为储层有效孔隙，根据电阻率计算孔隙度关系，有效孔隙度与浅侧向电阻率之间应存在如下关系： (4) 式中： 为有效孔隙度，%；ρ mf 为泥浆滤液电阻率，Ω×m；ρ lls 为浅侧向电阻率，Ω×m；S wmf 为泥浆滤液饱和度，%。 浅地层完全被冲洗时，有效孔隙中的泥浆滤液饱和度为1，则： (5) 利用岩心刻度测井，根据孔隙度曲线拟合计算总孔隙度 ，则 (其中， 为束缚流体孔隙度，%) 。当束缚流体孔隙度越小时，岩性越纯，孔隙结构相对简单；当束缚流体孔隙度越大时，孔隙结构越复杂。图6为建立的束缚流体孔隙度与储层品质系数(K/f)之间的关系(K为岩心分析渗透率，mD)可以看出，束缚流体孔隙度与储层品质系数呈负相关关系，且相关性较好。该方法有效地解决了复杂孔隙结构下同孔、不同渗的问题。 (6) 图6. 束缚流体孔隙度与储层品质系数交会图 胜北洼陷为低幅度构造，岩石亲水特性明显。由于油藏幅度低，油藏成藏压力小，岩石表面束缚水含量较高 [ 4 ] ，则岩石主体导电特性可以分为泥质、束缚水和纯砂岩条件下的自由水并联导电。理论模型如下： (7) 式中：σ为岩石电导率，S/m；φ(sh)为泥质体积分数，% ； 为束缚水饱和度，% ；σ w 为地层水电导率，S/m；σ sand 为含水纯砂岩电导率，S/m。 根据电阻率与电导率之间关系，利用式(7)可以推演出： (8) 式中：ρ a 为地层视电阻率，Ω×m；ρ sh 为泥岩视电阻率，Ω×m；ρ sand 为砂岩视电阻率，Ω×m。 排除泥岩和束缚水导电，岩石导电模型中砂岩导电特性更加符合胜北洼陷纯砂岩阿尔奇公式岩电关系(a = 0.9844、b = 0.916、m = 1.987、n = 1.791)，利用纯砂岩模型计算含水饱和度近似为可动水饱和度S wm ： (9) 则含水饱和度为： (10)
system: [SYSTEM]

user: 5. 应用效果 利用上述方法对胜北洼陷A井含水饱和度进行计算处理，结果如图7所示，可以看出，在相对高部位、 图7. 胜北洼陷A井测井处理综合图 岩性较纯的井段，新模型计算的含水饱和度(S w，c )与阿尔奇公式计算的含水饱和度(S w，a )相差较小；在位置较低的井段，S w，c 小于S w，a ，S w，c 的曲线形态与录井信息更为相符。胜北洼陷A井于 2013 年 11 月 6 日 对2979~2991 m、2994.5~2999 m和3003.4~3005.4 m井段射孔，4 mm油嘴，日产液13 m 3 ，日产油10.1 t，日产气15000 m 3 ，与解释结论相符。综合统计，该方法在胜北洼陷的解释符合率达到83.4%，效果较好。
system: [SYSTEM]

user: 6. 结论 1) 胜北洼陷油藏为低幅度构造油藏，由于受储层岩石成熟度较低、粒径分布宽、岩石分选差等特征影响，储层孔隙结构差异较大。 2) 研究发现低幅度油藏的储层品质系数与束缚流体孔隙度之间呈较好的负相关关系。 3) 利用束缚流体孔隙度导电模型，推导出低幅度构造油藏的储层含水饱和度模型，实际应用效果较好。
system: [SYSTEM]

user: 摘 要 Karamkas油田油井供液能力强，挖潜空间大，目前的工作制度无法充分发挥油井的供液能力。将产液量看作冲程、冲次、泵径、下泵深度的函数，建立以产液量最大为目标的优化模型。根据现场设备的工作参数，该模型的解集是有限的。利用节点分析原理，计算每一种工作制度对应的产液量、悬点载荷、泵效等，然后选择最优方案。根据Karamkas油田的特点，通过动态控制图分析，并考虑采液指数和含水率，建立抽油机井挖潜选井原则。最后选择该油田22口潜力井进行设计，现场实践证明，设计结果合理可行。 关键词 :抽油机井，节点分析，动态控制图，采液指数 Copyright © 2017 by authors， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 Karamkas油田位于哈萨克斯坦的曼格什套(Mangistau)州，该油田储层埋深浅(平均油层深度864 m)，储层物性好，供液能力强，采液指数变化范围大(最大可以达到1000 m 3 /(d∙MPa))，地层平均静压为9.2 MPa，平均饱和压力为6.5 MPa。目前以游梁式有杆泵为主要的举升方式，平均抽汲参数小，最大理论排液量为135 m 3 /d，无法充分发挥地层的供液能力。为了提高产液量，有必要对抽油机井系统进行挖潜设计。 有杆泵采油系统具有结构简单、适应性强、成本低等特点，目前是国内外应用最广泛的采油方式 [ 1 ] 。在以往的研究中，抽油机井系统设计大多以系统效率为目标，产液量为约束条件，以追求系统能耗最低 [ 2 ] [ 3 ] [ 4 ] [ 5 ] ；文献 [ 6 ] 综合考虑井下效率和地面系统合理性得到一个抽汲参数组合下的经济指标，并以此作为目标函数。有杆泵采油系统挖潜设计可以定义为通过调整工作制度，包括冲程、冲次、泵径、下泵深度等最大程度地提高抽油机井产量，同时要保证抽油机系统安全高效地运行。所以挖潜设计数学模型是一个以抽油机井产量最大为目标，以抽油机井的工作制度为设计参数的最优化模型。该模型是一个非线性优化模型，且产量方程无法用一个显式表达。对于这类方程一般采用智能算法求解，比如遗传算法 [ 7 ] 、群算法等。但由于模型解集是有限的、非连续的，采用枚举法求解速度更快，而且不会陷于局部最优。节点分析原理是抽油系统设计 [ 8 ] [ 9 ] 、参数敏感性分析 [ 10 ] 及优化 [ 11 ] 的重要理论基础。求解抽油机井挖潜数学模型同样要以生产系统的供排协调为基础，根据节点分析的原理进行求解。
system: [SYSTEM]

user: 2. 数学模型 在抽油机系统中，对一口固定的井，其产液量可以看成是冲程、冲次、泵径和下泵深度的函数，即： (1) 式中：Q l 为产液量m 3 /d；r p 为泵径，m；s为抽油机冲程，m；n为抽油机冲次，min −1 ；h p 为下泵深度，m。 其中冲程、冲次和泵径主要影响泵的理论排量。高冲次下，杆柱的工作状况恶化，特别容易造成抽油杆的损坏，同时高冲次容易引起抽油机振动，所以冲次不能过高。另外，更换抽油机的工作量大，在现场一般不采取更换抽油机的方式进行挖潜。 建立以产液量最大为目标的最优化模型，以目前现场设备所能提供的抽汲参数以及考虑安全经济生产作为约束条件： (2) 式中：S为抽油机能提供的冲程集合；N为抽油机能提供的冲次集合；R p 为现场提供的泵径集合；H为油层中深，m； 为泵效，%； 是允许最小泵效，%； 为应力范围比，%； 、 分别是抽油杆允许最小和最大应力范围比，%； 为抽油机载荷利用率，%； 为抽油机允许最大载荷利用率，%。
system: [SYSTEM]

user: 3. 求解方法 求解该最优化模型，最主要的是求解在某种工作制度下的产量，二者的关系无法用一个显式方程表达。可以采用节点分析的原理求解。 图1所示为抽油机井生产系统，以泵吸入口为求解点进行分析。流入过程为从油层到泵吸入口，流出过程为从泵吸入口到泵筒内。泵排出口以上不能和泵以下作为一个系统进行分析，主要是因为在该过程中，泵对流体做功，系统不满足能量守恒。 流入曲线函数可以用下式表示： (3) 式中：J为采液指数，m 3 /(d∙MPa)；p r 为地层平均压力，MPa；p in 为泵吸入口压力，MPa；Δp为泵挂深度到油层中深的油管流体压差(利用气液两相流压力降公式计算)， ，MPa；p wf 为井底流压，MPa。 流出过程反映的是泵吸入口压力与进泵液量之间的关系，可以按下式计算： 图1. 抽油机井生产系统 (4) 式中：Q t 为泵的理论排量，m 3 /d；f p 为柱塞截面积，m 2 。对于给定的工作制度，在同一坐标系里，根据上述计算方法绘制流入、流出曲线(图2)，其交点即为协调产量和对应的泵吸入口压力。 图2. 系统分析曲线 在该基础上，计算泵效、抽油机载荷利用率和抽油杆的应力范围比，根据约束条件，判断该工作制度是否满足要求。 冲程集合S、冲次集合N以及泵径集合R p 都是有限的几个点且不连续，下泵深度h p 也可以取有限的一系列非连续的点，因此模型的解集是有限的。对于该模型如果采用智能算法求解，不但工作量大，还会因为陷入局部最优而得不到最优解。笔者采用枚举法求解，即求出每一个解对应的产液量、载荷利用率、泵效以及抽油杆应力范围比，然后选择最优工作制度。
system: [SYSTEM]

user: 4. Karamkas油田挖潜设计 4.1. 挖潜选井原则 挖潜设计必须要选有潜力而且抽汲参数偏小的井。根据Karamkas油田实际情况，制定以下挖潜选井原则。 1) 选择供液能力强、目前抽汲参数偏小的油井。这类油井可以根据抽油机井动态控制图选择。 2) Karamkas油田普遍含水率高，当含水率大于95%时，则认为没有挖潜的必要。 3) 采液指数对产液量、泵效和载荷都有很大影响。该油田采液指数范围较大，通过计算分析可以得到不同泵径下采液指数和悬点载荷、泵效之间的相互关系，如图3、图4所示。从图3和图4可以看出，随着采液指数增加，悬点载荷逐渐减小，泵效逐渐增加，且变化幅度逐渐变小。当采液指数大于80 m 3 /(d∙MPa)时，载荷小、泵效高，满足约束条件，可以采用大泵快速抽汲，不需要进行设计。 根据以上原则，对含水率小于95%、采液指数小于80 m 3 /(d∙MPa)的53口油井进行动态控制分析，结果如图5所示。其中潜力井有22口。由于该油田地层压力高，部分井属于连喷带抽状态，校核泵效大于100%。 图3. 采液指数与悬点载荷关系曲线 图4. 采液指数与泵效关系曲线 图5. 抽油机井动态控制图 4.2. 挖潜设计 现场设备能提供的工作参数见表1。典型井基本参数见表2。假设一个工作制度，冲程为2.1 m，冲次为6 min − 1 ，下泵深度为800 m，以泵吸入口为求解点，绘制不同泵径下的系统分析图，从而得到对应的产液量和泵吸入口压力，然后计算抽油机载荷和抽油杆应力范围比，结果见图6和表3。 采用上述方法，对不同的工作制度进行求解，然后根据约束条件优选产液量最高的几种方案。约束条件为：泵效大于80%，应力范围比上限取90%，不考虑下限；载荷利用率小于80%，同时考虑到冲次过大会影响抽油机的稳定性，在选择方案的时候，冲次最大不超过7 min − 1 。挖潜方案设计结果见表4。 Table 1 The list of operating parameter 抽油机参数 抽油泵系列/mm 抽油杆系列/mm 额定载荷/kN 冲程/m 冲次/min -1 60 0.9、1.2、1.5、1.8、2.1 3.5、4、4.5、5、5.5、5、6.5、7、7.5 38、44、57、70、83、95 19、22 表1. 工作参数列表 Table 2 The basic parameters of a typical wel 原油相对密度/1 天然气相对密度/1 含水率/% 生产气油比/(m 3 ∙m −3 ) 饱和压力/MPa 下泵深度/m 0.9 0.85 85 27.5 7 800 表2. 典型井基本参数 图6. 典型井系统分析曲线(冲程2.1 m，冲次6 min −1 ，泵深800 m)
system: [SYSTEM]

user: 致密岩心带压渗吸实验及多尺度模拟研究 李帅 1，2 ，丁云宏 2 ，杨艳明 3 ，卢拥军 2 ，才博 2 1 中国石油勘探开发研究院，北京 2 中国石油勘探开发研究院廊坊分院，河北 廊坊 3 渤海钻探工程有限公司井下作业分公司，河北 任丘 作者简介：李帅(1987-)，男，博士生，主要从事储层改造与油藏数值模拟。 收稿日期：2016年4月8日；录用日期： 2016 年 9 月 30 日 ；发表日期：2017年4月15日
system: [SYSTEM]

user: 摘 要 致密储层孔喉半径小，毛细管压力高，渗吸动力强，压裂后若不立即放喷，可依靠压后焖井过程中的压差驱替和自发渗吸提高原油采出程度。为探索该过程机理，进行了带压渗吸实验和多尺度模拟研究：首先，将岩心单面开启，进行带压渗吸实验，模拟压裂完成后裂缝壁面在停泵压力下的渗吸行为。其次，建立基于CT扫描的孔喉尺度模型，分别获得驱替和渗吸的相对渗透率曲线和毛细管压力曲线。最后，以相渗曲线和毛管压力曲线为纽带建立矿场尺度模型，进行油藏尺度模拟。结果表明：①岩心带压渗吸条件下的采收率比自发渗吸提高约10%~15%；②通过调整孔喉尺度模型的微观孔隙参数进行渗吸采收率的拟合，可分别计算得到驱替和渗吸的相对渗透率及毛细管压力曲线；③油藏尺度下对基质赋予渗吸属性、对裂缝赋予驱替属性，可对开采初期的产量变化有较好描述。 关键词 :致密岩心，体积压裂，自发渗吸，强制渗吸，孔隙尺度模型 Copyright © 2017 by authors， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/