CN104747163B - 一种在致密砂岩中识别储层裂缝的方法及装置 - Google Patents
一种在致密砂岩中识别储层裂缝的方法及装置 Download PDFInfo
- Publication number
- CN104747163B CN104747163B CN201310750713.5A CN201310750713A CN104747163B CN 104747163 B CN104747163 B CN 104747163B CN 201310750713 A CN201310750713 A CN 201310750713A CN 104747163 B CN104747163 B CN 104747163B
- Authority
- CN
- China
- Prior art keywords
- log
- normalization
- weight coefficient
- lateral resistivity
- well
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 206010017076 Fracture Diseases 0.000 claims abstract description 82
- 208000010392 Bone Fractures Diseases 0.000 claims abstract description 75
- 238000011161 development Methods 0.000 claims abstract description 55
- 239000002343 natural gas well Substances 0.000 claims abstract description 46
- 238000010606 normalization Methods 0.000 claims description 43
- 239000004576 sand Substances 0.000 claims description 33
- 238000005070 sampling Methods 0.000 claims description 18
- 230000005611 electricity Effects 0.000 claims description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 42
- 239000003345 natural gas Substances 0.000 abstract description 21
- 239000007789 gas Substances 0.000 abstract description 4
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 239000011435 rock Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种在致密砂岩中识别储层裂缝的方法及装置，属于天然气开发领域。所述方法包括：采集在致密砂岩中钻开的天然气井的测井数据，所述测井数据包括声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率；获取每种测井数据对应的测井曲线，并计算每条测井曲线的权值系数；根据所述每条测井曲线和所述每条测井曲线的权值系数，计算所述天然气井包括的每个井段的裂缝发育概率；如果存在裂缝发育概率大于或等于预设阈值的井段，则确定所述井段的储层存在裂缝。所述装置包括：采集模块、第一计算模块、第二计算模块和确定模块。本发明提高了识别储层裂缝的准确率，进而提高在该井段处开采天然气的产量。
Description
技术领域
本发明涉及天然气开发领域，特别涉及一种在致密砂岩中识别储层裂缝的方法及装置。
背景技术
在天然气开发中，钻开一个天然气井后，技术人员可以识别出开采天然气的储层在该天然气井中的位置，在该位置处射开储层来开采天然气。
目前，发现两种砂岩中都存在储藏天然气的储层，该两种砂岩为普通砂岩和致密砂岩。对于普通砂岩，技术人员可以通过常规测井解释方式识别出开采天然气的储层。而对于致密砂岩，由于致密砂岩与普通砂岩不同，致密砂岩的储层孔隙度小于5%，渗透率小于0.2mD，储层的渗流能力差，从致密砂岩的储层中开采天然气的产量会特别低；但是如果在致密砂岩的储层裂缝发育段开采天然气，可以大大提高天然气的产量，所以现有技术从致密砂岩中开采天然气之前，技术人员首先通过常规测井解释方法从该致密砂岩中识别储层裂缝。
由于通过常规测井解释方法在致密砂岩中识别储层裂缝的能力有限，一旦识别错误，射开储层后发现天然气的产量较低或者不产天然气，技术人员只能采取酸化、压裂和改层等措施进行补救，从而导致大量物力和财力的浪费，增加了生产成本。
发明内容
为了解决现有技术的问题，本发明实施例提供了一种在致密砂岩中识别储层裂缝的方法及装置。所述技术方案如下：
一方面，提供了一种在致密砂岩中识别储层裂缝的方法，所述方法包括：
采集在致密砂岩中钻开的天然气井的测井数据，所述测井数据包括声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率；
获取每种测井数据对应的测井曲线，并计算每条测井曲线的权值系数；
根据所述每条测井曲线和所述每条测井曲线的权值系数，计算所述天然气井包括的每个井段的裂缝发育概率；
如果存在裂缝发育概率大于或等于预设阈值的井段，则确定所述井段的储层存在裂缝。
其中，所述计算每条测井曲线的权值系数，包括：
根据预设的井段位置，对每条测井曲线分别进行采样，得到多个采样点；
根据所述多个采样点的测井数据，计算所述每条测井曲线的权值系数。
进一步地，所述根据预设的井段位置，对每条测井曲线分别进行采样，得到多个采样点，包括：
对每条测井曲线分别进行归一化处理，得到所述每条测井曲线的归一化测井曲线；
根据预设的井段位置，对所述每条测井曲线的归一化测井曲线分别进行采样，得到多个采样点。
进一步地，所述根据所述多个采样点的测井数据，计算所述每条测井曲线的权值系数，包括：
从所述多个采样点中，获取位于同一井段位置的采样点；
根据所述获取的采样点的测井数据，分别计算所述获取的每个采样点的权值系数；
根据所述每条测井曲线包括的每个采样点的权值系数，分别计算所述每条测井曲线的权值系数。
其中，所述根据所述每条测井曲线和所述每条测井曲线的权值系数，计算所述天然气井包括的每个井段的裂缝发育概率，包括：
根据所述每条测井曲线和所述每条测井曲线的权值系数，按照如下公式（1）计算所述天然气井包括的每个井段的裂缝发育概率；
其中，在上述公式（1）中，G为所述天然气井包括的每个井段的裂缝发育概率，ai为第i条测井曲线的校正值，m为测井曲线的条数，Wi为所述第i条测井曲线的权值系数，Xi为所述第i条测井曲线的测井数据。
另一方面，提供了一种在致密砂岩中识别储层裂缝的装置，所述装置包括：
采集模块，用于采集在致密砂岩中钻开的天然气井的测井数据，所述测井数据包括声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率；
第一计算模块，用于获取每种测井数据对应的测井曲线，并计算每条测井曲线的权值系数；
第二计算模块，用于根据所述每条测井曲线和所述每条测井曲线的权值系数，计算所述天然气井包括的每个井段的裂缝发育概率；
确定模块，用于如果存在裂缝发育概率大于或等于预设阈值的井段，则确定所述井段的储层存在裂缝。
其中，所述第一计算模块包括：
获取单元，用于获取每种测井数据对应的测井曲线；
采样单元，用于根据预设的井段位置，对每条测井曲线分别进行采样，得到多个采样点；
第一计算单元，用于根据所述多个采样点的测井数据，计算所述每条测井曲线的权值系数。
进一步地，所述采样单元包括：
归一化子单元，用于对每条测井曲线分别进行归一化处理，得到所述每条测井曲线的归一化测井曲线；
采样子单元，用于根据预设的井段位置，对所述每条测井曲线的归一化测井曲线分别进行采样，得到多个采样点。
进一步地，所述第一计算单元包括：
获取子单元，用于从所述多个采样点中，获取位于同一井段位置的采样点；
第一计算子单元，用于根据所述获取的采样点的测井数据，分别计算所述获取的每个采样点的权值系数；
第二计算子单元，用于根据所述每条测井曲线包括的每个采样点的权值系数，分别计算所述每条测井曲线的权值系数。
其中，所述第二计算模块包括：
第二计算单元，用于根据所述每条测井曲线和所述每条测井曲线的权值系数，按照如下公式（1）计算所述天然气井包括的每个井段的裂缝发育概率；
其中，在上述公式（1）中，G为所述天然气井包括的每个井段的裂缝发育概率，ai为第i条测井曲线的校正值，m为测井曲线的条数，Wi为所述第i条测井曲线的权值系数，Xi为所述第i条测井曲线的测井数据。
在本发明实施例中，计算声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率的测井曲线的权值系数，根据每条测井曲线和每条测井曲线的权值系数，计算天然气井包括的每个井段的裂缝发育概率；如果该天然气井包括的井段中存在裂缝发育概率大于或等于预设阈值的井段，则确定该井段的储层存在裂缝。由于根据声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率等多个测井曲线识别储层裂缝，提高了识别储层裂缝的准确率，进而提高在该井段处开采天然气的产量。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例一提供的一种在致密砂岩中识别储层裂缝的方法流程图；
图2是本发明实施例二提供的一种在致密砂岩中识别储层裂缝的方法流程图；
图3是本发明实施例三提供的一种在致密砂岩中识别储层裂缝的装置结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
实施例一
本发明实施例提供了一种在致密砂岩中识别储层裂缝的方法，参见图1，该方法包括：
步骤101：采集在致密砂岩中钻开的天然气井的测井数据，该测井数据包括声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率；
步骤102：获取每种测井数据对应的测井曲线，并计算每条测井曲线的权值系数；
步骤103：根据每条测井曲线和每条测井曲线的权值系数，计算该天然气井包括的每个井段的裂缝发育概率；
步骤104：如果存在裂缝发育概率大于或等于预设阈值的井段，则确定该井段的储层存在裂缝。
其中，计算每条测井曲线的权值系数，包括：
根据预设的井段位置，对每条测井曲线分别进行采样，得到多个采样点；
根据多个采样点的测井数据，计算每条测井曲线的权值系数。
进一步地，根据预设的井段位置，对每条测井曲线分别进行采样，得到多个采样点，包括：
对每条测井曲线分别进行归一化处理，得到每条测井曲线的归一化测井曲线；
根据预设的井段位置，对每条测井曲线的归一化测井曲线分别进行采样，得到多个采样点。
进一步地，根据多个采样点的测井数据，计算每条测井曲线的权值系数，包括：
从多个采样点中，获取位于同一井段位置的采样点；
根据获取的采样点的测井数据，分别计算获取的每个采样点的权值系数；
根据每条测井曲线包括的每个采样点的权值系数，分别计算每条测井曲线的权值系数。
其中，根据每条测井曲线和每条测井曲线的权值系数，计算天然气井包括的每个井段的裂缝发育概率，包括：
根据每条测井曲线和每条测井曲线的权值系数，按照如下公式（1）计算该天然气井包括的每个井段的裂缝发育概率；
其中，在上述公式（1）中，G为该天然气井包括的每个井段的裂缝发育概率，ai为第i条测井曲线的校正值，m为测井曲线的条数，Wi为第i条测井曲线的权值系数，Xi为第i条测井曲线的测井数据。
在本发明实施例中，计算声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率的测井曲线的权值系数，根据每条测井曲线和每条测井曲线的权值系数，计算天然气井包括的每个井段的裂缝发育概率；如果该天然气井包括的井段中存在裂缝发育概率大于或等于预设阈值的井段，则确定该井段的储层存在裂缝。由于根据声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率等多个测井曲线识别储层裂缝，提高了识别储层裂缝的准确率，进而提高在该井段处开采天然气的产量。
实施例二
本发明实施例提供了一种在致密砂岩中识别储层裂缝的方法，参见图2，该方法包括：
步骤201：采集在致密砂岩中钻开的天然气井的测井数据，该测井数据包括声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率；
具体地，采用测井仪器对致密砂岩中钻开的天然气井进行数据采集，获取该天然气井的测井数据，该测井数据包括声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率。
其中，声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率分别对应一个测井曲线，且测井曲线是测井数据与井段的对应关系。
可选地，由于采集天然气井的测井数据的测井仪器是通过电缆，将该测井仪器从该天然气井的井底匀速地向井口移动，所以测井曲线还可以是测井数据与时间的对应关系。
优选地，当测井仪器采集的测井数据不准确时，还可以通过自然伽马对测井仪器采集的测井数据进行补充。
步骤202：获取每种测井数据对应的测井曲线，并根据预设的井段位置，对每条测井曲线分别进行采样，得到多个采样点；
具体地，对每条测井曲线分别进行归一化处理，得到每条测井曲线的归一化测井曲线；根据预设的井段位置，对每条测井曲线的归一化测井曲线分别进行采样，得到多个采样点。
需要补充说明的是，由于每条测井曲线对不同的地质因素、自身的响应机理、曲线的量纲和采集装置不同等因素存在差异，所以根据每条测井曲线进行储层裂缝识别时，需要对每条测井曲线进行归一化处理。
其中，对声波和补偿中子的测井曲线进行归一化处理的具体操作可以为：根据声波的测井曲线，按照如下公式（1）计算归一化后的声波的测井曲线，以及根据补偿中子的测井曲线，按照如下公式（1）计算归一化后的补偿中子的测井曲线；
其中，在上述公式（1）中，X'为归一化后的测井曲线中的测井数据，X为归一化之前的测井曲线中的测井数据，Xmin为储层裂缝不发育层段的测井数据，Xmax为储层裂缝发育层段的测井数据。
其中，对密度、深侧向电阻率、浅侧向电阻率和微侧向电阻率的测井曲线进行归一化处理的具体操作可以为：根据密度的测井曲线，按照如下公式（2）计算归一化后的密度的测井曲线；根据深侧向电阻率的测井数据，按照如下公式（2）计算归一化后的深侧向电阻率的测井曲线；根据浅侧向电阻率的测井数据，按照如下公式（2）计算归一化后的浅侧向电阻率的测井曲线；根据微侧向电阻率的测井数据，按照如下公式（2）计算归一化后的微侧向电阻率的测井曲线；
其中，在上述公式（2）中，Y'为归一化后的测井曲线中的测井数据，Y为归一化之前的测井曲线中的测井数据，Ymin为储层裂缝不发育层段的测井数据，Ymax为储层裂缝发育层段的测井数据。
由于裂缝发育的地方，密度、深侧向电阻率、浅侧向电阻率和微侧向电阻率都是向偏小的方向发展，所以对密度、深侧向电阻率、浅侧向电阻率和微侧向电阻率的测井曲线进行归一化的公式和对声波和补偿中子的测井曲线进行归一化的公式不同。
步骤203：根据多个采样点的测井数据，计算每条测井曲线的权值系数；
具体地，从采样的多个采样点中，获取位于同一井段位置的采样点；根据获取的采样点的测井数据，分别计算获取的每个采样点的权值系数；根据每条测井曲线包括的每个采样点的权值系数，分别计算每条测井曲线的权值系数。
其中，根据获取的采样点的测井数据，分别计算获取的每个采样点的权值系数的具体操作可以为：根据获取的采样点的测井数据，按照如下公式（3）计算获取的每个采样点的权值系数；
其中，在上述公式（3）中，wij为第i条测井曲线的第j个采样点的权值系数，xij为第i条测井曲线的第j个采样点的测井数据，m为测井曲线的条数。
其中，计算每条测井曲线的权值系数之前，可以对每条测井曲线进行编号，以及对每条测井曲线上的采样点进行编号，且每条测井曲线中的同一井段位置的采样点的编号相同，在上述公式（3）中每条测井曲线的第j个采样点都是同一井段位置的采样点。
其中，根据每条测井曲线包括的每个采样点的权值系数，分别计算每条测井曲线的权值系数的具体操作可以为：根据每条测井曲线包括的每个采样点的权值系数，按照如下公式（4）计算每条测井曲线的权值系数；
其中，在上述公式（4）中，Wi为第i条测井曲线的权值系数，ni为第i条测井曲线与裂缝响应的符合率。
步骤204：根据每条测井曲线和每条测井曲线的权值系数，计算该天然气井包括的每个井段的裂缝发育概率；
具体地，根据每条测井曲线和每条测井曲线的权值系数，按照如下公式（5）计算该天然气井包括的每个井段的裂缝发育概率；
其中，在上述公式（5）中，G为该天然气井包括的每个井段的裂缝发育概率，ai为第i条测井曲线的校正值，m为测井曲线的条数，Wi为第i条测井曲线的权值系数，Xi为第i条测井曲线的测井数据。
步骤205：如果存在裂缝发育概率大于或等于预设阈值的井段，则确定该井段的储层存在裂缝。
具体地，将计算的每个井段的裂缝发育概率与预设阈值进行比较，如果存在裂缝发育概率大于或等于预设阈值的井段，则确定该井段的储层存在裂缝。
进一步地，确定储层存在裂缝的井段之后，可以射开该井段的储层来开采天然气。
其中，如果该天然气井包括的井段中不存在裂缝发育概率大于或等于预设阈值的井段，则确定该天然气井的储层中不存在裂缝，所以不从该天然气井中开采天然气。
进一步地，可以事先从在致密砂岩中钻开的天然气井中取出一块岩心，观察取出的岩心的裂缝发育概率，当确定该井段的储层存在裂缝之后，还可以将岩心的裂缝发育概率与该井段的裂缝发育概率进行比较，如果两者之间的差值小于预设数值时，可以确定该井段的储层确实存在裂缝。
在本发明实施例中，计算声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率的测井曲线的权值系数，根据每条测井曲线和每条测井曲线的权值系数，计算天然气井包括的每个井段的裂缝发育概率；如果该天然气井包括的井段中存在裂缝发育概率大于或等于预设阈值的井段，则确定该井段的储层存在裂缝。由于根据声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率等多个测井曲线识别储层裂缝，提高了识别储层裂缝的准确率，进而提高在该井段处开采天然气的产量。
实施例三
参见图3，本发明实施例提供了一种在致密砂岩中识别储层裂缝的装置，该装置包括：
采集模块301，用于采集在致密砂岩中钻开的天然气井的测井数据，该测井数据包括声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率；
第一计算模块302，用于获取每种测井数据对应的测井曲线，并计算每条测井曲线的权值系数；
第二计算模块303，用于根据每条测井曲线和每条测井曲线的权值系数，计算该天然气井包括的每个井段的裂缝发育概率；
确定模块304，用于如果存在裂缝发育概率大于或等于预设阈值的井段，则确定该井段的储层存在裂缝。
其中，第一计算模块302包括：
获取单元，用于获取每种测井数据对应的测井曲线；
采样单元，用于根据预设的井段位置，对每条测井曲线分别进行采样，得到多个采样点；
第一计算单元，用于根据多个采样点的测井数据，计算每条测井曲线的权值系数。
进一步地，采样单元包括：
归一化子单元，用于对每条测井曲线分别进行归一化处理，得到每条测井曲线的归一化测井曲线；
采样子单元，用于根据预设的井段位置，对每条测井曲线的归一化测井曲线分别进行采样，得到多个采样点。
进一步地，第一计算单元包括：
获取子单元，用于从多个采样点中，获取位于同一井段位置的采样点；
第一计算子单元，用于根据获取的采样点的测井数据，分别计算获取的每个采样点的权值系数；
第二计算子单元，用于根据每条测井曲线包括的每个采样点的权值系数，分别计算每条测井曲线的权值系数。
其中，第二计算模块303包括：
第二计算单元，用于根据每条测井曲线和每条测井曲线的权值系数，按照如下公式（1）计算天然气井包括的每个井段的裂缝发育概率；
其中，在上述公式（1）中，G为该天然气井包括的每个井段的裂缝发育概率，ai为第i条测井曲线的校正值，m为测井曲线的条数，Wi为第i条测井曲线的权值系数，Xi为第i条测井曲线的测井数据。
在本发明实施例中，计算声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率的测井曲线的权值系数，根据每条测井曲线和每条测井曲线的权值系数，计算天然气井包括的每个井段的裂缝发育概率；如果该天然气井包括的井段中存在裂缝发育概率大于或等于预设阈值的井段，则确定该井段的储层存在裂缝。由于根据声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率等多个测井曲线识别储层裂缝，提高了识别储层裂缝的准确率，进而提高在该井段处开采天然气的产量。
需要说明的是：上述实施例提供的在致密砂岩中识别储层裂缝的装置在致密砂岩中识别储层裂缝时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的在致密砂岩中识别储层裂缝的装置与在致密砂岩中识别储层裂缝的方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种在致密砂岩中识别储层裂缝的方法，其特征在于，所述方法包括：
采集在致密砂岩中钻开的天然气井的测井数据，所述测井数据包括声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率；
获取每种测井数据对应的测井曲线，对每条测井曲线分别进行归一化处理，得到所述每条测井曲线的归一化测井曲线，并计算每条测井曲线的权值系数；
根据所述每条测井曲线和所述每条测井曲线的权值系数，计算所述天然气井包括的每个井段的裂缝发育概率；
如果存在裂缝发育概率大于或等于预设阈值的井段，则确定所述井段的储层存在裂缝，
其中，对声波和补偿中子的测井曲线进行归一化处理的操作为：根据声波的测井曲线，按照公式计算归一化后的声波的测井曲线，以及根据补偿中子的测井曲线，按照公式计算归一化后的补偿中子的测井曲线；
其中，X'为归一化后的测井曲线中的测井数据，X为归一化之前的测井曲线中的测井数据，Xmin为储层裂缝不发育层段的测井数据，Xmax为储层裂缝发育层段的测井数据；
其中，对密度、深侧向电阻率、浅侧向电阻率和微侧向电阻率的测井曲线进行归一化处理的操作为：根据密度的测井曲线，按照公式计算归一化后的密度的测井曲线；根据深侧向电阻率的测井数据，按照公式 计算归一化后的深侧向电阻率的测井曲线；根据浅侧向电阻率的测井数据，按照公式计算归一化后的浅侧向电阻率的测井曲线；根据微侧向电阻率的测井数据，按照公式计算归一化后的微侧向电阻率的测井曲线；
其中，Y'为归一化后的测井曲线中的测井数据，Y为归一化之前的测井曲线中的测井数据，Ymin为储层裂缝不发育层段的测井数据，Ymax为储层裂缝发育层段的测井数据。
2.如权利要求1所述的方法，其特征在于，所述计算每条测井曲线的权值系数，包括：
根据预设的井段位置，对每条测井曲线分别进行采样，得到多个采样点；
根据所述多个采样点的测井数据，计算所述每条测井曲线的权值系数。
3.如权利要求2所述的方法，其特征在于，所述根据预设的井段位置，对每条测井曲线分别进行采样，得到多个采样点，包括：
根据预设的井段位置，对所述每条测井曲线的归一化测井曲线分别进行采样，得到多个采样点。
4.如权利要求2所述的方法，其特征在于，所述根据所述多个采样点的测井数据，计算所述每条测井曲线的权值系数，包括：
从所述多个采样点中，获取位于同一井段位置的采样点；
根据所述获取的采样点的测井数据，分别计算所述获取的每个采样点的权值系数；
根据所述每条测井曲线包括的每个采样点的权值系数，分别计算所述每条测井曲线的权值系数。
5.如权利要求1所述的方法，其特征在于，所述根据所述每条测井曲线和所述每条测井曲线的权值系数，计算所述天然气井包括的每个井段的裂缝发育概率，包括：
根据所述每条测井曲线和所述每条测井曲线的权值系数，按照如下公式(1)计算所述天然气井包括的每个井段的裂缝发育概率；
其中，在上述公式(1)中，G为所述天然气井包括的每个井段的裂缝发育概率，ai为第i条测井曲线的校正值，m为测井曲线的条数，Wi为所述第i条测 井曲线的权值系数，Xi为所述第i条测井曲线的测井数据。
6.一种在致密砂岩中识别储层裂缝的装置，其特征在于，所述装置包括：
采集模块，用于采集在致密砂岩中钻开的天然气井的测井数据，所述测井数据包括声波、密度、补偿中子、深侧向电阻率、浅侧向电阻率和微侧向电阻率；
第一计算模块，用于获取每种测井数据对应的测井曲线，用于对每条测井曲线分别进行归一化处理，得到所述每条测井曲线的归一化测井曲线，并计算每条测井曲线的权值系数；
第二计算模块，用于根据所述每条测井曲线和所述每条测井曲线的权值系数，计算所述天然气井包括的每个井段的裂缝发育概率；
确定模块，用于如果存在裂缝发育概率大于或等于预设阈值的井段，则确定所述井段的储层存在裂缝；
其中，对声波和补偿中子的测井曲线进行归一化处理的操作为：根据声波的测井曲线，按照公式计算归一化后的声波的测井曲线，以及根据补偿中子的测井曲线，按照公式计算归一化后的补偿中子的测井曲线；
其中，X'为归一化后的测井曲线中的测井数据，X为归一化之前的测井曲线中的测井数据，Xmin为储层裂缝不发育层段的测井数据，Xmax为储层裂缝发育层段的测井数据；
其中，对密度、深侧向电阻率、浅侧向电阻率和微侧向电阻率的测井曲线进行归一化处理的操作为：根据密度的测井曲线，按照公式计算归一化后的密度的测井曲线；根据深侧向电阻率的测井数据，按照公式 计算归一化后的深侧向电阻率的测井曲线；根据浅侧向电阻率的测井数据，按照公式计算归一化后的浅侧向电阻率的测井曲线； 根据微侧向电阻率的测井数据，按照公式计算归一化后的微侧向电阻率的测井曲线；
其中，Y'为归一化后的测井曲线中的测井数据，Y为归一化之前的测井曲线中的测井数据，Ymin为储层裂缝不发育层段的测井数据，Ymax为储层裂缝发育层段的测井数据。
7.如权利要求6所述的装置，其特征在于，所述第一计算模块包括：
获取单元，用于获取每种测井数据对应的测井曲线；
采样单元，用于根据预设的井段位置，对每条测井曲线分别进行采样，得到多个采样点；
第一计算单元，用于根据所述多个采样点的测井数据，计算所述每条测井曲线的权值系数。
8.如权利要求7所述的装置，其特征在于，所述采样单元包括：
采样子单元，用于根据预设的井段位置，对所述每条测井曲线的归一化测井曲线分别进行采样，得到多个采样点。
9.如权利要求7所述的装置，其特征在于，所述第一计算单元包括：
获取子单元，用于从所述多个采样点中，获取位于同一井段位置的采样点；
第一计算子单元，用于根据所述获取的采样点的测井数据，分别计算所述获取的每个采样点的权值系数；
第二计算子单元，用于根据所述每条测井曲线包括的每个采样点的权值系数，分别计算所述每条测井曲线的权值系数。
10.如权利要求6所述的装置，其特征在于，所述第二计算模块包括：
第二计算单元，用于根据所述每条测井曲线和所述每条测井曲线的权值系数，按照如下公式(1)计算所述天然气井包括的每个井段的裂缝发育概率；
其中，在上述公式(1)中，G为所述天然气井包括的每个井段的裂缝发育 概率，ai为第i条测井曲线的校正值，m为测井曲线的条数，Wi为所述第i条测井曲线的权值系数，Xi为所述第i条测井曲线的测井数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310750713.5A CN104747163B (zh) | 2013-12-31 | 2013-12-31 | 一种在致密砂岩中识别储层裂缝的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310750713.5A CN104747163B (zh) | 2013-12-31 | 2013-12-31 | 一种在致密砂岩中识别储层裂缝的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104747163A CN104747163A (zh) | 2015-07-01 |
| CN104747163B true CN104747163B (zh) | 2017-10-17 |
Family
ID=53587342
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310750713.5A Active CN104747163B (zh) | 2013-12-31 | 2013-12-31 | 一种在致密砂岩中识别储层裂缝的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104747163B (zh) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106094052A (zh) * | 2016-06-01 | 2016-11-09 | 中国地质大学（武汉） | 一种致密白云岩储层的裂缝发育程度识别方法 |
| CN109254330B (zh) * | 2017-07-12 | 2020-11-03 | 中国石油天然气股份有限公司 | 火成岩储层的裂缝段的识别方法及装置 |
| CN109425912A (zh) * | 2017-09-05 | 2019-03-05 | 中国石油化工股份有限公司 | 一种碳酸盐岩油藏储层有效裂缝的识别方法 |
| CN107956472B (zh) * | 2017-11-10 | 2021-01-01 | 中国石油天然气股份有限公司 | 裂缝性碳酸盐岩储集层完井方式确定方法、装置及介质 |
| CN110847898B (zh) * | 2018-07-24 | 2023-05-16 | 中国石油化工股份有限公司 | 一种建立双介质储层分类标准的方法 |
| CN110939428B (zh) * | 2018-09-21 | 2023-08-01 | 中国石油化工股份有限公司 | 一种致密砂岩油气藏储层裂缝的识别方法 |
| CN109614584B (zh) * | 2018-11-16 | 2020-10-27 | 中国科学院计算技术研究所 | 一种资源测井数据的重建方法 |
| CN109870720B (zh) * | 2019-01-25 | 2021-01-01 | 中国石油天然气集团有限公司 | 一种页岩气微裂缝测井识别方法 |
| CN111561266B (zh) * | 2019-01-29 | 2022-05-17 | 中国石油化工股份有限公司 | 一种地层裂缝识别方法及系统 |
| CN112709566B (zh) * | 2019-10-25 | 2024-03-19 | 中国石油化工股份有限公司 | 一种碳酸盐岩气藏产层的预测方法 |
| CN110847887B (zh) * | 2019-10-28 | 2023-01-10 | 中国石油天然气股份有限公司大港油田分公司 | 一种细粒沉积陆相页岩裂缝识别评价方法 |
| CN114592857A (zh) * | 2020-12-07 | 2022-06-07 | 中国石油天然气股份有限公司 | 测井曲线填充处理方法及装置 |
| CN115061219B (zh) * | 2022-08-17 | 2022-11-08 | 北京派特杰奥科技有限公司 | 基于石油天然气探测的裂缝型储层预测识别方法及系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1040245A (zh) * | 1988-07-25 | 1990-03-07 | 石油工业部石油勘探开发科学研究院石油地质研究所 | 一种测井信息的解释处理方法 |
| CN102454399A (zh) * | 2010-10-26 | 2012-05-16 | 中国石油化工股份有限公司 | 测井声波时差信号校正方法 |
| CN103422852A (zh) * | 2012-05-24 | 2013-12-04 | 中国石油化工股份有限公司 | 一种不同井间气测值转换对比方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120199345A1 (en) * | 2011-02-04 | 2012-08-09 | Mark Anderson | Unconventional Gas Fracture Logging Method and Apparatus |
-
2013
- 2013-12-31 CN CN201310750713.5A patent/CN104747163B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1040245A (zh) * | 1988-07-25 | 1990-03-07 | 石油工业部石油勘探开发科学研究院石油地质研究所 | 一种测井信息的解释处理方法 |
| CN102454399A (zh) * | 2010-10-26 | 2012-05-16 | 中国石油化工股份有限公司 | 测井声波时差信号校正方法 |
| CN103422852A (zh) * | 2012-05-24 | 2013-12-04 | 中国石油化工股份有限公司 | 一种不同井间气测值转换对比方法 |
Non-Patent Citations (1)
| Title |
|---|
| 碳酸盐岩测井解释方法在塔里木北部地区储量计算中的应用;魏历灵等;《石油实验地质》;19951231;第17卷(第4期);367-371 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104747163A (zh) | 2015-07-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104747163B (zh) | 一种在致密砂岩中识别储层裂缝的方法及装置 | |
| CN103592690B (zh) | 基于电成像测井孔隙度谱信息自动识别储层裂缝的方法 | |
| CN104360415B (zh) | 一种致密砂岩储层裂缝识别的方法 | |
| CN104948176B (zh) | 一种基于渗透增大率识别碳酸盐岩储层裂缝的方法 | |
| CN109085663A (zh) | 一种致密砂岩储层层理缝识别方法 | |
| CN104977617B (zh) | 储层裂缝识别方法和成像测井储层裂缝识别方法 | |
| US9176255B2 (en) | Permeability prediction systems and methods using quadratic discriminant analysis | |
| CN106353792A (zh) | 一种适用于水力压裂微震震源定位的方法 | |
| CN105986815A (zh) | 一种用于识别页岩地层地质甜点的方法 | |
| US20190113639A1 (en) | Method and Device for Estimating Sonic Slowness in A Subterranean Formation | |
| CN104615843B (zh) | 一种识别隐伏火山岩型铀矿化信息的物化探综合方法 | |
| CN103649781A (zh) | 方位角脆性测井系统和方法 | |
| CN105863628B (zh) | 一种油田开发阶段储层精细预测方法 | |
| CN103485758A (zh) | 一种利用多种资料拟合反演复杂岩性解释方法及装置 | |
| CN104614764B (zh) | 一种基于扰动弹性阻抗反演的裂缝储层流体识别方法 | |
| CN103590828B (zh) | 一种录井dck指数法评价地层压力的方法 | |
| CN106979006A (zh) | 地层压力的确定方法和装置 | |
| CN108828661B (zh) | 基于地震脉冲响应谱测定场地卓越周期的方法 | |
| CN104678455A (zh) | 一种陆相缝洞型储层识别方法 | |
| CN104182651A (zh) | 用于三分量检波器接收的微地震事件方位角自动质控方法 | |
| CN104181610B (zh) | 一种射线路径弹性反演方法以及系统 | |
| CN104314558B (zh) | 用斯通利波能量损失度判别储层有效性的方法 | |
| CN116413831A (zh) | 多井合成地震记录自动标定及构造解释和储层预测方法 | |
| CN107330103A (zh) | 一种岩土工程钻探信息采集及处理方法和系统 | |
| CN109598049B (zh) | 钻孔岩石裂隙发育程度及区域岩石裂隙发育规律的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |