一种支撑剂有效沉降粒径的测量方法.pdf

上传人:a2 文档编号:4669194 上传时间:2018-10-26 格式:PDF 页数:14 大小:1.89MB
返回 下载 相关 举报
摘要
申请专利号:

CN201410471855.2

申请日:

2014.09.16

公开号:

CN104198345A

公开日:

2014.12.10

当前法律状态:

终止

有效性:

无权

法律详情:

未缴年费专利权终止IPC(主分类):G01N 15/04申请日:20140916授权公告日:20160629终止日期:20170916|||授权|||实质审查的生效IPC(主分类):G01N 15/04申请日:20140916|||公开

IPC分类号:

G01N15/04

主分类号:

G01N15/04

申请人:

中国石油大学(华东)

发明人:

张国栋; 李明忠; 谢建勇; 石彦; 李吉斌

地址:

266580 山东省青岛市经济技术开发区长江西路66号

优先权:

专利代理机构:

济南金迪知识产权代理有限公司 37219

代理人:

杨磊

PDF下载: PDF下载
内容摘要

本发明涉及一种支撑剂有效沉降粒径的测量方法,步骤如下:S1、测定支撑剂颗粒沉降速度;S2、给定一支撑剂颗粒有效沉降粒径初值,计算支撑剂雷诺数;S3、计算支撑剂有效沉降粒径;S4、将步骤S3计算得到的支撑剂有效沉降粒径de与初值de0比较,重复步骤S1~S4直到满足精度要求。本发明建立了支撑剂在牛顿流体、幂律流体以及粘弹性流体中沉降时的有效沉降粒径测量方法,弥补了不规则颗粒沉降速度计算时,形状因子修正系数测量、计算复杂的不足,能够精确计算支撑剂沉降速度进而精确预测裂缝中砂堤展布形态,提高压裂施工效果预测精度。

权利要求书

权利要求书
1.  一种支撑剂有效沉降粒径的测量方法,步骤如下:
S1、测定支撑剂颗粒在流体中的自由沉降速度v;
S2、给出一支撑剂颗粒有效沉降粒径初值,输入相应的各参数值,计算支撑剂雷诺数;
当流体为牛顿流体时,雷诺数按公式(1)计算,
Rep0=ρlvde0μ---(1)]]>
当流体为幂律流体时,雷诺数按公式(2)计算,
Rep0=ρlv2-nde0nK---(2)]]>
当流体为粘弹性流体时,雷诺数按公式(3)计算,
Rep0=ρlv2-nde0nK---(3)]]>
公式(1)~(3)中,Rep0表示支撑剂雷诺数,ρl表示流体密度,v表示支撑剂自由沉降速度,de0表示支撑剂有效沉降粒径初值,μ表示流体粘度,K表示流体稠度系数,n表示流体流性指数;
S3、当流体为牛顿流体时,根据步骤S2得到的Rep0判断流动类型并根据公式(4a)或(4b)计算支撑剂有效沉降粒径de,
当Rep0≤1则de=1.3553μ0.5v0.5(ρp-ρl)-0.5     (4a)
当1<Rep0≤500,则de=1.2427μ0.375v0.875ρl0.25(ρp-ρl)-0.625---(4b);]]>
当流体为幂律流体时,根据公式(5)计算
CDnRep0n=43K1ngv-n+2nρl-n+1n(ρp-ρl)---(5)]]>
然后根据公式(6)求Rep,
CDnRep0n=24XRep-n+1n---(6)]]>
公式(6)中,X表示修正系数,根据步骤S2得到的Rep0判断流动类型并按公式(7a)、(7b)或(7c)计算得到X,
当Rep0≤1×10-5,X=32n-3n2-n+3n3n---(7a)]]>
当1×10-5<Rep0≤1×10-3,X=32n-3n2-n+3n3n+1-n23n+1log(103Rep0)---(7b)]]>
当1×10-3<Rep0≤1×103,X=32n-3n2-n+3n3n+4n424Rep0n-33---(7c);]]>
进而根据公式(2)反求支撑剂有效沉降粒径de;
当流体为粘弹性流体时,根据公式(8)计算
CDnRep0n=43K1ngv-n+2nρl-n+1n(ρp-ρl)---(8)]]>
然后根据公式(9)求Rep,
CDnRep0n=24XRep-n+1n[1-α(RepWi)β]---(9)]]>
公式(9)中,X表示修正系数,根据步骤S2得到的Rep0判断流动类型并按公式(7a)、(7b)或(7c)计算得到X;Wi表示韦森伯格数,按公式(10)计算得到,
Wi=A(v/de0)b-nK---(10)]]>
进而根据公式(3)反求支撑剂有效沉降粒径de;
公式(4a)/(4b)~(10)中,de表示支撑剂有效沉降粒径,de0表示支撑剂有效沉降粒径初值,μ表示流体粘度,v表示支撑剂自由沉降速度,ρl表示流体密度,ρp表示支撑剂密度,Rep(Rep0)表示支撑剂雷诺数,K表示流体稠度系数,n表示流体流性指数,g为重力加速度,CD表示支撑剂沉降时的拖曳系数,α、β表示常数,A、b表示材料参数;
S4、将步骤S3计算得到的支撑剂有效沉降粒径de与初值de0比较,
如果︱de-de0︱≤ε,则计算结束,输出计算结果de;
如果︱de-de0︱>ε,则更新de0=(de+de0)/2,重复步骤S2~S4,直到︱de-de0︱≤ε为止;步骤S4中,ε表示精度要求。

2.  根据权利要求1所述的支撑剂有效沉降粒径的测量方法,其特征在于,步骤S1中,测定支撑剂颗粒在流体中的自由沉降速度的方法为:
首先配置实验所需溶液并置入可视化温控沉降管道中,静置使溶液温度达到平衡,通过测定支撑剂颗粒在溶液中以平衡速度自由沉降固定距离所需时间,计算得到支撑剂沉降速 度。

3.  根据权利要求2所述的支撑剂有效沉降粒径的测量方法,其特征在于,步骤S1中,实验前溶液在恒温下静置时间为24~30小时。

4.  根据权利要求2所述的支撑剂有效沉降粒径的测量方法,其特征在于,步骤S1中,实验前支撑剂在恒温下置于溶液中的时间为24~30小时。

5.  根据权利要求2所述的支撑剂有效沉降粒径的测量方法,其特征在于,步骤S1中,支撑剂粒径与管道内径之比≤0.05。

6.  根据权利要求2所述的支撑剂有效沉降粒径的测量方法,其特征在于,步骤S1中,测定沉降速度的管段固定距离位于管道中部,距管口和管底的距离为0.2~0.3m。

说明书

说明书一种支撑剂有效沉降粒径的测量方法
技术领域
本发明涉及油气田开发水力压裂技术领域,具体涉及一种支撑剂有效沉降粒径的测量方法。
背景技术
水力压裂技术自20世纪30年代诞生以来推动了石油工业的快速发展。在水力压裂设计过程中,支撑剂的沉降规律始终是石油工作者研究的主要问题之一,它直接决定裂缝中砂堤的展布形态,影响压裂施工效果。准确计算支撑剂沉降速度,精确预测裂缝中砂堤展布形态,对压裂设计和压后产能分析具有重要的意义。近年来,随着全球能源需求的增加以及常规油气资源产量的大幅降低,世界能源结构发生了巨大变化,页岩气、页岩油以及致密油气等非常规油气藏成为开发热点。与常规油气藏相比,该类储层渗透率极低,通常需要进行水平井+滑溜水压裂技术,而由于滑溜水的低粘度,支撑剂在其中的沉降行为更加显著,使得精确计算支撑剂在滑溜水压裂液中的沉降速度显得尤为重要。
在水力压裂过程中,支撑剂的沉降速度受压裂液性质、裂缝壁面、浓度以及颗粒形状等多重因素影响,其中支撑剂的不规则性是主要影响因素之一。针对非球形颗粒的沉降速度,通常引入表征颗粒不规则性的形状因子进行修正。Wadell提出了球形度的概念,它表征了与固体颗粒等体积球的表面积与固体颗粒表面积的比值,然而对于规则程度较差的颗粒,其表面积难以精确计算。Corey和Krumbein分别提出了形状系数和Krumbein球形度的概念,该两种方法需要测量不规则颗粒三维方向上的轴径,对于小颗粒这显然是很难做到的。佐藤博引入固体颗粒等价球直径,提出了面积指数的概念,虽然随着高精度电子天平的出现,不规则颗粒等价球直径相对容易测量,但对于小粒径支撑剂,该方法操作相对复杂,仍较难实现。
支撑剂有效沉降粒径表征的是与支撑剂颗粒具有相同沉降速度的等价球直径。在进行压裂设计之前,预先通过室内实验测定支撑剂颗粒在所需压裂液中的有效沉降粒径,进而在压裂设计中以此计算支撑剂沉降速度。该方法有效消除支撑剂颗粒的不规则性对其沉降速度的影响,能够精确预测裂缝中砂堤展布形态,提高压裂施工效果预测精度。
发明内容
针对现有技术的不足,尤其是不规则颗粒沉降速度计算时,球形度等修正系数测量、计算困难,本发明提供一种支撑剂在不同性质流体中沉降时的有效沉降粒径的测量方法,该方法可有效消除支撑剂不规则性对沉降速度的影响,满足现场压裂施工的设计要求。
术语说明:
牛顿流体:是指在受力后极易变形,且剪切应力与剪切速率成正比的低粘性流体。
幂律流体:是指符合τ=Kγn流变规律的流体。式中:τ--剪切应力,K--稠度系数,单位(Pa·sn),γ—剪切速率,单位(s-1),n--流性指数,K值是粘度的度量。
粘弹性流体:是指介于黏性流体和弹性固体之间,同时表现出黏性和弹性。在不超过屈服强度的条件下,剪切应力除去以后,其变形能部分复原。
本发明的技术方案如下:
一种支撑剂有效沉降粒径的测量方法,步骤如下:
S1、测定支撑剂颗粒在流体中的自由沉降速度v;
S2、给出一支撑剂颗粒有效沉降粒径初值de0,输入相应的各参数值,计算支撑剂雷诺数;
当流体为牛顿流体时,雷诺数按公式(1)计算,
Rep0=ρlvde0μ---(1)]]>
当流体为幂律流体时,雷诺数按公式(2)计算,
Rep0=ρlv2-nde0nK---(2)]]>
当流体为粘弹性流体时,雷诺数按公式(3)计算,
Rep0=ρlv2-nde0nK---(3)]]>
公式(1)~(3)中,Rep0表示支撑剂雷诺数,ρl表示流体密度(单位kg/m3),v表示支撑剂自由沉降速度(单位m/s),de0表示支撑剂有效沉降粒径初值(单位m),μ表示流体粘度(单位Pa·s),K表示流体稠度系数(单位Pa·sn),n表示流体流性指数;
S3、当流体为牛顿流体时,根据步骤S2得到的Rep0判断流动类型并根据公式(4a)或(4b)计算支撑剂有效沉降粒径de,
当Rep0≤1,则de=1.3553μ0.5v0.5(ρp-ρl)-0.5   (4a)
当1<Rep0≤500,则de=1.2427μ0.375v0.875ρl0.25(ρp-ρl)-0.625---(4b);]]>
当流体为幂律流体时,根据公式(5)计算
CDnRep0n=43K1ngv-n+2nρl-n+1n(ρp-ρl)---(5)]]>
然后根据公式(6)求Rep,
CDnRep0n=24XRep-n+1n---(6)]]>
公式(6)中,X表示修正系数,根据步骤S2得到的Rep0判断流动类型并按公式(7a)、(7b)或(7c)计算得到X,
当Rep0≤1×10-5,X=32n-3n2-n+3n3n---(7a)]]>
当1×10-5<Rep0≤1×10-3,X=32n-3n2-n+3n3n+1-n23n+1log(103Rep0)---(7b)]]>
当1×10-3<Rep0≤1×103,X=32n-3n2-n+3n3n+4n424Rep0n-33---(7c);]]>
进而根据公式(2)反求支撑剂有效沉降粒径de;
当流体为粘弹性流体时,根据公式(8)计算
CDnRep0n=43K1ngv-n+2nρl-n+1n(ρp-ρl)---(8)]]>
然后根据公式(9)求Rep,
CDnRep0n=24XRep-n+1n[1-α(RepWi)β]---(9)]]>
公式(9)中,X表示修正系数,根据步骤S2得到的Rep0判断流动类型并按公式(7a)、(7b)或(7c)计算得到X;Wi表示韦森伯格数,按公式(10)计算得到,
Wi=A(v/de0)b-nK---(10)]]>
进而根据公式(3)反求支撑剂有效沉降粒径de;
公式(4a)/(4b)~(10)中,de表示支撑剂有效沉降粒径(单位m),de0表示支撑剂有效沉降粒径初值(单位m),μ表示流体粘度(单位Pa·s),v表示支撑剂自由沉降速度(单位m/s),ρl表示流体密度(单位kg/m3),ρp表示支撑剂密度(单位kg/m3),Rep(Rep0)表示支撑剂雷诺数,K表示流体稠度系数(单位Pa·sn),n表示流体流性指数,g为重力加速度(单位m·s-2),CD表示支撑剂沉降时的拖曳系数,α、β表示常数,A、b表示材料参数;
S4、将步骤S3计算得到的支撑剂有效沉降粒径de与初值de0比较,
如果︱de-de0︱≤ε,则计算结束,输出计算结果de;
如果︱de-de0︱>ε,则更新de0=(de+de0)/2,重复步骤S2~S4,直到︱de-de0︱≤ε为止;步骤S4中,ε表示精度要求。
根据本发明,优选的,步骤S1中,测定支撑剂颗粒在流体中的自由沉降速度的方法为:
首先配置实验所需溶液并置入可视化温控沉降管道中,静置使溶液温度达到平衡,通过 测定支撑剂颗粒在溶液中以平衡速度自由沉降固定距离所需时间,计算得到支撑剂沉降速度;
进一步优选的,实验前溶液在恒温下静置时间为24~30小时,以使溶液中气体完全逃逸出,并达到温度平衡;
进一步优选的,实验前支撑剂在恒温下置于溶液中的时间为24~30小时,以使支撑剂颗粒表面气体完全逃逸出,并达到温度平衡;
进一步优选的,支撑剂粒径与管道内径之比≤0.05,以消除管道壁面对沉降速度的影响;
进一步优选的,测定沉降速度的管段固定距离位于管道中部,距管口和管底的距离为0.2~0.3m,以使支撑剂颗粒在进入测试段之前达到平衡速度,且消除管底对支撑剂沉降速度的影响。
本发明各参数单位,如无特殊说明,均采用国际单位制。
本发明的有益效果:
1、本发明提出的支撑剂有效沉降粒径的测量方法实用性强,简单易操作,克服了不规则颗粒沉降速度计算时,形状因子等参数测量、计算困难的不足。
2、本发明提出的支撑剂有效沉降粒径的测量方法有效消除了支撑剂颗粒的不规则性对其沉降速度的影响,能够精确预测裂缝中砂堤展布形态,提高压裂施工效果预测精度。
附图说明
图1为本发明实施例1中20/40目支撑剂有效沉降粒径计算结果;其中图1-a为支撑剂有效沉降粒径频率分布直方图,图1-b为支撑剂有效沉降粒径累计分布曲线。
图2为本发明实施例1中30/50目支撑剂有效沉降粒径计算结果;其中图2-a为支撑剂有效沉降粒径频率分布直方图,图2-b为支撑剂有效沉降粒径累计分布曲线。
图3为本发明实施例1中40/70目支撑剂有效沉降粒径计算结果;其中图3-a为支撑剂有效沉降粒径频率分布直方图,图3-b为支撑剂有效沉降粒径累计分布曲线。
具体实施方式
下面通过具体实施例并结合附图对本发明作进一步说明,但不限于此。
实施例1
一种支撑剂有效沉降粒径的测量方法,步骤如下:
S1、配置90wt%浓度的甘油溶液并置入可视化温控沉降管道中,测定其在25℃下的粘度、密度以及支撑剂颗粒密度;溶液静置25小时,以使溶液中气体完全逃逸出,达到温度平衡;实验前支撑剂置于溶液中25小时,以使支撑剂颗粒表面气体完全逃逸出,达到温度平衡;可视化温控沉降管道内径为0.05m,管道长1m;
在25℃恒温下,测定支撑剂颗粒在90wt%浓度的甘油溶液中的自由沉降速度;
90wt%浓度的甘油溶液为牛顿流体;
S2、假定一支撑剂颗粒有效沉降粒径初值,输入相应的各参数值,计算支撑剂雷诺数;
雷诺数按公式(1)计算,
Rep0=ρlvde0μ---(1)]]>
公式(1)中,Rep0表示支撑剂雷诺数,ρl表示流体密度(单位kg/m3),v表示支撑剂自由沉降速度(单位m/s),de0表示支撑剂有效沉降粒径初值(单位m),μ表示流体粘度(单位Pa·s);
S3、根据步骤S2得到的Rep0判断流动类型并根据公式(4a)或(4b)计算支撑剂有效沉降粒径de,
当Rep0≤1,则de=1.3553μ0.5v0.5(ρp-ρl)-0.5   (4a)
当1<Rep0<500,则de=1.2427μ0.375v0.875ρl0.25(ρp-ρl)-0.625---(4b);]]>
公式(4a)或(4b)中,de表示支撑剂有效沉降粒径(单位m),μ表示流体粘度(单位Pa·s),v表示支撑剂自由沉降速度(单位m/s),ρl表示流体密度(单位kg/m3),ρp表示支撑剂密度(单位kg/m3),Rep(Rep0)表示支撑剂雷诺数;
S4、将步骤S3计算得到的支撑剂有效沉降粒径de与初值de0比较,
如果︱de-de0︱≤ε,则计算结束,输出计算结果de;
如果︱de-de0︱>ε,则更新de0=(de+de0)/2,重复步骤S2~S4,直到︱de-de0︱≤ε为止;步骤S4中,ε表示精度要求,本实施例中ε=10-6。
重复上述步骤S1~S4,对大量支撑剂颗粒有效沉降粒径进行测定,根据Roisn-Rammler曲线对计算结果进行统计分析。
本实施例中,所述Roisn-Rammler曲线函数为:
P(di)=1-e-(didm)n---(11)]]>
其中,P(di)为小于粒径di的百分比;dm为特征粒径。n为不均匀系数,一般0.75≤n≤0.9。
如图1、图2和图3所示,支撑剂有效沉降粒径服从Roisn-Rammler分布,用Roisn-Rammler曲线对支撑剂有效沉降粒径累计分布曲线进行拟合时,拟合度较高。根据拟合结果,通过公式(11)对支撑剂有效沉降平均粒径、中值粒径等参数进行计算,以满足现场应用要求。
表1列出了20/40目、30/50目以及40/70目支撑剂有效沉降粒径统计分析结果,分别对其进行了182组、200组以及197组实验,统计计算分析得到了不同尺寸支撑剂的中值直径(de50)、平均直径(de)、de25、de75、分选系数以及偏度系数。
表1  支撑剂有效沉降粒径统计分析结果

实施例2
本实施例的流体类型为幂律流体。
一种支撑剂有效沉降粒径的测量方法,步骤如下:
S1、测定支撑剂颗粒在流体中的自由沉降速度v,步骤同实施例1;
S2、给出一支撑剂颗粒有效沉降粒径初值,输入相应的各参数值,计算支撑剂雷诺数;
雷诺数按公式(2)计算,
Rep0=ρlv2-nde0nK---(2)]]>
公式(2)中,Rep0表示支撑剂雷诺数,ρl表示流体密度(单位kg/m3),v表示支撑剂自由沉降速度(单位m/s),de0表示支撑剂有效沉降粒径初值(单位m),K表示流体稠度系数(单位Pa·sn),n表示流体流性指数;
S3、根据公式(5)计算
CDnRep0n=43K1ngv-n+2nρl-n+1n(ρp-ρl)---(5)]]>
然后根据公式(6)求Rep,
CDnRep0n=24XRep-n+1n---(6)]]>
公式(6)中,X表示修正系数,根据步骤S2得到的Rep0判断流动类型并按公式(7a)、(7b)或(7c)计算得到X,
当Rep0≤1×10-5,X=32n-3n2-n+3n3n---(7a)]]>
当1×10-5<Rep0≤1×10-3,X=32n-3n2-n+3n3n+1-n23n+1log(103Rep0)---(7b)]]>
当1×10-3<Rep0≤1×103,X=32n-3n2-n+3n3n+4n424Rep0n-33---(7c);]]>
进而根据公式(2)反求支撑剂有效沉降粒径de;
上述公式中,de表示支撑剂有效沉降粒径(单位m),v表示支撑剂自由沉降速度(单位m/s),ρl表示流体密度(单位kg/m3),ρp表示支撑剂密度(单位kg/m3),Rep(Rep0)表示支撑剂雷诺数,K表示流体稠度系数(单位Pa·sn),n表示流体流性指数,g为重力加速度(单位m·s-2),CD表示支撑剂沉降时的拖曳系数;
S4、将步骤S3计算得到的支撑剂有效沉降粒径de与初值de0比较,
如果︱de-de0︱≤ε,则计算结束,输出计算结果de;
如果︱de-de0︱>ε,则更新de0=(de+de0)/2,重复步骤S2~S4,直到︱de-de0︱≤ε为止;步骤S4中,ε表示精度要求。
实施例3
本实施例的流体类型为粘弹性流体。
一种支撑剂有效沉降粒径的测量方法,步骤如下:
S1、测定支撑剂颗粒在流体中的自由沉降速度v,步骤同实施例1;
S2、给出一支撑剂颗粒有效沉降粒径初值,输入相应的各参数值,计算支撑剂雷诺数;
雷诺数按公式(3)计算,
Rep0=ρlv2-nde0nK---(3)]]>
公式(3)中,Rep0表示支撑剂雷诺数,ρl表示流体密度(单位kg/m3),v表示支撑剂自由沉降速度(单位m/s),de0表示支撑剂有效沉降粒径初值(单位m),K表示流体稠度系数(单位Pa·sn),n表示流体流性指数;
S3、根据公式(8)计算
CDnRep0n=43K1ngv-n+2nρl-n+1n(ρp-ρl)---(8)]]>
然后根据公式(9)求Rep,
CDnRep0n=24XRep-n+1n[1-α(RepWi)β]---(9)]]>
公式(9)中,X表示修正系数,根据步骤S2得到的Rep0判断流动类型并按公式(7a)、(7b)或(7c)计算得到X;
当Rep0≤1×10-5,X=32n-3n2-n+3n3n---(7a)]]>
当1×10-5<Rep0≤1×10-3,X=32n-3n2-n+3n3n+1-n23n+1log(103Rep0)---(7b)]]>
当1×10-3<Rep0≤1×103,X=32n-3n2-n+3n3n+4n424Rep0n-33---(7c);]]>
Wi表示韦森伯格数,按公式(10)计算得到,
Wi=A(v/de0)b-nK---(10)]]>
进而根据公式(3)反求支撑剂有效沉降粒径de;
上述公式中,de表示支撑剂有效沉降粒径(单位m),v表示支撑剂自由沉降速度(单位m/s),ρl表示流体密度(单位kg/m3),ρp表示支撑剂密度(单位kg/m3),Rep(Rep0)表示支撑剂雷诺数,K表示流体稠度系数(单位Pa·sn),n表示流体流性指数,g为重力加速度(单位m·s-2),CD表示支撑剂沉降时的拖曳系数,α、β表示常数,A、b表示材料参数;
S4、将步骤S3计算得到的支撑剂有效沉降粒径de与初值de0比较,
如果︱de-de0︱≤ε,则计算结束,输出计算结果de;
如果︱de-de0︱>ε,则更新de0=(de+de0)/2,重复步骤S2~S4,直到︱de-de0︱≤ε为止;步骤S4中,ε表示精度要求。
本发明未详细阐述的技术内容属于本领域技术人员的公知技术。
尽管对本发明优选实施例进行了说明和描述,本领域技术人员应该清楚,本发明不限于优选实施例的范围,对于本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。

一种支撑剂有效沉降粒径的测量方法.pdf_第1页
第1页 / 共14页
一种支撑剂有效沉降粒径的测量方法.pdf_第2页
第2页 / 共14页
一种支撑剂有效沉降粒径的测量方法.pdf_第3页
第3页 / 共14页
点击查看更多>>
资源描述

《一种支撑剂有效沉降粒径的测量方法.pdf》由会员分享,可在线阅读,更多相关《一种支撑剂有效沉降粒径的测量方法.pdf(14页珍藏版)》请在专利查询网上搜索。

1、(10)申请公布号 CN 104198345 A (43)申请公布日 2014.12.10 CN 104198345 A (21)申请号 201410471855.2 (22)申请日 2014.09.16 G01N 15/04(2006.01) (71)申请人 中国石油大学 (华东) 地址 266580 山东省青岛市经济技术开发区 长江西路 66 号 (72)发明人 张国栋 李明忠 谢建勇 石彦 李吉斌 (74)专利代理机构 济南金迪知识产权代理有限 公司 37219 代理人 杨磊 (54) 发明名称 一种支撑剂有效沉降粒径的测量方法 (57) 摘要 本发明涉及一种支撑剂有效沉降粒径的测 量方。

2、法, 步骤如下 : S1、 测定支撑剂颗粒沉降速度 ; S2、 给定一支撑剂颗粒有效沉降粒径初值, 计算支 撑剂雷诺数 ; S3、 计算支撑剂有效沉降粒径 ; S4、 将 步骤S3计算得到的支撑剂有效沉降粒径de与初值 de0比较, 重复步骤 S1 S4直到满足精度要求。本 发明建立了支撑剂在牛顿流体、 幂律流体以及粘 弹性流体中沉降时的有效沉降粒径测量方法, 弥 补了不规则颗粒沉降速度计算时, 形状因子修正 系数测量、 计算复杂的不足, 能够精确计算支撑剂 沉降速度进而精确预测裂缝中砂堤展布形态, 提 高压裂施工效果预测精度。 (51)Int.Cl. 权利要求书 2 页 说明书 8 页 附图。

3、 3 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书8页 附图3页 (10)申请公布号 CN 104198345 A CN 104198345 A 1/2 页 2 1. 一种支撑剂有效沉降粒径的测量方法, 步骤如下 : S1、 测定支撑剂颗粒在流体中的自由沉降速度 v ; S2、 给出一支撑剂颗粒有效沉降粒径初值, 输入相应的各参数值, 计算支撑剂雷诺数 ; 当流体为牛顿流体时, 雷诺数按公式 (1) 计算, 当流体为幂律流体时, 雷诺数按公式 (2) 计算, 当流体为粘弹性流体时, 雷诺数按公式 (3) 计算, 公式 (1) (3) 中, Rep0表示。

4、支撑剂雷诺数, l表示流体密度, v 表示支撑剂自由沉降 速度, de0表示支撑剂有效沉降粒径初值, 表示流体粘度, K 表示流体稠度系数, n 表示流 体流性指数 ; S3、 当流体为牛顿流体时, 根据步骤 S2得到的 Rep0判断流动类型并根据公式 (4a) 或 (4b) 计算支撑剂有效沉降粒径 de, 当 Rep0 1 则 de 1.35530.5v0.5(p-l)-0.5 (4a) 当 1 Rep0 500, 则 当流体为幂律流体时, 根据公式 (5) 计算 然后根据公式 (6) 求 Rep, 公式 (6) 中, X 表示修正系数, 根据步骤 S2得到的 Rep0判断流动类型并按公式 。

5、(7a)、 (7b) 或 (7c) 计算得到 X, 当 Rep0 110-5, 当 110-5 Rep0 110-3, 权 利 要 求 书 CN 104198345 A 2 2/2 页 3 当 110-3 Rep0 1103, 进而根据公式 (2) 反求支撑剂有效沉降粒径 de; 当流体为粘弹性流体时, 根据公式 (8) 计算 然后根据公式 (9) 求 Rep, 公式 (9) 中, X 表示修正系数, 根据步骤 S2得到的 Rep0判断流动类型并按公式 (7a)、 (7b) 或 (7c) 计算得到 X ; Wi表示韦森伯格数, 按公式 (10) 计算得到, 进而根据公式 (3) 反求支撑剂有效。

6、沉降粒径 de; 公式 (4a)/(4b) (10) 中, de表示支撑剂有效沉降粒径, de0表示支撑剂有效沉降粒径 初值, 表示流体粘度, v 表示支撑剂自由沉降速度, l表示流体密度, p表示支撑剂密 度, Rep(Rep0) 表示支撑剂雷诺数, K 表示流体稠度系数, n 表示流体流性指数, g 为重力加速 度, CD 表示支撑剂沉降时的拖曳系数, 、 表示常数, A、 b 表示材料参数 ; S4、 将步骤 S3计算得到的支撑剂有效沉降粒径 de与初值 de0比较, 如果 de-de0 , 则计算结束, 输出计算结果 de; 如果 de-de0 , 则更新 de0 (de+de0)/2。

7、, 重复步骤 S2 S4, 直到 de-de0 为止 ; 步骤 S4中, 表示精度要求。 2.根据权利要求1所述的支撑剂有效沉降粒径的测量方法, 其特征在于, 步骤S1中, 测 定支撑剂颗粒在流体中的自由沉降速度的方法为 : 首先配置实验所需溶液并置入可视化温控沉降管道中, 静置使溶液温度达到平衡, 通 过测定支撑剂颗粒在溶液中以平衡速度自由沉降固定距离所需时间, 计算得到支撑剂沉降 速度。 3.根据权利要求2所述的支撑剂有效沉降粒径的测量方法, 其特征在于, 步骤S1中, 实 验前溶液在恒温下静置时间为 24 30 小时。 4.根据权利要求2所述的支撑剂有效沉降粒径的测量方法, 其特征在于,。

8、 步骤S1中, 实 验前支撑剂在恒温下置于溶液中的时间为 24 30 小时。 5.根据权利要求2所述的支撑剂有效沉降粒径的测量方法, 其特征在于, 步骤S1中, 支 撑剂粒径与管道内径之比 0.05。 6.根据权利要求2所述的支撑剂有效沉降粒径的测量方法, 其特征在于, 步骤S1中, 测 定沉降速度的管段固定距离位于管道中部, 距管口和管底的距离为 0.2 0.3m。 权 利 要 求 书 CN 104198345 A 3 1/8 页 4 一种支撑剂有效沉降粒径的测量方法 技术领域 0001 本发明涉及油气田开发水力压裂技术领域, 具体涉及一种支撑剂有效沉降粒径的 测量方法。 背景技术 0002。

9、 水力压裂技术自 20 世纪 30 年代诞生以来推动了石油工业的快速发展。在水力压 裂设计过程中, 支撑剂的沉降规律始终是石油工作者研究的主要问题之一, 它直接决定裂 缝中砂堤的展布形态, 影响压裂施工效果。 准确计算支撑剂沉降速度, 精确预测裂缝中砂堤 展布形态, 对压裂设计和压后产能分析具有重要的意义。 近年来, 随着全球能源需求的增加 以及常规油气资源产量的大幅降低, 世界能源结构发生了巨大变化, 页岩气、 页岩油以及致 密油气等非常规油气藏成为开发热点。 与常规油气藏相比, 该类储层渗透率极低, 通常需要 进行水平井 + 滑溜水压裂技术, 而由于滑溜水的低粘度, 支撑剂在其中的沉降行为。

10、更加显 著, 使得精确计算支撑剂在滑溜水压裂液中的沉降速度显得尤为重要。 0003 在水力压裂过程中, 支撑剂的沉降速度受压裂液性质、 裂缝壁面、 浓度以及颗粒形 状等多重因素影响, 其中支撑剂的不规则性是主要影响因素之一。针对非球形颗粒的沉降 速度, 通常引入表征颗粒不规则性的形状因子进行修正。Wadell 提出了球形度的概念, 它 表征了与固体颗粒等体积球的表面积与固体颗粒表面积的比值, 然而对于规则程度较差的 颗粒, 其表面积难以精确计算。Corey 和 Krumbein 分别提出了形状系数和 Krumbein 球形 度的概念, 该两种方法需要测量不规则颗粒三维方向上的轴径, 对于小颗粒。

11、这显然是很难 做到的。 佐藤博引入固体颗粒等价球直径, 提出了面积指数的概念, 虽然随着高精度电子天 平的出现, 不规则颗粒等价球直径相对容易测量, 但对于小粒径支撑剂, 该方法操作相对复 杂, 仍较难实现。 0004 支撑剂有效沉降粒径表征的是与支撑剂颗粒具有相同沉降速度的等价球直径。 在 进行压裂设计之前, 预先通过室内实验测定支撑剂颗粒在所需压裂液中的有效沉降粒径, 进而在压裂设计中以此计算支撑剂沉降速度。 该方法有效消除支撑剂颗粒的不规则性对其 沉降速度的影响, 能够精确预测裂缝中砂堤展布形态, 提高压裂施工效果预测精度。 发明内容 0005 针对现有技术的不足, 尤其是不规则颗粒沉降。

12、速度计算时, 球形度等修正系数测 量、 计算困难, 本发明提供一种支撑剂在不同性质流体中沉降时的有效沉降粒径的测量方 法, 该方法可有效消除支撑剂不规则性对沉降速度的影响, 满足现场压裂施工的设计要求。 0006 术语说明 : 0007 牛顿流体 : 是指在受力后极易变形, 且剪切应力与剪切速率成正比的低粘性流体。 0008 幂律流体 : 是指符合 Kn流变规律的流体。式中 : - 剪切应力, K- 稠度 系数, 单位 (Pasn), 剪切速率, 单位 (s-1), n- 流性指数, K 值是粘度的度量。 0009 粘弹性流体 : 是指介于黏性流体和弹性固体之间, 同时表现出黏性和弹性。 在不。

13、超 说 明 书 CN 104198345 A 4 2/8 页 5 过屈服强度的条件下, 剪切应力除去以后, 其变形能部分复原。 0010 本发明的技术方案如下 : 0011 一种支撑剂有效沉降粒径的测量方法, 步骤如下 : 0012 S1、 测定支撑剂颗粒在流体中的自由沉降速度 v ; 0013 S2、 给出一支撑剂颗粒有效沉降粒径初值 de0, 输入相应的各参数值, 计算支撑剂雷 诺数 ; 0014 当流体为牛顿流体时, 雷诺数按公式 (1) 计算, 0015 0016 当流体为幂律流体时, 雷诺数按公式 (2) 计算, 0017 0018 当流体为粘弹性流体时, 雷诺数按公式 (3) 计算。

14、, 0019 0020 公式 (1) (3) 中, Rep0表示支撑剂雷诺数, l表示流体密度 ( 单位 kg/m3), v 表 示支撑剂自由沉降速度 ( 单位 m/s), de0表示支撑剂有效沉降粒径初值 ( 单位 m), 表示 流体粘度 ( 单位 Pas), K 表示流体稠度系数 ( 单位 Pasn), n 表示流体流性指数 ; 0021 S3、 当流体为牛顿流体时, 根据步骤 S2得到的 Rep0判断流动类型并根据公式 (4a) 或 (4b) 计算支撑剂有效沉降粒径 de, 0022 当 Rep0 1, 则 de 1.35530.5v0.5(p-l)-0.5 (4a) 0023 当 1 。

15、Rep0 500, 则 0024 当流体为幂律流体时, 根据公式 (5) 计算 0025 0026 然后根据公式 (6) 求 Rep, 0027 0028 公式 (6) 中, X 表示修正系数, 根据步骤 S2得到的 Rep0判断流动类型并按公式 (7a)、 (7b) 或 (7c) 计算得到 X, 0029 当 Rep0 110-5, 0030 当 110-5 Rep0 110-3, 说 明 书 CN 104198345 A 5 3/8 页 6 0031 当 110-3 Rep0 1103, 0032 进而根据公式 (2) 反求支撑剂有效沉降粒径 de; 0033 当流体为粘弹性流体时, 根据。

16、公式 (8) 计算 0034 0035 然后根据公式 (9) 求 Rep, 0036 0037 公式 (9) 中, X 表示修正系数, 根据步骤 S2得到的 Rep0判断流动类型并按公式 (7a)、 (7b) 或 (7c) 计算得到 X ; Wi表示韦森伯格数, 按公式 (10) 计算得到, 0038 0039 进而根据公式 (3) 反求支撑剂有效沉降粒径 de; 0040 公式 (4a)/(4b) (10) 中, de表示支撑剂有效沉降粒径 ( 单位 m), de0表示支撑剂 有效沉降粒径初值 ( 单位 m), 表示流体粘度 ( 单位 Pas), v 表示支撑剂自由沉降速度 ( 单位 m/s。

17、), l表示流体密度 ( 单位 kg/m3), p表示支撑剂密度 ( 单位 kg/m3), Rep(Rep0) 表示支撑剂雷诺数, K 表示流体稠度系数 ( 单位 Pasn), n 表示流体流性指数, g 为重力加 速度 ( 单位 ms-2), CD 表示支撑剂沉降时的拖曳系数, 、 表示常数, A、 b 表示材料参 数 ; 0041 S4、 将步骤 S3计算得到的支撑剂有效沉降粒径 de与初值 de0比较, 0042 如果 de-de0 , 则计算结束, 输出计算结果 de; 0043 如果 de-de0 , 则更新 de0 (de+de0)/2, 重复步骤 S2 S4, 直到 de-de0。

18、 为止 ; 步骤 S4中, 表示精度要求。 0044 根据本发明, 优选的, 步骤 S1中, 测定支撑剂颗粒在流体中的自由沉降速度的方法 为 : 0045 首先配置实验所需溶液并置入可视化温控沉降管道中, 静置使溶液温度达到平 衡, 通过测定支撑剂颗粒在溶液中以平衡速度自由沉降固定距离所需时间, 计算得到支撑 剂沉降速度 ; 0046 进一步优选的, 实验前溶液在恒温下静置时间为2430小时, 以使溶液中气体完 全逃逸出, 并达到温度平衡 ; 0047 进一步优选的, 实验前支撑剂在恒温下置于溶液中的时间为2430小时, 以使支 撑剂颗粒表面气体完全逃逸出, 并达到温度平衡 ; 0048 进一。

19、步优选的, 支撑剂粒径与管道内径之比 0.05, 以消除管道壁面对沉降速度 的影响 ; 0049 进一步优选的, 测定沉降速度的管段固定距离位于管道中部, 距管口和管底的距 说 明 书 CN 104198345 A 6 4/8 页 7 离为 0.2 0.3m, 以使支撑剂颗粒在进入测试段之前达到平衡速度, 且消除管底对支撑剂 沉降速度的影响。 0050 本发明各参数单位, 如无特殊说明, 均采用国际单位制。 0051 本发明的有益效果 : 0052 1、 本发明提出的支撑剂有效沉降粒径的测量方法实用性强, 简单易操作, 克服了 不规则颗粒沉降速度计算时, 形状因子等参数测量、 计算困难的不足。。

20、 0053 2、 本发明提出的支撑剂有效沉降粒径的测量方法有效消除了支撑剂颗粒的不规 则性对其沉降速度的影响, 能够精确预测裂缝中砂堤展布形态, 提高压裂施工效果预测精 度。 附图说明 0054 图 1 为本发明实施例 1 中 20/40 目支撑剂有效沉降粒径计算结果 ; 其中图 1-a 为 支撑剂有效沉降粒径频率分布直方图, 图 1-b 为支撑剂有效沉降粒径累计分布曲线。 0055 图 2 为本发明实施例 1 中 30/50 目支撑剂有效沉降粒径计算结果 ; 其中图 2-a 为 支撑剂有效沉降粒径频率分布直方图, 图 2-b 为支撑剂有效沉降粒径累计分布曲线。 0056 图 3 为本发明实施。

21、例 1 中 40/70 目支撑剂有效沉降粒径计算结果 ; 其中图 3-a 为 支撑剂有效沉降粒径频率分布直方图, 图 3-b 为支撑剂有效沉降粒径累计分布曲线。 具体实施方式 0057 下面通过具体实施例并结合附图对本发明作进一步说明, 但不限于此。 0058 实施例 1 0059 一种支撑剂有效沉降粒径的测量方法, 步骤如下 : 0060 S1、 配置 90wt浓度的甘油溶液并置入可视化温控沉降管道中, 测定其在 25下 的粘度、 密度以及支撑剂颗粒密度 ; 溶液静置 25 小时, 以使溶液中气体完全逃逸出, 达到温 度平衡 ; 实验前支撑剂置于溶液中 25 小时, 以使支撑剂颗粒表面气体完。

22、全逃逸出, 达到温 度平衡 ; 可视化温控沉降管道内径为 0.05m, 管道长 1m ; 0061 在 25恒温下, 测定支撑剂颗粒在 90wt浓度的甘油溶液中的自由沉降速度 ; 0062 90wt浓度的甘油溶液为牛顿流体 ; 0063 S2、 假定一支撑剂颗粒有效沉降粒径初值, 输入相应的各参数值, 计算支撑剂雷诺 数 ; 0064 雷诺数按公式 (1) 计算, 0065 0066 公式 (1) 中, Rep0表示支撑剂雷诺数, l表示流体密度 ( 单位 kg/m3), v 表示支撑 剂自由沉降速度 ( 单位 m/s), de0表示支撑剂有效沉降粒径初值 ( 单位 m), 表示流体粘 度 (。

23、 单位 Pas) ; 0067 S3、 根据步骤 S2得到的 Rep0判断流动类型并根据公式 (4a) 或 (4b) 计算支撑剂有 效沉降粒径 de, 0068 当 Rep0 1, 则 de 1.35530.5v0.5(p-l)-0.5 (4a) 说 明 书 CN 104198345 A 7 5/8 页 8 0069 当 1 Rep0 500, 则 0070 公式(4a)或(4b)中, de表示支撑剂有效沉降粒径(单位m), 表示流体粘度(单 位 Pas), v 表示支撑剂自由沉降速度 ( 单位 m/s), l表示流体密度 ( 单位 kg/m3), p表 示支撑剂密度 ( 单位 kg/m3),。

24、 Rep(Rep0) 表示支撑剂雷诺数 ; 0071 S4、 将步骤 S3计算得到的支撑剂有效沉降粒径 de与初值 de0比较, 0072 如果 de-de0 , 则计算结束, 输出计算结果 de; 0073 如果 de-de0 , 则更新 de0 (de+de0)/2, 重复步骤 S2 S4, 直到 de-de0 为止 ; 步骤 S4中, 表示精度要求, 本实施例中 10-6。 0074 重复上述步骤 S1 S4, 对大量支撑剂颗粒有效沉降粒径进行测定, 根据 Roisn-Rammler 曲线对计算结果进行统计分析。 0075 本实施例中, 所述 Roisn-Rammler 曲线函数为 : 。

25、0076 0077 其中, P(di) 为小于粒径 di的百分比 ; dm为特征粒径。n 为不均匀系数, 一般 0.75 n 0.9。 0078 如图 1、 图 2 和图 3 所示, 支撑剂有效沉降粒径服从 Roisn-Rammler 分布, 用 Roisn-Rammler曲线对支撑剂有效沉降粒径累计分布曲线进行拟合时, 拟合度较高。 根据拟 合结果, 通过公式 (11) 对支撑剂有效沉降平均粒径、 中值粒径等参数进行计算, 以满足现 场应用要求。 0079 表 1 列出了 20/40 目、 30/50 目以及 40/70 目支撑剂有效沉降粒径统计分析结果, 分别对其进行了 182 组、 20。

26、0 组以及 197 组实验, 统计计算分析得到了不同尺寸支撑剂的中 值直径 (de50)、 平均直径 (de)、 de25、 de75、 分选系数以及偏度系数。 0080 表 1 支撑剂有效沉降粒径统计分析结果 0081 0082 实施例 2 0083 本实施例的流体类型为幂律流体。 0084 一种支撑剂有效沉降粒径的测量方法, 步骤如下 : 0085 S1、 测定支撑剂颗粒在流体中的自由沉降速度 v, 步骤同实施例 1 ; 0086 S2、 给出一支撑剂颗粒有效沉降粒径初值, 输入相应的各参数值, 计算支撑剂雷诺 数 ; 0087 雷诺数按公式 (2) 计算, 说 明 书 CN 104198。

27、345 A 8 6/8 页 9 0088 0089 公式 (2) 中, Rep0表示支撑剂雷诺数, l表示流体密度 ( 单位 kg/m3), v 表示支撑 剂自由沉降速度 ( 单位 m/s), de0表示支撑剂有效沉降粒径初值 ( 单位 m), K 表示流体稠度 系数 ( 单位 Pasn), n 表示流体流性指数 ; 0090 S3、 根据公式 (5) 计算 0091 0092 然后根据公式 (6) 求 Rep, 0093 0094 公式 (6) 中, X 表示修正系数, 根据步骤 S2得到的 Rep0判断流动类型并按公式 (7a)、 (7b) 或 (7c) 计算得到 X, 0095 当 Re。

28、p0 110-5, 0096 当 110-5 Rep0 110-3, 0097 当 110-3 Rep0 1103, 0098 进而根据公式 (2) 反求支撑剂有效沉降粒径 de; 0099 上述公式中, de表示支撑剂有效沉降粒径 ( 单位 m), v 表示支撑剂自由沉降速度 ( 单位 m/s), l表示流体密度 ( 单位 kg/m3), p表示支撑剂密度 ( 单位 kg/m3), Rep(Rep0) 表示支撑剂雷诺数, K 表示流体稠度系数 ( 单位 Pasn), n 表示流体流性指数, g 为重力加 速度 ( 单位 ms-2), CD 表示支撑剂沉降时的拖曳系数 ; 0100 S4、 将。

29、步骤 S3计算得到的支撑剂有效沉降粒径 de与初值 de0比较, 0101 如果 de-de0 , 则计算结束, 输出计算结果 de; 0102 如果 de-de0 , 则更新 de0 (de+de0)/2, 重复步骤 S2 S4, 直到 de-de0 为止 ; 步骤 S4中, 表示精度要求。 0103 实施例 3 0104 本实施例的流体类型为粘弹性流体。 0105 一种支撑剂有效沉降粒径的测量方法, 步骤如下 : 0106 S1、 测定支撑剂颗粒在流体中的自由沉降速度 v, 步骤同实施例 1 ; 0107 S2、 给出一支撑剂颗粒有效沉降粒径初值, 输入相应的各参数值, 计算支撑剂雷诺 数。

30、 ; 说 明 书 CN 104198345 A 9 7/8 页 10 0108 雷诺数按公式 (3) 计算, 0109 0110 公式 (3) 中, Rep0表示支撑剂雷诺数, l表示流体密度 ( 单位 kg/m3), v 表示支撑 剂自由沉降速度 ( 单位 m/s), de0表示支撑剂有效沉降粒径初值 ( 单位 m), K 表示流体稠度 系数 ( 单位 Pasn), n 表示流体流性指数 ; 0111 S3、 根据公式 (8) 计算 0112 0113 然后根据公式 (9) 求 Rep, 0114 0115 公式 (9) 中, X 表示修正系数, 根据步骤 S2得到的 Rep0判断流动类型并。

31、按公式 (7a)、 (7b) 或 (7c) 计算得到 X ; 0116 当 Rep0 110-5, 0117 当 110-5 Rep0 110-3, 0118 当 110-3 Rep0 1103, 0119 Wi表示韦森伯格数, 按公式 (10) 计算得到, 0120 0121 进而根据公式 (3) 反求支撑剂有效沉降粒径 de; 0122 上述公式中, de表示支撑剂有效沉降粒径 ( 单位 m), v 表示支撑剂自由沉降速度 ( 单位 m/s), l表示流体密度 ( 单位 kg/m3), p表示支撑剂密度 ( 单位 kg/m3), Rep(Rep0) 表示支撑剂雷诺数, K 表示流体稠度系数。

32、 ( 单位 Pasn), n 表示流体流性指数, g 为重力加 速度 ( 单位 ms-2), CD 表示支撑剂沉降时的拖曳系数, 、 表示常数, A、 b 表示材料参 数 ; 0123 S4、 将步骤 S3计算得到的支撑剂有效沉降粒径 de与初值 de0比较, 0124 如果 de-de0 , 则计算结束, 输出计算结果 de; 0125 如果 de-de0 , 则更新 de0 (de+de0)/2, 重复步骤 S2 S4, 直到 de-de0 为止 ; 步骤 S4中, 表示精度要求。 0126 本发明未详细阐述的技术内容属于本领域技术人员的公知技术。 说 明 书 CN 104198345 A。

33、 10 8/8 页 11 0127 尽管对本发明优选实施例进行了说明和描述, 本领域技术人员应该清楚, 本发明 不限于优选实施例的范围, 对于本技术领域的普通技术人员来讲, 只要各种变化在所附的 权利要求限定和确定的本发明的精神和范围内, 这些变化是显而易见的, 一切利用本发明 构思的发明创造均在保护之列。 说 明 书 CN 104198345 A 11 1/3 页 12 图 1-a 图 1-b 说 明 书 附 图 CN 104198345 A 12 2/3 页 13 图 2-a 图 2-b 说 明 书 附 图 CN 104198345 A 13 3/3 页 14 图 3-a 图 3-b 说 明 书 附 图 CN 104198345 A 14 。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 物理 > 测量;测试


copyright@ 2017-2020 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备2021068784号-1