Variability is one of the major characteristics of Active Galactic Nuclei (AGN), and it is used for understanding the energy generation mechanism in the center of AGN and/or related physical phenomena. It it known that there exists a time lag between AGN light curves simultaneously observed at different wavelengths, which can be used as a tool to estimate the size of the area that produce the radiation. In this paper, We present long term near-infrared variability of optically bright type 1 AGN using the Wide-field Infrared Survey Explorer data. From the Milliquas catalogue v6.4, 73 type 1 QSOs/AGN and 140 quasar candidates are selected that are brighter than 18 mag in optical and located within 5 degree around the ecliptic poles. Light curves in the W1 band (3.4 ㎛) and W2 band (4.6 ㎛) during the period of 2010-2019 were constructed for these objects by extracting multi-epoch photometry data from WISE and NEOWISE all sky survey database. Variability was analyzed based on the excess variance and the probability Pvar . Applying both criteria, the numbers of variable objects are 19 (i.e., 26%) for confirmed AGN and 12 (i.e., 9%) for AGN candidates. The characteristic time scale of the variability (τ) and the variability amplitude (σ) were derived by fitting the DRW model to W1 and W2 light curves. No significant correlation is found between the W1/W2 magnitude and the derived variability parameters. Based on the subsample that are identified in the X-ray source catalog, there exists little correlation between the X-ray luminosity and the variability parameters. We also found four AGN with changing W1-W2 color.