CN103926267B - 一种定量评价应力敏感过程中孔喉变化程度的方法 - Google Patents

一种定量评价应力敏感过程中孔喉变化程度的方法 Download PDF

Info

Publication number
CN103926267B
CN103926267B CN201410174667.3A CN201410174667A CN103926267B CN 103926267 B CN103926267 B CN 103926267B CN 201410174667 A CN201410174667 A CN 201410174667A CN 103926267 B CN103926267 B CN 103926267B
Authority
CN
China
Prior art keywords
pore throat
core
formation water
rock core
variation
Prior art date
Application number
CN201410174667.3A
Other languages
English (en)
Other versions
CN103926267A (zh
Inventor
高辉
王雅楠
蔡文斌
高永利
张荣军
Original Assignee
西安石油大学
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 西安石油大学 filed Critical 西安石油大学
Priority to CN201410174667.3A priority Critical patent/CN103926267B/zh
Publication of CN103926267A publication Critical patent/CN103926267A/zh
Application granted granted Critical
Publication of CN103926267B publication Critical patent/CN103926267B/zh

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Abstract

本发明公开了一种定量评价应力敏感过程中孔喉变化程度的方法,包括:取实验用岩心,洗油烘干,测岩心孔隙度和气测渗透率;配制实验用模拟地层水,达到地层水矿化度;选取模拟地层水驱替速度,将岩心置于岩心夹持器中,加围压,用模拟地层水驱替岩心;根据岩心入口压力和驱替速度计算水测渗透率,并测该围压下的核磁共振T2谱;逐渐升高围压,测不同围压下的核磁共振T2谱;不同围压下的核磁共振T2谱换算成孔喉半径r;绘制不同围压下的孔喉半径分布曲线图,定量得出不同围压下的孔喉变化程度。该方法可实时在线测试,无所多次卸装样品,避免卸装过程中围压变化对孔喉变化的影响;从而实现对应力敏感过程中孔喉变化程度定量评价。

Description

一种定量评价应力敏感过程中孔喉变化程度的方法
技术领域
[0001 ]本发明涉及油气开发实验技术领域,特别涉及一种定量评价应力敏感过程中孔喉 变化程度的方法。
背景技术
[0002] 油气藏开发过程中,孔隙压力不断减小,而上覆岩石压力没有变化,导致有效应力 增大,孔喉不断减小,这一特征在低渗透油气藏表现的尤为明显,明确应力敏感过程中孔喉 的变化程度,对于减小应力敏感程度,采取相应措施,进而指导油田合理开发具有重要的参 考价值。现有研究中,专利CN103528934A,一种测量超低渗岩石渗透率应力敏感性的互相关 技术。CN202066847U,疏松砂岩储层应力敏感模拟测试装置及其专用岩心夹持器。 CN102721636A,一种基于尺度模型确定裂缝性低渗储层应力敏感性的方法。CN102011581A, 应力敏感性地层钻井堵漏模拟评价装置。2002年第29卷第6期,成都理工学院学报,向阳等 人在《致密砂岩气藏应力敏感的全模拟试验研究》一文中利用高温高压全模拟全直径致密 砂岩气藏岩石降压开采试验数据揭示了致密砂岩气藏高速开采可对储层造成渗透率损害。 2007年第27卷第3期,天然气工业,康毅力等人在《致密砂岩微观孔隙结构参数对有效应力 变化的响应》一文中采用有效应力增加单面进汞测试技术分析致密砂岩微观孔隙结构特 征。2007年第28卷第4期,石油学报,于忠良等人在《致密储层应力敏感性及其对油田开发的 影响》一文中对致密储层的渗透率及应力敏感性进行了实验研究,建立了以原始储层有效 应力为有效应力起点的压敏评价方法。2011年第32卷第3期,焦春艳等人在《超低渗透砂岩 储层应力敏感性实验》一文中采用全自动岩心驱替系统,应用双重有效应力理论对鄂尔多 斯盆三叠系超低渗露头全直径岩心进行应力敏感性评价。2012年第19卷第2期,油气地质与 采收率,刘丽在《低渗透油藏启动压力梯度的应力敏感性实验研究》一文中借助室内岩心流 动实验,模拟再现地层压力的下降过程,研究了低渗透油藏有效上覆压力对岩石启动压力 梯度的影响。2013年第40卷第2期,石油勘探与开发,赵伦等人在《异常高压碳酸盐岩油藏应 力敏感实验评价一以滨里海盆地肯基亚克裂缝-孔隙型低渗透碳酸盐岩油藏为例》一文中 采取气测法,先逐级增加人造岩心样品围压,然后逐级降压,稳定后测定每个压力点的孔隙 度和渗透率,分析岩心的应力敏感性。2013年第12卷第3期,岩石力学与工程学报,孙军昌等 人在《特低渗储层不同渗流介质应力敏感特征及其评价方法研究》一文中选取大庆油田外 围特低渗油藏及长庆油田某露头储层砂岩岩石进行不同渗流介质(氮气、盐水、煤油)渗透 率应力性试验,分析岩样渗透率、孔喉变形、流体压缩性及流固耦合作用机制等因素对应力 敏感性的影响。
[0003] 现有研究主要是通过应力敏感过程中渗透率的变化或结合压汞、铸体薄片、扫描 电镜等手段去分析应力敏感前后的孔喉变化,存在的主要问题是:渗透率只是微观孔喉的 宏观反应,推测只能是定性评价;而应力敏感前后的压汞、铸体薄片、扫描电镜分析并不是 在相同的样品上完成,所得结论有待商榷。而核磁共振能够实现无损测试,可以在实验过程 中任何阶段进行测试,能够保证整个过程在相同样品上完成,而且核磁共振T2谱与孔喉半 径存在-对应关系,能够实现定量评价。
发明内容
[0004] 本发明的目的是提供一种定量评价应力敏感过程中孔喉变化程度的方法,该方法 通过实验,可始终保持地层条件,实时在线测试,无所多次卸装样品,避免卸装过程中围压 变化对孔喉变化的影响;从而实现对应力敏感过程中孔喉变化程度定量评价,克服定性分 析的不足。
[0005] 本发明的目的是通过下述技术方案来实现的。
[0006] -种定量评价应力敏感过程中孔喉变化程度的方法,包括下述步骤:
[0007] 步骤一、从全直径岩心上钻取直径2.5cm、长度介于4-6cm的岩心,洗油烘干,测岩 心孔隙度和气测渗透率;
[0008] 步骤二、配制实验用模拟地层水,达到地层水矿化度;
[0009] 步骤三、根据岩心气测渗透率大小选取模拟地层水驱替速度,将岩心置于岩心夹 持器中,加围压至2MPa,用模拟地层水驱替岩心;
[0010]步骤四、驱替稳定后,记录岩心入口压力(一般稳定等待时间需2小时左右),根据 岩心入口压力和驱替速度计算水测渗透率,并测该围压下的核磁共振T2谱;
[0011] 步骤五、保持岩样入口驱替速度不变,逐渐升高围压,重复步骤四,测不同围压下 的核磁共振Τ2谱;
[0012] 步骤六、根据下列公式将不同围压下的核磁共振Τ2谱换算成孔喉半径; ΙχΤ,χρ
[0013] r.=~-- !000
[0014] 式中:T2为核磁共振测得的横向弛豫时间,ms;r为孔喉半径,μπι;Ρ为表面弛豫率,μ m/s;
[0015] 步骤七、将不同围压下的孔喉半径分布曲线绘制在同一张图上,对比可定量得出 不同围压下的孔喉变化程度。
[0016] 进一步地,所述地层水矿化度达到所取岩心所在油藏的地层水矿化度。
[0017] 进一步地,所述根据岩心气测渗透率大小选取驱替速度满足条件为小于发生速敏 的临界流速;所述速敏为流体在岩心中流动时,因流体流动速度变化引起岩心中微粒运移、 堵塞喉道,导致岩心渗透率下降的现象。
[0018] 进一步地,所述根据岩心入口压力和水驱速度计算水测渗透率通过下式实现:
[0019] [,Ly"Xl()0 AAP
[0020] 式中:K为水测渗透率,10-3μπι2; Q为驱替速度,ml/s; L为岩心长度,cm; μ为模拟地层 水粘度,mPa· s;Α为岩心横截面积,cm2; Δ Ρ为岩心入口压力,MPa。
[0021 ]与现有技术相比,本发明具有以下优点:
[0022] (1)能够保证整个实验过程都在同一样品上完成,克服由于样品不同而带来的分 析结果缺陷;
[0023] (2)实验过程可始终保持地层条件,实时在线测试,无所多次卸装样品,避免卸装 过程中围压变化对孔喉变化的影响;
[0024] (3)核磁共振!^谱与孔喉半径存在一一对应关系,可以实现定量评价,克服定性分 析的不足。
[0025] 本方法主要利用核磁共振无损测试,定量、实时、在线评价应力敏感实验过程中, 围压变化对同一样品孔喉变化的影响。能够明确应力敏感过程中孔喉的变化程度,对于减 小应力敏感程度,采取相应措施,进而指导油田合理开发具有重要的参考价值。
附图说明
[0026]图1为本发明方法流程框图。
[0027]图2为核磁共振T2谱图。
[0028]图3为孔喉半径分布图。
具体实施方式
[0029] 下面选取某油田样品结合附图对本发明做详细叙述。
[0030] 如图1所示,本发明定量评价应力敏感过程中孔喉变化程度的方法,包括下述步 骤:
[0031] 步骤一、从全直径岩心上钻取直径2.5cm、长度介于5cm的岩心,洗油烘干,测孔隙 度为10.70%,气测渗透率为0.33 X 10-3μπι2;
[0032] 步骤二、配制实验用模拟地层水,达到地层水矿化度为65670mg/L;
[0033] 步骤三、根据岩心气测渗透率大小选取模拟地层水驱替速度为O.llml/min,根据 岩心气测渗透率大小选取驱替速度满足条件应小于发生速敏的临界流速;速敏是指流体在 岩心中流动时,因流体流动速度变化引起岩心中微粒运移、堵塞喉道,导致岩心渗透率下降 的现象。然后将岩心置于岩心夹持器中,加围压至2MPa,用模拟地层水驱替岩心;
[0034] 步骤四、驱替稳定后,记录岩心入口压力为2.31MPa,计算水测渗透率为0.081 X 10-3μπι2,测核磁共振T2谱;
[0035] 计算水测渗透率通过下式实现:
[0036] …遞 ΑΛΡ
[0037] 式中:Κ为水测渗透率,10-3μπι2; Q为驱替速度,ml/s; L为岩心长度,cm; μ为模拟地层 水粘度,mPa· s;Α为岩心横截面积,cm2; Δ Ρ为岩心入口压力,MPa。
[0038] 步骤五、保持岩样入口驱替速度为O.llml/min不变,逐渐升高围压至6MPa、10MPa、 15MPa、20MPa,驱替稳定后岩心入口压力分别为2 · 8 IMPa、3 · 06MPa、3 · 17MPa、3 · 25MPa,计算 水测渗透率分别为Ο · 066 X 10-3μπι2、0 · 061 X 10-3μπι2、0 · 058 X 10-3μπι2、0 · 057 X 10-3ym2、测核 磁共振Τ2谱见图2;
[0039] 步骤六、根据下列公式将不同围压下的核磁共振Τ2谱换算成孔喉半径; Γ Ί 2x7:, χ 2.39 0040] r 二-- 1000
[0041]步骤七、将不同围压下的孔喉半径分布曲线绘制在同一张图上,参照图3,根据图3 对比可知当围压从2MPa增加至6MPa,分布于0.14μπι~5.22μπι范围内的孔喉含量减小;当围 压增加至lOMPa,分布于0·30μπι~2·61μπι范围内的孔喉含量减小;当围压增加至15MPa,主要 是分布于ο. 49μπι~1.13μπι和0.05μπι~0.09μπι范围内的孔喉含量减小;当围压增加至20MPa, 主要是分布于〇. 〇2μπι~0.25μπι范围内的孔喉含量减小。
[0042]实验方法的原理说明
[0043] 岩心饱和流体后置于均匀分布的静磁场中,流体中的氢核(¾)会被磁场极化,产 生磁化矢量。核磁共振T2谱主要反映了孔隙介质中流体的流动特性,而这一流动特性受孔 喉制约,核磁共振T2谱与孔喉半径存在对应关系,T2谱的变化反映了孔喉的分布情况。实验 过程中模拟地层温度和压力,流体采用模拟地层水,岩心夹持器用不含磁PEEK材料制成,用 不含氢信号的氟油加围压来模拟地层压力,实验室在恒温箱中进行,可在线实时测试T2谱。 不同有效应力条件下,核磁共振T2谱分布不同,将核磁共振T2谱转换成孔喉半径,即可反映 孔喉的变化,可定量评价应力敏感过程中孔喉的变化程度。
[0044] 以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定 本发明的具体实施方式仅限于此,对于本发明所属技术领域的普通技术人员来说,在不脱 离本发明构思的前提下,还可以做出若干简单的推演或替换,都应当视为属于本发明由所 提交的权利要求书确定专利保护范围。

Claims (4)

1. 一种定量评价应力敏感过程中孔喉变化程度的方法,其特征在于,包括下述步骤: 步骤一、从全直径岩屯、上钻取直径2.5cm和长度介于4-6cm的岩屯、,洗油烘干,测岩屯、孔 隙度和气测渗透率; 步骤二、配制实验用模拟地层水,达到地层水矿化度; 步骤Ξ、根据岩屯、气测渗透率大小选取模拟地层水驱替速度,将岩屯、置于岩屯、夹持器 中,加围压至2MPa,用模拟地层水驱替岩屯、; 所述根据岩屯、气测渗透率大小选取模拟地层水驱替速度满足条件为小于发生速敏的 临界流速;所述速敏是为流体在岩屯、中流动时,因流体流动速度变化引起岩屯、中微粒运移 和堵塞喉道,导致岩屯、渗透率下降的现象; 步骤四、驱替稳定后,记录岩屯、入口压力,根据岩屯、入口压力和驱替速度计算水测渗透 率,并测该围压下的核磁共振T2谱; 步骤五、保持岩屯、入口驱替速度不变,逐渐升高围压,重复步骤四,测不同围压下的核 磁共振T2谱; 步骤六、根据下列公式将不同围压下的核磁共振T2谱换算成孔喉半径。
Figure CN103926267BC00021
式中:T2为核磁共振测得的横向弛豫时间,ms; r为孔喉半径,μπι; Ρ为表面弛豫率,ym/s; 步骤屯、绘制不同围压下的孔喉半径分布曲线图,定量得出不同围压下的孔喉的尺寸 大小及含量变化程度。
2. 根据权利要求1所述的定量评价应力敏感过程中孔喉变化程度的方法,其特征在于, 所述地层水矿化度达到所取岩屯、所在油藏的地层水矿化度。
3. 根据权利要求1所述的定量评价应力敏感过程中孔喉变化程度的方法,其特征在于, 所述模拟地层水驱替岩屯、稳定后等待时间需2-3小时。
4. 根据权利要求1所述的定量评价应力敏感过程中孔喉变化程度的方法,其特征在于, 所述根据岩屯、入口压力和驱替速度计算水测渗透率通过下式实现:
Figure CN103926267BC00022
式中:K为水测渗透率,10-3μπι2; Q为驱替速度,ml/s; L为岩屯、长度,cm; μ为模拟地层水粘度,mPa. S; A为岩屯、横截面积,cm2; Δ P为岩屯、入口压力,MPa。
CN201410174667.3A 2014-04-28 2014-04-28 一种定量评价应力敏感过程中孔喉变化程度的方法 CN103926267B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410174667.3A CN103926267B (zh) 2014-04-28 2014-04-28 一种定量评价应力敏感过程中孔喉变化程度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410174667.3A CN103926267B (zh) 2014-04-28 2014-04-28 一种定量评价应力敏感过程中孔喉变化程度的方法

Publications (2)

Publication Number Publication Date
CN103926267A CN103926267A (zh) 2014-07-16
CN103926267B true CN103926267B (zh) 2016-05-04

Family

ID=51144564

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410174667.3A CN103926267B (zh) 2014-04-28 2014-04-28 一种定量评价应力敏感过程中孔喉变化程度的方法

Country Status (1)

Country Link
CN (1) CN103926267B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104405374B (zh) * 2014-10-16 2017-12-05 中国石油天然气股份有限公司 一种致密气藏储层应力敏感性的测量方法
CN106324008B (zh) * 2015-07-02 2018-10-16 中国石油天然气股份有限公司 核磁共振测井波谱的校正方法
CN105158287B (zh) * 2015-07-08 2017-05-10 上海大学 一种凝胶在岩石微通道中的运移性的可视化评价方法
CN105067794A (zh) * 2015-07-24 2015-11-18 成都理工大学 测试页岩储层水敏感性、盐敏感性及碱敏感性的方法
CN106522931A (zh) * 2015-09-14 2017-03-22 中国石油化工股份有限公司 模拟地层条件下的注入水伤害评价方法
CN105606513B (zh) * 2016-01-29 2018-07-24 河海大学 确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法
CN106153662A (zh) * 2016-06-17 2016-11-23 北京大学 岩心应力敏感性的测量方法
CN106526079B (zh) * 2016-10-27 2019-04-16 中国石油大学(北京) 一种研究致密砂岩孔喉结构动态变化的方法
CN106770415B (zh) * 2016-12-13 2019-03-15 中国石油天然气股份有限公司 一种高粘度低流度致密油在不同温压条件下的微观可流动性的评价方法
CN106884634B (zh) * 2016-12-30 2019-09-17 浙江海洋大学 一种研究高矿化度地层水对co2驱影响的实验方法
CN106950242B (zh) * 2017-03-24 2018-11-20 西安石油大学 一种定量评价矿化度对页岩油藏自吸水驱油能力影响的方法
CN107607451A (zh) * 2017-08-08 2018-01-19 广州海洋地质调查局 一种天然气水合物岩心流速敏感性测定方法
CN107688037B (zh) * 2017-08-17 2020-04-07 中国海洋石油集团有限公司 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法
CN109612897A (zh) * 2018-09-25 2019-04-12 西安石油大学 污水回注对致密砂岩物性及孔喉分布影响的定量评价方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4542343A (en) * 1983-06-17 1985-09-17 Chevron Research Company Computer-controlled, portable pulsed NMR instrument and method of use wherein depolarization is ultrashortened using pseudo randomly spaced and oriented RF depolarization pulses
CN201902206U (zh) * 2010-11-24 2011-07-20 中国石油大学(北京) 应力敏感性地层钻井堵漏模拟评价装置
CN102353625A (zh) * 2011-08-31 2012-02-15 中国石油大学(北京) 渗流力学实验中水测覆压孔隙度的测定方法
CN103257151A (zh) * 2013-04-24 2013-08-21 西安石油大学 一种定量评价油气二次运移过程中孔喉动用规律的方法
WO2013148508A1 (en) * 2012-03-29 2013-10-03 Baker Hughes Incorporated Carbonate permeability by pore typing
CN103528934A (zh) * 2013-10-29 2014-01-22 中国石油大学(北京) 一种测量超低渗岩石渗透率应力敏感性的互相关技术

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4542343A (en) * 1983-06-17 1985-09-17 Chevron Research Company Computer-controlled, portable pulsed NMR instrument and method of use wherein depolarization is ultrashortened using pseudo randomly spaced and oriented RF depolarization pulses
CN201902206U (zh) * 2010-11-24 2011-07-20 中国石油大学(北京) 应力敏感性地层钻井堵漏模拟评价装置
CN102353625A (zh) * 2011-08-31 2012-02-15 中国石油大学(北京) 渗流力学实验中水测覆压孔隙度的测定方法
WO2013148508A1 (en) * 2012-03-29 2013-10-03 Baker Hughes Incorporated Carbonate permeability by pore typing
CN103257151A (zh) * 2013-04-24 2013-08-21 西安石油大学 一种定量评价油气二次运移过程中孔喉动用规律的方法
CN103528934A (zh) * 2013-10-29 2014-01-22 中国石油大学(北京) 一种测量超低渗岩石渗透率应力敏感性的互相关技术

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
低渗储层物性压力敏感性研究;刘建军 程林松;《新疆石油科技》;20050630;第15卷(第2期);16-18 *
低渗透油藏两相渗流参数及驱油效率的应力敏感性研究;刘丽;《科学技术与工程》;20120831;第12卷(第23期);5869-5872 *

Also Published As

Publication number Publication date
CN103926267A (zh) 2014-07-16

Similar Documents

Publication Publication Date Title
Chong et al. Experimental investigations on energy recovery from water-saturated hydrate bearing sediments via depressurization approach
Zeinijahromi et al. Effect of fines migration on oil–water relative permeability during two-phase flow in porous media
Pini et al. Capillary pressure and heterogeneity for the CO2/water system in sandstone rocks at reservoir conditions
Wang et al. Effect of low-velocity non-Darcy flow on well production performance in shale and tight oil reservoirs
Ghanizadeh et al. A comparison of shale permeability coefficients derived using multiple non-steady-state measurement techniques: Examples from the Duvernay Formation, Alberta (Canada)
CN104990851B (zh) 一种新的页岩敏感性实验研究方法
Liu et al. Experimental evaluation of CO2 enhanced recovery of adsorbed-gas from shale
Delli et al. Experimental determination of permeability of porous media in the presence of gas hydrates
Gao et al. Pore structure characterization, permeability evaluation and enhanced gas recovery techniques of tight gas sandstones
Vavra et al. Geological applications of capillary pressure: a review
Tanikawa et al. Klinkenberg effect for gas permeability and its comparison to water permeability for porous sedimentary rocks
Tan et al. Experimental study of permeability and its anisotropy for shale fracture supported with proppant
Al‐Wardy et al. Effective stress law for the permeability of clay‐rich sandstones
Huang et al. Characterization of pore structure, gas adsorption, and spontaneous imbibition in shale gas reservoirs
Meng et al. Monitor the process of shale spontaneous imbibition in co-current and counter-current displacing gas by using low field nuclear magnetic resonance method
Yang et al. Effect of depressurization pressure on methane recovery from hydrate–gas–water bearing sediments
Li et al. Characterization of spontaneous water imbibition into gas-saturated rocks
McLatchie et al. The effective compressibility of reservoir rock and its effects on permeability
Zemke et al. Petrophysical analysis to investigate the effects of carbon dioxide storage in a subsurface saline aquifer at Ketzin, Germany (CO2SINK)
Wang et al. Characteristics of oil distributions in forced and spontaneous imbibition of tight oil reservoir
Jin et al. Growth of methane clathrate hydrates in porous media
Roychaudhuri et al. An experimental and numerical investigation of spontaneous imbibition in gas shales
Mason et al. Correlation for the effect of fluid viscosities on counter-current spontaneous imbibition
US2618151A (en) Cell for measuring relative permeability
Cao et al. A multiscale-multiphase simulation model for the evaluation of shale gas recovery coupled the effect of water flowback

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model