Research article

Three-dimensional inversion analysis of transient electromagnetic response signals of water-bearing abnormal bodies in tunnels based on numerical characteristic parameters


  • Received: 10 August 2022 Revised: 09 October 2022 Accepted: 21 October 2022 Published: 25 October 2022
  • The transient electromagnetic inversion of detection signals mainly depends on fast inversion in the half-space state. However, the interpretation results have a certain degree of uncertainty and blindness, so the accuracy and applicability of the three-dimensional full-space inversion need to be investigated. Two different three-dimensional full-space inversions were carried out. First, the numerical characteristic parameters of the response signals were extracted. Then, the correlations between the numerical characteristic parameters and physical parameters of the water-bearing abnormal bodies were judged, and the judgment criterion of the iterative direction was proposed. Finally, the inversion methods of the iterative algorithm and the BP neural network were utilized based on the virtual example samples. The results illustrate that the proposed numerical characteristic parameters can accurately reflect the response curve of the full-space surrounding rock. The difference in the numerical characteristic parameters was used to determine the update direction and correction value. Both inversion methods have their advantages and disadvantages. A single inversion method cannot realize the three-dimensional inversion of the physical parameters of water-bearing abnormal bodies quickly, effectively and intelligently. Therefore, the benefits of different inversion methods need to be considered to comprehensively select a reasonable inversion method. The results can provide essential ideas for the subsequent interpretation of the three-dimensional spatial response signals of water-bearing abnormal bodies.

    Citation: Yikang Xu, Zhaohua Sun, Wei Gu, Wangping Qian, Qiangru Shen, Jian Gong. Three-dimensional inversion analysis of transient electromagnetic response signals of water-bearing abnormal bodies in tunnels based on numerical characteristic parameters[J]. Mathematical Biosciences and Engineering, 2023, 20(1): 1106-1121. doi: 10.3934/mbe.2023051

    Related Papers:

    [1] Xiaoming Su, Jiahui Wang, Adiya Bao . Stability analysis and chaos control in a discrete predator-prey system with Allee effect, fear effect, and refuge. AIMS Mathematics, 2024, 9(5): 13462-13491. doi: 10.3934/math.2024656
    [2] Kottakkaran Sooppy Nisar, G Ranjith Kumar, K Ramesh . The study on the complex nature of a predator-prey model with fractional-order derivatives incorporating refuge and nonlinear prey harvesting. AIMS Mathematics, 2024, 9(5): 13492-13507. doi: 10.3934/math.2024657
    [3] Nehad Ali Shah, Iftikhar Ahmed, Kanayo K. Asogwa, Azhar Ali Zafar, Wajaree Weera, Ali Akgül . Numerical study of a nonlinear fractional chaotic Chua's circuit. AIMS Mathematics, 2023, 8(1): 1636-1655. doi: 10.3934/math.2023083
    [4] A. Q. Khan, Ibraheem M. Alsulami . Complicate dynamical analysis of a discrete predator-prey model with a prey refuge. AIMS Mathematics, 2023, 8(7): 15035-15057. doi: 10.3934/math.2023768
    [5] Xiao-Long Gao, Hao-Lu Zhang, Xiao-Yu Li . Research on pattern dynamics of a class of predator-prey model with interval biological coefficients for capture. AIMS Mathematics, 2024, 9(7): 18506-18527. doi: 10.3934/math.2024901
    [6] Weili Kong, Yuanfu Shao . The effects of fear and delay on a predator-prey model with Crowley-Martin functional response and stage structure for predator. AIMS Mathematics, 2023, 8(12): 29260-29289. doi: 10.3934/math.20231498
    [7] Asharani J. Rangappa, Chandrali Baishya, Reny George, Sina Etemad, Zaher Mundher Yaseen . On the existence, stability and chaos analysis of a novel 4D atmospheric dynamical system in the context of the Caputo fractional derivatives. AIMS Mathematics, 2024, 9(10): 28560-28588. doi: 10.3934/math.20241386
    [8] Yao Shi, Zhenyu Wang . Bifurcation analysis and chaos control of a discrete fractional-order Leslie-Gower model with fear factor. AIMS Mathematics, 2024, 9(11): 30298-30319. doi: 10.3934/math.20241462
    [9] Guilin Tang, Ning Li . Chaotic behavior and controlling chaos in a fast-slow plankton-fish model. AIMS Mathematics, 2024, 9(6): 14376-14404. doi: 10.3934/math.2024699
    [10] Xuyang Cao, Qinglong Wang, Jie Liu . Hopf bifurcation in a predator-prey model under fuzzy parameters involving prey refuge and fear effects. AIMS Mathematics, 2024, 9(9): 23945-23970. doi: 10.3934/math.20241164
  • The transient electromagnetic inversion of detection signals mainly depends on fast inversion in the half-space state. However, the interpretation results have a certain degree of uncertainty and blindness, so the accuracy and applicability of the three-dimensional full-space inversion need to be investigated. Two different three-dimensional full-space inversions were carried out. First, the numerical characteristic parameters of the response signals were extracted. Then, the correlations between the numerical characteristic parameters and physical parameters of the water-bearing abnormal bodies were judged, and the judgment criterion of the iterative direction was proposed. Finally, the inversion methods of the iterative algorithm and the BP neural network were utilized based on the virtual example samples. The results illustrate that the proposed numerical characteristic parameters can accurately reflect the response curve of the full-space surrounding rock. The difference in the numerical characteristic parameters was used to determine the update direction and correction value. Both inversion methods have their advantages and disadvantages. A single inversion method cannot realize the three-dimensional inversion of the physical parameters of water-bearing abnormal bodies quickly, effectively and intelligently. Therefore, the benefits of different inversion methods need to be considered to comprehensively select a reasonable inversion method. The results can provide essential ideas for the subsequent interpretation of the three-dimensional spatial response signals of water-bearing abnormal bodies.



    Throughout the paper, we work over an algebraically closed field k of characteristic zero. Let C be a nonsingular projective curve of genus g0, and L be a very ample line bundle on C. The complete linear system |L| embeds C into a projective space Pr:=P(H0(C,L)). For an integer k0, the k-th secant variety

    Σk=Σk(C,L)Pr

    of C in Pr is the Zariski closure of the union of (k+1)-secant k-planes to C.

    Assume that degL2g+2k+1. Then the k-th secant variety Σk can be defined by using the secant sheaf Ek+1,L and the secant bundle Bk(L) as follows. Denote by Cm the m-th symmetric product of C. Let

    σk+1:Ck×CCk+1

    be the morphism sending (ξ,x) to ξ+x, and p:Ck×CC the projection to C. The secant sheaf Ek+1,L on Ck+1 associated to L is defined by

    Ek+1,L:=σk+1,pL,

    which is a locally free sheaf of rank k+1. Notice that the fiber of Ek+1,L over ξCk+1 can be identified with H0(ξ,L|ξ). The secant bundle of k-planes over Ck+1 is

    Bk(L):=P(Ek+1,L)

    equipped with the natural projection πk:Bk(L)Ck+1. We say that a line bundle L on a variety X separates m+1 points if the natural restriction map H0(X,L)H0(ξ,L|ξ) is surjective for any effective zero-cycle ξX with length(ξ)=m+1. Notice that a line bundle L is globally generated if and only if L separates 1 point, and L is very ample if and only if L separates 2 points. Since degL2g+k, it follows from Riemann–Roch that L separates k+1 points. Then the tautological bundle OBk(L)(1) is globally generated. We have natural identifications

    H0(Bk(L),OBk(L)(1))=H0(Ck+1,Ek+1,)=H0(C,L),

    and therefore, the complete linear system |OBk(L)(1)| induces a morphism

    βk:Bk(L)Pr=P(H0(C,L)).

    The k-th secant variety Σk=Σk(C,L) of C in Pr can be defined to be the image βk(Bk(L)). Bertram proved that βk:Bk(L)Σk is a resolution of singularities (see [1,Section 1]).

    It is clear that there are natural inclusions

    C=Σ0Σ1Σk1ΣkPr.

    The preimage of Σk1 under the morphism βk is actually a divisor on Bk(L). Thus there exits a natural morphism from Bk(L) to the blowup of Σk along Σk1. Vermeire proved that B1(L) is indeed the blowup of Σ1 along Σ0=C ([3,Theorem 3.9]). In the recent work [2], we showed that Bk(L) is the normalization of the blowup of Σk along Σk1 ([2,Proposition 5.13]), and raised the problem asking whether Bk(L) is indeed the blowup itself ([2,Problem 6.1]). The purpose of this paper is to give an affirmative answer to this problem by proving the following:

    Theorem 1.1. Let C be a nonsingular projective curve of genus g, and L be a line bundle on C. If degL2g+2k+1 for an integer k1, then the morphism βk:Bk(L)Σk(C,L) is the blowup of Σk(C,L) along Σk1(C,L).

    To prove the theorem, we utilize several line bundles defined on symmetric products of the curve. Let us recall the definitions here and refer the reader to [2] for further details. Let

    Ck+1=C××Ck+1times

    be the (k+1)-fold ordinary product of the curve C, and pi:Ck+1C be the projection to the i-th component. The symmetric group Sk+1 acts on p1Lpk+1L in a natural way: a permutation μSk sends a local section s1sk+1 to sμ(1)sμ(k+1). Then p1Lpk+1L is invariant under the action of Sk+1, so it descends to a line bundle Tk+1(L) on the symmetric product Ck+1 via the quotient map q:Ck+1Ck+1. We have qTk+1(L)=p1Lpk+1L. Define a divisor δk+1 on Ck+1 such that the associated line bundle OCk+1(δk+1)=det(σk+1,(OCk×C)). Let

    Ak+1,L:=Tk+1(L)(2δk+1)

    be a line bundle on Ck+1. When k=0, we use the convention that T1(L)=E1,L=L and δ1=0.

    The main ingredient in the proof of Theorem 1.1 is to study the positivity of the line bundle Ak+1,L. Some partial results and their geometric consequences have been discussed in [2,Lemma 5.12 and Proposition 5.13]. Along this direction, we establish the following proposition to give a full picture in a general result describing the positivity of the line bundle Ak+1,L. This may be of independent interest.

    Proposition 1.2. Let C be a nonsingular projective curve of genus g, and L be a line bundle on C. If degL2g+2k+ for integers k,0, then the line bundle Ak+1,L on Ck+1 separates +1 points.

    In particular, if degL2g+2k, then Ak+1,L is globally generated, and if degL2g+2k+1, then Ak+1,L is very ample.

    In this section, we prove Theorem 1.1. We begin with showing Proposition 1.2.

    Proof of Proposition 1.2. We proceed by induction on k and . If k=0, then A1,L=L and degL2g+. It immediately follows from Riemann–Roch that L separates +1 points. If =0, then degL2g+2k. By [2,Lemma 5.12], Ak+1,L separates 1 point.

    Assume that k1 and 1. Let z be a length +1 zero-dimensional subscheme of Ck+1. We aim to show that the natural restriction map

    rz,k+1,L:H0(Ck+1,Ak+1,L)H0(z,Ak+1,L|z)

    is surjective. We can choose a point pC such that Xp contains a point in the support of z, where Xp is the divisor on Ck+1 defined by the image of the morphism CkCk+1 sending ξ to ξ+p. Let y:=zXp be the scheme-theoretic intersection, and Ix:=(Iz:IXp), which defines a subscheme x of z in Ck+1, where Iz and IXp are ideal sheaves of z and Xp in Ck+1, respectively. We have the following commutative diagram

    rz,k+1,L:H0(Ck+1,Ak+1,L)H0(z,Ak+1,L|z)

    where all rows and columns are short exact sequences. By tensoring with Ak+1,L and taking the global sections of last two rows, we obtain the commutative diagram with exact sequences

    rz,k+1,L:H0(Ck+1,Ak+1,L)H0(z,Ak+1,L|z)

    in which we use the fact that H1(Ak+1,L(Xp))=0 (see the proof of [2,Lemma 5.12]). Note that Ak+1,L(Xp)=Ak+1,L(p) and Ak+1,L|XpAk,L(2p), where we identify Xp=Ck.

    Since length(y)length(z)=+1 and degL(2p)2g+2(k1)+, the induction hypothesis on k implies that ry,k,L(2p) is surjective. On the other hand, if x=, which means that z is a subscheme of Xp, then trivially rx,k+1,L(p) is surjective. Otherwise, suppose that x. By the choice of Xp, we know that y is not empty, and therefore, we have length(x)length(z)1=. Now, degL(p)2g+2k+(1), so the induction hypothesis on implies that L(p) separates points. In particular, rx,k+1,L(p) is surjective. Hence rz,k+1,L is surjective as desired.

    Lemma 2.1. Let φ:XY be a finite surjective morphism between two varieties. If φ1(q) is scheme theoretically a reduced point for each closed point qY, then φ is an isomorphism.

    Proof. Note that φ is proper, injective, and unramifield. Then it is indeed a classical result that φ is an isomorphism. Here we give a short proof for reader's convenience. The problem is local. We may assume that X=SpecB and Y=SpecA for some rings A,B. We may regard A as a subring of B. For any qY, let p:=φ1(q)X. It is enough to show that the localizations A:=Amq and B:=Bmp are isomorphic. Let mq,mp be the maximal ideals of the local rings A,B, respectively. The assumption says that mqB=mp. We have

    B/AAA/mq=B/(mqB+A)=B/(mp+A)=0.

    By Nakayama lemma, we obtain B/A=0.

    We keep using the notations used in the introduction. Recall that C is a nonsingular projective curve of genus g0, and L is a very ample line bundle on C. Consider ξkCk and xC, and let ξ:=ξk+xCk+1. The divisor ξk spans a k-secant (k1)-plane P(H0(ξk,L|ξk)) to C in P(H0(C,L)), and it is naturally embedded in the (k+1)-secant k-plane P(H0(ξ,L|ξ)) spanned by ξ. This observation naturally induces a morphism

    αk,1:Bk1(L)×CBk(L).

    To see it in details, we refer to [1,p.432,line –5]. We define the relative secant variety Z=Zk1 of (k1)-planes in Bk(L) to be the image of the morphism αk,1. The relative secant variety Z is a divisor in the secant bundle Bk(L), and it is the preimage of (k1)-th secant variety Σk1 under the morphism βk. It plays the role of transferring the codimension two situation (Σk,Σk1) into the codimension one situation (Bk(L),Z). We collect several properties of Z.

    Proposition 2.2. ([2,Proposition 3.15,Theorem 5.2,and Proposition 5.13]) Recall the situation described in the diagram

    αk,1:Bk1(L)×CBk(L).

    Let H be the pull back of a hyperplane divisor of Pr by βk, and let IΣk1|Σk be the ideal sheaf on Σk defining the subvariety Σk1. Then one has

    1. OBk(L)((k+1)HZ)=πkAk+1,L.

    2. Riβk,OBk(L)(Z)={IΣk1|Σkifi=00ifi>0.

    3. IΣk1|ΣkOBk(L)=OBk(L)(Z).

    As a direct consequence of the above proposition, we have an identification

    H0(Ck+1,Ak+1,L)=H0(Σk,IΣk1|Σk(k+1)).

    We are now ready to give the proof of Theorem 1.1.

    Proof of Theorem 1.1. Let

    b:˜Σk:=BlΣk1ΣkΣk

    be the blowup of Σk along Σk1 with exceptional divisor E. As IΣk1|ΣkOBk(L)=OBk(L)(Z) (see Proposition 2.2), there exists a morphism α from Bk(L) to the blowup ˜Σk fitting into the following commutative diagram

    b:˜Σk:=BlΣk1ΣkΣk

    We shall show that α is an isomorphism.

    Write V:=H0(Σk,IΣk1|Σk(k+1)). As proved in [2,Theorem 5.2], IΣk1|Σk(k+1) is globally generated by V. This particularly implies that on the blowup ˜Σk one has a surjective morphism VO˜ΣkbOΣk(k+1)(E), which induces a morphism

    γ:˜ΣkP(V).

    On the other hand, one has an identification V=H0(Ck+1,Ak+1,L) by Proposition 2.2. Recall from Proposition 1.2 that Ak+1,L is very ample. So the complete linear system |V|=|Ak+1,L| on Ck+1 induces an embedding

    ψ:Ck+1P(V).

    Also note that α(bOΣk(k+1)(E))=βkOΣk(k+1)(Z)=πkAk+1,L by Proposition 2.2. Hence we obtain the following commutative diagram

    ψ:Ck+1P(V).

    Take an arbitrary closed point x˜Σk, and consider its image x:=b(x) on Σk. There is a nonnegative integer mk such that xΣmΣm1Σk. In addition, the point x uniquely determines a degree m+1 divisor ξm+1,x on C in such a way that ξm+1,x=ΛC, where Λ is a unique (m+1)-secant m-plane to C with xΛ (see [2,Definition 3.12]). By [2,Proposition 3.13], β1k(x)Ckm and πk(β1k(x))=ξm+1,x+CkmCk+1. Consider also x:=γ(x) which lies in the image of ψ. As ψ is an embedding, we may think x as a point of Ck+1. Now, through forming fiber products, we see scheme-theoretically

    α1(x)π1k(x)β1k(x).

    However, the restriction of the morphism πk on β1k(x) gives an embedding of Ckm into Ck+1. This suggests that π1k(x)β1k(x) is indeed a single reduced point, and so is α1(x). Finally by Lemma 2.1, α is an isomorphism as desired.



    [1] M. S. Zhdanov, Electromagnetic geophysics: Notes from the past and the road ahead, Geophys., 75 (2010), 557–560. https://doi.org/10.1190/1.3483901 doi: 10.1190/1.3483901
    [2] E. Auken, A. V. Christiansen, C. Kirkegaard, G. Fiandaca, C. Schamper, A. A. Behroozmand, et al., An overview of a highly versatile forward and stable inverse algorithm for airborne, ground-based and borehole electromagnetic and electric data, Explor. Geophys., 46 (2015), 223–235. https://doi.org/10.1071/EG13097 doi: 10.1071/EG13097
    [3] W. Qian, H. Li, J. Yu, Z. Gu, Theoretical and experimental investigation of vehicle-mounted transient electromagnetic method detection for internal defects of operational tunnels, Appl. Sci., 11 (2021), 6909. https://doi.org/10.3390/app11156906 doi: 10.3390/app11156906
    [4] Z. Li, T. Qi, S. Qin, W. Qian, The analysis of the early electromagnetic response of the receiving coil and its application at close-range TEM detection, J. Appl. Geophys., 193 (2021), 104409. https://doi.org/10.1016/j.jappgeo.2021.104409 doi: 10.1016/j.jappgeo.2021.104409
    [5] S. C. Constable, R. L. Parker, C. G. Constable, Occam's inversion: A practical algorithm for generating smooth models from electromagnetic sounding data, Geophysics, 52 (1987), 289–300. https://doi.org/10.1190/1.1442303 doi: 10.1190/1.1442303
    [6] A. V. Christiansen, N. B Christensen, A quantitative appraisal of airborne and ground-based transient electromagnetic (TEM) measurements in Denmark, Geophysics, 68 (2003), 523–534. https://doi.org/10.1190/1.1567220 doi: 10.1190/1.1567220
    [7] Y. Chang, M. Xiao, Y. Wu, Studies on initial parameter selection of one-dimensional inversion for transient electromagnetic Data, Geophys. Prospect. Petrol., 45 (2010), 295–298. https://doi.org/10.13810/j.cnki.issn.1000-7210.2010.02.016 doi: 10.13810/j.cnki.issn.1000-7210.2010.02.016
    [8] P. Yogeshwar, B. Tezkan, Analyzing two-dimensional effects in central loop transient electromagnetic sounding data using a semi-synthetic tipper approach, Geophys. Prospect., 66 (2018), 444–456. https://doi.org/10.1111/1365-2478.12520 doi: 10.1111/1365-2478.12520
    [9] E. Haber, D. W. Oldenburg, R. Shekhtman, Inversion of time domain three-dimensional electromagnetic data, Geophys. J. Int., 171 (2010), 550–564. https://doi.org/10.1111/j.1365-246X.2007.03365.x doi: 10.1111/j.1365-246X.2007.03365.x
    [10] W. Zhao, T. Yan, Z. Gao, Magnetotelluric nonlinear conjugate gradient inversion experiments: an example from data acquired in the Jarud Basin, Inner Mongolia, China, Prog. Geophys., 29 (2014), 2128–2135. https://doi.org/10.6038/pg20140520 doi: 10.6038/pg20140520
    [11] Z. Y. Zhou, The 3-D numerical implementation of the whole tunnel space with TEM, Southwest Jiaotong University, (2014), 89–112.
    [12] Y. Li, T. Qi, B. Lei, Z. Li, W. Qian, An iterative inversion method using transient electromagnetic data to predict water-filled caves during the excavation of a tunnel, Geophys., 84 (2019), E89–E103. https://doi.org/10.1190/geo2018-0253.1 doi: 10.1190/geo2018-0253.1
    [13] X. Sun, Y. Wang, X. Yang, Y. Wang, Three-dimensional transient electromagnetic inversion with optimal transport, J. Inverse Ill-posed Probl., 30 (2021), 549–565. https://doi.org/10.1515/jiip-2020-0159 doi: 10.1515/jiip-2020-0159
    [14] X. Wang, N. You, Q. Di, J. Deng, Y. Chang, 3-D parallel inversion of multichannel transient electromagnetic data using a moving footprint, Geophys. J. Int., 226 (2021), 1783–1799. https://doi.org/10.1093/gji/ggab187 doi: 10.1093/gji/ggab187
    [15] M. Liu, H. Cai, H. Yang, Y. Xiong, X. Hu, Three dimensional joint inversion of ground and semi-airborne transient electromagnetic method, Chin. J. Geophys., 65 (2022), 3997–4011. https://doi.org/10.6038/cjg2022P0897 doi: 10.6038/cjg2022P0897
    [16] D. Tan, T. Qi, Application of expert system in transient electromagnetic inversion, Technol. Highway Transp., 11 (2009), 124–127.
    [17] P. Ertan, Y. Türker, K. Yekta, Ö. Coşkun, Application of particle swarm optimization on self-potential data, J. Appl. Geophys., 75 (2011), 305–318. https://doi.org/10.1016/j.jappgeo.2011.07.013 doi: 10.1016/j.jappgeo.2011.07.013
    [18] M. Wang, G. Liu, D. Wang Y. Liu, Application of artificial bee colony algorithm in the version of transient electromagnetic sounding files, Prog. Geophys., 30 (2015), 133–139. https://doi.org/10.6038/pg20150120 doi: 10.6038/pg20150120
    [19] X. Liang, T. Qi, Z. Jin, W. Qian, Hybrid support vector machine optimization model for inversion of tunnel transient electromagnetic method, Math. Biosci. Eng., 17 (2020), 3998–4017. https://doi.org/10.3934/mbe.2020221 doi: 10.3934/mbe.2020221
    [20] H. Sun, Q. Wu, R. Chen, H. Li, Q. Fan, G. Liu, et al., Experimental study on transient electromagnetic responses to shallow karst, Chin. J. Rock Mech. Eng., 37 (2018), 652–661. https://doi.org/10.13722/j.cnki.jrme.2017.0038 doi: 10.13722/j.cnki.jrme.2017.0038
    [21] W. Qian, T. Qi, X. Liang, S. Qin, Z. Li, Y. Li, Vehicle-borne transient electromagnetic numerical characteristic parameter of water-bearing body behind tunnel linings, Math. Probl. Eng., 19 (2020), 8514913. https://doi.org/10.1155/2020/8514913 doi: 10.1155/2020/8514913
    [22] S. L. Butler, Z. Zhang, Forward modeling of geophysical electromagnetic methods using Comsol, Comput. Geosci., 87 (2016), 1–10. https://doi.org/10.1016/j.cageo.2015.11.004 doi: 10.1016/j.cageo.2015.11.004
    [23] Y. Qi, H. El-Kaliouby, A. Revil, A. Ahmed, A. Ghorbani, J. Li, Three-dimensional modeling of frequency-domain and time-domain electromagnetic methods with induced polarization effects, Comput. Geosci., 124 (2019), 85–92. https://doi.org/10.1016/j.cageo.2018.12.011 doi: 10.1016/j.cageo.2018.12.011
    [24] S. Li, Y. Liu, W. Sun, Intelligent Calculation and Parameter Inversion, Science Press, (2008), 29–38.
    [25] B. Liu, H. Guo, MATLAB Neural Network Super-Learning Manual, China Communications Press, (2014), 143–158.
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1869) PDF downloads(94) Cited by(0)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog