CN103674801B - 在压力脉冲衰减试验中减少不确定度的方法 - Google Patents

在压力脉冲衰减试验中减少不确定度的方法 Download PDF

Info

Publication number
CN103674801B
CN103674801B CN201310410037.7A CN201310410037A CN103674801B CN 103674801 B CN103674801 B CN 103674801B CN 201310410037 A CN201310410037 A CN 201310410037A CN 103674801 B CN103674801 B CN 103674801B
Authority
CN
China
Prior art keywords
coefficient
pressure pulse
rock sample
permeability
storage
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.)
Expired - Fee Related
Application number
CN201310410037.7A
Other languages
English (en)
Other versions
CN103674801A (zh
Inventor
宋寅先
金正赞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Korea Institute of Geoscience and Mineral Resources KIGAM
Original Assignee
Korea Institute of Geoscience and Mineral Resources KIGAM
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Korea Institute of Geoscience and Mineral Resources KIGAM filed Critical Korea Institute of Geoscience and Mineral Resources KIGAM
Publication of CN103674801A publication Critical patent/CN103674801A/zh
Application granted granted Critical
Publication of CN103674801B publication Critical patent/CN103674801B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/082Investigating permeability by forcing a fluid through a sample
    • G01N15/0826Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change

Landscapes

  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Dispersion Chemistry (AREA)
  • Fluid Mechanics (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明公开在压力脉冲衰减试验中减少不确定度的方法,能准确获得岩石样品的透水系数和单位储水系数。具体地,本发明的在压力脉冲衰减试验中减少不确定度的方法的特征在于,对岩石样品施加轴压和围压,在上述岩石样品连接上游段和下游段,在此状态下,从目标函数的图表读取表示等高线的最小值的坐标值,并将该坐标值决定为上述岩石样品的透水系数和单位储水系数,其中上述目标函数的图表中,在横轴和纵轴分别表示利用从外部施加压力脉冲的压力脉冲衰减试验法得出的透水系数和单位储水系数。并且,本发明的特征在于,将在不同的边界条件下重复进行上述压力脉冲衰减试验而获得的目标函数的图表重叠表示,能够减少岩石样品的透水系数和单位储水系数。

Description

在压力脉冲衰减试验中减少不确定度的方法
技术领域
本发明涉及一种在测定岩石样品的透水系数和单位储水系数时减少不确定度的方法。
背景技术
在开发地下水资源或石油资源时,为了防止各种填埋地的污染扩散,或者为了勘察深地层地下储藏用地的地下岩石环境,并且在进行其他各种土木施工等时,需要对地基进行勘察,尤其需要对基岩的透水性或储水性进行测定。
其中,基岩的透水性是指,通过基岩内的空隙,沿着相连接的流路产生的流体的易流动性。由于通过上述流路的流动受人力或压力差等的阻力,因此其透水性随着基岩或地层而不同。
并且,储水性是指,沿着通过基岩内的空隙来相连接的流路注入流体时,岩石对流体储藏能力的特性。如果没有施加于流体的势能或动能,流体将不会移动,保持静止状态。
一般而言,流体(代表性的是“水”)流入地下基岩层或地层(以下统称为“基岩层”)内或者从地下基岩层或地层内流出时,该流体的压力扩散主要受形成基岩的岩石或该地层的透水性或储水性的影响。
此时,岩石层的透水性或储水性是预测流体的流入/流出地点以及预测根据时间变化的流体变化(扩散)所必需的参数。
这些岩石层的透水性和储水性是在地下的自然状态下产生的,现实中我们不可能准确掌握实际透水性或储水性,一般都是通过实验室的水利试验来获取所需的数据。
测定这些透水性或储水性的方法有多种,形成岩石层的岩石致密、具有不透水性时,为了减少测定时间,多利用向测定用岩石样品施加压力脉冲的脉冲-衰减(过渡)法。
通过如上所述的实验室的水利试验,能够计算出透水系数(hydraulicconductivity)和储水系数(storage coefficient),这两个系数是表示形成目标岩石层的岩石样品的水利特性的代表性值。
其中,上述透水系数的定义为水通过的速度,即,通过的距离除以时间的值。该透水系数的值越小,就意味着流体(水)在地下基岩层的流动越难。
并且,上述储水系数的定义为,根据单位水头的变化,流入/流出单位面积的含水层的水量。
另一方面,单位储水系数(Specific storage,specific storativity)的定义为,根据含水层内的单位体积的单位水头的上升或下降,流入/流出该含水层的水量(体积),该单位储水系数乘以含水层厚度,即可求出储水系数。
以往的脉冲-衰减法将单位储水系数假设为0,只根据透水系数测定或者逆运算压力曲线,求出透水系数和单位储水系数。
利用现有技术来求出透水系数和单位储水系数的方法存在的第一个问题就是,完全没有提示各参数具有何种程度的不确定度(uncertainty)。
即,由于在实验室通过水利试验来测定的压力值肯定存在实验误差,因此最终求出的透水系数和单位储水系数均具有某种程度的不确定度,这是不可避免的,但由于以往的技术并未提供与此相关的数据,因此完全无法得知实验的可信度。
第二个问题是,因实验误差而产生的参数的不确定度因实验系统(上游段和下游段的大小)(上游段和下游段分别参照图1中的附图标记10和20)而异,无法获知上述参数是在何种条件下求出的。
第三个问题是,可能会出现根据水利试验的组成条件,各参数之间的不确定度相互产生影响的情况,在此情况下,不知道任一个参数的不确定度时,其他参数的不确定度也不得而知。
现有的透水系数和单位储水系数计算方法对如上所述的问题完全没有拿出解决方案。
另一方面,本发明相关的现有技术文献有以下的非专利文献(论文)。
(1)Brace,W.F.,J.B.Walsh,and W.T.Frangos,Permeability of graniteunder high pressure(花岗岩在高压下的透水性能),J.Geophys.(《地球物理学》)Res.1968732225-2236。
(2)Hsieh,P.A.,J.V.Tracy,C.E.Neuzil,J.D.Bredehoeft,and S.E.Silliman,A transient laboratory method for determining the hydraulicproperties of tight rocks(用于确定致密岩石的水硬性的瞬时实验方法)-I.Theory(理论I),Int.J.Rock Mech.Min.Sci.&Geomech.Abstr.198118245-252。
(3)Wang,H.F.,and D.J.Hart,Experimental error for permeability andspecific storage from pulse decay measurements(脉冲式衰变测量的透水性及特定存储的实验误差),Int.J.Mech.Min.Sci.Geomech.Abstr.199330,1173-1176。
(4)Zhang,M.,M.Takahashi,R.H.Morin,and T.Esaki,Evaluation andapplication of the transient-pulse technique for determining the hydraulicproperties of low-permeability rocks-Part2(确定低透水性岩石的水硬性的瞬时脉冲技术的演变及应用-第2部分):Experimental application(实验应用),Geotechnical Testing Journal(《地球技术测试期刊》)2000b;23,091-099。
发明内容
本发明的一个目的是,明确提示各参数具有的不确定度(uncertainty)相关数据,来提高实验的可靠度。
并且,本发明的再一个目的是,指出各参数是在何种实验条件下求出的,来明确因实验误差而导致的参数不确定度。
并且,本发明的另一个目的是,在各参数的不确定度相互产生影响时,同样明确指出各参数的不确定度。
根据旨在解决上述本发明目的的本发明一个优选实施例,在压力脉冲衰减试验中减少不确定度的方法,其特征在于,包括如下步骤:将岩石样品位于外罩内的套管的内部的步骤;在上述岩石样品的左右侧分别安装端塞和多孔盘来固定上述岩石样品的步骤;在上述固定的岩石样品的多孔盘的轴方向外侧分别连接上游段和下游段的步骤;在上述上游段和上述下游段分别附着测压传感器的步骤;对上述岩石样品施加轴压和围压的步骤,对位于上述套管的内部的上述岩石样品施加轴压和围压以形成封闭的压力系统;通过上述上游段,对上述封闭的压力系统的内部的上述岩石样品施加压力脉冲的步骤;以及求出透水系数和单位储水系数的步骤,通过上述压力脉冲被衰减而获得的压力脉冲的变动来求出透水系数和单位储水系数。
其中,在上述求出透水系数和单位储水系数的步骤中,从表示通过计算获得的透水系数和单位储水系数的计算值与通过上述测压传感器获得的上述压力脉冲的变动测定值之差的和/或其平均值的目标函数的图表,读取表示上述图表内的等高线的最小值的坐标值,将上述最小值作为上述岩石样品的透水系数和单位储水系数。
并且,上述套管的材质优选为橡胶。
并且,上述对岩石样品施加轴压和围压的步骤可通过向流体流入空间供给流体并施加压力的方法来实现。
并且,施加围压及轴压的上述流体优选为油(oil)。
根据旨在解决上述本发明目的的本发明的再一个优选实施例,在压力脉冲衰减试验中减少不确定度的方法,其特征在于,包括如下步骤:将岩石样品置于外罩的内部的套管内的步骤;在上述岩石样品的左右侧分别安装端塞和多孔盘来固定上述岩石样品的步骤;将固定的上述岩石样品与上游段和下游段相连接的步骤;在上述上游段和上述下游段分别附着测压传感器的步骤;分隔设置用于调节上述上游段的大小的两个阀的步骤;对岩石样品施加轴压和围压的步骤,对位于上述套管内的上述岩石样品施加轴压和围压来形成封闭的压力系统;施加压力脉冲的步骤,通过上述上游段,对上述封闭的压力系统内的上述岩石样品施加压力脉冲;以及求出透水系数和单位储水系数的步骤,通过上述压力脉冲被衰减而获得的压力脉冲的变动来求出透水系数和单位储水系数。
其中,优选地,在上述求出透水系数和单位储水系数的步骤中,从表示通过计算获得的透水系数和单位储水系数的计算值与通过上述测压传感器获得的上述压力脉冲的变动测定值之差的和或其平均值的目标函数的图表,读取表示上述图表内的等高线的最小值的坐标值,将上述最小值决定为上述岩石样品的透水系数和单位储水系数。
并且,优选地,在上述施加压力脉冲的步骤中,分别控制上述分隔的两个阀,在上述上游段的大小得到调节的状态下,通过上述上游段施加压力脉冲。
其中,通过重叠表示出多次重复上述施加压力脉冲的步骤而获得的目标函数的图表,可减少上述岩石样品的上述透水系数和单位储水系数的不确定度。
并且,上述套管的材质优选为橡胶。
并且,上述对岩石样品施加轴压和围压的步骤可通过向流体流入空间供给流体并施加压力的方法来实现。
并且,施加围压及轴压的上述流体优选为油。
由于本发明可提示因实验误差而导致的透水系数和单位储水系数的不确定度,因此可以评价实验的可靠度。
并且,由于本发明的实验条件明确,因此可以在进行实验前实施简单的模拟演练,以此来设计出最佳的实验系统。由此,可提示各参数间的不确定度相互不产生影响的边界条件,来减少不确定度。
并且,根据本发明,由于在水利试验中各参数之间的不确定度相互产生影响的情况下,也可以通过其他实验来限制各参数中一个参数的不确定度,因此可以获得更加准确的水利试验结果。
附图说明
图1是为了测定岩石样品的透水系数和单位储水系数而在本发明的压力脉冲-衰减试验中使用的三轴-岩芯夹持装置的简要剖视图。
图2是同时表示目标函数的图表和相关测定误差(Error)的图,其中测定误差表示与压力脉冲大小的比例。
图3是同时表示相互产生影响的目标函数的图表和相关测定误差的图。
图4是表示将总共经过三次水利试验获得的各个目标函数重叠时的目标函数的图表中的透水系数和单位储水系数的不确定度的图。
附图标记的说明
500:三轴-岩芯夹持装置
2:上侧部夹持器维持部  2':下侧部夹持器维持部
4:外罩                6:O型密封圈
8:上侧部夹持器        8':下侧部夹持器
10:上游段             20:下游段
12:上侧部端塞         22:下侧部端塞
30,35:多孔盘          40:套管
50:本体部流体供给口   60:上侧部流体供给口
70:下游侧测压器       80上游侧测压器
90、95:阀             100:岩石样品
P:压力脉冲            a、b:(流体流入)空间
具体实施方式
以下参照附图详细说明的实施例会让本发明的优点和特征以及实现这些优点和特征的方法更加明确。但是,本发明不局限于以下所公开的实施例,能够以互不相同的各种方式实施,本实施例只用于使本发明的公开内容更加完整,有助于本发明所属技术领域的普通技术人员完整地理解本发明的范畴,本发明根据权利要求书的范畴而定义。
以下,参照附图,对本发明的优选实施例的在压力脉冲衰减试验中减少不确定度的方法进行详细说明。
首先,对本发明的在压力脉冲衰减试验中减少不确定度的方法所使用的水利试验设备的结构进行说明。
图1是三轴-岩芯夹持装置的简要剖视图。
上述三轴-岩芯夹持装置500大分为上部、本体部及下部。在图中,左侧部分表示下部,右侧部分表示上部。
本发明的压力脉冲衰减试验中使用的岩石样品100位于装置500外罩4的套管40的内部,如图1所示,优选地,套管40以两侧稍微拉开的气缸形态形成。
上述套管40的材质优选为橡胶,但如果在施加后述的围压(confiningpressure)时,只要是该围压可以适当地传达到岩石样品100的结构,就可以使用其他材料。
上述岩石样品100在本体部的套管40内部,依次被按上下部轴方向配置的上侧部端塞(end plug)12和下侧部端塞22第一次固定;在上侧部端塞12和下侧部端塞22的上下部轴方向外侧,被多孔盘30、35第二次固定;在该多孔盘30、35的上下部轴方向外侧,被上侧部夹持器8和下侧部夹持器8'最终固定。
其中,上下部轴方向分别意味着图中的右侧及左侧方向。
并且,考虑到因流体而导致的腐蚀或耐久性问题,上侧部端塞12和下侧部端塞22的材质优选为不锈钢。
在三轴-岩芯夹持装置500的上部,相隔规定间距形成有从外部传达压力脉冲P的上游段10、用于测定上述上游段10的内部的压力脉冲P的上游侧测压器80、用于调节(控制)上述上游段10的大小的两个阀90、95。
并且,在装置的上部,分别形成有上侧部夹持器8、上侧部夹持器维持部2、上侧部流体供给口60,上述上侧部夹持器8在上侧部端塞12的轴方向外侧,重新固定用于固定本体部内的岩石样品100的外侧的多孔盘35的上侧部端塞12,上述上侧部流体供给口60围绕上述上侧部夹持器8,上述上侧部流体供给口60与上侧部夹持器8和上侧部夹持器维持部2之间的空间b相连接。
其中,上游段10和下游段20的位置可倒置,上述两个阀90、95可调节下游段20的大小,而不是调节上游段10的大小。
为了向三轴-岩芯夹持装置500施加轴压(axial pressure),可将岩石样品100位于本体部之后,通过上述上侧部流体供给口60,向流体流入空间b供给具有规定粘度的流体,该流体优选为油,在供给的油压作用下,上述上侧部夹持器8可向图中的左侧方向被推开的同时向岩石样品100施加所需的轴压。
此时,优选地,在上述上侧部流体供给口60的外部形成有能够对油进行加压来供给的油加压部(未图示)。
本发明的三轴-岩芯夹持装置500的本体部中配置有:套管40,内部置有岩石样品100;多孔盘30、35,在上述套管40的内部用于从两侧支撑上述岩石样品100;上侧部端塞12和下侧部端塞22,分别支撑上述多孔盘30、35;流体流入空间a,位于上述套管40的外侧;以及本体部流体供给口50、50,向上述流体流入空间a供给流体来赋予围压。
其中,优选地,用于形成上述围压的流体为与用于形成轴压的具有规定粘度的油相同的油。优选地,如同上侧部流体供给口60一样,其外侧形成有油加压部(未图示)。
本发明的三轴-岩芯夹持装置500的下部配置有:下游段20,具有与上部侧的上游段10相对应的结构;下游侧测压器70,附着于上述下游段20的下侧部的端部,用于测定三轴-岩芯夹持装置500内的压力;下侧部夹持器8',用于固定并维持形成于岩石样品100的下部的轴方向外侧的下侧部端塞22;以及下侧部夹持器维持部2',围绕下侧部夹持器8'。
值得注意的是,就三轴-岩芯夹持装置500的上部而言,在上侧部夹持器8和上侧部夹持器维持部2之间形成有连通的流体流入空间b,就三轴-岩芯夹持装置500的下部而言,在下侧部夹持器8'和下侧部夹持器维持部2'之间并没有形成有连通空间。
另一方面,优选地,上述三轴-岩芯夹持装置500的上侧部夹持器8和上侧部夹持器维持部2之间、上侧部夹持器维持部2和外罩4之间以及外罩4和下侧部夹持器维持部2'之间设置有适当数量的流体密封用O型密封圈(o-ring)6、6、6、6。
并且,优选地,上侧部夹持器维持部2和外罩4、外罩4和下侧部夹持器维持部2'以及下侧部夹持器8'和下侧部夹持器维持部2'分别进行螺丝结合。尤其是,下侧部夹持器8'和下侧部夹持器维持部2'的螺丝结合执行适当调节施加于岩石样品100的流体的轴压的功能。
即,如上所述,将岩石样品100位于套管40的内部,安装完三轴-岩芯夹持装置500之后,从上侧部流体供给口60供给流体,以此可形成规定的轴压,但为了形成所需的轴压,最终可通过下侧部夹持器8'和下侧部夹持器维持部2'之间的螺丝结合来进行调节。
其次,对向本发明的三轴-岩芯夹持装置500赋予围压(confiningpressure)的结构进行说明。
将岩石样品100位于套管40的内部,安装完三轴-岩芯夹持装置500之后,通过本体部流体供给口50、50供给流体,如此供给的流体在气缸形状的套管40的同心圆方向的外侧向套管40实施加压,并向岩石样品100供给最终的围压。
其中,优选地,围压和轴压的形成顺序几乎同时进行。
并且,优选地,用于形成上述围压的流体为与用于形成轴压的具有规定粘度的油相同的油。优选地,如同上侧部流体供给口60一样,其外侧形成有油加压部(未图示)。
如上所示,以下参照在图1所示的三轴-岩芯夹持装置500固定岩石样品100的结构,对本发明的在压力脉冲衰减试验中减少不确定度的方法进行说明。
在三轴-岩芯夹持装置500固定岩石样品100之后,在岩石样品100的上游侧和下游侧的两侧分别连接上游段10和下游段20,将上述上游段10与泵(未图示)相连接,来施加压力脉冲P。
接着,如上所述,向岩石样品100施加轴压和围压,来形成封闭的压力系统,并向上游段10施加由上述泵生成的压力脉冲P。
随着时间的经过,上游段10和下游段20的压力脉冲将出现变化(衰减),通过上游侧测压器80和下游侧测压器70来得出表示脉冲的压力变化的曲线之后,对其进行逆运算,来求出岩石样品100的透水系数和单位储水系数。
岩石样品100的透水系数和单位储水系数不是通过压力变化曲线来直接测定的,而是利用测定值和理论值的曲线拟合(curve fitting)法来求出,测定值和计算值的差异最小时,将预测的两个参数定为该岩石样品100的参数。
上游段和下游段的压力(pu和pd)的理论性变化曲线可依据以下数学式1来得出。
数学式1
p u P = 1 1 + α + αβ + 2 Σ m = 1 ∞ ( cos φ m - βφ m sin φ m ) exp ( - φ m 2 τ ) ( 1 + α + αβ - βφ m 2 ) cos φ m - φ m ( 1 + αβ + 2 β ) sin φ m
p d P = 1 1 + α + αβ + 2 Σ m = 1 ∞ exp ( - φ m 2 τ ) ( 1 + α + αβ - βφ m 2 ) cos φ m - φ m ( 1 + αβ + 2 β ) sin φ m
上述数学式1是以脉冲压力P正常化的公式。
其中,α=SsAL/Su,β=Sd/SsAL,τ=kt/(SsL)。
A和L分别为试样的截面积及长度,Su和Sd分别为上游段和下游段的压缩存储常数(compressive storage),Ss为单位储水系数,k为透水系数,t为从发生脉冲的时刻开始计算的时间,Φm为以下数学式2的解。
数学式2
tan φ = ( 1 + αβ ) φ β φ 2 - α
通过如上数学式求出的计算值和本发明的水利试验测定值之差的和(或平均值)被称为目标函数(Objective function),如图2、图3所示,可以用横轴显示的透水系数和纵轴显示的单位储水系数的函数(图表)来表示。
另一方面,求出透水系数和单位储水系数包含目标函数的最小值的范围内的目标函数之后,在等高线读取该目标函数最小时的坐标值(最小值),将其定为作为岩石样品100的测定值的透水系数和单位储水系数。
离上述最小值的坐标越远,目标函数的值越大,此时目标函数的值表示压力测定值的误差(Error)。
在图2的右侧,以条形图表示了该误差。
在这里,测定误差越大,即表示误差的条形图越靠近上侧,各参数的不确定度越大。
利用该目标函数的图表,就可以从压力测定值得知透水系数和单位储水系数的不确定度。
具体地,图2的测定误差的平均值为压力脉冲的2%时,在目标函数的图表中央可以用亮的灰色部分来表示。
在此情况下,透水系数(Permeability)的不确定度约为-20~60%,单位储水系数(Specific storage)的不确定度约为-30~60%。
简要整理图2内容,本发明的一个结构如下:向岩石样品施加轴压和围压,将上述岩石样品与上游段和下游段相连接,在此状态下,从目标函数的图表中读取表示等高线的最小值的坐标值,并将该坐标值决定为上述岩石样品的透水系数和单位储水系数,其中上述目标函数的图表中,在横轴和纵轴分别表示利用从外部施加压力脉冲的压力脉冲衰减试验法来得出的透水系数和单位储水系数。
其次,图3是表示显示相互产生影响的不确定度的目标函数的图。
图3表示不确定度相互产生影响的情况,如图所示,表示试验误差的等高线呈椭圆形、稍微倾斜的形状。
在此情况下,相比于已知其中一个参数的不确定度的情况,未知任一个参数的不确定度时,其他参数的不确定度将增大不少。
例如,单位储水系数的不确定度为0时,透水系数的不确定度约为20%,但未知单位储水系数时,透水系数的不确定度将大幅增大。
在此情况下,由于重叠表示图2和图3的目标函数,就会限制单位储水系数的不确定度,因此透水系数的不确定度也将被严格限制,可以减少两个参数的不确定度。
对此,将参照图4进行说明。
为了获得重叠图2和图3的目标函数的图表,需要两个以上的具有不同的边界条件的试验。
为此,利用形成于图1所示的三轴-岩芯夹持装置500的上游段10的两个阀90、95,将上述上游段10调成不同的大小,在不同的边界条件下多次重复进行实验之后,将该实验结果重叠起来进行分析,就可以显著减少不确定度。
即,相比于只进行一次实验的情况,在不同的边界条件下多次进行一系列实验,并将通过每次实验来得出的不确定度重叠之后表示交集合,就可以显著减少不确定度。
按如上所述的方法,相比于其他实验,在如图3所示的不确定度相互产生影响的情况下进行的实验中也可限制某一个参数的不确定度,由此可降低整个系统的测定不确定度。
如图4所示,相比于只使用一种水利试验结果的情况,使用图表中央用四边形表示总共经过三次的水利试验获得的目标函数重叠时的目标函数图表中的透水系数和单位储水系数的不确定度时,可大幅降低不确定度。
简要整理图3及图4的内容,本发明的再一个结构如下:将在不同的边界条件下重复进行压力脉冲衰减试验而获得的目标函数的图表重叠起来表示交集,由此就可以减少岩石样品的透水系数和单位储水系数的不确定度。
如上,举例对本发明的优选实施例的在压力脉冲衰减试验中减少不确定度的方法进行了说明,但这仅用于例示,对于本发明所属技术领域的普通技术人员来说,由此可实施本发明的各种变形,或者可实施与本发明等同的技术方案。

Claims (10)

1.一种在压力脉冲衰减试验中减少不确定度的方法,其特征在于,包括如下步骤:
将岩石样品置于外罩内的套管的内部的步骤;
在所述岩石样品的左右侧分别安装端塞和多孔盘来固定所述岩石样品的步骤;
在固定的所述岩石样品的多孔盘的轴方向外侧分别连接上游段和下游段的步骤;
在所述上游段和所述下游段分别附着测压传感器的步骤;
对位于所述套管的内部的所述岩石样品施加轴压和围压以形成封闭的压力系统的步骤;
通过所述上游段对所述封闭的压力系统的内部的所述岩石样品施加压力脉冲的步骤;以及
通过所述压力脉冲被衰减而获得的压力脉冲的变动来求出透水系数和单位储水系数的步骤,
其中,在所述求出透水系数和单位储水系数的步骤中,从表示通过计算获得的透水系数和单位储水系数的计算值与通过所述测压传感器获得的所述压力脉冲的变动测定值之差的和、或所述差的平均值的目标函数的图表,读取表示所述图表内的等高线的最小值的坐标值,将所述最小值决定为所述岩石样品的透水系数和单位储水系数。
2.根据权利要求1所述的在压力脉冲衰减试验中减少不确定度的方法,其特征在于,所述套管的材质为橡胶。
3.根据权利要求1所述的在压力脉冲衰减试验中减少不确定度的方法,其特征在于,在所述对岩石样品施加轴压和围压的步骤中,向流体流入空间供给流体并施加压力。
4.根据权利要求3所述的在压力脉冲衰减试验中减少不确定度的方法,其特征在于,所述流体为油。
5.一种在压力脉冲衰减试验中减少不确定度的方法,其特征在于,包括如下步骤:
将岩石样品置于外罩的内部的套管内的步骤;
在所述岩石样品的左右侧分别安装端塞和多孔盘来固定所述岩石样品的步骤;
将固定的所述岩石样品与上游段和下游段相连接的步骤;
在所述上游段和所述下游段分别附着测压传感器的步骤;
分隔设置用于调节所述上游段的大小的两个阀的步骤;
对位于所述套管内的所述岩石样品施加轴压和围压来形成封闭的压力系统的步骤;
通过所述上游段对所述封闭的压力系统内的所述岩石样品施加压力脉冲的步骤;以及
通过所述压力脉冲被衰减而获得的压力脉冲的变动来求出透水系数和单位储水系数的步骤,
其中,在所述求出透水系数和单位储水系数的步骤中,从表示通过计算获得的透水系数和单位储水系数的计算值与通过所述测压传感器获得的所述压力脉冲的变动测定值之差的和、或所述差的平均值的目标函数的图表,读取表示所述图表内的等高线的最小值的坐标值,将所述最小值决定为所述岩石样品的透水系数和单位储水系数。
6.根据权利要求5所述的在压力脉冲衰减试验中减少不确定度的方法,其特征在于,所述套管的材质为橡胶。
7.根据权利要求5所述的在压力脉冲衰减试验中减少不确定度的方法,其特征在于,所述施加压力脉冲的步骤中,分别控制所述分隔的两个阀,在所述上游段的大小得到调节的状态下,通过所述上游段施加压力脉冲。
8.根据权利要求7所述的在压力脉冲衰减试验中减少不确定度的方法,其特征在于,通过重叠表示出多次重复所述施加压力脉冲的步骤而获得的目标函数的图表,减少所述岩石样品的所述透水系数和单位储水系数的不确定度。
9.根据权利要求5所述的在压力脉冲衰减试验中减少不确定度的方法,其特征在于,在所述对岩石样品施加轴压和围压的步骤中,向流体流入空间供给流体来施加压力。
10.根据权利要求9所述的在压力脉冲衰减试验中减少不确定度的方法,其特征在于,所述流体为油。
CN201310410037.7A 2012-09-11 2013-09-10 在压力脉冲衰减试验中减少不确定度的方法 Expired - Fee Related CN103674801B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR10-2012-0100488 2012-09-11
KR1020120100488A KR101245333B1 (ko) 2012-09-11 2012-09-11 압력 펄스 붕괴 시험에서의 불확도 감소 기법

Publications (2)

Publication Number Publication Date
CN103674801A CN103674801A (zh) 2014-03-26
CN103674801B true CN103674801B (zh) 2015-10-07

Family

ID=48182202

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310410037.7A Expired - Fee Related CN103674801B (zh) 2012-09-11 2013-09-10 在压力脉冲衰减试验中减少不确定度的方法

Country Status (4)

Country Link
US (1) US9335245B2 (zh)
JP (1) JP5544443B2 (zh)
KR (1) KR101245333B1 (zh)
CN (1) CN103674801B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016179593A1 (en) * 2015-05-07 2016-11-10 The Uab Research Foundation Full immersion pressure-pulse decay
CN104914028B (zh) * 2015-05-08 2018-02-09 山东科技大学 假三轴高压渗透试验仪
CN105424575B (zh) * 2015-11-10 2019-06-11 中国石油天然气股份有限公司 一种确定渗透率的方法和装置
KR101800796B1 (ko) * 2016-08-16 2017-11-27 한국지질자원연구원 극저투수율암석의 투수율 측정 장치 및 그 방법
KR101999279B1 (ko) * 2017-06-13 2019-07-11 이성종 탄소섬유를 이용한 x-ray CT촬영용 압력장치
US10845292B2 (en) * 2017-06-19 2020-11-24 Saudi Arabian Oil Company Method for correcting low permeability laboratory measurements for leaks
US11079313B2 (en) * 2019-05-17 2021-08-03 Saudi Arabian Oil Company Methods and systems for determining core permeability pulse decay experiments
CN110823779B (zh) * 2019-11-19 2022-09-16 西南石油大学 一种岩心夹持器和组合式岩心夹持模型
US12050168B2 (en) 2022-04-26 2024-07-30 Saudi Arabian Oil Company Method for determining a matrix permeability of a subsurface formation

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2188205Y (zh) * 1993-10-22 1995-01-25 西南石油学院 岩石物性参数测试装置
JP2000009631A (ja) * 1998-06-19 2000-01-14 Agency Of Ind Science & Technol 岩石試験体の透水試験装置
US6178808B1 (en) * 1999-08-12 2001-01-30 Wisconsin Alumni Research Foundation Apparatus and method for testing the hydraulic conductivity of geologic materials
CN201532351U (zh) * 2009-11-13 2010-07-21 中国科学院武汉岩土力学研究所 一种利用变容积压力脉冲法测试岩石气体渗透系数的装置

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4552011A (en) * 1982-09-28 1985-11-12 Phillips Petroleum Company Establishing a differential pressure across a core sample or calibrating a differential pressure transducer
US4506542A (en) * 1983-04-22 1985-03-26 Chandler Engineering Company Apparatus and procedure for relative permeability measurements
US4571985A (en) * 1983-11-17 1986-02-25 The United States Army Corps Of Engineers As Represented By The Secretary Of The Army Method and apparatus for measuring the hydraulic conductivity of porous materials
US5417104A (en) * 1993-05-28 1995-05-23 Gas Research Institute Determination of permeability of porous media by streaming potential and electro-osmotic coefficients
US5813798A (en) * 1997-03-28 1998-09-29 Whiffen; Greg Piecewise continuous control of groundwater remediation
JP3503927B2 (ja) * 1998-09-07 2004-03-08 財団法人電力中央研究所 三軸セル、三軸試験装置および三軸試験方法
US6321317B1 (en) * 1998-12-16 2001-11-20 Hewlett-Packard Co Apparatus for and method of multi-dimensional constraint optimization in storage system configuration
US6415648B1 (en) * 1999-02-18 2002-07-09 Colorado School Of Mines Method for measuring reservoir permeability using slow compressional waves
WO2003036330A2 (en) * 2001-10-22 2003-05-01 Coppola Emery A Jr Neural network based predication and optimization for groundwater / surface water system
US6986281B1 (en) * 2004-06-30 2006-01-17 Battelle Energy Alliance, Llc Exfiltrometer apparatus and method for measuring unsaturated hydrologic properties in soil
FI20045353A (fi) 2004-09-24 2006-03-25 Sandvik Tamrock Oy Menetelmä kiven rikkomiseksi
US7882726B2 (en) * 2007-05-08 2011-02-08 Porous Materials, Inc. Compression vacuapore for determination of pore structure characteristics of hydrophobic materials under compressive stress
US20100089124A1 (en) * 2008-09-26 2010-04-15 North Dakota State University Integrated porous rigid wall and flexible wall permeability test device for soils
US20110067857A1 (en) * 2009-09-23 2011-03-24 Schlumberger Technology Corporation Determining properties of a subterranean structure during hydraulic fracturing
US8806954B1 (en) * 2010-11-16 2014-08-19 U.S. Department Of Energy Steady state method to determine unsaturated hydraulic conductivity at the ambient water potential

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2188205Y (zh) * 1993-10-22 1995-01-25 西南石油学院 岩石物性参数测试装置
JP2000009631A (ja) * 1998-06-19 2000-01-14 Agency Of Ind Science & Technol 岩石試験体の透水試験装置
US6178808B1 (en) * 1999-08-12 2001-01-30 Wisconsin Alumni Research Foundation Apparatus and method for testing the hydraulic conductivity of geologic materials
CN201532351U (zh) * 2009-11-13 2010-07-21 中国科学院武汉岩土力学研究所 一种利用变容积压力脉冲法测试岩石气体渗透系数的装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
徳永朋祥等.流体圧力&#12539 *
歪の同時計測による透水係数および貯留係数の 求め方にっいて.《応用地質》.2001,第42卷(第4期),第208-213页. *
瞬态压力脉冲法及其在岩石三轴试验中的应用;李小春等;《岩石力学与工程学报》;20011031;第20卷;1725-1733页 *

Also Published As

Publication number Publication date
US20140069183A1 (en) 2014-03-13
JP5544443B2 (ja) 2014-07-09
KR101245333B1 (ko) 2013-03-19
CN103674801A (zh) 2014-03-26
JP2014055500A (ja) 2014-03-27
US9335245B2 (en) 2016-05-10

Similar Documents

Publication Publication Date Title
CN103674801B (zh) 在压力脉冲衰减试验中减少不确定度的方法
US10571384B2 (en) Methods and systems for determining gas permeability of a subsurface formation
US10416064B2 (en) Methods and systems for determining gas permeability of a subsurface formation
Wang et al. Gas permeability and porosity evolution of a porous sandstone under repeated loading and unloading conditions
CN103558136B (zh) 温度应力环向渗流耦合作用下岩石损伤与渗透测试系统和测试方法
Huy et al. Carbon dioxide gas permeability of coal core samples and estimation of fracture aperture width
CN103616322B (zh) 低渗透岩石渗透率非稳态测定方法
US10732086B2 (en) Device and method for measuring magnitude of seepage force and its influence on effective stress of formation
CN104034644B (zh) 一种可快速测量孔隙率的多相渗流介质三轴应力渗流耦合试验装置
KR100953650B1 (ko) 지반의 퇴적특성과 하중이력을 고려한 역학적 성능평가를 위한 일체형 시료 조성장치
CN102012347B (zh) 一种水泥静胶凝强度测窜方法
CN108119132B (zh) 致密砂岩气藏近井带径向渗流含水饱和度模拟装置及方法
CN203630018U (zh) 低渗透岩石渗透率非稳态测定装置
Sato et al. Groundwater hydraulics
Zhou et al. Delayed drainage of aquitard in response to sudden change in groundwater level in adjacent confined aquifer: Analytical and experimental studies
CN107831103B (zh) 一种压力脉冲衰减气测渗透率测试装置的精度评估方法
CN103926184A (zh) 岩心气测孔隙度检测方法及其检测装置
CN106769684B (zh) 页岩气扩散能力测试系统
CN205484324U (zh) 一种固结渗透联合实验装置
Jang et al. Novel apparatus to measure the low-permeability and porosity in tight gas reservoir
Uehara et al. Evolution of permeability and microstructure of experimentally-created shear zones in Neogene siliceous mudstones from Horonobe, Japan
Yang et al. A new measurement method for radial permeability and porosity of shale
Skoczylas Determining the gas permeability coefficient of a porous medium by means of the bubble-counting flow meter
Sharp Jr et al. An electronic minipermeameter for use in the field and laboratory
Jang et al. Experimental apparatus and method to investigate permeability and porosity of shale matrix from Haenam Basin in Korea

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151007

Termination date: 20200910

CF01 Termination of patent right due to non-payment of annual fee