Caracterização geológico-geofísica do meio aqüífero fissural : uma contribuição aos modelos de fluxo e armazenamento de água subterrânea

Northeastern Brazil is mainly formed by crystalline terrains (around 60% in area). Moreover, this region presents a semi-arid climate so that it is periodically subject to drought seasons. Furthermore, ground water quality extracted fromwells usually presents poor quality because of their high salinity contents. Nevertheless, ground water is still a very important source of water for human and animal consumption in this region. Well sitting in hard rocks terrains in Northeastern Brazil offers a mean success index of aboul 60%, given that a successful siting is defined by a well producing at least 0.5 m³/h. This low index reveals lack of knowledga about the true conditions of storage and percolation of ground water in crystalline rocks. Two models for structures storing and producing ground water in crystalline rocks in Northeastem Brazil have been proposed in the literature. The first model,tradnionally used for well sitting since the sixties are controlled by faults or fractures zones. This model is commonly referred, in Brazilian hydrogeological literature, as the "creek-crack" model (riacho-fenda in Portuguese). Sites appearing to present dense drainage network are preferred for water well siting - particularly at points where the drainages cross-cul each other. Field follow up work is usually based only on geological criteria. The second model is the "eluvio-alluvial through" (calha eluvio-aluvionar in Portuguese); it is also described in the literature but it is not yet incorporated in well sitting practice. This model is based on the hypothesis that reclilinear drainages can also be controlled by the folietion of the rock. Eventually, depending upon the degree of weathering, a through-shaped structure filled with sediments (alluvium and regolith) can be developed which can store and water can be produced from. Using severalfield case studies, this Thesis presents a thorough analysis ofthe two above cited models and proposes a new model. The analysis is based on an integrated methodological approach using geophysics and structural geology. Both land (Resitiviy and Ground Penetrating Radar- GPR) and aerogeophysical (magnetics and frequency domain eletromagnetics) surveys were used. Slructural analysis emphasized neolectonic aspects; in general, itwas found that fractures in the E-W direction are relatively open, as compared to fracturas inthe N-S direction, probably because E-W fractures were opened by the neotectonic stress regime in Northeastern Brazil, which is controlled by E-W compression and N-S extension. The riacho-fenda model is valid where drainages are controlled by fractures. The degree of fracturing and associated weathering dictale the hydrogeological potential of the structure. Field work in structural analogues reveals that subvertical fractures show consistent directions both in outcrop and aerophotograph scales. Geophysical surveys reveal subvertical conductive anomalies associated to the fracture network controlling the drainage; one of the borders of the conductive anomaly usually coincide wih the drainage. An aspect of particular importance to the validation of fracture control are the possible presence of relalively deep conductive anomalies wihoul continuation or propagalion to the surface. The conductive nature of lhe anomaly is due to the presence of wealhered rock and sedirnenls (alluvium and/or regolilh) storing ground waler which occur associated to the fracture network. Magnetic surveys are not very sensisnive to these structures.lf soil or covering sedirnents are resislive (> 100 Ohm.m), GPR can ba used to image precisely lhe fracture network. A major limialion of riacho-fenda model, revealed by GPR images, is associated to the fact thal subhorizontal fractures do play a very important role in connecting the fracture network, besides connect shallow recharge zones to relalively deep subvertical frecture zones. Iffractures play just a secondary control on the drainage, however, r/acho-fenda model may have a very limiled validny; in these cases, large portions oflhe drainage do nol coincide wilh frectures and mosl oflhewells localed in lhe drainage surrounding would resull dry. Usually, a secondary conlrol on lhe drainage by Ihefraclure networkcan be revealed only wilh detailed geophysical survey. The calha elClv1o-aluvlonarmodel is valid where drainages are conlrolled by folialion. The degree 01 wealhering 01 lhe lolialion planes dictales lhe hydrogeological polenlial 01 lhe slruclure. Outcrop analysis reveals Ihal lolialion and drainage direclions are parallel and Ihal no Iraclures, orfraclures wilh diflerent directions 01 lhe drainage direclion occur. Geophysical surveys reveal conduclive anomalies in a slab lorm associaled 10 lhe Ihrough 01 lhe wealhered rock and sedimenls (alluvium and/or regolith). Magnelic surveys can ofler a very good conlrol on lolialion direclion. An importanl aspect 10 validale lolialion conlrol are lhe presence 01 conductive anomalies showing shallow and deep portions area which are linked. Illhere is an exlensive soil cover, r/acho-fenda and calha eIClv1o-aluv/onar conlrols can be easily misinlerpreled in lhe absence 01 geophysical conlrol. Certainly, Ihis lacl could explain at leasl a part of lhe failure index in well sitting. The model wealhering sack (bolsllo de Intempertsmo in Portuguese) is proposed to explain cases where a very inlensive wealhering occur over lhe crystalline rock so Ihal a secondary inlerslilial porosity is crealed. The waler is Ihen stored in lhe porous of lhe regolilh in a similar mannerlo sedimentary rocks. A possible example ofthis model was delecled by using land geophysical survey where a relalivelyvery deep isolaled conduclive anomaly, in a slab form, was delected. Iflhis structure does store ground waler, certainly Ihere must be a link 01 lhe deep slructure wilh lhe surface in orderlo provide walerfeeding. This model mighl explain anomalous waler yields as greal as 50 m³/h Ihalsomelimescan occur in crystalline rocks in Northeaslern Brazil

Reavaliação de critérios estruturais na hidrogeologia de terrenos cristalinos, com ênfase na neotectônica e sensoriamento remoto

Hydrogeological prospecting in Northeast Brazil and in other crystalline terrains has been developed on the basis of structural and regional geology concepts that date back to the 50-60 decades and, as such, demand a natural re-evaluation and update. In this kind of terrain, the percolation and accumulation of ground water are controlled by fractures and other types of discontinuities, such as foliations and geological contacts that, through weathering, impart porosity and permeability to the rocks, allowing water flow and storage. Several factors should be considered in the process of locating water wells, as discussed in the literature. Among these, the kind of structures, fracture geometry (including aperture and connectivity) and their geological and chronological context. It is important to correlate fracture systems with the regional neotectonic framework. Fractures at low angle (sub parallel) with the principal stress axis (s1) are those which tend to open (actually they work as tension joints) and, in principle, would present major hydric potential; in the opposite side, fractures at high angle to s1 would behave as closed by a compressional component. Fractures diagonal to the compression and tension axes correspond to shear fractures and, due to their connectivity with second fractures, are also important in terms of hydric potential. Uplift followed by terrain denudation leads to decompression and a general tendency to open (aided by weathering processes) fractures and other rock discontinuities, at different orientations. Low angle fractures, formed in this context, are equally important to increase connectivity, collection of water and recharge of the aquifer systems. In a general way, an opening component (neotectonic or by terrain decompression) and several models to increase fracture connectivity correlate with a greater hydric potential of these structures. Together with parallel research, this thesis addresses models of ground water occurrence in crystalline terrains, either improving well established concepts like the (Riacho-Fenda model), but also stressing other possibilities, like the role of alluvium and paleo-regoliths (the Calha Elúvio-Aluvionar model) and of strongly altered, permo-porous zones placed at variable depths below the present surface, flanking several types of discontinuities, especially interconnected fracture arrays (the Bolsões de Intemperismo model). Different methodological approaches are also discussed in order to improve success rates in the location of water wells in crystalline terrains. In this methodological review, a number of case studies were selected in the eastern domain of the State of Rio Grande do Norte, involving the localities of Santa Cruz, Santo Antônio, Serrinha, Nova Cruz, Montanhas, Lagoa de Pedras and Lagoa Salgada. Besides the neotectonic analysis of brittle structures, this Thesis addresses the validation of remote sensing as a tool for ground water prospecting. Several techniques were tested in order to detect and select areas with higher potential for ground water accumulation, using Landsat 5-TM and RADARSAT images, besides conventional aerial photos. A number of filters were tested to emphasize lineaments in the images, improving their discrimination, to identify areas with higher overburden humidity, which could reflect subsurface water accumulation, as well as alluvium and other sedimentary covers that might act as recharge zones. The work started with a regional analysis with the orbital images, followed by analysis of aerial photos, up to a detailed structural study of rock exposures in the terrain. This last step involved the analysis of outcrops surrounding wells (in a ray of approximately 10 to 100 m) with distinct productivities, including dry examples. At the level required for detail, it was not possible to accomplish a statistical approach using the available well data catalogs, which lack the desired specific information. The methodology worked out in this Thesis must undergo a testing phase through location of new water wells. An increase in the success rates as desired will led to a further consolidation step with wider divulgation of the methodology to private companies and governmental agencies involved in ground water prospecting in crystalline terrains

Estudo de Fractalidade e Evolução Dinâmica de Sistemas Complexos

In this work, the study of some complex systems is done with use of two distinct procedures. In the first part, we have studied the usage of Wavelet transform on analysis and characterization of (multi)fractal time series. We have test the reliability of Wavelet Transform Modulus Maxima method (WTMM) in respect to the multifractal formalism, trough the calculation of the singularity spectrum of time series whose fractality is well known a priori. Next, we have use the Wavelet Transform Modulus Maxima method to study the fractality of lungs crackles sounds, a biological time series. Since the crackles sounds are due to the opening of a pulmonary airway bronchi, bronchioles and alveoli which was initially closed, we can get information on the phenomenon of the airway opening cascade of the whole lung. Once this phenomenon is associated with the pulmonar tree architecture, which displays fractal geometry, the analysis and fractal characterization of this noise may provide us with important parameters for comparison between healthy lungs and those affected by disorders that affect the geometry of the tree lung, such as the obstructive and parenchymal degenerative diseases, which occurs, for example, in pulmonary emphysema. In the second part, we study a site percolation model for square lattices, where the percolating cluster grows governed by a control rule, corresponding to a method of automatic search. In this model of percolation, which have characteristics of self-organized criticality, the method does not use the automated search on Leaths algorithm. It uses the following control rule: pt+1 = pt + k(Rc − Rt), where p is the probability of percolation, k is a kinetic parameter where 0 < k < 1 and R is the fraction of percolating finite square lattices with side L, LxL. This rule provides a time series corresponding to the dynamical evolution of the system, in particular the likelihood of percolation p. We proceed an analysis of scaling of the signal obtained in this way. The model used here enables the study of the automatic search method used for site percolation in square lattices, evaluating the dynamics of their parameters when the system goes to the critical point. It shows that the scaling of , the time elapsed until the system reaches the critical point, and tcor, the time required for the system loses its correlations, are both inversely proportional to k, the kinetic parameter of the control rule. We verify yet that the system has two different time scales after: one in which the system shows noise of type 1 f , indicating to be strongly correlated. Another in which it shows white noise, indicating that the correlation is lost. For large intervals of time the dynamics of the system shows ergodicity

Universaidade em sistemas da mecâniaestatística de não equilíbrio com estados absorventes e percolação geográfica

Complex systems have stimulated much interest in the scientific community in the last twenty years. Examples this area are the Domany-Kinzel cellular automaton and Contact Process that are studied in the first chapter this tesis. We determine the critical behavior of these systems using the spontaneous-search method and short-time dynamics (STD). Ours results confirm that the DKCA e CP belong to universality class of Directed Percolation. In the second chapter, we study the particle difusion in two models of stochastic sandpiles. We characterize the difusion through diffusion constant D, definite through in the relation h(x)2i = 2Dt. The results of our simulations, using finite size scalling and STD, show that the diffusion constant can be used to study critical properties. Both models belong to universality class of Conserved Directed Percolation. We also study that the mean-square particle displacement in time, and characterize its dependence on the initial configuration and particle density. In the third chapter, we introduce a computacional model, called Geographic Percolation, to study watersheds, fractals with aplications in various areas of science. In this model, sites of a network are assigned values between 0 and 1 following a given probability distribution, we order this values, keeping always its localization, and search pk site that percolate network. Once we find this site, we remove it from the network, and search for the next that has the network to percole newly. We repeat these steps until the complete occupation of the network. We study the model in 2 and 3 dimension, and compare the bidimensional case with networks form at start real data (Alps e Himalayas)

Estimativa de expoentes críticos em Percolação

In Percolation Theory, functions like the probability that a given site belongs to the infinite cluster, average size of clusters, etc. are described through power laws and critical exponents. This dissertation uses a method called Finite Size Scaling to provide a estimative of those exponents. The dissertation is divided in four parts. The first one briefly presents the main results for Site Percolation Theory for d = 2 dimension. Besides, some important quantities for the determination of the critical exponents and for the phase transistions understanding are defined. The second shows an introduction to the fractal concept, dimension and classification. Concluded the base of our study, in the third part the Scale Theory is mentioned, wich relates critical exponents and the quantities described in Chapter 2. In the last part, through the Finite Size Scaling method, we determine the critical exponents fi and. Based on them, we used the previous Chapter scale relations in order to determine the remaining critical exponents

Petrologia do plúton granitico serra verde, porção lestre do domínio seridó

The study area is located at the eastern-central portion of the Seridó Belt, on the interface between the Seridó Group Metasediments and the crystalline basement rocks of the Caicó Complex (RN). Petrographic and geochemical data allow us to define aspects related to the genesis and evolution of the Serra Verde Pluton magmas, which composes the goal of this dissertation The Serra Verde Pluton is a stock with outcropping area of about 25 km², which is intrusive into metasedimentary sequence and the basement gneisses. The pluton intrusion is sintectonic to the Brasiliano event, elongated along the NE direction, developing a cornue geometry. The rock is a monzogranite mainly composed by K-feldspar, plagioclase and quartz, which usually compose more than 85% of the modal analisys. The main mafic mineral is the biotite, while amphibole, sphene, epidote, opaque minerals, allanite, zircon and apatite occur as accessory minerals. It features still a latemagmatic paragenesis composed by chlorite, granular epidote, carbonates and muscovite, developed through the percolation of late CO2 and H2O rich fluids. Chemically, the Serra Verde Pluton rocks may be classified as metaluminous, of calc-alkaline affiliation, sometimes showing trondhjemític characteristics, with high Na2O (>4,5%), Sr (>400ppm) and Ba (>800ppm) and low K2O (≤3,0%), MgO (<1,0%), TiO2 (<0,5%), Rb (<90ppm), Y (≤16ppm) and Zr (≤13ppm). Micropetrographic evidences (mineral assembly and microtextures) indicate that the magma evolution occurred in moderated to high fO2 conditions, above the FMQ buffer. Thermo-barometric data obtained by minor elements geochemistry and the CIPW data, suggest a final/minimal pressure crystallization for the Serra Verde Pluton samples of about 3 to 5 kbar, liquidus temperature around 800o C, solidus temperature between 680o and 660o C. This data is compatible with those observed by many authors for the Neoproterozoic granites of the Seridó Belt. The group of analyzed data (Petrographic, microtextural and geochemical), suggests that the dominant process of the generation and evolution of the Serra Verde Granite magma was the fractional crystallization, probably from basement quartz-dioritic and tonalitic orthogneisses source

Produção de estrutura com gradiente de porosidade a partir da sinterização por plasma da mistura de pós Ti-Nb-Sn

Porous structures are being widely investigated for use in biomedical implants, aiming to mechanically integrate and functionally the implant inside the bone tissue. Moreover, this structure is also important for drugs that can be stored and can induce and accelerate the process of osseointegration. With the purpose to investigate this effect, Ti, Nb and Sn metal powders, were sintered by plasma using a hollow cathode discharge. Sintering was performed in argon plasma set at 4 mbar pressure and temperatures of 500 ° C, 600 ° C and 700 ° C. Samples were also sintered in the electrical resistance furnace at 1200 ° C in order to compare plasma sintering with the conventional method. It was observed that plasma samples sintered with the hollow cathode configuration showed a gradient in porosity, while the samples sintered in the resistive furnace did not. Furthermore, differences in the microstructure of the samples were found, were a surface with higher porosity and ales porous core were obtained at different temperatures. The percolation profile of distilled water and the chemical compositions of the porous layers of the plasma treated samples were the main results obtained. Based on these results, we can conclude that this structure is particularly important for application in the biomedical field such as scaffolds for drug delivery and implants

Desempenho de uma lagoa de infiltração na absorção de cheias e na recarga de aquífero

A type of macro drainage solution widely used in urban areas with predomi-nance of closed catchments (basins without outlet) is the implementation of detention and infiltration reservoirs (DIR). This type of solution has the main function of storing surface runoff and to promote soil infiltration and, consequently, aquifer recharge. The practice is to avoid floods in the drainage basin low-lying areas. The catchment waterproofing reduces the distributed groundwater recharge in urban areas, as is the case of Natal city, RN. However, the advantage of DIR is to concentrate the runoff and to promote aquifer recharge to an amount that can surpass the distributed natu-ral recharge. In this paper, we proposed studying a small urban drainage catchment, named Experimental Mirassol Watershed (EMW) in Natal, RN, whose outlet is a DIR. The rainfall-runoff transformation processes, water accumulation in DIR and the pro-cess of infiltration and percolation in the soil profile until the free aquifer were mod-eled and, from rainfall event observations, water levels in DIR and free aquifer water level measurements, and also, parameter values determination, it is was enabled to calibrate and modeling these combined processes. The mathematical modeling was carried out from two numerical models. We used the rainfall-runoff model developed by RIGHETTO (2014), and besides, we developed a one-dimensional model to simu-late the soil infiltration, percolation, redistribution soil water and groundwater in a combined system to the reservoir water balance. Continuous simulation was run over a period of eighteen months in time intervals of one minute. The drainage basin was discretized in blocks units as well as street reaches and the soil profile in vertical cells of 2 cm deep to a total depth of 30 m. The generated hydrographs were transformed into inlet volumes to the DIR and then, it was carried out water balance in these time intervals, considering infiltration and percolation of water in the soil profile. As a re-sult, we get to evaluate the storage water process in DIR as well as the infiltration of water, redistribution into the soil and the groundwater aquifer recharge, in continuous temporal simulation. We found that the DIR has good performance to storage excess water drainage and to contribute to the local aquifer recharge process (Aquifer Dunas / Barreiras).

Estudo da transição de fase da percolação através da entropia da informação

Various physical systems have dynamics that can be modeled by percolation processes. Percolation is used to study issues ranging from fluid diffusion through disordered media to fragmentation of a computer network caused by hacker attacks. A common feature of all of these systems is the presence of two non-coexistent regimes associated to certain properties of the system. For example: the disordered media can allow or not allow the flow of the fluid depending on its porosity. The change from one regime to another characterizes the percolation phase transition. The standard way of analyzing this transition uses the order parameter, a variable related to some characteristic of the system that exhibits zero value in one of the regimes and a nonzero value in the other. The proposal introduced in this thesis is that this phase transition can be investigated without the explicit use of the order parameter, but rather through the Shannon entropy. This entropy is a measure of the uncertainty degree in the information content of a probability distribution. The proposal is evaluated in the context of cluster formation in random graphs, and we apply the method to both classical percolation (Erd¨os- R´enyi) and explosive percolation. It is based in the computation of the entropy contained in the cluster size probability distribution and the results show that the transition critical point relates to the derivatives of the entropy. Furthermore, the difference between the smooth and abrupt aspects of the classical and explosive percolation transitions, respectively, is reinforced by the observation that the entropy has a maximum value in the classical transition critical point, while that correspondence does not occurs during the explosive percolation.

Estudo da transição de fase da percolação através da entropia da informação

Various physical systems have dynamics that can be modeled by percolation processes. Percolation is used to study issues ranging from fluid diffusion through disordered media to fragmentation of a computer network caused by hacker attacks. A common feature of all of these systems is the presence of two non-coexistent regimes associated to certain properties of the system. For example: the disordered media can allow or not allow the flow of the fluid depending on its porosity. The change from one regime to another characterizes the percolation phase transition. The standard way of analyzing this transition uses the order parameter, a variable related to some characteristic of the system that exhibits zero value in one of the regimes and a nonzero value in the other. The proposal introduced in this thesis is that this phase transition can be investigated without the explicit use of the order parameter, but rather through the Shannon entropy. This entropy is a measure of the uncertainty degree in the information content of a probability distribution. The proposal is evaluated in the context of cluster formation in random graphs, and we apply the method to both classical percolation (Erd¨os- R´enyi) and explosive percolation. It is based in the computation of the entropy contained in the cluster size probability distribution and the results show that the transition critical point relates to the derivatives of the entropy. Furthermore, the difference between the smooth and abrupt aspects of the classical and explosive percolation transitions, respectively, is reinforced by the observation that the entropy has a maximum value in the classical transition critical point, while that correspondence does not occurs during the explosive percolation.

Análise de modelos geológicos utilizando percolação dinâmica

In the present study we elaborated algorithms by using concepts from percolation theory which analyze the connectivity conditions in geological models of petroleum reservoirs. From the petrophysical parameters such as permeability, porosity, transmittivity and others, which may be generated by any statistical process, it is possible to determine the portion of the model with more connected cells, what the interconnected wells are, and the critical path between injector and source wells. This allows to classify the reservoir according to the modeled petrophysical parameters. This also make it possible to determine the percentage of the reservoir to which each well is connected. Generally, the connected regions and the respective minima and/or maxima in the occurrence of the petrophysical parameters studied constitute a good manner to characterize a reservoir volumetrically. Therefore, the algorithms allow to optimize the positioning of wells, offering a preview of the general conditions of the given model s connectivity. The intent is not to evaluate geological models, but to show how to interpret the deposits, how their petrophysical characteristics are spatially distributed, and how the connections between the several parts of the system are resolved, showing their critical paths and backbones. The execution of these algorithms allows us to know the properties of the model s connectivity before the work on reservoir flux simulation is started

Estudo de adsorvente obtido de terras diatomáceas para remoção de metais pesados em águas produzidas: processamento, caracterização e dinâmica em leito fixo

The acceleration of industrial growth in recent decades on all continents aroused the interest of the companies to counter the impacts produced on the environment, spurred primarily by major disasters in the petroleum industry. In this context, the water produced is responsible for the largest volume of effluent from the production and extraction of oil and natural gas. This effluent has in its composition some critical components such as inorganic salts, heavy metals (Fe, Cu, Zn, Pb, Cd, ), presence of oil and chemicals added in the various production processes. In response to impact, have been triggered by research alternative adsorbent materials for water treatment and water produced, in order to removing oils and acids and heavy metals. Many surveys of diatomaceous earth (diatomite) in Brazil involve studies on the physico-chemical, mineral deposits, extraction, processing and applications. The official estimated Jazi are around 2.5 million tonnes, the main located in the states of Bahia (44%) and Rio Grande do Norte (37,4%). Moreover, these two states appear as large offshore producers, earning a prominent role in research of adsorbents such as diatomite for treatment of water produced. Its main applications are as an agent of filtration, adsorption of oils and greases, industrial load and thermal insulator. The objective of this work was the processing and characterization of diatomite diatomaceous earth obtained from the municipality of Macaíba-RN (known locally as tabatinga) as a low cost regenerative adsorbent for removal of heavy metals in the application of water produced treatment. In this work we adopted a methodology for batch processing, practiced by small businesses located in producing regions of Brazil. The characterization was made by X-ray diffraction (XRD), scanning electron microscopy (SEM) and specific surface area (BET). Research conducted showed that the improvement process used was effective for small volume production of diatomite concentrated. The diatomite obtained was treated by calcination at temperature of 900 oC for 2 hours, with and without fluxing Na2CO3 (4%), according to optimal results in the literature. Column adsorption experiments were conducted to percolation of the in nature, calcined and calcined fluxing diatomites. Effluent was used as a saline solution containing ions of Cu, Zn, Na, Ca and Mg simulating the composition of produced waters in the state of Rio Grande do Norte, Brazil. The breakthrough curves for simultaneous removal of copper ions and zinc as a result, 84.3% for calcined diatomite and diatomite with 97.3 % for fluxing. The calcined fluxing diatomite was more efficient permeability through the bed and removal of copper and zinc ions. The fresh diatomite had trouble with the permeability through the bed under the conditions tested, compared with the other obtained diatomite. The results are presented as promising for application in the petroleum industry

Estudo das propriedades críticas do processo epidêmico por par com difusão de pares

The pair contact process - PCP is a nonequilibrium stochastic model which, like the basic contact process - CP, exhibits a phase transition to an absorbing state. While the absorbing state CP corresponds to a unique configuration (empty lattice), the PCP process infinitely many. Numerical and theoretical studies, nevertheless, indicate that the PCP belongs to the same universality class as the CP (direct percolation class), but with anomalies in the critical spreading dynamics. An infinite number of absorbing configurations arise in the PCP because all process (creation and annihilation) require a nearest-neighbor pair of particles. The diffusive pair contact process - PCPD) was proposed by Grassberger in 1982. But the interest in the problem follows its rediscovery by the Langevin description. On the basis of numerical results and renormalization group arguments, Carlon, Henkel and Schollwöck (2001), suggested that certain critical exponents in the PCPD had values similar to those of the party-conserving - PC class. On the other hand, Hinrichsen (2001), reported simulation results inconsistent with the PC class, and proposed that the PCPD belongs to a new universality class. The controversy regarding the universality of the PCPD remains unresolved. In the PCPD, a nearest-neighbor pair of particles is necessary for the process of creation and annihilation, but the particles to diffuse individually. In this work we study the PCPD with diffusion of pair, in which isolated particles cannot move; a nearest-neighbor pair diffuses as a unit. Using quasistationary simulation, we determined with good precision the critical point and critical exponents for three values of the diffusive probability: D=0.5 and D=0.1. For D=0.5: PC=0.89007(3), β/v=0.252(9), z=1.573(1), =1.10(2), m=1.1758(24). For D=0.1: PC=0.9172(1), β/v=0.252(9), z=1.579(11), =1.11(4), m=1.173(4)