1、第九章中子测井第九章 中子测井(Neutron log)利用中子与地层相互作用的各种效应,来研究钻井地质剖面的一类测井方法统称中子测井。 它是利用岩石的另一种特性,即岩石中的含氢量来研究岩石性质和孔隙度等地质问题。 这种测井方法在于将装有中子源和探测器的井下仪器下入井中,由中子源中子进入岩层, 同物质的原子核发生碰撞将产生减速、扩散和被俘获几个过程,到达探测器。 在这些过程中,探测器周围的中子分布状况,以及中子被俘获后所放出的伽马射线强度,与仪器周围的岩石性质,特别是岩石的含氢量有关。 而储集层的含氢量又取决于它的孔隙度,因此,中子测井是目前广泛使用的一种孔隙度测井。 根据中子测井的记录内容:
2、可以将它分为中子-中子测井和中子-伽马测井。根据仪器的结构特点,中子中子测井又可分为 中子-超热中测井(SNP)井壁中子测井 中子-热中子测井(CNL)补偿中子测井一、中子测井的核物理基础 1 中子和中子源 中子是组成原子核的一种不带电荷的中性粒子,其质量与氢核的质量相近。中子与物质作用时,能穿过原子的电子壳层而与原子核相碰撞,所以它对物质的穿透能力较强。 通常中子与质子以很强的核力结合在一起,形成稳定的原子核。要使中子从原子核里释放出来,就必须供给一定的能量。如果使原子核获得的能量大于中子结合能,中子就可能从核中发射出来。 可以用粒子、氘核d、质子p或光子轰击原子核,引起各种核反应,使中子从
3、核内释放出来。这种产生中子的装置称中子源。一、中子测井的核物理基础 因为不同能量的中子与原子核作用时有着不同的特点,所以通常根据中子的能量大小,可以把它分成几类:v高能快中子:能量大于10万电子伏特;v中能中子:能量在100电子伏特10万电子伏特之间;v慢中子:能量小于100电子伏特;其中0.1100电子伏特的中子为超热中子;能量等于0.025电子伏特的中子为热中子。一、 中子测井的核物理基础 1 中子和中子源 中子测井所用的中子源有两类: 即同位素中子源和加速器中子源。v 同位素中子源:如镅铍(Am-Be)中子源,利用镅衰变产生的粒子去轰击铍原子核,发生核反应而放出中子。产生的中子的平均能量
4、约5MeV。 该类中子源的特点是连续发射中子。v 加速器中子源:(亦称脉冲中子源),如D-T加速器中子源,用加速器加速氘核(D)去轰击氚核(T)产生快中子,其能量是14MeV。 该类中子源的特点是人为控制脉冲式发射中子。二、中子与物质作用几种作用形式:(1)非弹性作用:高能快中子与原子核碰撞(2)弹性散射:高能快中子经一、二次非弹性散射后,能量降低,继续碰撞原子,降低能量和运动速度,而总能量不变,经多次碰撞,能量损失,速度降低,最后变为热中子。(3)辐射俘获:能量低的热中子在其他物质附近漫游,很容易被其他物质俘获而被吸收,其他物质由于俘获中子后则处于激发态,在由激发态向稳定态转变时,则易放出r
5、射线。二、中子与物质的相互作用 由中子源发射出来的快中子与组成物质的原子核发生作用,可以分为以下几个阶段:1.快中子的减速过程 平均能量约4Mev的高能快中子碰撞原子核-发生弹性散射中子一部分能量传给原子核,成为原子核动能,中子本身的能量减少,运动速度降低-继续碰撞其它原子核.反复多次,能量不断损失,速度不断减慢,最后中子成为热中子,此过程为快中子的减速过程。 岩石中不同元素对中子产生弹性散射几率(散射截面)不同,H元素弹性散射截面最大。不同元素减速能力不同,轻原子核对中子减速起主要作用,特别是氢原子核与H碰撞,减速成热中子过程最快,因此,高含H岩石中,快中子将很快减速成热中子。1.快中子的减
6、速过程 在减速过程中,中子与原子核正面碰撞一次可损失的最大能量E为式中El中子碰撞前的能量 式中A为原子量 对于氢元素,质量A1,因而a0,EEl,即中子与氢核发生碰撞时,中子就失去全部能量。对于碳元素,A12,a0.716,中子与碳核碰撞时,中子损失的最大能量为0.28E1。显然A越大的元素,中子与它碰撞时能量损失越小。 在实际的弹性散射过程中,中子与靶核并不总是正面碰撞,因此,每次碰撞后,中子损失的能量并不相同,这与散射角有关。当快中子与原子核碰撞多次,使中子能量降低为0.025电子伏特时,这时的中子为热中子。中子变为热中子之后,就象分子热运动一样在物质中进行扩散,当它再与原子核发生碰撞时
7、,失去和得到的能量几乎相等。1.快中子的减速过程 一个中子与一个原子核发生弹性散射的几率称为微观弹性散射截面s,单位为巴(10-24cm2)。 1cm3物质的原子核的微观弹性散射截面之和叫宏观弹性散射截面s。 通常可以利用宏观弹性散射截面来描述这个减速过程。 此外,还可用“减速长度Ls”来描述快中子变为热中子的减速过程。减速长度定义为由快中子减速成热中子所经过的直线距离的平均值,单位为厘米。 下表为沉积岩中常见元素的散射截面和每次碰撞的最大能量损失以及中子能量由2百万电子伏特减速为热中子所需的平均碰撞次数。 1.快中子的减速过程 从表中可以看出,沉积岩中不同元素对中子产生弹性散射的几率(散射截
8、面)不同,氢元素的弹性散射截面最大。另外,不同元素对中子的减速能力也不同,和氢核相碰撞,能量损失最大,减速成热中子的过程也最快。因此,在含氢量高的岩石中,快中子将很快减速成热中子。 介质含氢越多,减速长度越短,这也说明氢元素对快中子的减速能力最大。 氢是所有元素中最强的中子减速剂,这是中子测井法测定地层含氢量及解决与含氢量有关的各种地质问题的依据。2.热中子的扩散及被俘获 快中子减速成热中子之后,同气体分子的扩散类似,便从密度大的地方向密度小的地方扩散。热中子扩散时,由于速度较慢,在原子核周围停留的时间相对较长,因而很容易被原子核俘获。标准热中子能量为0.025MeV,速度为2.2105cm/
9、s。 热中子被元素原子核俘获的几率取决于元素的俘获能力,通常用“宏观俘获截面a”来量度。单位为巴。 下表给出了沉积岩中常见的几种元素的微观俘获截面。 氯元素特别是硼的俘获截面很大。在油、气井中,氯元素是常见的,因此,它的存在将使热中子被俘获的几率显著增加,热中子扩散的过程或扩散距离将缩短。所以含有高矿化度水的岩石比含油的同类岩石宏观俘获截面大。2.热中子的扩散及被俘获 描述热中子扩散及俘获特性的参数有扩散长度Ld、宏观俘获截面a和热中子寿命t参数。 扩散长度:从产生热中子起到其被俘获吸收为止,热中子移动的距离。物质对热中子俘获吸收能力越强,扩散长度Ld就越短。 微观俘获截面:一个原子核俘获热中
10、子的几率称之。 宏观俘获截面a:1cm3物质中所有原子核的微观俘获截面之和称之为宏观俘获截面。 热中子寿命t:从热中子生成开始到它被俘获吸收为止所经过的平均时间叫热中子寿命,它和宏观俘获截面的关系是: 式中v为热中子移动速度,常温下,v0.22cm/s,所以上式可写成: 当地层中含有俘获截面高的元素时,t就大大减小。高矿化度水的t要比油层小的多,因此可以确定油水界面和区分油水层。2.热中子的扩散及被俘获 元素的原子核俘获热中子之后,处于激发状态,当它回到稳定的基态时,多余的能量便以伽马射线的形式释放出来。该射线称为俘获伽马射线,或次生伽马射线。不同元素俘获热中子后放出的俘获伽马射线的能量存在一
11、定的差别,特别是氯元素释放出的俘获伽马射线能量要比一般元素高一些,且伽马射线的数目也相对多些。因此,当岩石中有氯元素存在时,测得的热中子数将显著减少,但测得的俘获伽马射线却又会普遍增高。3.中子探测器 中子探测器探测的是超热中子和热中子。 热中子探测器通常由普通的闪烁计数器在其外壁上涂上锂或硼构成。由于锂和硼对热中子有强吸收性,并在吸收热中子后发生核反应而放射出粒子,该粒子能使闪烁计数器中萤光体发光,从而在记数管中的阳极产生电脉冲,然后送入地面记录仪便可对其记录。 超热中子探测器是热中子计数器在其外壁上加一层石蜡和一层镉构成。镉的作用是吸收探测器周围的热中子,而只让超热中子通过,并进入石蜡层,
12、然后再经石蜡减速为热中子,便可被热中子计数管对其记录。4 地层的含氢指数 前面已经讲过,地层对快中子的减速能力主要决定于地层的含氢量。含氢量高的地层宏观减速能力强、减速长度小。 为了方便,在中子测井中把淡水的含氢量规定为一个单位,用它来衡量地层中所有其它岩石或矿物的含氢量。单位体积的任何岩石或矿物中氢核数与同样体积的淡水中氢核数的比值,称为该岩石或矿物的含氢指数,用H表示。H与单位体积介质里的氢核数成正比,因而可用下式表示:式中, 是介质密度,单位为gcm3;M是该化合物的克分子量;x是该化合物每个分子中的氢原子数; K是比例常数。 纯水的含氢指数 按定义,淡水的含氢指数为1,由此确定出上式中
13、的K值。因水的分子式为H2O,所以x2,M18,而水的密度lgcm3,由此求出K9代入上式得则用该式可求出任何密度为,分子量为M且每个分子中有x个氢核的单一分子组成的物质的含氢指数。盐水的含氢指数 NaCl溶于水后占据了空间,而使盐水中氢的密度减小。计算盐水含氢指数的一般公式为 Hww(1-p)式中w为盐水的密度;p为NaCl的浓度。 在测裸眼井时,渗透性地层一般都有侵入,中子测井探测范围内的水的矿化度,可以认为是与泥浆滤液的矿化度接近的。油、气的含氢指数 液态烃的含氢指数与水接近,然而天然气的氢浓度很低,并且随温度和压力而变化。因而当天然气很靠近井眼而处于中子测井探测范围内时,中子测井测出的
14、含氢指数较小。 烃的含氢指数可根据其组分和密度来估算。分子式为CHx(其分子量为12十x)和密度为h的烃的含氢指数为 用此式可算出甲烷(CH4)的含氢指数为2.25甲烷,而石油(nCH2)的含氢指数为1.28油。若油的密度为油0.85g/cm3,则含氢指数为1.09。与有效孔隙度无关的含氢指数泥质:泥质伴生有化学结晶水和束缚水,所以它具有很大的含氢指数,一般可达0.150.30,因而在含泥质的地层中,含氢指数大于地层的有效孔隙度。石膏:石膏的分子式是CaSO42H2O,其含氢指数约为0.49,与孔隙度为49的灰岩相当。岩性影响:当仪器以纯石灰岩为标准进行刻度时,其它岩性的岩石骨架显示为一定数值
15、的等效含氢量。孔隙度等于零的砂岩,显示为负含氢指数(-3%),而白云岩显示为正的含氢指数(5%)。二、中子测井中子测井包括两种方法:v一种是记录探测器周围超热中子密度的中子超热中子测井(SNP),亦称井壁中子测井;v另一种是记录探测器周围热中子密度的中子热中子测井(CNL),亦称补偿中子测井。1、中子超热中子测井 超热中子测井是探测探测器周围中子变为热中子之前的超热中子密度,以反映地层的中子减速特性,进而计算储层孔隙度和对储集层进行评价。右图是一种超热中子测井仪的井下仪器示意图,也叫井壁超热中子测井仪(SNP)。超热中子探测器和中子源贴靠井壁测量以减小井眼的影响。由中子源发出的快中子与地层中的
16、原子核发生弹性散射,能量逐渐降低,而成为超热中子,其过程主要取决于上述所述的快中子的减速过程。1、中子超热中子测井 在组成地层的所有元素中,氢是减速能力最强的元素,远远超过其它元素,它的存在和含量就决定着地层的减速长度的大小。因此,当孔隙中100充满水时,孔隙度越大则地层减速长度就越短。 右图描述了这种关系,给出了充满水的砂岩、石灰岩和白云岩等三种岩性的岩石减速长度和孔隙度的关系曲线。由图可以看出Ls随增大而缩短,而且孔隙度相同、岩性不同的地层减速长度也不同。 因此,孔隙度不同,岩性不同,超热中子在中子源周围的分布不同。孔隙度越大,含氢量越多,减速长度Ls越小,则在源附近的超热中子越多。相反,
17、孔隙度越小,减速长度Ls越大,则在较远的空间形成有较多的超热中子。如果把探测器放在较远的地方,超热中子探测器的计数率,则对于孔隙度大的计数率低,孔隙度小的计数率高。如图1所示。 如果把探测器放在较近的地方接收记录超热中子的计数率则有相反的情况,即孔隙度大,计数率高,孔隙度小,计数率低。当探测器放在某一个位置时,计数率与孔隙度的大小无关。如图2所示探测器到源之间的距离叫源距,第一种情况的源距叫长源距,第二种情况的源距叫短源距,第三种情况叫零源距。约40cm在实际工作中,通常用的是长源距探测器。 所以,测井记录的超热中子的计数率越大,反映岩层的孔隙度越小,反之计数率越小,反映岩层的孔隙度越大。 在
18、不含有氢元素的地层中,超热中子读数随含氢量增高呈指数规律降低,在岩石中,含H量直接反映着孔隙度的大小,因此 lgN=a+b 其中b-仪器常数 a-与井径、源距等有关的参数 N-中子计数率 这就是利用超热中子测井可以测量岩层孔隙度的原理。 由于超热中子被元素俘获的截面非常小,所以超热中子的空间分布不受岩层含氯量的影响(即地层水矿化度的影响),所以能够较好地反映氢含量的多少,即较好地反映岩层孔隙度的大小。 由中子源发出的快中子在周围介质中减速成热中子,探测热中子密度的测井方法叫热中子测井。 补偿中子测井CNL是较好的一种热中子测井方法。 超热中子测井,仅反映地层对快中子的减速性质,而与地层对热中子
19、的俘获能力无关,因而有利于测定地层的含氢指数,贴靠井壁测量减少了井对测量结果的影响,能很好地测定地层的孔隙度。但也有不足之处,超热中子分布范围小,探测深度浅,源距小,井条件和贴井壁状态的变化都会影响测量的结果,例如当井径不规则,井壁不平整,探头推靠状态不合适,在井壁和探测器之间有泥饼或泥浆薄膜时,测得的孔隙度就偏大。因而在继续改进这仪方法的同时,还必须研究其它测井方法。 热中子的分布范围比超热中子大的多,探测范围大,其空间分布规律与超热中子的空间分布规律是一样的,即在长源距的情况下饱含流体的岩层的孔隙度越大,热中子的计数率越低;孔隙度越小计数率越高。2 中子热中子测井 快中子与地层作用减速成热
20、中子,探测器周围热中子的密度不仅与地层的减速特性有关,而且亦与地层的俘获特性有关。这就决定了热中子的空间分布既与岩层的含氢量有关,又与含氯量有关。这对于用热中子计数率大小反映岩层含氢量,进而反映岩层孔隙度值来说,氯含量就是个干扰因素。 在均匀无限介质中,对点状快中子源造成的热中子分布进行了理论推导,得到下列关系:式中 Nt热中子计数率; r探测器到中子源的眨离(源距); D扩散系数; Ls、Ld分别为减速长度和扩散长度; K与仪器有关的系数。 由上式可见计数率的大小不仅决定于岩层的减速性质(反映含氢量),还与岩层的俘获性质有关(反映含氯量)。2 中子热中子测井 若采用源距不同的两个探测器,记录
21、两个计数率Nt(r1)和Nt(r2),取这两个计数率比值,当源距r足够大时,则有 因热中子的扩散长度Ld比快中子的减速长度Ls小很多,所以当源距r足够大时,含有Ld的指数项与含有Ls的指数项相比可以忽略不计,如r=70cm,对=15%的砂岩 。 故上式可简化为: 当源距r1、r2选定后,这个比值只与地层的减速性质有关,所以比值能很好地反映地层的含氢量。该式即为双源距热中子测井的理论依据。 此外这种方法能减小井参数及岩石对热中子俘获性质对测量结果的影响,所以通常称为补偿中子测井。2 中子热中子测井 左图是补偿中子测井示意图,用长、短源距两个探测器接收热中子,得到两个计数率Nt(r1),Nt(r2
22、),根据用石灰岩刻度的仪器得到的计数率比值Nt(r1)Nt(r2)(r1r2)与岩石石灰岩孔隙度N的关系,补偿中子测井直接给出石灰岩孔隙度值曲线。 实际工作中,补偿中子测井仪通常都在刻度井内已知孔隙度的纯石灰岩地层上进行刻度,由此获得的孔隙度单位称为“石灰岩孔隙度”。 显然,它在纯石灰岩地层上等于地层的真孔隙度。 但在非纯石灰岩地层,用这种方式刻度的仪器测得的孔隙度将与地层的真孔隙度不同,称为“视石灰岩孔隙度”。3、中子测井的应用 确定地层孔隙度 中子测井(SNP、CNL)主要的用途是确定地层孔隙度。 右图SNP孔隙度岩性影响校正图版中子测井仪器通常以石灰岩孔隙度为标准刻度的,所以它所记录的孔
23、隙度是石灰岩孔隙度。因此,我们通常称中子测井值为视石灰岩孔隙度。 对于除石灰岩以外的其它岩性的岩石,石灰岩孔隙度包含有由于岩性不同,岩石骨架造成的附加孔隙度,例如孔隙度为零的纯砂岩和白云岩,用以石灰岩为标准刻度的仪器测量得到的石灰岩孔隙度,前者是3,后者为十5。所以要求砂岩或白云岩的孔隙度必须做岩性校正,用图所示的岩性影响校正图版进行校正。3、中子测井的应用 FDC与CNL石灰岩孔隙度曲线重叠定性判断气层 天然气使FDC测井计算孔隙度增大,而使CNL测井计算孔隙度偏小。故二者在气层上有一定的幅度差,而且DN。 右图是FDC与CNL的实测石灰岩孔隙度曲线重叠图。图中在三个深度上明显显示出含气层。
24、3、中子测井的应用 CNL与FDC测井交会求孔隙度、确定岩性 由密度测井(FDC)的体积密度值和CNL的石灰岩孔隙度值的交会点,可确定地层的孔隙度ND的大小和岩性。若是双矿物岩石,可以确定双矿物的比例。如图。 三、 中子伽马测井 1 基本原理 这种方法在于测量岩石中元素原子核俘获热中子之后所放出的俘获伽马射线的强度。 这一强度与两个因素有关,即岩石对中子的减速能力和对热中子的俘获能力。在沉积岩的元素中,对这两种特性起决定作用的是氢和氯,因此,中子伽马测井结果主要与岩石中的含氢量和含氯量有关。右图是实验得到的中子伽马计数率Nn-r与源距L的关系曲线。图中,曲线1淡水;曲线2盐水;曲线3孔隙度分别
25、为0、10、20,含淡水的模拟地层中的关系曲线。三、 中子伽马测井 Nn-r随L增大而按指数规律下降;L35cm为零源距。在此源距时,Nn-r与地层含氢量(孔隙度)无关,但仍能反映含氯量的变化。含氯量增大,Nn-r也增大;在长源距(即源距大于零源距)的情况下,含氢量增大(孔隙度增大),计数率Nn-r减小。 实际测井时,中子伽马测井采用的是长源距(即正源距)。国内通常采用L6065cm。所以中子伽马测井的计数率Nn-r大,说明地层的含氢量小,孔隙度小。若Nn-r小,则说明地层的含氢大,孔隙度大。因此,也可用中子伽马测井来确定岩石的孔隙度。 三、 中子伽马测井 对相同含氢量的岩石而言,如果含氯量不
26、同时,在含氯量高的岩石中,无论采用的源距如何,测得的中子伽马射线强度均有所增高。这是因为氯元素俘获热中子之后放出的伽马射线能量较高,且数量较多(能放出三个伽马量子,而一般元素俘获热中子后仅放出一个伽马量子)的缘故。从这一点看出,氯元素对中子伽马测井结果的影响,正好与氢对热中子测井的影响相反。 根据以上的分析看出,用中子伽马测井结果来定量研究岩石的孔隙度时,岩石中氯元素的存在,必将使计算结果产生较大的误差。因此,用这种方法来研究岩石的孔隙度不如上述两种方法理想。 中子伽马测井的下井仪器包括中子源和射线探测器,在源和探测器之间放有屏蔽体铅,防止中子源伴生射线由仪器内部直接进入探测器。中子伽马测井探
27、测深度略大于热中子(CNL)和超热中子测井(SNP)。2 中子-伽马测井曲线的应用中子伽马测井曲线的主要应用有:划分气层。中子伽马测井曲线可以用来划分气层。 天然气中的含氢量比水和油中的含氢量低得多,因此气层上得中子伽马测井显示出很高的计数率值。右图是中子伽马测井显示气层的实例。 在利用中子伽马曲线划分气层时,要注意和岩性致密造成的高幅度值相区别。为此,应参照声波和微电极等曲线综合分析。2 中子-伽马测井曲线的应用确定油水界面。 因为油水层的含氢量基本上是相同的,只有地层水的矿化度高时,水层的含氯量显著大于油层,油层和水层的中子伽马测井曲线的计数率值才有明显的差别(水层的中子伽马测井计数率值大
28、于油层的中子伽马测井计数率值),所以只有在地层水矿化度比较高的情况下,才能利用中子伽马测井曲线划分油水界面,区分油水层。 右图是利用中子伽马测井曲线划分油水界面的实例。v天然气中含氢量比水和油中含氢量低得多,因此气层中子测井显示比油、水层幅度高,计算孔隙度偏低。v高记数率与致密岩性高值区别在于声波(周波跳跃),微电极(是否分开等)。 v低矿化度地层中,油水在中子-r曲线上特征相同难以分开。 v高矿化度地层水中,由于cl作用,使中子-r能量高、吸收大、区分油水层、划分界面。仪器的刻度-视石灰岩刻度 岩性:几种常见岩石的中子曲线反映如下: 泥岩:泥岩总孔隙度大,含有大量吸附水和结晶水,在中子-中子
29、及中子-r测井曲线上均显示为低值。特别是在泥岩段井径常扩大,热中子密度和次生r射线强度进一步降低。泥岩在自然r曲线上显示为高值,两者结合,易把这和大孔隙度的渗透层区分开。 胶结很好致密砂岩、致密灰岩、白云岩和硬石膏,中子-中子和中子-r曲线上显示均为高值,而自然r曲线上为低值。 含有低矿化度地层水的渗透性砂岩和裂隙石灰岩:中子-中子和中子-r测井曲线上显示为中等数值。自然r曲线显示为低值,随孔隙度增大或含泥质增多,中子测井读数降低,自然r测井随泥质增多而升高。 孔隙中充满高矿化度地层水砂岩和石灰岩:中子-中子测井读数低中子-r读数高cl,自然r测井曲线显示为低值。 孔隙中充满天然气的砂岩和石灰
30、岩:中子-中子和中子-r测井曲线显高值,自然r测井显低值。 石膏层:中子-中子和中子-r测井曲线上显低值,自然r上为低值(含结晶水)岩盐层:中子-中子低值,中子-r高值。 中子测井确定地层的孔隙度 v中子测井读数是和岩石中含氢的总量有关,所以它反映的是总孔隙度。 、ma、f v中子孔隙度与含氢指数呈正比 N=(1- )(N)ma+ (N)f =(N-(N)ma)/(N)f-(N)ma)v流体(水)的中子孔隙度(N)f=1 v(N)ma对于不同岩石,不同中子测井法、有不同的值。 v若地层为含油气纯地层 N=(1- )(N)ma+(N)mf(1-Shr)+Shr(N)hr =(N-(N)ma)/(N)mf-(N)ma) -(Shr (N)hr-(N)mf)/(N)mf-(N)ma) =N/(1+Shr(N)hr-(N)mf)/(N)mf-(N)ma) N=(N-(N)ma)/(N)mf-(N)ma)为中子测井视孔隙度v有油气影响对,N减小,计算孔隙度低。 v利用上求时,对于油(N)hr=油-1.03 v对于气(N)hr=2.25气 v油和气分别表示油和气的密度 vv中子测井孔隙度 v*用中子测井计算的孔隙度是地层的含氢孔隙度或总孔隙度,当地层中有含氢量很高的石膏存在时,计算的孔隙度比实际值偏高,此时,需求石膏含量加以校正。二、中子-中子测井原理 1.中子-热中子测井 热中子密
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1