1.一种裂缝各向异性油藏注水开发可预测物理模型建立方法，其特征在于，该方法包括：(a)根据裂缝性油藏水驱油开发过程的特点，利用渗流力学理论和相似性分析，建立裂缝性油藏开发模拟的相似性准则，所述相似性准则包括外形与空间相似、井筒几何相似、岩石物性相似、油水粘度相似、重力-压力相似、基质与裂缝可动油量比相似、基质内含油量分布的相似、基质渗吸与裂缝驱替特征时间相似、渗吸强度分布相似及时间过程相似；(b)相似性准则的实现及模型参数设计方法；(c)建立满足多重相似性的油藏宏观物理模型，以全面模拟预测实际裂缝性油藏的渗流特征和开发过程。
2.如权利要求1所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特
征在于，所述步骤(a)中所述岩石物性分布相似包括裂缝渗透率分布相似及裂缝孔隙度相似，且所述相似性准则还包括饱和度分布相似及位势分布相似，该相似性准则共包括下面的裂缝各向异性油藏水驱油模拟相似性准则表中所列的22个相似准数。
裂缝各向异性油藏水驱油模拟相似性准则表
3.如权利要求2所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特征在于，所述步骤(a)中包括(a1)确定相似性准则的建模条件，所述建模条件包括：
(a11)油藏介质为双孔单渗，即基质和裂缝均为流体存储空间，裂缝系统是渗流通道；
(a12)考虑基质-裂缝间的渗吸作用；
(a13)考虑重力及油水重度差的影响；
(a14)考虑裂缝渗透率的各向异性；
(a15)忽略裂缝中的毛管力；
(a16)忽略流体及岩石的压缩性。
4.如权利要求3所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特征在于，所述步骤(a)中，还包括(a2)建立无量纲化渗流数学模型：
裂缝中的油水运动方程：
裂缝中的物质平衡方程：
自然限制条件：
So+Sw＝1，qw+qo＝0 (3)
动态渗吸方程：
初始条件：
Φo(x，y，z，t＝0)＝Φi，Φw(x，y，z，t＝0)＝0，Sw(x，y，z，t＝0)＝0 (5)边界条件：假设油藏边界Γ为封闭边界，n为边界法向，则
对于井筒边界，假设为定压注采，则
上述各式中，v、q、S分别表示渗流速度、渗吸强度、饱和度，下标o、w分别表示油和水。
x、y、z为直角坐标系的三个坐标， 表示油藏空间点， 和 分别表示注水井和生
产井井筒上任意一点。Aw、Ao、A分别为张量形式的水相流度、油相流度及流体总流度：
K为各向异性渗透率张量，φ为孔隙度，Kro和Krw分别是油相和水相的相对渗透率，Φo，Φw分别为油相和水相的位势，Φi为初始位势，p为压力，Δp为注采压差。记γ为重度，G为油水重力势差，则有如下表达式：
Φo＝p+γoz，Φw＝p+γwz，G＝(γw-γo)z (8)
*
R为单位体积基岩所含可动油体积，λ表示渗吸强度。记T 为渗吸半周期，则
*
λln2/T (9)
把运动方程带入物质平衡方程，得：
(10)和(11)两式相加，得：
(10)可化为：
故渗流数学模型可写为：
将(14)式在以渗透率主方向为坐标方向的直角坐标系中展开，并设渗透率主值分别
为Kx、Ky、Kz，可得：
裂缝中的油水相渗曲线为相互交叉的对角线，其表达式如下：
Kw＝K·Sw (16)
Ko＝K·(1-Sw) (17)
故(15)式可化为
下面将数学模型(18)式无量纲化：
5个自变量无量纲化：
4个孔渗参数无量纲化：
3个流体参数无量纲化：
2个渗吸常数无量纲化：
3个 因 变 量 无 量 纲 化 ： SwD＝Sw (23)其中，Lx、Ly、Lz——x、y、z方向的特征长度，可取最大或平均长度；
——水驱特征时间；
——平均孔隙度；
——x、y、z方向渗透率主值的平均值；
——油相、水相流体平均粘度；因忽略流体压缩性，
——平均油水重度差；因忽略流体压缩性，
——R，λ的平均值，
把(19)～(23)带入(18)式，得
考虑到(21)式，上式变为
5.如权利要求1所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特征在于，所述步骤(b)中相似性准则的实现方法包括：
(b11)各向异性渗透率相似性(π7～π11)的实现方法：
相似准数π7～π11要求裂缝性油藏物理模型与实际油藏保持渗透率分布的非均质
性及各向异性相似，即在任意区域满足：
(32)
按照(32)式要求，物理模型内部渗透率的相对分布必须与实际油藏分布成正比，但不限定其绝对大小；
实现方法：根据模型渗透率主值的平均值 利用(32)中各式计算模型内所
有区域、所有方向的渗透率主值，依此设计制作模型便能满足相似准数π7～π11的要求。
6.如权利要求5所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特征在于，所述步骤(b)中相似性准则的实现方法还包括(b12)裂缝孔隙度相似性(π12)的实现方法，具体为：
相似准数π12要求模型内裂缝孔隙度分布与实际油藏保持相似，即在任意区域满足：
根据裂缝渗流理论，当模型内裂缝分布满足渗透率的相似性时，其孔隙度相似性会同时得到满足。
7.如权利要求6所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特征在于，所述步骤(b)中相似性准则的实现方法还包括(b13)油水粘度比相似性(π13)的实现方法：
因为不计流体的压缩性，所以可认为驱替液和被驱替液的粘度均为常数；根据相似准则π13，选择粘度合适的流体作为驱替液和被驱替液，使得油藏和模型的油水粘度比满足下述关系：
(μo/μw)|模型＝(μo/μw)|油藏 (34)
考虑到安全性，具有较强挥发性和毒性的轻质烃组分不适合在实验室内使用，因此选择柴油或更重的原油成分进行调合后作为被驱替液，以模拟油藏内的原油；选择水和适当的增粘剂进行调合后作为驱替液，以模拟油藏内的水
如果实际油藏的原油粘度大于或等于柴油的粘度，则根据(34)式，物理模型可用纯水作为驱替液，可用柴油和重质油调合后做被驱替液，便可使之满足油水粘度比的相似性要求；
如果实际油藏的原油粘度小于柴油的粘度，则根据(34)式，物理模型用柴油做被驱替液，其常温下粘度最小约为3.0mP·s，用纯水和适当的增粘剂调合后作为驱替液，便可使之满足油水粘度比的相似性要求；反复试验证明，用蔗糖作为增粘剂可以满足一般油藏模拟对驱替液粘度的要求，且糖水使用安全，对物理模型伤害小，因此最终确定用糖水做驱替液。
8.如权利要求7所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特征在于，所述步骤(b)中相似性准则的实现方法还包括(b14)利用有限真空饱和技术和天然砂岩选择方法实现模型和油藏的裂缝-基质可动原油储量的相似性(π15和π16)。
根据相似性准则π15和π16，实际油藏和实验模型的基质-裂缝可动油量比以及基质含可动油量分布应满足(38)式：
(38)式中，实际油藏单位体积基质中可动油量R油藏和裂缝孔隙度φ油藏可通过现场测试和文献资料调研得到，裂缝孔隙度φ油藏由模型测试得到；由(38)式可确定实验模型内任意区域单位体积基质中可动油量R模型的值；
实现基质岩块R模型合理取值的有限真空饱和技术和天然砂岩选择方法：
(1)往干燥、空气饱和的小岩块中饱和驱替液，由于单纯依靠毛管力驱替液很难进入基质岩块，因此先用真空机将岩块中的大部分空气抽取出来，使岩块内的孔隙处于“有限真空”状态，然后从小岩块周围充入驱替液，使小岩块内压力恢复原态；进入小岩块内的驱替液与残留的空气形成液包气的气液分布形态，即残留的气体位于岩块中心圆球形区域内，圆球以外的区域由驱替液占据；
(2)往上述小岩块中饱和被驱替液，从外部对小岩块抽真空，使之再次处于“有限真空”状态；利用小岩块中空气的膨胀性，向外驱出部分驱替液，然后从周围表面同时往模型中补充被驱替液，在基质小岩块中由内而外形成空气、驱替液和被驱替液三相流体依次分布的形态，即空气位于岩块中心的圆球形区域内，驱替液位于圆球以外的环形区域，环形以外的区域由被驱替液占据；
(3)把前述含有三相流体的小岩块浸入驱替液中，观测其渗吸过程及被驱替液的最终渗吸量Q，设小岩块的体积为V，则R模型＝Q/V；
(4)R模型值主要决定于饱和过程真空度和小岩块物理性质，利用上述步骤(1)～(3)进行多种试验，就可以确定真空度、岩石类型和R模型之间的关系；然后采用相同的真空度，对多种岩石进行试验选择，找出所有符合R模型取值要求的岩石类型。
9.如权利要求8所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特
征在于，所述步骤(b)中相似性准则的实现方法还包括(b15)裂缝-基质供液能力相似性(π17和π18)的实现方法：
根据相似准则π17和π18要求实际油藏和实验模型的渗吸半周期-裂缝驱替特征时
间比应满足(39)式：
(39)式中，实际油藏的渗吸半周期 可通过现场测试或文献资料调研得到，实际油
藏水驱特征时间T油藏和实验模型的水驱特征时间T模型可由水驱特征时间的定义式求得：
由(39)式确定实验模型内任意区域渗吸半周期 的取值； 值的实现及岩石类型
的确定方法：对符合R横型取值要求的所有类型的岩石材料，进行渗吸半周期测试实验，找到渗吸半周期 符合(39)式要求的天然砂岩；用这些砂岩制作物理模型，即可实现模型与油藏之间的裂缝-基质原油储量和供油能力的相似性(π15～π18)。
10.如权利要求9所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特
征在于，所述步骤(b)中模型参数设计方法包括：
(b21)根据实际油藏的尺度和形状、实验室空间条件及关系式(25)，计算确定模型的几何尺度(Lx，Ly，Lz)和形状，以及模型中小岩块的大小、数量；
(b22)根据实际油藏及其井筒的几何参数确定模型内的井筒半径rw模型：首先利用(25)式得到
rw1＝Lr模型·rw油藏/Lr油藏
如果rw1≥6.0mm，则取rw模型＝rw1，此时(30)式中rw1＝rw2；如果rw1＜6.0mm，则取rw模型＝rw2＝6.0mm，此时(30)式中rw1≠rw2；
(b23)根据小岩块加工和粘接工艺过程，确定 和 的值；
(b24)根据实际油藏的裂缝渗透率与孔隙度分布及相似关系(32)、(33)式，计算确定模型中渗透率分布、孔隙度分布和裂缝分布，确定模型中每个小岩块的粘结方式；
(b25)根据实际油藏的油水粘度和(34)式，利用(b13)所述方法，试验设计具有合适粘度的驱替液和被驱替液；
(b26)根据实际油藏的注采井底压力和油水密度，以及物理模型驱替液和被驱替液的密度，并考虑扩大井径和表皮系数的影响，利用重力压差与注采压差之比的相似性计算确定模型的注采压力；
(b27)根据实际油藏参数及物理模型的裂缝孔隙度和水驱特征时间，利用(38)、(39)式计算确定模型内每个区域的单位体积基质可动油储量R模型和基质-裂缝渗吸半周期
11.如权利要求10所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特征在于，所述步骤(c)中包括：(c1)天然砂岩的选择及小岩块的制备，具体为：
(c11)针对物理模型的每个区域，利用所述裂缝-基质原油储量和供液能力的相似性
(π15～π18)实现方法，在相同的操作条件，其两次抽真空的压力分别为 和 对小岩块进行饱和及渗吸半周期测试，选择合适的小岩块砂岩品种，使之同时满足单位体积基质中可动油量R模型和渗吸半周期 两方面的要求；
(c12)使用步骤(c11)中选定的天然砂岩加工制作正方形小岩块；小岩块的边长一般
可取25mm～50mm，同一物理模型中所有小岩块的尺寸必须严格相等；
(c13)在设计井筒穿过的小岩块上钻孔，形成预设的井眼；井径取步骤(b22)中的设计值；
(c14)对需要预设各种测试管线及流体饱和通道的小岩块进行加工处理。
12.如权利要求11所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特征在于，所述步骤(c)中还包括(c2)物理模型制作过程，具体为：
(c21)根据步骤(b24)的设计方案，将小岩块顺序粘结形成大尺度的物理模型岩体；
(c22)在模型岩体表面均匀涂刷环氧树脂胶，待其凝固形成封闭的模型边界；
(c23)连接各井筒及测试点的管线，在模型底部和顶部设置流体饱和通道。
13.如权利要求12所述的裂缝各向异性油藏注水开发可预测物理模型建立方法，其特征在于，所述步骤(c)中还包括(c3)流体饱和过程，具体为：
(c31)模型饱和过程采用有限真空技术，可采用与步骤(c11)相同的操作时间和同样
的真空度；
(c32)饱和过程开始，首先利用真空泵从顶部饱和通道将模型内压力降至 然后保
持压力 不变，向模型底部饱和通道注入驱替液，直到裂缝系统全部充满驱替液，形成底部注、顶部采的循环流动，使得模型内每个小岩块所处的流体环境相同，然后关闭顶部通道，模型内压力逐渐上升至初始状态，驱替液进入所有小岩块；
(c33)利用真空泵从模型底部饱和通道抽取驱替液，将模型内压力降至 然后保持
压力 不变，向模型顶部饱和通道注入被驱替液，直到裂缝系统全部充满被驱替液，形成顶部注入、底部采出的循环流动，使得模型内每个小岩块所处的流体环境相同，然后关闭顶部通道，模型内压力逐渐上升至初始状态，被驱替液进入所有小岩块，此时模型内每个小岩块内的流体分布达到模型饱和要求。