技术领域Technical Field
本发明涉及油田开发技术领域，特别是涉及到一种均衡平面驱替的确定径向水射流分支长度的优化设计方法。The invention relates to the technical field of oilfield development, and in particular to an optimization design method for determining radial water jet branch lengths in a balanced plane displacement.
背景技术Background Art
径向水射流钻孔技术可以连通地层到井筒的通道，穿透近井筒地带的污染区，提高地层的渗透率，减小井筒周围地带的渗流阻力，最终使生产井产量增加、注水井注水量增加。Radial water jet drilling technology can connect the channel from the formation to the wellbore, penetrate the contaminated area near the wellbore, improve the permeability of the formation, reduce the seepage resistance in the area around the wellbore, and ultimately increase the output of production wells and the water injection volume of water injection wells.
低渗透油田中有一些井的产量很低，采用压裂酸化等常规措施的效果不明显，治理潜力比较小。而径向水射流钻孔技术可以在地层定向钻孔，提高地层渗透性，降低井筒周围区域的渗流阻力，从而能够提高油井产量和注水井的注水量，所以研究径向水射流钻孔技术与井网适配优化具有重要的意义。李坤采用数值模拟技术研究径向水射流的渗流特征，在此基础上运用Green函数、Newman乘积、势叠加原理推导出单井情况下的产能公式，并分析增产增注机理。研究地质因素、开发因素、径向水射流设计参数等对开发效果的影响因素及政策界限。确定技术极限井距和经济极限井距，并根据径向水射流与井网形式的匹配模式，优化确定了适配井网。Some wells in low permeability oil fields have very low production, and the effects of conventional measures such as fracturing and acidizing are not obvious, and the potential for governance is relatively small. However, radial water jet drilling technology can directional drill holes in the formation, improve the permeability of the formation, and reduce the seepage resistance in the area around the wellbore, thereby increasing the production of oil wells and the injection volume of water injection wells. Therefore, it is of great significance to study the adaptation and optimization of radial water jet drilling technology and well network. Li Kun used numerical simulation technology to study the seepage characteristics of radial water jets. On this basis, he used Green function, Newman product, and potential superposition principle to derive the production capacity formula under the condition of a single well, and analyzed the mechanism of increased production and injection. The influencing factors and policy boundaries of geological factors, development factors, radial water jet design parameters, etc. on the development effect were studied. The technical limit well spacing and economic limit well spacing were determined, and the adaptive well network was optimized and determined according to the matching mode of radial water jets and well network forms.
通过研究发现，径向钻孔可以降低压降损失，变径向流为线性流，增加了地层有效渗透率和井网动用程度。采用五点法井网形式开发，钻孔方向与井排平行；采用九点法井网形式，钻孔方向与井排夹角45°，研究成果为今后径向水射流在低渗透油藏的推广应用给予了可靠的技术支持。Through research, it is found that radial drilling can reduce pressure drop loss, change radial flow into linear flow, and increase the effective permeability of the formation and the utilization of the well pattern. The five-point well pattern is used for development, and the drilling direction is parallel to the well pattern; the nine-point well pattern is used, and the drilling direction is at an angle of 45° to the well pattern. The research results provide reliable technical support for the promotion and application of radial water jets in low permeability reservoirs in the future.
胜利油区低渗透油藏资源丰富，但是油藏品质差，平面非均质性强，导致平面驱替不均衡现象严重，针对多油层油藏平面驱替不均衡的现象，现场主要应用径向水射流来进行流线调配，但既没有流线分层显示的方法，也没有对径向水射流分支长度的优化设计展开更为系统的研究。The Shengli Oilfield is rich in low-permeability reservoir resources, but the reservoir quality is poor and the planar heterogeneity is strong, resulting in serious unbalanced planar displacement. To address the unbalanced planar displacement of multi-layer reservoirs, radial water jets are mainly used on site to adjust streamlines, but there is no method for displaying streamline stratification, nor is there a more systematic study on the optimal design of radial water jet branch lengths.
发明内容Summary of the invention
本发明的目的是提供一种流线分层显示的方法，以此为基础下开展均衡平面驱替的径向水射流分支长度的优化设计方法。The purpose of the present invention is to provide a method for displaying streamline layers, based on which a method for optimizing the design of radial water jet branch lengths for balanced plane displacement is developed.
为实现上述目的，本发明采取以下技术方案：一种均衡平面驱替的确定径向水射流分支长度的优化设计方法，包括以下步骤：To achieve the above object, the present invention adopts the following technical solution: an optimization design method for determining the radial water jet branch length for balanced plane displacement, comprising the following steps:
1)将工区的精细地质模型拆分为与小层数相同的多个地质模型，并得到各小层的流量分配系数；1) Split the fine geological model of the work area into multiple geological models with the same number of small layers, and obtain the flow distribution coefficient of each small layer;
2)根据储层非均质状况，计算非均质条件下的理论采液强度；2) According to the heterogeneous condition of the reservoir, calculate the theoretical liquid production intensity under heterogeneous conditions;
3)筛选径向水射流单井；3) Screening of radial water jet single wells;
4)优化最佳的均衡平面驱替的径向水射流分支长度。4) Optimize the radial water jet branch length for the best balanced plane displacement.
优选地，步骤1)具体方法为：Preferably, the specific method of step 1) is:
①依据模型的分层情况，以隔层所在的纵向网格为界，拆分为与小层数相同的地质模型；① According to the stratification of the model, the vertical grid where the interlayer is located is used as the boundary to split it into geological models with the same number of small layers;
②将工区数值模型的历史油水井数据，按照网格净流出情况导出，产油量和产水量均为正值，注水量为负值，按照纵向小层网格划分情况整理历史数据，形成与小层数数量相同的历史数据文件；② The historical oil and water well data of the numerical model of the work area are exported according to the net outflow of the grid. The oil production and water production are both positive values, and the water injection volume is negative. The historical data are sorted according to the vertical small layer grid division to form a historical data file with the same number of small layers;
③将工区数值模型的射孔数据，以隔层所在纵向网格为界，拆分为与小层数数量相同的射孔数据文件；③ Split the perforation data of the numerical model of the work area into perforation data files with the same number of small layers, based on the longitudinal grid where the interlayer is located;
④将整理的各小层历史数据文件和射孔数据文件分别导入相应的小层地质模型，形成各小层的流线模型，从而得到各小层与水井对应油井的流量分配系数。④ Import the sorted historical data files and perforation data files of each sub-layer into the corresponding sub-layer geological model to form the streamline model of each sub-layer, so as to obtain the flow distribution coefficient of each sub-layer and the corresponding oil well of the water well.
优选地，模型分为4小层，1小层纵向1-9网格，2小层纵向11-21网格，3小层纵向23-27网格，4小层纵向29-34网格；Preferably, the model is divided into 4 small layers, 1 small layer has 1-9 grids in the longitudinal direction, 2 small layers have 11-21 grids in the longitudinal direction, 3 small layers have 23-27 grids in the longitudinal direction, and 4 small layers have 29-34 grids in the longitudinal direction;
优选地，隔层所在的纵向网格为纵向10、22、28网格。Preferably, the longitudinal grid where the partition is located is a longitudinal 10, 22, or 28 grid.
优选地，步骤2)非均质条件下的理论采液强度的计算方法为：Preferably, the calculation method of the theoretical liquid extraction intensity under heterogeneous conditions in step 2) is:
(1)利用工区的油水相对渗透率曲线绘制无因次采液指数随含水变化理论曲线，得到对应当前含水条件下的理论无因次采液指数；(1) Using the oil-water relative permeability curve of the work area, a theoretical curve of the dimensionless liquid production index changing with water content is drawn to obtain the theoretical dimensionless liquid production index corresponding to the current water content condition;
(2)统计工区单井初期平均日产液量，与当前含水条件下的理论无因次采液指数乘积，即为当前含水条件下的理论单井液量；理论单井液量与平均射开有效厚度的比值，即为理论采液强度。(2) The initial average daily liquid production of a single well in the work area is calculated and multiplied by the theoretical dimensionless liquid production index under the current water-cut conditions to obtain the theoretical single-well liquid production under the current water-cut conditions. The ratio of the theoretical single-well liquid production to the average effective perforation thickness is the theoretical liquid production intensity.
(3)建立非均质性理论模型，并对物性进行加权平均建立均质理论模型；定义非均质理论模型与均质理论模型的日产液量比值为干扰系数；(3) Establishing a heterogeneous theoretical model and establishing a homogeneous theoretical model by weighted average of the physical properties; defining the ratio of the daily liquid production of the heterogeneous theoretical model to that of the homogeneous theoretical model as the interference coefficient;
(4)当前含水条件下的理论无因次采液指数与干扰系数的乘积即为非均质条件下的理论采液强度。(4) The product of the theoretical dimensionless production index under the current water content conditions and the interference coefficient is the theoretical production intensity under heterogeneous conditions.
优选地，计算无因次采油指数αo的公式为：Preferably, the formula for calculating the dimensionless oil recovery indexαo is:
式中，Kro(Sw)—不同含水饱和度Sw下的油相相对渗透率；Kromax—束缚水饱和度Swi下的油相相对渗透率；K—fw＝0时的油层绝对渗透率；Kw—含水为fw时的油层绝对渗透率；Where, Kro (Sw )—relative permeability of oil phase at different water saturations Sw ; Kromax —relative permeability of oil phase at irreducible water saturation Swi ; K—absolute permeability of oil layer when fw = 0; Kw —absolute permeability of oil layer when water content is fw ;
优选地，令K＝Kw，Preferably, let K = Kw ,
无因次采液指数αl的计算公式为：The calculation formula of dimensionless liquid extraction index αl is:
其中，in,
式中，Kro—不同含水饱和度Sw下的油相相对渗透率；Krw—不同含水饱和度Sw下的水相相对渗透率；μo—地层条件下的原油粘度；μw—地层条件下的水粘度。Wherein, Kro — relative permeability of oil phase at different water saturationsSw ; Krw — relative permeability of water phase at different water saturationsSw ; μo — viscosity of crude oil under formation conditions; μw — viscosity of water under formation conditions.
优选地，步骤3)径向水射流单井的筛选方法为：将实际采液强度与理论采液强度对比，实际采液强度小于理论采液强度的油井即为筛选出的径向水射流油井；Preferably, the radial water jet single well screening method in step 3) is: comparing the actual liquid production intensity with the theoretical liquid production intensity, and the oil wells whose actual liquid production intensity is less than the theoretical liquid production intensity are the screened radial water jet oil wells;
优选地，实际采液强度为小层模型油井的产液量与该油井在该模型中的射开有效厚度的比值。Preferably, the actual liquid production intensity is the ratio of the liquid production of the oil well in the small layer model to the effective perforation thickness of the oil well in the model.
优选地，步骤4)优化最佳的均衡平面驱替的径向水射流分支长度，包括以下步骤：Preferably, step 4) optimizing the radial water jet branch length for optimal balanced plane displacement comprises the following steps:
S1.根据分层流线的形成可以显示在该层流线的分布状况以及水淹推进状况，沿平行水线推进方向确定径向水射流分支方向；S1. According to the formation of layered streamlines, the distribution of streamlines in the layer and the flooding advancement status can be displayed, and the radial water jet branching direction can be determined along the advancement direction parallel to the waterline;
S2.根据工艺可实现的最大径向水射流分支100m为限，按照分支20、40、60、80、100m分别得到流量分配系数；S2. According to the maximum radial water jet branch 100m that can be achieved by the process, the flow distribution coefficient is obtained according to the branches 20, 40, 60, 80, and 100m respectively;
S3.计算分流量曲线；S3. Calculate the flow rate distribution curve;
S4.计算各油井可采剩余油饱和度；S4. Calculate the recoverable remaining oil saturation of each oil well;
S5.以各油井地层系数与可采剩余油饱和度的乘积(Kh·Sor)和流量分配系数的比值作为参数，当级差最小时，对应的径向水射流分支长度即为优化得到的平面均衡驱替的径向水射流分支最优长度。S5. Taking the product of the formation coefficient of each oil well and the recoverable remaining oil saturation (Kh·Sor ) and the ratio of the flow distribution coefficient as parameters, when the step difference is the smallest, the corresponding radial water jet branch length is the optimal length of the radial water jet branch for the optimized plane balanced displacement.
优选地，水相分流公式为：Preferably, the water phase split formula is:
式中Kro—不同含水饱和度Sw下的油相相对渗透率；Krw—不同含水饱和度Sw下的水相相对渗透率；μo—地层条件下的原油粘度；μw—地层条件下的水粘度；a—回归系数b—回归系数Sw—含水饱和度。WhereKro is the relative permeability of the oil phase at different water saturationsSw ;Krw is the relative permeability of the water phase at different water saturationsSw ;μo is the viscosity of the crude oil at the formation conditions;μw is the viscosity of water at the formation conditions; a is the regression coefficient b—Regression coefficient Sw — water saturation.
优选地，剩余油饱和度计算公式为：Preferably, the residual oil saturation calculation formula is:
Syor＝(1-Sor)-(1-Sfw)＝Sfw-SorSyor =(1-Sor )-(1-Sfw )=Sfw -Sor
式中，Syor—可采剩余油饱和度；Sor—残余油饱和度；Sfw—目前含水条件下含水饱和度；其中Sfw可根据目前各井的含水状况采用分流量曲线获得。In the formula,Syor is the recoverable remaining oil saturation;Sor is the residual oil saturation;Sfw is the water saturation under the current water-containing conditions;Sfw can be obtained by using the fractional flow curve according to the current water-containing conditions of each well.
与现有技术相比，本发明具有以下优势：Compared with the prior art, the present invention has the following advantages:
本发明实现了分层流线的显示，分小层确定流量分配系数，并考虑非均质对产能带来的影响，以同时驱替完毕所有可动剩余油为目标，最大程度地减少无效水循环，实现最终的驱替均衡。The present invention realizes the display of stratified streamlines, determines the flow distribution coefficient by small layers, and considers the influence of heterogeneity on production capacity, with the goal of simultaneously displacing all movable residual oil, minimizing invalid water circulation, and achieving final displacement balance.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1为本发明实施例的分层流线形成的技术路线图。FIG. 1 is a technical roadmap for forming layered streamlines according to an embodiment of the present invention.
图2为本发明实施例的理论采液强度技术路线图。FIG. 2 is a technical roadmap of theoretical liquid extraction intensity according to an embodiment of the present invention.
图3为本发明实施例的工区整体流线分布示意图。FIG3 is a schematic diagram of the overall streamline distribution of the work area according to an embodiment of the present invention.
图4为本发明实施例的工区分层流线分布示意图。FIG. 4 is a schematic diagram of the distribution of stratified streamlines in a work area according to an embodiment of the present invention.
图5为本发明实施例的油水相对渗透率曲线图。FIG. 5 is a graph showing the oil-water relative permeability according to an embodiment of the present invention.
图6为本发明实施例的无因次采液指数随含水变化理论曲线图。FIG6 is a theoretical curve diagram showing the dimensionless liquid production index changing with water content according to an embodiment of the present invention.
图7为本发明实施例的非均质模型示意图。FIG. 7 is a schematic diagram of a heterogeneous model according to an embodiment of the present invention.
图8为本发明实施例的非流量曲线图。FIG. 8 is a non-flow curve diagram of an embodiment of the present invention.
图9为本发明实施例的径向水射流实施前后流线变化示意图。FIG. 9 is a schematic diagram of streamline changes before and after the implementation of the radial water jet according to an embodiment of the present invention.
具体实施方式DETAILED DESCRIPTION
应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。It should be noted that the following detailed descriptions are exemplary and are intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meanings as those commonly understood by those skilled in the art to which the present invention belongs.
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作和/或它们的组合。It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprising" and/or "including" are used in this specification, it indicates the presence of features, steps, operations and/or combinations thereof.
为了使得本领域技术人员能够更加清楚地了解本发明的技术方案，以下将结合具体的实施例详细说明本发明的技术方案。In order to enable those skilled in the art to more clearly understand the technical solution of the present invention, the technical solution of the present invention will be described in detail below in conjunction with specific embodiments.
实施例Example
一种均衡平面驱替的确定径向水射流分支长度的优化设计方法，包括如下步骤：An optimization design method for determining radial water jet branch lengths for balanced plane displacement includes the following steps:
1)将工区的精细地质模型拆分为与小层数相同的多个地质模型，并得到各小层的流量分配系数，包括以下步骤：1) Split the fine geological model of the work area into multiple geological models with the same number of small layers, and obtain the flow distribution coefficient of each small layer, including the following steps:
①依据模型的分层情况(1小层纵向1-9网格；2小层纵向11-21网格；3小层纵向23-27网格；4小层纵向29-34网格)，以隔层(纵向10、22、28网格)所在的纵向网格为界，拆分为4个小层的地质模型；① According to the stratification of the model (1st layer has 1-9 vertical grids; 2nd layer has 11-21 vertical grids; 3rd layer has 23-27 vertical grids; 4th layer has 29-34 vertical grids), the geological model is divided into 4 small layers with the vertical grids where the interlayer (10, 22, 28 vertical grids) as the boundary;
②将工区数值模型的历史油水井数据，按照网格净流出情况导出，产油量和产水量均为正值，注水量为负值，按照纵向小层网格划分情况整理历史数据，形成4个小层的历史数据文件；② The historical oil and water well data of the numerical model of the work area are exported according to the net outflow of the grid. The oil production and water production are both positive values, and the water injection volume is negative. The historical data are sorted according to the vertical small layer grid division to form a historical data file of 4 small layers;
③将工区数值模型的射孔数据，以隔层所在纵向网格为界，拆分为4个小层的射孔数据文件；③ Split the perforation data of the numerical model of the work area into four small layers of perforation data files based on the longitudinal grid where the interlayer is located;
④将整理的4个历史数据文件和4个射孔数据文件分别导入4个小层地质模型，形成4个小层的流线模型(图4)，从而得到各小层与水井对应油井的流量分配系数(表1-表4)。④ Import the sorted 4 historical data files and 4 perforation data files into the 4 small-layer geological models respectively to form the streamline model of the 4 small layers (Figure 4), so as to obtain the flow distribution coefficient of each small layer and the corresponding oil well of the water well (Table 1-Table 4).
表1 1小层流量分配系数Table 1 1 Small layer flow distribution coefficient
表2 2小层流量分配系数Table 2 2 Small layer flow distribution coefficient
表3 3小层流量分配系数Table 3 3-layer flow distribution coefficient
表4 4小层流量分配系数Table 4 4-layer flow distribution coefficient
2)根据储层非均质状况，计算理论采液强度，包括以下步骤：2) Calculate the theoretical liquid production intensity according to the reservoir heterogeneity, including the following steps:
①利用工区的油水相对渗透率曲线(图5)和地层条件下的油藏参数绘制无因次采液指数随含水变化理论曲线(图6)。① Using the oil-water relative permeability curve of the work area (Figure 5) and the reservoir parameters under formation conditions, a theoretical curve of dimensionless liquid production index changing with water content was drawn (Figure 6).
计算无因次采油指数αo的公式为：The formula for calculating the dimensionless oil recovery indexαo is:
式中，Kro(Sw)—不同含水饱和度Sw下的油相相对渗透率；Where, Kro (Sw )—relative permeability of oil phase at different water saturations Sw ;
Kromax—束缚水饱和度Swi下的油相相对渗透率；Kromax — relative permeability of oil phase at irreducible water saturation Swi ;
K—fw＝0时的油层绝对渗透率；K—absolute permeability of the oil layer whenfw =0;
Kw—含水为fw时的油层绝对渗透率。Kw — absolute permeability of the oil layer when the water content is fw .
在此不考虑注水开发过程中绝对渗透率的变化，令K＝Kw，则上式变为：The change of absolute permeability during water injection development is not considered here, and K =Kw , then the above formula becomes:
无因次采液指数αl的计算公式为：The calculation formula of dimensionless liquid extraction index αl is:
其中，in,
式中，Kro—不同含水饱和度Sw下的油相相对渗透率；Where, Kro — relative permeability of oil phase at different water saturation Sw ;
Krw—不同含水饱和度Sw下的水相相对渗透率；Krw — relative permeability of water phase at different water saturation Sw ;
μo—地层条件下的原油粘度；本实施例中μo为0.96mPa·s；μo —crude oil viscosity under formation conditions; in this embodiment, μo is 0.96 mPa·s;
μw—地层条件下的水粘度，本实施例中μw为0.30mPa·s。μw —water viscosity under formation conditions. In this embodiment, μw is 0.30 mPa·s.
②统计工区单井初期平均日产液量，与目前含水18.3％条件下的理论无因次采液指数0.64相乘，即为目前含水条件下的理论单井液量8.1m3/d。理论单井液量与平均射开有效厚度的比值，即为理论采液强度。② The average daily liquid production of a single well in the work area at the beginning is calculated and multiplied by the theoretical dimensionless liquid production index of 0.64 under the current water content of 18.3%, which is8.1m3 /d of theoretical single well liquid production under the current water content. The ratio of theoretical single well liquid production to average effective perforation thickness is the theoretical liquid production intensity.
计算目前含水下油井理论采液强度Lst的公式为：The formula for calculating the theoretical liquid production intensityLst of the current water-bearing oil well is:
式中，ql—油井平均日产液量，本实施例ql大小为8.1m3/d；Wherein, ql is the average daily liquid production of the oil well. In this embodiment, ql is 8.1 m3 /d;
h—油井平均射开有效厚度，本实施例油井平均射开有效厚度为6.3m。h—average effective perforation thickness of the oil well. In this embodiment, the average effective perforation thickness of the oil well is 6.3 m.
目前含水条件下油井理论采液强度为1.29m3/(d·m)。At present, the theoretical liquid production intensity of oil wells under water-containing conditions is 1.29m3 /(d·m).
③按照层间非均质性建立平面均质纵向非均质的理论模型(表5和图7)；以建立非均质性理论模型为基础，对各小层物性加权平均值建立均质理论模型，其中均质模型有效厚度h＝h1+h2+h3+h4＝10.1m；渗透率K＝(K1h1+K2h2+K3h3+K4h4)/h＝28.5md。③ According to the interlayer heterogeneity, a theoretical model of planar homogeneity and longitudinal heterogeneity is established (Table 5 and Figure 7); based on the establishment of the heterogeneity theoretical model, a homogeneous theoretical model is established for the weighted average value of the physical properties of each small layer, where the effective thickness of the homogeneous model is h = h1 +h2 +h3 +h4 = 10.1m; and the permeability K = (K1 h1 +K2 h2 +K3 h3 +K4 h4 )/h = 28.5md.
表5模型物性和有效厚度统计表Table 5 Statistics of model properties and effective thickness
利用数值模拟软件ecllipse分别计算非均质模型与均质模型的日液能力分别8.2t/d和6.4t/d，则干扰系数为0.78。The numerical simulation software Eclipse was used to calculate the daily liquid capacities of the heterogeneous model and the homogeneous model to be 8.2 t/d and 6.4 t/d respectively, and the interference coefficient was 0.78.
④目前含水条件下的油井理论采液强度1.29m3/(d·m)与干扰系数0.78的乘积即为理论采液强度1.0m3/(d·m)。④ The product of the theoretical production intensity of oil wells under current water-bearing conditions, 1.29 m3 /(d·m), and the interference coefficient, 0.78, is the theoretical production intensity of 1.0 m3 /(d·m).
3)筛选径向水射流单井，包括以下步骤：3) Screening radial water jet single wells, including the following steps:
①小层模型油井的产液量与该油井在该模型中的射开有效厚度的比值，得到油井在各小层的实际采液强度(表6)；① The ratio of the liquid production of the oil well in the small layer model to the effective perforation thickness of the oil well in the model is used to obtain the actual liquid production intensity of the oil well in each small layer (Table 6);
表6油井在各小层的实际采液强度统计表Table 6 Statistics of actual fluid production intensity of oil wells in each sublayer
②将实际采液强度与理论采液强度对比，实际采液强度小于理论采液强度1.0m3/(d·m)的油井层位即为筛选出的径向水射流油井层位。② Compare the actual production intensity with the theoretical production intensity. The oil well layers where the actual production intensity is less than the theoretical production intensity by 1.0m3 /(d·m) are the selected radial water jet oil well layers.
4)优化最佳的均衡平面驱替的径向水射流分支长度，包括以下步骤：4) Optimizing the radial water jet branch length for optimal balanced plane displacement, including the following steps:
①以1小层水井L87井所在井组为例，井区涉及油井3口，分别为L87-X11、L87-X12和L87-X2井。其中L87-X2井实际采液强度仅为0.35m3/(t·d)，需要利用径向水射流提高该井采液强度，均衡平面驱替。① Taking the well group where the L87 well in the first layer is located as an example, the well area involves 3 oil wells, namely L87-X11, L87-X12 and L87-X2. Among them, the actual liquid production intensity of the L87-X2 well is only0.35m3 /(t·d), and radial water jets are needed to increase the liquid production intensity of the well and balance the plane displacement.
②根据分层流线的形成可以显示在该层流线的分布状况以及水淹推进状况，沿平行水线推进方向确定径向水射流分支方向NE30°；② According to the formation of stratified streamlines, the distribution of streamlines in this layer and the flooding advancement can be displayed, and the radial water jet branch direction NE30° can be determined along the advancement direction parallel to the waterline;
③根据工艺可实现的最大径向水射流分支100m为限，按照分支20、40、60、80、100m分别得到流量分配系数(表7)；③ According to the maximum radial water jet branch 100m that can be achieved by the process, the flow distribution coefficients are obtained according to the branches of 20, 40, 60, 80, and 100m respectively (Table 7);
表7不同长度分支流量分配系数统计表Table 7 Statistics of flow distribution coefficients of branches with different lengths
④计算分流量曲线④Calculate the flow rate distribution curve
含水率计算公式为：The moisture content calculation formula is:
由于油水两相相对渗透率比值表示为含水饱和度的函数，即：Since the relative permeability ratio of oil and water phases is expressed as a function of water saturation, that is:
得到：称为水相分流量公式，结合相对渗透率数据，可以得到不同含水饱和度Sw下的含水fw,即水相分流量曲线。get: It is called the water phase fraction flow formula. Combined with the relative permeability data, the water contentfw under different water saturationsSw can be obtained, that is, the water phase fraction flow curve.
⑤计算各采油井的可采剩余油饱和度，即：⑤ Calculate the recoverable remaining oil saturation of each oil production well, that is:
Syor＝(1-Sor)-(1-Sfw)＝Sfw-SorSyor =(1-Sor )-(1-Sfw )=Sfw -Sor
式中，Syor—可采剩余油饱和度，小数；Where,Syor —recoverable remaining oil saturation, decimal;
Sor—残余油饱和度，小数；Sor — residual oil saturation, decimal;
Sfw—目前含水条件下含水饱和度，小数。Sfw — water saturation under current water conditions, decimal.
其中Sfw可根据目前各井的含水状况采用分流量曲线获得。Among them, Sfw can be obtained by using the flow rate distribution curve according to the current water content of each well.
⑥以各油井地层系数与可采剩余油饱和度的乘积(Kh·Syor)(表8)和流量分配系数的比值(定义为无因次均衡系数)作为参数(表9)，当分支长度80m时，级差最小，仅为1.51，此时平面驱替最为均衡(图8)。⑥ Taking the product of the formation coefficient of each oil well and the recoverable remaining oil saturation (Kh·Syor ) (Table 8) and the ratio of the flow distribution coefficient (defined as the dimensionless equilibrium coefficient) as parameters (Table 9), when the branch length is 80m, the step difference is the smallest, only 1.51, at this time the plane displacement is most balanced (Figure 8).
表8油井地层系数与可采剩余油饱和度统计表Table 8 Statistics of oil well formation coefficient and recoverable remaining oil saturation
表9各油井无因次均衡系数统计表Table 9 Statistics of dimensionless equilibrium coefficients of various oil wells
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。The above embodiments are preferred implementation modes of the present invention, but the implementation modes of the present invention are not limited to the above embodiments. Any other changes, modifications, substitutions, combinations, and simplifications that do not deviate from the spirit and principles of the present invention should be equivalent replacement methods and are included in the protection scope of the present invention.
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010008572.XACN111209666B (en) | 2020-01-03 | 2020-01-03 | Design method for determining radial water jet branch length for balanced plane displacement |
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010008572.XACN111209666B (en) | 2020-01-03 | 2020-01-03 | Design method for determining radial water jet branch length for balanced plane displacement |
| Publication Number | Publication Date |
|---|---|
| CN111209666A CN111209666A (en) | 2020-05-29 |
| CN111209666Btrue CN111209666B (en) | 2024-11-08 |
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010008572.XAActiveCN111209666B (en) | 2020-01-03 | 2020-01-03 | Design method for determining radial water jet branch length for balanced plane displacement |
| Country | Link |
|---|---|
| CN (1) | CN111209666B (en) |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102606129A (en)* | 2012-04-10 | 2012-07-25 | 中国海洋石油总公司 | Method and system for thin interbed oilfield development |
| CN103472484A (en)* | 2013-09-18 | 2013-12-25 | 西南石油大学 | Horizontal well track optimization method based on RS three-dimensional sensitivity seismic attribution analysis |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104060985B (en)* | 2014-06-16 | 2016-05-25 | 中国地质大学（北京） | A kind of stratified deposit profile control and water plugging blocking agent penetration depth method of testing and system |
| CN110552691B (en)* | 2018-06-04 | 2020-12-04 | 中国石油大学(北京) | A method and system for identifying dominant channels in reservoir fractures |
| CN109002574B (en)* | 2018-06-06 | 2022-11-22 | 西安石油大学 | A method for predicting development indexes of multi-layer reservoirs with pulse periodic water injection |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102606129A (en)* | 2012-04-10 | 2012-07-25 | 中国海洋石油总公司 | Method and system for thin interbed oilfield development |
| CN103472484A (en)* | 2013-09-18 | 2013-12-25 | 西南石油大学 | Horizontal well track optimization method based on RS three-dimensional sensitivity seismic attribution analysis |
| Publication number | Publication date |
|---|---|
| CN111209666A (en) | 2020-05-29 |
| Publication | Publication Date | Title |
|---|---|---|
| CN111322050B (en) | Shale horizontal well section internal osculating temporary plugging fracturing construction optimization method | |
| CN107622328B (en) | A kind of more cluster pressure break horizontal well production prediction methods of shale gas reservoir segmentation | |
| CN108301813B (en) | Method for developing and adjusting multilayer sandstone oil reservoir for development of scattered sand bodies | |
| CN103527163B (en) | A kind of compact reservoir horizontal well volume fracturing technique | |
| Deng et al. | Integrated optimization of fracture parameters for subdivision cutting fractured horizontal wells in shale oil reservoirs | |
| CN109209333B (en) | Shale gas multi-well group efficient mining interval optimization method | |
| CN109958413B (en) | A dynamic flow unit division method for ultra-high water-cut reservoirs | |
| CN102287172A (en) | Heavy oil reservoir oil extraction method | |
| CN102052067B (en) | Equal pressure drop gradient step-by-step deep control and flooding method | |
| CN106437674A (en) | Imitation horizontal well waterflooding well pattern adapting method | |
| CN103899285A (en) | Similar limit water cut period take-turn water drive method for multi-layer sandstone reservoir | |
| CN112593901A (en) | Shale gas well full production cycle production system and drainage and extraction process intelligent optimization method | |
| CN111523279A (en) | Oil reservoir layered water injection efficiency evaluation method, medium, terminal and device | |
| CN110952976B (en) | Single-well exploitation stable yield potential evaluation method in gas reservoir development mode | |
| CN112943230B (en) | Residual oil distribution prediction method for common heavy oil reservoir | |
| CN106321051A (en) | Method for optimizing multi-section fractured horizontal well network crack parameter | |
| CN107832540A (en) | A kind of compact oil reservoir technical limit well space determines method | |
| CN106651158B (en) | Quantitative evaluation method for water injection development effectiveness degree of ultra-low permeability tight reservoir horizontal well | |
| CN116025331A (en) | Invalid water identification and quantitative characterization method | |
| CN111209666B (en) | Design method for determining radial water jet branch length for balanced plane displacement | |
| CN116108572A (en) | Shale gas condensate well volume fracturing outer zone productivity contribution analysis method | |
| CN114912281A (en) | An evaluation method based on the influencing factors of shale gas aerodynamics | |
| CN103628848B (en) | Multidirectional interlayer displacement in flooding oil production method and system | |
| CN114547748B (en) | Method and system for optimally designing construction scheme of horizontal well staged fracturing technology | |
| CN118669110A (en) | A three-dimensional acid fracturing design method for carbonate reservoirs based on integrated flow |
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |