Laboratory Creep Tests For Time-Dependent Properties of A Marble in Jinping II Hydropower Station
Laboratory Creep Tests For Time-Dependent Properties of A Marble in Jinping II Hydropower Station
        Abstract: In order to investigate the time-dependent behaviors of deep hard rocks in the diversion tunnel of Jinping II
        hydropower station, uniaxial creep tests were carried out by using the triaxial testing machine RC-2000. The axial compressive
        load was applied step by step and each creep stage was kept for over several days. Test results show that: (1) The lateral
        deformation of rock specimens is 2–3 times the axial compressive deformation and accelerates drastically before damage, which
        may be employed as an indicator to predict the excavation-induced instability of rocks. (2) The resultant deformation changes
        from compression to expansion when the Poisson’s ratio is larger than 0.5, indicating the starting point of damage. (3) In the
        step-loading stages, the Poisson’s ratio approximately remains constant; under constantly imposed load, the Poisson’s ratio
        changes with elapsed time, growing continuously before the specimen is damaged. (4) When the applied load reaches a certain
        threshold value, the rock deteriorates with time, and the strength of rocks approximately has a negative exponent relation with
        time. (5) The failure modes of the deep marble are different in long- and short-term loading conditions. Under the condition of
        short-term loading, the specimen presents a mode of tensile failure; while under the condition of long-term loading, the specimen
        presents a mode of shear failure, followed by tensile failure.
        Key words: time-dependent mechanical behaviors; marble; long-term strength; the Poisson’s ratio of rocks; rock creep
mudstone, and considered that the lateral acceleration                2.2 Rock samples
creep stage occurred earlier than the axial                              The rock samples taken from the auxiliary tunnel
acceleration creep stage and the Poisson’s ratio had a                of Jinping II hydropower station are gray-white,
nonlinear relationship with creep. Cui and Fu (2006)                  medium-fine crystalline marbles of Baishan group
conducted uniaxial compression rheological tests on                   (T2b), the sampling location is shown in Fig. 2. The
red freestone, and observed that in the failure stage,                rocks are hard, brittle and compact. Elasticity
the lateral creep was evident than the axial creep, and               modulus of the marble varies from 25 to 40 GPa, and
the lateral acceleration creep stage appeared earlier                 the deformation modulus is 8–16 GPa. No palpable
than the axial stage. Fujii et al. (1999) believed that               infilling materials, filling belt or structural plane in
                                                                      marbles can be observed with naked eyes.
the lateral deformation was larger than the axial
                                                                      Homogeneous lithology shows that the components
deformation. In this regard, the lateral deformation
                                                                      of minerals are mainly composed of quartz, potash
can be used as an index of rock damage in creep tests
                                                                      feldspar, calcite, mica, chlorite and smectite. The
of rocks.
                                                                      cylinder rock samples, with a diameter of 50 mm and
   Most rocks in the diversion tunnel of Jinping II
                                                                      height of 100 mm, are obtained by drill-and-blast
hydropower station are composed of marble (Chen et
                                                                      method. In addition, the two end faces of rock
al., 2010; Li et al., 2011; Wu et al., 2011). The in-situ
                                                                      samples are polished smoothly. Processing accuracy
stress of 60-70 MPa after excavation is observed in
                                                                      of the test is controlled properly in accordance with
field. Rock falling caused by time-dependency and
                                                                      specification (2007).
time-delayed rockbursts are the two main
engineering disasters in Jinping II hydropower
station. In this paper, by using deep marble samples
taken 2 500 m below the surface, axial compression
                                                                                    Sampling location
creep test is conducted to analyze the time-dependent
property, deformation and damage of marble after
tunnel excavation.
5.1 MPa should be considered till rocks fail in the                                               When a certain stress level is reached, stable creep
form of rheological acceleration. A loading rate of                                            rate increases remarkably. When the deformation
0.1 MPa/s is adopted for all loading stages. The axial                                         increases rapidly, rock samples present creep failure.
stress keeps constant after load reaches the preset                                            The three stages of attenuation, stabilization and
threshold value, and then the next loading stage can                                           acceleration can be evidently captured on rock
be considered when the lateral deformation of rocks                                            samples. The results of axial and lateral creep strains
becomes stable under previous loading stage.                                                   of rocks at different stress levels are shown in Table 1.
Recording is stopped when rock sample is damaged.                                              In Table 1, the instantaneous elastic strain is not
                                                                                               included in either the axial strain or the lateral strain.
3 Time-dependent mechanical pro-
perties of marble                                                                              Table 1 Axial and lateral creep strains of rocks at different
                                                                                               stress levels.
   Time-dependent properties of marble are strongly                                                1         1          2          v
                                                                                                                                                 1 / 1    2 /  2    2 / 1
affected by various stress levels on the axial or lateral                                        (MPa)          6
                                                                                                            (10 )           6
                                                                                                                        (10 )       (106)
deformation. Under low stress levels, the two stages
                                                                                                  84.7        33         101         169        0.031        0.04         3.06
of attenuation and stabilization are reflected to some
                                                                                                 105.1        55         126         186        0.052        0.05         2.29
extent on axial and lateral creeps of rock samples.
                                                                                                 115.3        66         162         258        0.062        0.063        2.45
The duration of attenuated stage is very short and the
                                                                                                 120.4        53         143         206        0.049        0.056         2.7
stable stage of rock creep comes quickly. With the
                                                                                                 125.5        88         190         292        0.083        0.074        2.16
increase in stress level, instantaneous elastic strain
comes up observably and rock creep decreases                                                     130.6       175         405         635        0.165        0.159        2.31
rapidly with time, which is called the attenuated 135.7 588 1 425 3 309 0.556 0.558 2.42
stage of rock creep. In addition, there is an increasing Note:  1 is the total axial creep strain,  2 is the total lateral creep
trend for the duration of attenuated stage with                                                strain, 1 is the axial creep strain, 2 is the lateral creep strain, and v is the
                                                                                               volumetric strain.
increase in stress level, which can be observed on the
representative results shown in Figs. 3 and 4. As
                                                                                               3.1 Variation in axial and lateral strains of time-
shown in Fig. 3, compressive strain is positive and
                                                                                               dependent properties of rocks
expansive strain is negative. Duration of attenuated
                                                                                                  At different stress levels, the axial and lateral
creep is defined from the end of instantaneous elastic
                                                                                               strains of rocks exhibit different time-dependent
strain to the start of constant creep rate.
                                                                                               properties. The increment in lateral creep is larger
                                                                                               than that of axial creep under the same stress level,
                                           Axial creep strain
                     0.4                                                                       i.e. the lateral creep strain is 2–3 times the axial creep
                                            Lateral creep strain             135.7 MPa
                     0.3
                                         115.3 MPa 120.4 MPa125.5 MPa
                                                                      130.6 MPa                strain. From this point of view, with the increase in
                                 105.1 MPa
                     0.2    84.7 MPa                                                           stress level, the increment in lateral creep becomes
      Strain (%)
                     0.1                                                                       larger than that in axial creep at the same stress level.
                     0.0
                           084.7 MPa50     100      150      200      250        300     350
                                                                                               The reason may be interpreted that with the increase
                   0.1         105.1 MPa 115.3 MPa        Time (hour)
                                                   120.4 MPa
                                                                                               in stress level, rocks begin to appear damage and
                                                           125.5 MPa 130.6 MPa
                   0.2                                                                        dilatation, which will result in the gradually
                                                                            135.7 MPa
                   0.3                                                                        dominant role of the lateral strain. And it can be
Fig. 3 The axial and lateral creep strain curves of marble                                     attributed to the increase in the volumetric strain. In
under uniaxial compression test (UCT).
                                                                                               the failure stage, the lateral steady creep rate is
                                                                                               evidently larger than the axial one, the lateral strain
                   35
                   30                                                                          increases rapidly, and its increment is much larger
Time (hour)
                   25
                   20
                                                                                               than that of axial strain.
                   15                                                                             From Table 1, it can also be observed that with the
                   10
                    5                                                                          increase in stress level, both the axial and lateral
                    0                                                                          creep strains exhibit an increasing trend. For
                     80          90        100        110      120            130        140
                                                  Stress (MPa)                                 instances, at the stress level of 84.7 MPa and after
Fig. 4 Curve between the duration of attenuated stage under                                    38.45 hours from the instantaneous strain, the
uniaxial compression creep and axial stress levels.                                            increment in axial creep is 33×106, and that in
Xiaojun Zhao et al. / J Rock Mech Geotech Eng. 2012, 4 (2): 168–176                                                                                                      171
                                                                                                                       Poisson’s ratio
                       0.25
  Lateral strain (%)
                                                                                               Poisson’s ratio
                                                                                         0.6                                             0.6
                       0.20
                                            120.4 MPa       130.6 MPa                                                                    0.5                                      0.1
                       0.15                             125.5 MPa                                                                        0.4
                                  Creep                                                  0.4                                             0.3                                      0.2
                                               115.3 MPa
                       0.10              105.1 MPa                                                                                       0.2
                              84.7 MPa                                                                                                   0.1                                      0.3
                                                                                         0.2
                       0.05                                                                                                              0.0
                                                                                                                                            0   50   100   150 200 250     300   350
                       0.00                                                              0.0                                                               Time (hour)
                          0.10           0.15        0.20           0.25        0.30   0.35
                                                                                                                 Fig. 7 Evolution law of Poisson’s ratio and volumetric strain
                                                     Axial strain (%)
Fig. 6 The relation between the lateral strain, the Poisson’s                                                    with time.
ratio and the axial strain.
                                                                                                                 loading step, which can reflect that the growth rate of
                                                                                                                 lateral creep deformation is larger than that of axial
(Wang, 2007) similar to the Poisson’s ratio needs to
                                                                                                                 creep deformation. The Poisson’s ratio becomes
be introduced to describe the lateral deformation of
rock sample. Thus, the Poisson’s ratio herein does                                                               larger than 0.5 when the applied load on the rock
not have the same meaning as that in the elastic range.                                                          specimen lasts for 225 hours, and when the
It includes the range of viscoelastic-plasticity in                                                              volumetric deformation  v is zero and the stress of
consideration of the time effect.                                                                                the rock is 130 MPa. It is basically consistent with
   From Fig. 6, it is observed that the relation curve                                                           the time point when the volumetric deformation
of lateral strain and axial strain increases stepwise                                                            shifts from compression to expansion, which shows
and the curves at each step consist of two parts:                                                                the initial damage development of rock specimen.
approximately straight line segment and slope-                                                                   The micro-fissures begin to reproduce and expand
increasing curve segment, i.e. the loading and creep                                                             when damage occurs in the rock specimen. With the
stages. The Poisson’s ratio in the two stages is                                                                 elapsed time, the micro-fissures increase and
defined as loading Poisson’s ratio and creep                                                                     propagate, resulting in the Poisson’s ratio larger than
Poisson’s ratio, respectively. Evidently, the lateral                                                            0.5. Sometimes, the Poisson’s ratio could exceed 1.0
creep has approximately a linear relation with the                                                               at the failure of the rock specimen.
axial creep, which shows that the loading Poisson’s                                                                 The tendency of the variation in the Poisson’s ratio
ratio can be roughly regarded as a constant during the                                                           with stress, under uniaxial compression and uniaxial
loading process. When the axially applied load keeps                                                             compression creep, is shown in Fig. 8. It can be seen
constant and the rock specimen reaches the creep                                                                 that:
stage, the relation curve of lateral and axial creep                                                                (1) Under uniaxial compression, the Poisson’s
gradually transfers from straight line segment to                                                                ratio decreases slightly with the increase in stress and
curve segment. The slope of the curve increases                                                                  it can be regarded as a constant since the decreasing
gradually and keeps steady finally, which is related                                                             amplitude is very small. When the applied stress
to the creep loading process, i.e. the rock specimen                                                             reaches 68 MPa, the Poisson’s ratio begins to
first reaches the decayed creep stage and then                                                                   increase. When the imposed stress reaches 83 MPa,
develops to the steady creep. The slopes of the curves                                                           the Poisson’s ratio increases to 0.66 drastically and
and the Poisson’s ratio increase gradually at the stage                                                          the rock specimen is damaged. Under the condition
of decayed creep, while the slopes of the curves and                                                             of low stress level, the cracks inside the rock
the Poisson’s ratio keep constant at the stage of                                                                specimen are coalesced and the axial deformation is
steady creep. It indicates that the creep Poisson’s                                                              larger than the lateral deformation, which results in
ratio has experienced the evolution from fluctuant to                                                            the decrease tendency of the Poisson’s ratio. When
steady state, not constant.                                                                                      the stress reaches the threshold that the damage
   The evolution law of the Poisson’s ratio is shown                                                             occurs in the specimen, the crack inside the specimen
in Fig. 7. In a global sense, the variation in loading                                                           begins to develop and the lateral deformation
Poisson’s ratio is small and keeps nearly constant.                                                              becomes larger than the axial deformation, which
The loading Poisson’s ratio varies under different                                                               results in increase in the Poisson’s ratio. When the
stress levels. At the stage of creep, the creep                                                                  stress reaches the critical failure stress, the internal
Poisson’s ratio increases continuously at each                                                                   crack propagates and expands sharply. The lateral
Xiaojun Zhao et al. / J Rock Mech Geotech Eng. 2012, 4 (2): 168–176                                                                                                                    173
                                                                                                                             2.0
                            1.0                                                                                             1.8                  130.6 MPa
                                                                                                                            1.6                  125.5 MPa
                            0.8                                                                                                                  120.4 MPa
                                                                                                                            1.4
                                                                                 135.7 MPa
                                                                                                                            1.2                  115.3 MPa
                            0.6                                                130.6 MPa    Creep
                                                                           125.5 MPa                                         1.0
                                                                                           Loading                                               105.1 MPa
                            0.4                                         120.4 MPa                                           0.8
                                                                    115.3 MPa                                                                    84.7 MPa
                                                        105.1 MPa                                                           0.6
                            0.2            84.7 MPa
                                                                                                                            0.4
                                                           Creep
                                                   Loading                                                                  0.2
                           0.0
                              50           60   70 80 90 100 110 120 130 140                                                 0.0
                                                    Axial stress (MPa)                                                             0    10       20    30         40   50    60   70
                                 (a) The relationship between Poisson’s ratio and stress.                                                              Time (hour)
                                 0.8
                                                                                                     Fig. 9 The evolution laws of deformation rate and time under
                                                                                                     different loads.
                                 0.6
               Poisson’s ratio
shows that in the low stress levels, rock’s creep only                                                                      0.000
                                                                                                                                    0        5         10        15         20    25
shows two stages: attenuation creep and steady creep.                                                                                                   Time (hour)
Under various high load levels, the creep rate quickly                                               Fig. 10 The relationship between creep rate and time in failure
                                                                                                     stage of marble.
reduces to a certain value, almost constant, with time,
which means that the rock’s creep is in a steady creep
stage, but the deformation of rocks still grows. The                                                    Assuming that the transferring point from steady
duration for decayed creep of rocks continuously                                                     creep stage to accelerated creep is point A, and that
increases with the applied stress. The increasing                                                    from axial creep to steady creep stage is point B, it
applied stress of the next stress level has a                                                        can be seen from Fig. 10 that point A grows ahead of
considerable effect on the initial creep rate of rock                                                point B, which means that the lateral deformation
sample, but the effect diminishes with time and                                                      accelerates earlier than the axial deformation before
finally disappears.                                                                                  rock failure. Therefore, if we pre-arrange
174                                                                                                        Xiaojun Zhao et al. / J Rock Mech Geotech Eng. 2012, 4 (2): 168–176
                                                                                                   Strength (MPa)
                                                                                                                                                    Fitting curve
                                                                                                                    130.7
prevention can be saved. In Fig. 11, after about 30
                                                                                                                    130.6
hours, the creep rate becomes stable.                                                                               130.5
                        4.0                                                                                         130.4
                                                                                                                    130.3
                        3.5
                                       Axial creep rate                                                             130.2
                        3.0            Lateral creep rate
                                                                                                                    130.1
  Creep rate (105)
                        2.5                                                                                                 2   4     6     8       10    12    14   16   18
                                                                                                                                                Time (hour)
                        2.0                                                                                                 Fig. 13 Strength-time curves of marble.
                        1.5
                        1.0                                                                         When t   , the long-term strength of hard
                        0.5                                                                      rock can be obtained:  t =130.1 MPa.
                        0.0
                                                                                                 3.6 Failure mechanism
                              0   10      20      30     40        50       60        70            Fig. 4 indicates that the accelerating stage of
                                                 Time (hour)                                     rock’s creep occurs at stress level of failure stage,
Fig. 11 The evolution laws of deformation rate and time under                                    and basically the creep failure of rocks occurs as long
applied stress of 130.6 MPa.                                                                     as the applied stress level is greater than the rock’s
                                                                                                 long-term strength. Wang et al. (2010) proposed that
3.5 Evolution law of rock strength vs. time                                                      the accelerating stage of rock creep would occur
   With the creep testing data, the stress-strain                                                when the applied stress was greater than the yield
isochronism curves at different times can be plotted                                             stress of rocks. As the applied stress is greater than
(Fig. 12). It shows that each isochronism curve is                                               the long-term strength of rocks, damage is observed
composed of linear and nonlinear segments with a                                                 on rock interior part, and micro-fissures begin to
turning point. As a whole, the turning point decreases                                           propagate. With the elapsed time, damages further
slowly with time and approximates to a limit value,                                              accumulate on interior of rocks and micro-fissures
i.e. long-term strength of rocks. Xu (1997) argued                                               further develop, thus the rock strength begins to
that the turning point is a symbol of rocks                                                      decrease. With the development of the micro-fissures
                                                                                                 in rocks, the rock strength reduces to a certain stress
transforming from the viscoelastic stage to
                                                                                                 level and the rock is finally damaged. The
viscoplastic stage, where rock interior structure
                                                                                                 acceleration of rock creep stage is believed to be
damage occurs and rock strength begins to decrease.
                                                                                                 rapid growing and connection of micro-fissures. In
                                                     2 h4 h           12 h                       this moment, the lateral creep rate and the axial creep
                        136                                6 h8 h10 h     14 h 18 h
                        134
                                                                                                 rate will both increase rapidly, leading to the ultimate
                                                                                                 damage of rock specimen. In rock projects, the start
                        132
         Stress (MPa)
changing from compression to expansion, and at the                                 Tempo-spatial characteristics and influential factors of rockburst: a
same time, the Poisson’s ratio is larger than 0.5.                                 case study of transportation and drainage tunnels in Jinping II
                                                                                   hydropower station. Journal of Rock Mechanics and Geotechnical
   (3) When the applied stress is beyond a certain
                                                                                   Engineering, 2011, 3 (2): 179–185.
threshold value, the rock strength decreases with time.
                                                                               Ma L. Experimental investigation of time dependent behavior of welded
The strength of the rock has approximately a                                       Tonopah spring tuff. Ann Arbor, USA: [s.n.], 2004.
negative exponential relation with time (Eq. (2)).                             Maranini E, Brignoli M. Creep behaviour of a weak rock: experimental
   (4) The failure modes of the deep marble are                                    characterization. International Journal of Rock Mechanics and
different under different loading conditions, such as                              Mining Sciences, 1999, 36 (1):127–138.
long- and short-term loadings. Under the short-term                            Munson D E. Constiutive model of creep in rock salt applied to
axial loading condition, rock sample presents tensile                              underground room closure. International Journal of Rock Mechanics
                                                                                   and Mining Sciences, 1997, 34 (2): 233–247.
failure when the compressive stress is imposed.
                                                                               Sun Jun. Rheology and engineering application of the geotechnical
Under the condition of long-term loading, rock
                                                                                   materials. Beijing: China Architecture and Building Press, 1999 (in
sample presents shear failure, followed by some                                    Chinese).
tensile failures. However, further study on the                                Sun Jun. Rock rheological mechanics and its advance in engineering
relationship between the failure mechanism and time                                applications. Chinese Journal of Rock Mechanics and Engineering,
is needed.                                                                         2007, 26 (6): 1 082–1 106 (in Chinese).
                                                                               The Professional Standards Compilation Group of People’s Republic of
References                                                                         China. DL/T 5368-2007 Code for rock tests of hydroelectric and
                                                                                   water conservancy engineering. Beijing: China Electric Power Press,
Chen Bingrui, Feng Xiating, Zeng Xionghui, Xiao Yaxun, Zhang Taizhao,              2007.
      Ming Huajun, Feng Guangliang. Real-time microseismic monitoring          Vyalov S S. Rheological fundamentals of soil mechanics. New York:
      and its characteristic analysis during TBM tunneling in deep-buried          Elsevier, 1986: 147–154.
      tunnel. Chinese Journal of Rock Mechanics and Engineering, 2010,         Wang Dengke, Liu Jian, Yin Guangzhi, Wei Xiaoji. Test study of creep
      29 (2): 349–357 (in Chinese).                                                properties of gas-bearing coal specimens under triaxial compression.
Cruden D M. The static fatigue of brittle rock under uniaxial compression.         Chinese Journal of Rock Mechanics and Engineering, 2010, 29 (2):
      International Journal of Rock Mechanics and Mining Sciences and              349–357 (in Chinese).
      Geomechanics Abstracts, 1974, 11 (2): 67–73.                             Wang Xuebin. Influence of the Poisson’s ratio on the failure mode and
Cui Xihai, Fu Zhiliang. Experimental study on rheology properties and              total deformation characteristics of the rock sample. Journal of
      long term strength of rocks. Chinese Journal of Rock Mechanics and           Geomechanics, 2007, 13 (3): 220–226 (in Chinese).
      Engineering, 2006, 25 (5): 1 021–1 024 (in Chinese).                     Wu Shiyong, Wang Ge. Rock mechanical problems and optimization for
Fabre G, Pellet F. Circumferential strain behavior during creep tests of
                                                                                   the long and deep diversion tunnels at Jinping II hydropower station.
      brittle rocks. International Journal of Rock Mechanics and Mining
                                                                                   Journal of Rock Mechanics and Geotechnical Engineering, 2011, 3
      Sciences, 2006, 43 (4): 950–960.
                                                                                   (4): 314–328.
Fan Qingzhong, Li Shucai, Gao Yanfa. Experimental study on creep
                                                                               Xu Hongfa. Time dependent behaviours of strength and elasticity modulus
      properties of soft rock under triaxial compression. Chinese Journal of
                                                                                   of weak rock. Chinese Journal of Rock Mechanics and Engineering,
      Rock Mechanics and Engineering, 2007, 26 (7): 1 381–1 385 (in
                                                                                   1997, 16 (3): 246–251 (in Chinese).
      Chinese).
Fu Zhiliang, Gao Yanfa, Ning Wei, Xu Jinpeng. Creep of anisotropic oil         Xu Weiya, Yang Shengqi, Yang Songlin, Xie Shouyi, Shao Jianfu, Wang
      shale. Journal of Mining and Safety Engineering, 2007, 24 (3):               Yifeng. Investigation on triaxial rheological mechanical properties of
      353–356 (in Chinese).                                                        green schist specimen (I): experimental results. Rock and Soil
Fujii Y, Kiyama T, Ishijima Y, Kodama J. Creep and time-dependent                  Mechanics, 2005, 26 (4): 531–537 (in Chinese).
      damage in argillaceous rocks. International Journal of Rock
                                                                               Yang Jianhui. A study of creep test phenomenon of sandstone under
      Mechanics and Mining Sciences, 1999, 36 (3): 323–337.
                                                                                   uniaxial compression. Journal of Shijiazhuang Railway Institute,
Hou Zhesheng, Gong Qiuming, Sun Zhuoheng. Primary failure types and
      their failure mechanisms of deep buried and intact marble at Jinping         1995, 8 (2): 77–80 (in Chinese).
      II hydropower station. Chinese Journal of Rock Mechanics and             Zhang Zhiliang, Xu Weiya, Wang Wei. Study of triaxial creep tests and its
      Engineering, 2011, 30 (4): 727–732 (in Chinese).                             nonlinear visco-elastoplastic creep model of rock from compressive
Li Huamin, Li Zhenhua, Su Chengdong. Testing study on creep                        zone of dam foundation in Xiangjiaba hydropower station. Chinese
      characteristics of marble. Chinese Journal of Rock Mechanics and
                                                                                   Journal of Rock Mechanics and Engineering, 2011, 30 (1): 132–140
      Engineering, 2004, 23 (22): 3 745–3 749 (in Chinese).
                                                                                   (in Chinese).
Li Yongsheng. Creep and relaxation of four kinds of rocks under uniaxial
      compression tests. Chinese Journal of Rock Mechanics and                 Zhao Yonghui, He Zhimin, Shen Mingrong. Testing study on rock
      Engineering, 1995, 14 (1): 39–47 (in Chinese).                               rheological characteristic of north anchor block of Runyang Bridge.
Liu Lipeng, Wang Xiaogang, Zhang Yizhong, Jia Zhixin, Duan Qingwei.                Rock and Soil Mechanics, 2003, 24 (4): 583–586 (in Chinese).