Une nouvelle approche de modélisation des milieux géologiques fracturés représentés par un modèle conceptuel de fractures discrètes et déterministes est présentée dans cette thèse. L'objectif principal de l'étude est de reproduire l'hétérogénéité et la complexité des milieux poreux fracturés dans un modèle géométrique tridimensionnel afin d'effectuer des simulations numériques dans le but d'améliorer les capacités de modélisation en hydrogéologie. Ceci est réalisé à travers le couplage entre une plateforme de modélisation géologique (GOCAD) et un code numérique (HydroGeoSphere). Les principaux défis à relever sont: la représentation géométrique du réseau de fractures, la sélection d'un maillage approprié pour la discrétisation spatiale du domaine de simulation et l'adaptation du code numérique à ce maillage. La nouvelle approche est basée sur une première phase de modélisation géologique 3D, suivie par la génération d'un maillage tétraédrique 3D et par la simulation numérique de l'écoulement souterrain en conditions saturées et du transport de solutés. En général, le maillage tétraédrique s'avère plus adéquat que les maillages de blocs ou de prismes pour discrétiser les geometries complexes telles que les milieux fracturés. De plus, une définition alternative du maillage "dual", qui est essentiel pour appliquer la méthode numérique élément finis - volume de contrôle utilisée par HydroGeoSphere, est analysée et intégrée dans le code numérique. Le code numérique proposé est d'abord vérifié par l'intermédiaire de simples scénarios de simulation dont les solutions, analytiques et numériques, sont déjà connues. La complexité des simulations est augmentée de façon graduelle. L'approche de modélisation est finalement appliquée au site Olkiluoto (Finlande) où un laboratoire de recherche souterrain est en construction afin d'évaluer la faisabilité du stockage géologique profond de déchets nucléaires à haute activité. Les techniques de modélisation géologique mises au point permettent de modéliser facilement la géométrie des fractures identifiées à travers la caractérisation géologique in situ. De plus, le modèle numérique s'avère adéquat pour la simulation de l'écoulement et du transport de solutés dans ce site complexe. Ce travail de recherche présente une contribution au développement des techniques de modélisation hydrogéologique des milieux fracturés.

La région de l’Outaouais (13760 km²), au Québec, connait une forte croissance démographique et les eaux souterraines sont davantage sollicitées. La perméabilité du socle rocheux, essentiellement cristallin, est contrôlée par des fractures entrecroisées et donc hétérogène. Pour identifier des zones-cibles de forage, nous avons confronté des données de capacité spécifique du Système d’Information Hydrogéologique du Québec, validés à l’aide d’essais de pompage réalisés par l’Université Laval, avec des caractéristiques du roc, parmi lesquelles la lithologie, l’épaisseur de roc exploitée par les puits et la présence de linéaments. Il apparait que la texture et la minéralogie des roches ont une influence sur le développement de fractures perméables. On constate aussi que ces fractures se raréfient en profondeur et qu’il n’est donc pas toujours opportun de forer des puits à des profondeurs importantes (< 80 m). En revanche, aucune corrélation entre la présence de linéaments et le rendement des captage situés à proximité n’a pu être observée. / The Outaouais Region (13760 km²) is one of the most rapidly growing communities in Quebec, and where groundwater is increasingly being relied upon as a water supply. The permeability of the mostly crystalline bedrock is controlled by interconnected fractures and then heterogeneous. In order to identify target areas for installing new pumping wells, specific capacity data from the SIH database, which were validated using pumping tests carried out by Université Laval, were compared with rock characteristics including lithology, thickness of rock exploited by each well and proximity to lineaments. The texture and mineralogy of the rock seemed to have an important influence on the development of permeable fractures. We also found that the fracture density decreases with depth and therefore it is not always profitable to drill to deeper depths (> 80 m) to try to increase productivity. In contrast, proximity to lineaments didn’t seem to have any impact on the productivity of wells.

La caractérisation des discontinuités mineures à l’échelle d’une galerie est cruciale en ingénierie minière, car elles contribuent à augmenter localement la perméabilité d’un massif rocheux et engendrer des instabilités structurales. L’utilisation de modèles discrets de joints (DFN) permet une représentation fidèle de ces conditions structurales. Cependant, les cartographies linéaires traditionnellement utilisées pour acquérir les intrants sont coûteuses en temps et entrainent des biais de mesures. Des logiciels disponibles commercialement reposant sur la photogrammétrie permettent d’opter pour la cartographie numérique surfacique. Les travaux réalisés permettent de confirmer que cette méthodologie réduit les biais engendrés par la cartographie linéaire. Lorsqu’intégrée à la calibration des modèles DFN de Baecher et de Veneziano, l’approche surfacique permet une représentativité statistique des conditions structurales d’orientation, de longueurs traces et d’intensité de fracturation observées sur le terrain. L’utilisation de DFN calibrés permet d’évaluer la conductivité hydraulique et le potentiel de génération d’instabilités structurales à l’échelle d’une galerie minière.

Les coûts de broyage dans l'industrie minière représentent une grande dépense pour une usine de traitement de minerai. En effet, ceux-ci peuvent atteindre plus de 50% des coûts énergétique du traitement. Le broyage pourrait être optimisé par une meilleure connaissance de la distribution spatiale de la broyabilité des massifs rocheux. Sur les sites miniers, on caractérise habituellement la broyabilité d'un gisement par des tests d(^ broyage sur des échantillons prélevés à plusieurs endroits dans la mine. Cependant, puisque ces tests sont coûteux, relativement peu de roc est ainsi caractérisé. Par ailleurs, certaines informations géologiques et géotechniques sont abondantes et facilement accessibles sur les sites miniers. Ce mémoire vise d'une part à montrer les relations statistiques entre les différentes méthodes usuelles de caractérisation de la roche intacte dans les opérations minières, et d'autre part à démontrer comment le modèle de blocs contenant l'information géologique référencée peut être utilisée pour maximiser notre connaissance de la broyabilité d'un gisement. Pratiquement, ce mémoire étudie les liens entre les résultats des méthodes de caractérisation de la roche en broyage, en mécanique des roches (propriétés du roc intact) et en géologie de 13 échantillons de minerai d'or qui proviennent de deux fosses de la mine Troilus : la fosse 87 et la fosse J4. Les échantillons, qui ont été caractérisés dans le cadre de ce projet, ont été prélevés de 3 trous de forage et de 2 niveaux d'exploitation de la fosse 87, de 3 trous de forage de la fosse ,J4, et de l’alimentation du broyeur en provenance de la fosse 87 (!t de la fosse J4. Ces échantillons ont été prélevés avant la définition du projet de maîtrise. Ils ont été prélevés dans le cadre d'une caractérisation métallurgique de certaines zones de la mine Troilus. Le premier travail réalisé dans le cadre de ce projet a été la caractérisation géologique des échantillons. Cette caractérisation des échantillons a été réalisée sur la base de la description visuelle du type de roche, de la quantité des minéraux majeurs et sur la base des valeurs de facteur de qualité de la roche RQD. Les résultats de la caractérisation géologique ont démontré que les échantillons testés provenaient de deux types de roches : la diorite et la roche intrusive felsique. Les valeurs de RQD ont montré que le massif de Troilus est de bonne qualité. Le deuxième travail a été la caractérisation de la broyabilité des échantillons. Les tests de broyabilité qui ont été choisis sont le test de Bond pour le broyage à boulet, le test de Bond pour le broyage à barres, le test de Bond à basse énergie qui permet de calculer l'indice Wic, le test d'impact qui permet d'évaluer la résistance à l'impact A*b et le test d'abrasion qui permet d'évaluer la résistance à l'abrasion ta- Les résultats de ces tests ont démontré que : - Le minerai de Troilus présente une résistance faible au broyage fin, une valeur élevée de Wic, une résistance à l'impact et à l'abrasion élevée. - La roche diorite est plus résistante au broyage fin et l'impact à faible énergie que le mélange de roche diorite et roche felsique. - La roche diorite est moins résistante à l'impact (A*b) et à l'abrasion {ta) que le mélange de diorite et de roche felsique. - Pour chaque type de roche, les indices de broyages présentent des corrélations acceptables. Le troisième travail a été la caractérisation de la roche par les tests de mécanique de roche. Les tests de mécanique des roches qui ont été choisis sont le test de la résistance en compression uniaxiale, le test de la résistance en tension et le test du double poinçonnement. Les résultats de ces tests ont montré que : - la résistance à la compression est variable pour un même échantillon et pour un même type de roche. - la roche diorite présente un indice de double poinçonnement inférieur au mélange de roche diorite et roche felsique. - la résistance en tension de la roche diorite est inférieure à celle du mélange diorite et roche felsique. - La résistance en tension et l'indice de double poinçonnement ne présentent pas une forte corrélation avec la résistance en compression uniaxiale. - La résistance en tension est fortement corrélée avec l'indice de double poinçonnement dans le cas du mélange. - La résistance en tension est fortement corrélée avec la compression uniaxiale pour les deux types de roches. L'analyse globale des résultats des tests de la broyabilité et des tests de mécanique des roches montre que : - La compression uniaxiale et la résistance en tension de la roche diorite sont liées à sa résistance à l'impact. - L'indice Wi BM est lié avec la résistance en compression uniaxiale et avec l'indice de double poinçonnement du mélange diorite et roche flelsique. - La broyabilité (A*b et i„ ) est liée au RQD moyen dans le cas où le RQD est défini par les 4 valeurs localisées à moins de 30m. - Les roches qui proviennent d'un massif rocheux de bonne qualité résistent mieux à l'impact et à l'abrasion. L'intégration de la relation entre le facteur RQD et les deux résistances de la roche (à l'impact et à l'abrasion) dans le modèle de blocs a permis d'attribuer la broyabilité aux différents blocs d'une fosse conceptuelle. Cette dernière application a montré qu'il est possible d'utiliser le modèle de blocs et l'information géologique disponible pour représenter la distribution dans l'espace de la broyabilité d'une mine à ciel ouvert.

La déshydratation, l’oxydation et la fracturation sont les trois principaux processus couplés subis par l’argilite callovo-oxfordienne recoupée par les galeries du Laboratoire Souterrain de Meuse/Haute-Marne de l'Andra lors de l'excavation de galeries. Une série de forages à différents temps après l'ouverture de la galerie a permis d'étudier l'évolution des perturbations en fonction du temps, jusqu'à 18 mois. Les profils de teneur en eau et les analyses isotopiques montrent une déshydratation en surface des parois d'environ 50%. Le front de déshydratation est rapidement sub-stationnaire car il résulte d'un équilibre entre la recharge par l'eau de la formation et l'évaporation imposée par la pression partielle de l'eau dans l'air de la galerie. La zone perturbée affecte principalement les premiers 70 cm, et ne se fait plus sentir à partir de 1 m 50. L’oxydation des pyrites provoque la formation de fluides sulfatés acides qui dissolvent les minéraux environnants, libèrent dans la solution des ions Ca2+, K+ et Na+ puis lors de la déshydratation de l’argilite, l'évaporation de l'eau porale induit la précipitation de minéraux sulfatés principalement. La dégradation de la matière organique dans les conditions des galeries montre une dégradation limitée, due soit à une oxydation très discrète, soit à une dégradation lente non détectable aux temps courts. Une série d'expérience de mise en équilibre d'eau avec une roche plus ou moins oxydée a permis d'évaluer l'impact de la présence des phases minérales néoformées sur les caractéristiques chimiques des eaux de réhydratation. Cette approche expérimentale a été confrontée aux résultats de modélisation géochimique. Les résultats obtenus permettent de mieux préciser les caractéristiques géométriques, minéralogiques et physico-chimiques de cette zone perturbée qui est susceptible d'évoluer par la suite en condition de réhydratation et d'élévation de température lors du stockage / Dehydration, oxidation, and fracturing are the three main and coupled processes underwent by Callovian-Oxfordian argillite cut by the galleries of the Laboratoire Souterrain of Meuse / Haute-Marne of Andra during the excavation. Several boreholes drilled at different times after the opening of the gallery were used to study the evolution of disturbances with time up to 18 months. Profiles of water content and isotopic analysis show dehydration at the surface of the walls of about 50%. The area of dehydration is rapidly sub-stationary because it results from a balance between the suction by the pore water and the evaporation imposed by the partial pressure of water in the air of the gallery. The disturbed area mainly affects the first 70 cm, and not exceeds 1 m 50. The oxidation of pyrite causes the formation of acid sulphate fluids which dissolve the surrounding minerals in the solution, release Ca2+, K+ and Na+ in the pore water then during the dehydration of argillite, the evaporation of the pore water leads to the precipitation of sulphate minerals mainly. The degradation of organic matter in the galleries shows a limited deterioration due to a very moderate oxidation or the slow degradation is not detectable at short times. Several experiments of equilibration between more or less oxidized samples and water were used to evaluate the impact of the presence of new-formed mineral phases (sulphates) on the chemical composition of the rehydration water. This experimental approach has been compared to results of geochemical modelling. The results obtained allow us to define geometrical, mineralogical and physico-chemical properties of the disturbed zone which could evolve thereafter during rehydration and temperature rise in storage condition

La caractérisation de l'agencement spatial des propriétés hydrauliques est essentielle pour prédire les écoulements et le transport des solutés dans les milieux hétérogènes. Les méthodes de tomographie hydraulique, principalement développées pour estimer les propriétés des milieux poreux, n'ont qu'une faible résolution spatiale qui ne reflète pas la vraie hétérogénéité des distributions de fractures des milieux fracturés. Le principal objectif de cette thèse est de développer une nouvelle méthode d'inversion spécifique pour imager les propriétés hydrauliques et de transport des milieux fracturés à l'échelle du site. Pour atteindre ces objectifs, des expériences in situ ainsi qu'une nouvelle approche de modélisation inverse sont proposées, notamment en utilisant la température comme marqueur des écoulements. Nous proposons tout d'abord la tomographie d'écoulement basée sur des tests séquentiels de débimétrie entre puits, comme une nouvelle approche pour caractériser la connectivité des fractures ainsi que leur transmissivité. À partir de simulations numériques reproduisant des cas d'études synthétiques, nous montrons que l'approche par tomographie réduit significativement l'incertitude sur les paramètres estimés, et fournit une caractérisation détaillée du réseau de fracture sans requérir à l'utilisation d'obturateurs hydrauliques. Nous montrons ensuite comment les mesures de température peuvent être utilisées pour quantifier les écoulements dans les milieux fracturés. Le grand intérêt d'utiliser la température est d'obtenir facilement et de façon continue en puits des profils de température. En utilisant un modèle numérique d'écoulement et de transfert de chaleur à l'échelle du puits, une méthode d'inversion pour estimer les vitesses d'écoulement dans le puits à partir des données de température est proposée. Nous couplons ensuite les deux approches présentées précédemment dans une nouvelle approche expérimentale consistant en des enregistrements séquentiels de température dans un puits dans des conditions de pompage entre puits. L'application de cette approche de tomographie en température sur le site de Stang er Brune montre des résultats encourageants pour l'identification du réseau global de connectivité et des zones d'écoulement principales. Enfin, nous discutons de l'intérêt d'utiliser la chaleur comme traceur par rapport à l'utilisation de traceurs classiques. Nous montrons que réaliser des tests de traçage thermiques en milieu fracturé fournit des contraintes supplémentaires importantes sur les propriétés de transport du milieu. / The accurate characterization of distribution of hydraulic properties and connectivity distribution is essential to predict flow and transport in fractured media. Classical approaches were developed for homogeneous aquifers and result in smooth tomograms that often do not match true heterogeneity distribution of fractured media. The main goal of this thesis is to develop new inverse approaches specifically for imaging hydraulic and transport properties in fractured media at the field-scale. To attain this objective new in situ measurement methods as well as new inverse modelling frameworks are proposed.We first propose flow tomography (i.e., sequential cross-borehole flowmeter tests) as a new approach for characterizing fracture connectivity and transmissivities. Based on a discrete fracture network approach, we present a general method to invert flow tomography data. From synthetic case studies, we show that the tomographic approach reduces significantly the uncertainty on the parameter estimation. Flow tomography approach provides detailed characterization on fracture networks without the necessity of using packers. We then study the contribution of temperature measurements for quantifying flow in fractured media. The advantage of using temperature data is that temperature profiles can be obtained more easily and continuously in space, compared to flowmeter profiles. Using a numerical model of flow and heat transfer at the borehole scale, a method to invert temperature measurements to derive borehole flow velocities was proposed. We then couple the two previously proposed approaches in a new experimental approach which we call temperature tomography. This experiment consists of sequential borehole temperature logging under cross-borehole flow conditions. The full inverse framework is then presented to interpret temperature tomography experiments. Application of the temperature tomography approach to Stang er Brune field site showed encouraging results for the identification of general connectivity patterns and transmissivities of the main flowpaths. Finally, we explore the interest of using push-pull thermal tracer tests. Through field experiments and numerical modelling, we demonstrate that conducting push-pull heat tracer tests provide important constraints on the effective transport behavior.

La gestion d'un aquifère repose sur l'étude de ses propriétés hydrauliques. Nous proposons une approche computationnelle afin de caractériser la structure physique d'un aquifère et les propriétés hydrauliques associées, en focalisant sur l'étude des milieux hétérogènes fracturés où les écoulements sont fortement chenalisés. Ainsi, une première partie cherche à définir le modèle de fracturation sous-jacent à deux sites naturels suédois à partir d'observations corrigées des biais d'échantillonnage. La modélisation choisie repose sur une distribution des longueurs en loi de puissance et sur un terme de densité dépendant des orientations. Les propriétés hydrauliques des milieux hétérogènes, poreux et fracturés, sont ensuite caractérisées grâce à deux indicateurs statistiques décrivant le degré de chenalisation des écoulements à partir de la distance entre deux chenaux principaux et de la longueur efficace des chenaux. Ces indicateurs permettent de relier les propriétés physiques du milieu (hétérogénéité, organisation) à ses propriétés hydrauliques. Finalement, nous tentons d'identifier les composantes principales de la structure de perméabilité à partir de données de charge hydraulique en résolvant le problème inverse associé. On montre ainsi qu'il est possible d'identifier les principales structures contrôlant les écoulements à condition d'utiliser des données contenant une information sur la perméabilité pertinente à l'échelle du site, une paramétrisation représentative des structures principales et une méthodologie appropriée. Ces différentes approches permettent de définir un modèle pertinent des propriétés hydrauliques d'un milieu complexe.

L'étude de l'écoulement, du transport de masse et des échanges isotopiques en milieu fracturé est nécessaire à la compréhension des gîtes minéraux d'affiliation hydrothermale. L'application du modèle de simulation numérique HydroSphère, intégrant le transport de masse par advection, diffusion, et dispersion hydrodynamique, ainsi que les échanges isotopiques, à l'écoulement tridimensionnel dans une matrice discrètement fracturée, permet d'approfondir les connaissances fondamentales sur ces sujets. Des simulations étudiant différentes configurations de fracture et conditions limites hydrogéologiques ont été faites. L'étude des patrons d'isoplèthes de O pour différents types de fractures et de systèmes fracturés ayant un degré croissant de complexité, révèle que les fractures dominent les mécanismes de transport et d'échange isotopique. C'est donc leur nombre, leur orientation, leurs propriétés hydrauliques et leur connectivité qui déterminent la forme des patrons d'isoplèthes modélisés. Il est possible d'interpréter certaines asymétries dans la géométrie des patrons d'isoplèthes pour déterminer la direction d'écoulement des fluides hydro thermaux. / An understanding of fluid flow, mass transport and isotopic exchange in fractured rock is required to understand the origin of hydrothermal deposits. The numerical model HydroSphere simulates advection, molecular diffusion, mechanical dispersion and isotopic exchange in a 3D discretely fractured media. HydroSphere was used to better understand the processes of transport and exchange in fractured rocks. Study of 18O isopleth patterns for different types of fractures and fracture networks with a range of complexity, shows that fracture properties and geometry control mass transport and isotopic exchange. The hydraulic properties, as well as the number, spacing and connectivity of fractures determine the isotopic patterns. Asymmetries in the geometry of isotopic patterns could be used to determine the direction of hydrothermal fluid flow.

La déformation due au fluage d'un massif rocheux autour d'un tunnel a été rencontrée fréquemment. Ce phénomène est évident où il y a des tunnels creusés dans la roche tendre, des masses rocheuses faible et fortement cisaillées, ou des massifs rocheux soumis à des contraintes in-situ élevées. La déformation due au fluage se produit fréquemment au moment d’excavation des tunnels longs où il y a des failles et des zones fracturées et cisaillées. Ce phénomène peut causer différents dommages sur des systèmes de soutènement en raison de la déformation excessive et des effondrements. La déformation excessive impose une ré-excavation de la section du tunnel, qui monte le coût supplémentaire, la durée de la réalisation du projet et le risque de la sécurité sur le projet. En plus, comme la stabilité de terrain est dans un état critique durant la ré-excavation, une petite négligence peut conduire à une grande caverne. Bien que la déformation de fluage est commune dans un massif rocheux à une faible résistance dans un tunnel très profond, mais ce phénomène a été observé dans des tunnels peu profonds. Une bonne compréhension des déformations causées par une excavation souterraine requiert la connaissance de l'interaction roche-support et l'interprétation des données de terrain. Auparavant, l’objet principal de la surveillance effectuée durant la construction du tunnel était des mesures de la pression au terrain imposé sur le revêtement du tunnel. Mais aujourd’hui, les méthodes modernes de construction de tunnel se concentrent sur la surveillance des déplacements pendant et après la construction. Afin de déterminer des déformations dans les tunnels, Panet et Sulem ont supposé que "Le tunnel a une section transversale circulaire et le milieu est homogène et isotrope, aussi le tunnel est suffisamment profond pour considérer que la distribution des contraintes est homogène". Mais dans le cas quasi réel, la distribution de la contrainte autour du tunnel est hétérogène et anisotrope. Dans cette étude, pour la modification des équations Panet et Sulem, certaines équations sont proposées en cas de matériau hétérogène et anisotrope pour généraliser le problème. La galerie de force motrice Seymareh a été considérée comme l’étude de cas. Celle-ci est une partie du conduit d’eau dans le projet de centrale électrique du barrage Seymareh. Ce projet est situé à l'ouest de l'Iran. Les données de surveillance de la galerie de force motrice sont collectées au moment de l’excavation du tunnel, et sont comparées avec les résultats de la modélisation numérique et de la solution analytique. Cette comparaison montre que les résultats des données expérimentales obtenues par la surveillance sont très proches des résultats de la solution analytique, mais il y a une différence entre les deux et la modélisation numérique. Il était prévisible, car l’effet d’autres activités comme l’excavation des tunnels verticaux n’est pas prise en compte dans l’analyse numérique et aussi dans la solution analytique. Il est évident que les autres activités comme l’excavation des tunnels verticaux et l'excavation du tunnel principal vers deux directions opposées, peuvent affecter sur les résultats de la surveillance. D'autre part, les données initiales utilisées dans l'analyse numérique et la solution analytique ne sont pas tout à fait exactes, car elles sont obtenues en tant que représentatives du massif rocheux de la région, mais pas pour une section particulière. Toutefois, le but de cette étude est le développement d'une solution analytique de la déformation dans les tunnels sur les conditions générales et la poursuite de cette étude pourra être plus développée. / The creep deformation of a rock mass around a tunnel has been encountered frequently. It is particularly common in tunnels excavated in soft rock, heavily sheared weak rock masses or rock masses subjected to high in-situ stresses. Creep deformation in fault and shear fractured zones are one of the frequently encountered difficulties in long tunnel construction, which tend to cause failure of supporting systems due to excessive deformation and cavern. Excessive deformation would necessitate re-mining of the tunnel cross section, thus imposing impacts such as extra cost, extended time schedule and safety risk on the project. Furthermore, as the ground stability is in critical condition during re-mining, the slightest negligence would lead to major cavern. Although creep deformation is common to extremely poor rock mass under high overburden in a tunnel alignment, but however this phenomenon is not limited to tunnels with high overburden. A good understanding of the deformations caused by an underground excavation requires simultaneously knowledge of the rock-support interaction and interpretation of field data. Formerly, the main purpose of the monitoring carried out during tunnel construction was to measure the ground pressures acting on the tunnel lining. Modern tunneling practice emphasizes the monitoring of the displacements occurring during and after the construction. Panet and Sulem for determining of deformations in tunnels have assumed that "The tunnel has a circular cross section and around the tunnel, the rock is homogeneous and isotropic and also the tunnel is deep enough to consider that the stress distribution is homogenous". But in almost real cases, the stresses distribution around the tunnel is not homogeneous and isotropic. In this study, for modification of the Panet and Sulem equations, some equations are proposed in case of nonhomogeneous and anisotropic for generalizing of the problem. Seymareh power tunnel which is considered as a case study is a part of the powerhouse waterways system of the Seymareh dam and hydroelectric power plant project. The project is located in west of Iran. The monitoring data of power tunnel which are collected during excavation of tunnel is compared with the results of numerical modelling and analytical solution results as well as. The results obtained from comparison show although the field data, which are collected through the monitoring, are very close to the analytical solution results (approximately), but there is a significant difference between both of them and numerical modelling results. It was predictable; because the influence of the other activities such as excavation of shaft and surge tank in the numerical analysis and also analytical solution are not considered. It is obvious that other activities such as excavation of shaft and surge tank and also excavation of mean tunnel from other direction which were under operation at the same time can effect on the results of monitoring. On the other hand, the initial data which are used in numerical analysis and analytical solution are not quite accurate; because they are extracted as a representative of the rock mass of region, not for a particular section. However the goal of this study is development of analytical solution of deformation in tunnels on general conditions and pursuit of the study could be leaded to more development in this field.